• hdu 1853 Cyclic Tour KM


    #include <stdio.h>
    #include <string.h>
    #define M 310
    #define inf 0x3f3f3f3f
    
    int n,nx,ny;
    int link[M],lx[M],ly[M],slack[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;
                    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|| w[link[i]][i]==-inf)
                return 0;
            else res += w[link[i]][i];;
        }
    
        return res;
    }
    
    int main()
    {
        int n,m;
        while(scanf("%d%d",&n,&m)!=EOF)
        {
            int i,j,x,y,z;
            for(i=1;i<=n;i++)
                for(j=1;j<=n;j++)
                    w[i][j]=-inf;
            for(i=1;i<=m;i++)
            {
                scanf("%d%d%d",&x,&y,&z);
                if(-z>w[x][y])
                    w[x][y]=-z;
            }
            nx=ny=n;
            int ans=KM();
            if(ans==0)
                printf("-1
    ");
            else printf("%d
    ",-ans);
        }
        return 0;
    }
    


     

  • 相关阅读:
    android OnActivityResult()的调用时机问题
    把汉化的eclipse还原为英语版的小绝招
    sunPKCS11加载动态库(转)
    android 源码分析
    android notification应用
    期待快点过去
    Linux 的Root 帐号
    C#之CMD
    更新Citrix Merchandising Server 2.2的默认Updater
    CRL的验证
  • 原文地址:https://www.cnblogs.com/vermouth/p/3710157.html
Copyright © 2020-2023  润新知