hdu4307

好题，详细题解在这里http://blog.csdn.net/weiguang_123/article/details/8077385

这里回顾一下：

当i和j都在一个集合里会产生新的收益，是经典题直接建

当i和j不在同一个集合产生新的收益，需要点是二分图XY并把X或Y点集s-t反建

这里是类似的（以前的总结：http://www.cnblogs.com/phile/p/4473226.html）

#include<bits/stdc++.h>

using namespace std;
typedef long long ll;
const int inf=1e9+7;
struct way{int po,flow,next;} e[2400010];
int t,n,p[1010],pre[1010],numh[1010],cur[1010],h[1010],d[1010],m,len;

void add(int x,int y,int f)
{
     e[++len].po=y;
     e[len].flow=f;
     e[len].next=p[x];
     p[x]=len;
}
void build(int x, int y, int f)
{
     add(x,y,f);
     add(y,x,0);
}
ll sap()
{
    memset(numh,0,sizeof(numh));
    memset(h,0,sizeof(h));
    numh[0]=t+1;
    for (int i=0; i<=t; i++) cur[i]=p[i];
    int j,u=0,neck=inf; ll s=0;
    while (h[0]<t+1)
    {
          d[u]=neck;
          bool ch=1;
          for (int i=cur[u]; i!=-1; i=e[i].next)
          {
              j=e[i].po;
              if (e[i].flow>0&&h[u]==h[j]+1)
              {
                 neck=min(neck,e[i].flow);
                 cur[u]=i;
                 pre[j]=u;  u=j;
                 if (u==t)
                 {
                    s+=neck;
                    while (u)
                    {
                          u=pre[u];
                          j=cur[u];
                          e[j].flow-=neck;
                          e[j^1].flow+=neck;
                    }
                    neck=inf;
                 }
                 ch=0;
                 break;
              }
          }
          if (ch)
          {
             if (--numh[h[u]]==0) return s;
             int q=-1,tmp=t;
             for (int i=p[u]; i!=-1; i=e[i].next)
             {
                 j=e[i].po;
                 if (e[i].flow&&h[j]<tmp)
                 {
                    tmp=h[j];
                    q=i;
                 }
             }
             cur[u]=q; h[u]=tmp+1;
             numh[h[u]]++;
             if (u)
             {
                u=pre[u];
                neck=d[u];
             }
          }
    }
    return s;
}

int main()
{
    int cas;
    scanf("%d",&cas);
    while (cas--)
    {
        len=-1;
        memset(p,255,sizeof(p));
        scanf("%d",&n);
        ll ans=0;
        for (int i=1; i<=n; i++)
        {
            int s=0;
            for (int j=1; j<=n; j++)
            {
                int x;
                scanf("%d",&x);
                s+=x; ans+=x;
                build(i,j,x);
            }
            build(0,i,s);
        }
        t=n+1;
        for (int i=1; i<=n; i++)
        {
            int x;
            scanf("%d",&x);
            build(i,t,x);
        }
        printf("%lld
",ans-sap());
    }
}