首页 > 代码库 > Acdream 1227 Beloved Sons KM
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; }
Acdream 1227 Beloved Sons KM
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。