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+;
 struct way{int po,flow,next;} e[];
 int t,n,p[],pre[],numh[],cur[],h[],d[],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,);
 }
 ll sap()
 {
     memset(numh,,sizeof(numh));
     memset(h,,sizeof(h));
     numh[]=t+;
     for (int i=; i<=t; i++) cur[i]=p[i];
     int j,u=,neck=inf; ll s=;
     while (h[]<t+)
     {
           d[u]=neck;
           bool ch=;
           for (int i=cur[u]; i!=-; i=e[i].next)
           {
               j=e[i].po;
               if (e[i].flow>&&h[u]==h[j]+)
               {
                  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^].flow+=neck;
                     }
                     neck=inf;
                  }
                  ch=;
                  break;
               }
           }
           if (ch)
           {
              if (--numh[h[u]]==) return s;
              int q=-,tmp=t;
              for (int i=p[u]; i!=-; 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+;
              numh[h[u]]++;
              if (u)
              {
                 u=pre[u];
                 neck=d[u];
              }
           }
     }
     return s;
 }

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