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;}