用“道”的思想解决费用流问题---取/不取皆是取 （有下界->有上界） / ACdreamoj 1171

题意： 给一个矩阵，给出约束：i（0<i<n）行至少去ai个数，j行至少取bi个数,要求取的数值之和最小。

开始一见，就直接建了二分图，但是，发现这是有下界无上界最小费用流问题，肿么办。。。问题转化：所谓正难则反！现在某行/列要至少取k个，总和最小，不就是那行/列最多留下K个，使留下的和最大？其实也就是最多取k个，使值最大，转化为下界为0，有上界的最大费用问题（普通问题）。“取”，“不取”，本质都是一样的，正是“无为”的思想！取，则最小；不取，最大。道也。道之道非常道，名可名非常名~~

还有一点，转化之后，最大费用时未必最大流。解决方法有二：

其一：每次增广后，加判断，若费用开始递减，则跳出，此时取最大。（据说二分图费用是先增后减函数：每次增广，当费用最大的时候，但是这时候流量不是最大，所以减小费用来增大流量，不知道一般图是不是。。。）

其二：释放法，X部所有点直接向汇点连边，费用0，流量Inf,我感觉这样，当X部还有流量的时候，直接就向汇点释放了，所有必是最大费用（不会再减少了）。

俩种方法我都试过，AC。

#include<cstdio>
#include<iostream>
#include<queue>
#include<cstring>
#include<string>
using namespace std;
const int maxv=200;
const int maxe=200*200*2+800;
const int inf=0x3f3f3f3f;
int nume=0;int e[maxe][4];int head[maxv];
int n,m;int ss,tt;
int val[105][105];
void inline adde(int i,int j,int c,int w)
{
    e[nume][0]=j;e[nume][1]=head[i];head[i]=nume;
    e[nume][2]=c;e[nume++][3]=w;
    e[nume][0]=i;e[nume][1]=head[j];head[j]=nume;
    e[nume][2]=0;e[nume++][3]=-w;
}
int inq[maxv];int pre[maxv];int prv[maxv];
int d[maxv];
bool spfa(int &sum,int &flow)
{
    for(int i=0;i<=tt;i++)
          {
              inq[i]=0;
              d[i]=inf;
          }
    queue<int>q;
    q.push(ss);
    inq[ss]=1;
    d[ss]=0;
    while(!q.empty())
    {
        int cur=q.front();
        q.pop();
        inq[cur]=0;
        for(int i=head[cur];i!=-1;i=e[i][1])
        {
            int v=e[i][0];
            if(e[i][2]>0&&d[cur]+e[i][3]<d[v])
            {
                d[v]=d[cur]+e[i][3];
                pre[v]=i;
                prv[v]=cur;
                if(!inq[v])
                {
                    q.push(v);
                    inq[v]=1;
                }
            }
        }
    }
    if(d[tt]==inf)return 0;
    int cur=tt;
    int minf=inf;
    while(cur!=ss)
    {
        int fe=pre[cur];
        minf=e[fe][2]<minf?e[fe][2]:minf;
        cur=prv[cur];
    }
     cur=tt;
    while(cur!=ss)
    {
        e[pre[cur]][2]-=minf;
        e[pre[cur]^1][2]+=minf;
        cur=prv[cur];
    }
    flow+=minf;
    sum+=d[tt]*minf;
    return 1;
}
int mincost(int &flow)
{
    int sum=0;
   // int lastsum=0;
   while(spfa(sum,flow))
   {
       ;
     //  if(-lastsum>-sum)return lastsum;                  //取最值法
      //  lastsum=sum; // cout<<sum<<endl;
   }
    return sum;
}
int sum_all=0;
void init()
{
    nume=0; sum_all=0;
    ss=n+m; tt=n+m+1;
    for(int i=0;i<=tt;i++)
       head[i]=-1;
}
void read_build()
{
    for(int i=0;i<n;i++)
     for(int j=0;j<m;j++)
         {
             scanf("%d",&val[i][j]);
             sum_all+=val[i][j];
             adde(i,j+n,1,-val[i][j]);
         }
      int aa;
     for(int i=0;i<n;i++)
     {
        scanf("%d",&aa);
        adde(ss,i,m-aa,0);
        adde(i,tt,m-aa,0);                            //X部直接向汇点连边（容量够释放就行）
     }
       for(int i=0;i<m;i++)
     {
        scanf("%d",&aa);
        adde(i+n,tt,n-aa,0);
     }
   /*  for(int i=0;i<=m+n+1;i++)
       for(int j=head[i];j!=-1;j=e[j][1])
       {
           printf("%d->%d:f %dw %d\n",i,e[j][0],e[j][2],e[j][3]);
       }*/
}
int main()
{
    int T;
    cin>>T;
    while(T--)
    {
        scanf("%d%d",&n,&m);
        init();
        read_build();
        int flow=0;
        int ans=sum_all+mincost(flow);
         printf("%d\n",ans);
    }
    return 0;
}