SGU 210 Acdream 1227 Beloved Sons KM

题目链接:点击打开链接

题意:

给定n个人

每一个人的点权

以下n行i行表示第i个人能够获得哪些数(数字从1-n。且不能反复分配)

若这个人获得了数字则你能够获得他的权值。

要你能获得的权值和最大。

问:

输出每一个人应该获得哪个数字。若没有获得到数字则输出0.


思路:

KM,每一个人给每一个数字连一条边,边权是这个人的权值。

对于这个人不能获得的数字连一条边权为0的边。


#include <stdio.h>
#include <string.h>
#define M 410
#define inf 0x3f3f3f3f
 
int n,nx,ny;
int link[M],lx[M],ly[M],slack[M], linkx[M];    //lx,ly为顶标。nx,ny分别为x点集y点集的个数
int visx[M],visy[M],w[M][M];
 
int DFS(int x)
{
    visx[x] = 1;
    for (int y = 1;y <= ny;y ++)
    {
        if (visy[y])
            continue;
        int t = lx[x] + ly[y] - w[x][y];
        if (t == 0)       //
        {
            visy[y] = 1;
            if (link[y] == -1||DFS(link[y]))
            {
                link[y] = x;
                linkx[x] = y;
                return 1;
            }
        }
        else if (slack[y] > t)  //不在相等子图中slack 取最小的
            slack[y] = t;
    }
    return 0;
}
int KM()
{
    int i,j;
    memset (link,-1,sizeof(link));
    memset (ly,0,sizeof(ly));
    for (i = 1;i <= nx;i ++)            //lx初始化为与它关联边中最大的
        for (j = 1,lx[i] = -inf;j <= ny;j ++)
            if (w[i][j] > lx[i])
                lx[i] = w[i][j];
 
    for (int x = 1;x <= nx;x ++)
    {
        for (i = 1;i <= ny;i ++)
            slack[i] = inf;
        while (1)
        {
            memset (visx,0,sizeof(visx));
            memset (visy,0,sizeof(visy));
            if (DFS(x))     //若成功(找到了增广轨),则该点增广完毕。进入下一个点的增广
                break;  //若失败(没有找到增广轨),则须要改变一些点的标号,使得图中可行边的数量添加。

//方法为:将全部在增广轨中(就是在增广过程中遍历到)的X方点的标号全部减去一个常数d, //全部在增广轨中的Y方点的标号全部加上一个常数d int d = inf; for (i = 1;i <= ny;i ++) if (!visy[i]&&d > slack[i]) d = slack[i]; for (i = 1;i <= nx;i ++) if (visx[i]) lx[i] -= d; for (i = 1;i <= ny;i ++) //改动顶标后。要把全部不在交错树中的Y顶点的slack值都减去d if (visy[i]) ly[i] += d; else slack[i] -= d; } } int res = 0; for (i = 1;i <= ny;i ++) if (link[i] > -1) res += w[link[i]][i]; return res; } template <class T> inline bool rd(T &ret) { char c; int sgn; if(c=getchar(),c==EOF) return 0; while(c!='-'&&(c<'0'||c>'9')) c=getchar(); sgn=(c=='-')?-1:1; ret=(c=='-')?

0:(c-'0'); while(c=getchar(),c>='0'&&c<='9') ret=ret*10+(c-'0'); ret*=sgn; return 1; } template <class T> inline void pt(T x) { if (x <0) { putchar('-'); x = -x; } if(x>9) pt(x/10); putchar(x%10+'0'); } using namespace std; int a[401]; void build(){ for(int i = 1; i <= n; i++) rd(a[i]), a[i] = a[i]*a[i]; for(int i = 1, siz, u; i <= n; i++) { memset(w[i], 0, sizeof w[i]); rd(siz); while(siz--) { rd(u); w[i][u] = a[i]; } } nx = ny = n; KM(); } int main(){ while(rd(n)){ build(); for(int i = 1; i <= n; i++) { if(w[i][linkx[i]]) pt(linkx[i]); else putchar('0'); i == n?puts("") : putchar(' '); } } return 0; }



原文地址:https://www.cnblogs.com/gavanwanggw/p/6920692.html