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



  • 相关阅读:
    自动构建部署
    EF 性能调优
    断点续传
    gis 相关资料
    easyui 特殊操作
    KJ面试
    前端面试题汇总
    es6之扩展运算符 三个点(...)
    vue.js开发环境搭建
    gulp 环境搭建
  • 原文地址:https://www.cnblogs.com/gavanwanggw/p/6920692.html
Copyright © 2020-2023  润新知