hdu 3435 A new Graph Game

//邻接表 2859MS
#include<cstdio>
#include<cstring>
using namespace std;
const int maxn=1010*2;
const int maxm=10010*10;//最大顶点数和边数
const int maxl=999999999;
inline int Min(int a,int b)
{
    return a<b?a:b;
}
inline int Max(int a,int b)
{
    return a>b?a:b;
}
struct st
{
    int y,d;
    int ne;
    int bro;
    int f;
} e[maxm];
int ee;
int st[maxn];
int n,m;
int F;
void addedge(int x,int y,int d,int f)
{
    //给顶点x和y间添加一条费用d,流量f的边
    e[ee].y=y;
    e[ee].d=d;
    e[ee].ne=st[x];
    e[ee].f=f;
    st[x]=ee++;
    e[ee].y=x;
    e[ee].d=-1*d;
    e[ee].ne=st[y];
    e[ee].f=0;
    st[y]=ee++;
    e[ee-2].bro=ee-1;
    e[ee-1].bro=ee-2;
}
int d[maxn],p[maxn];
//spfa所用到起点的最短距离(这里距离相当于cost)和路径记录之前的一个节点
int c[maxn];//spfa所用数组:是否在队列中
int que[maxn],head,tail;//spfa专用队列
int spfa(int sx,int ex)//求sx到ex的一次费用增广
{
    //如果没有增广路就返回maxl 否则返回费用
    int i,j,k;
    for (i=0; i<maxn; i++) d[i]=maxl;
    memset(c,0,sizeof(c));//初始化都没进
    d[sx]=0;
    que[head=0]=sx;
    tail=1;
    c[sx]=1;
    while (head!=tail)
    {
        k=que[head++];
        head%=maxn;
        c[k]=0;
        for (i=st[k]; i!=-1; i=e[i].ne) if (e[i].f)
                if (d[k]+e[i].d<d[e[i].y])
                {
                    d[e[i].y]=d[k]+e[i].d;
                    p[e[i].y]=i;
                    if (c[e[i].y]==0)
                    {
                        c[e[i].y]=1;
                        if (e[i].d<0)
                        {
                            head=(head-1+maxn)%maxn;
                            que[head]=e[i].y;
                        }
                        else
                        {
                            que[tail++]=e[i].y;
                            tail%=maxn;
                        }
                    }
                }
    }
    if (d[ex]==maxl) return maxl;
//如果无法到达终点返回maxl
    k=maxl;
    for (i=ex; i!=sx; i=e[e[p[i]].bro].y)
        k=Min(k,e[p[i]].f); //计算流
    //sum+=k; //sum计算流量
    F+=k;
    for (i=ex; i!=sx; i=e[e[p[i]].bro].y) //增加反向边
    {
        e[p[i]].f-=k;
        e[e[p[i]].bro].f+=k;
    }
    return d[ex]*k;
//返回费用为流大小*路径长度(cost累加)
}
int main()
{
    int T;
    int cas=1;
    int m;
    scanf("%d",&T);
    while(T--)
    {
        int nn,mm;
        scanf("%d%d",&nn,&mm);
        memset(st,-1,sizeof(st));
        ee=0;

        for(int i=1; i<=mm; i++)
        {
            int x,y,c;
            scanf("%d%d%d",&x,&y,&c);
            addedge(x,y+nn,c,1);
            addedge(y,x+nn,c,1);
        }

        for(int i=1; i<=nn; i++)
        {
            addedge(0,i,0,1);
            addedge(i+nn,nn+nn+1,0,1);
        }

        n=nn*2+2;
        F=0;
        int tot=0,k;
        while ((k=spfa(0,n-1))!=maxl) tot+=k;

        //printf("%")
        if(F!=nn)  printf("Case %d: NO\n",cas++);
        else printf("Case %d: %d\n",cas++,tot);
    }
    return 0;
}

//c++ 2390ms
//完美匹配
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int inf=1<<28;
int g[1005][1005],lx[1005],ly[1005];//顶点标号
bool sx[1005],sy[1005];//是否已经搜索过
int link[1005],n,stack[1005];
bool path(int k)
{
    sx[k]=true;
    for(int i=1; i<=n; i++)
    {
        if(sy[i])  continue;
        int t=ly[i]+lx[k]-g[k][i];
        if(t==0)
        {
            sy[i]=1;
            if(link[i]==-1||path(link[i]))
            {
                link[i]=k;
                return true;
            }
        }
        else if(stack[i]>t) stack[i]=t;
    }
    return false;
}
int BestMatch()
{
    int d,sum;
    memset(ly,0,sizeof(ly));
    memset(link,-1,sizeof(link));
    for(int i=1; i<=n; i++)
    {
        lx[i]=-inf;
        for(int j=1; j<=n; j++)
            if(lx[i]<g[i][j]&&g[i][j]!=0)  lx[i]=g[i][j];
    }
    for(int k=1; k<=n; k++)
    {
        for (int i=1; i<=n; i++) stack[i]=inf;
        while(1)
        {
            memset(sx,0,sizeof(sx));
            memset(sy,0,sizeof(sy));
            if(path(k))  break;
            d=inf;
            for(int i=1; i<=n; i++)
                if (!sy[i]&&stack[i]<d) d=stack[i];
            for(int i=1; i<=n; i++)
                if(sx[i]) lx[i]-=d;
            for(int i=1; i<=n; i++)
                if(sy[i]) ly[i]+=d;
                else  stack[i]-=d;
        }
    }
    sum=0;
    for(int i=1; i<=n; i++)
    {
        if(link[i]==-1||g[link[i]][i]==-inf)  return -1;
        sum+=g[link[i]][i];
    }
    return -sum;
}
int main()
{
    int T;
    int cas=1;
    int m;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&n,&m);
        for(int i=1; i<=n; i++)
            for(int j=1; j<=n; j++)
                g[i][j]=inf;
        for(int i=1; i<=m; i++)
        {
            int x,y,c;
            scanf("%d%d%d",&x,&y,&c);
            g[x][y]=g[y][x]=min(g[x][y],c);
        }
        for(int i=1; i<=n; i++)
            for(int j=1; j<=n; j++) g[i][j]=-g[i][j];

        int t=BestMatch();
        if(t==-1)  printf("Case %d: NO\n",cas++);
        else printf("Case %d: %d\n",cas++,t);
    }
    return 0;
}