uva11082 Matrix Decompressing 最大流

题意：已知一个n*m的矩阵每行没列之和，得出一个满足条件的矩阵。

分析：建一个二分图，每一行对应一个X节点，每一列对应一个Y节点，增加一个源点s，s到X节点容量为对应行的和，Y到汇点t的容量为对应列的和，Xi到Yj的容量为19.由于任意元素[1,20]，为了用容量表示值，所有元素减1，就变成[0,19]，所以把容量上限设为19；

#include<iostream>
#include<string>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<iomanip>
#include<map>
#include<algorithm>
#include<queue>
#include<set>
#define inf 10000000
#define pi acos(-1.0)
#define eps 1e-8
#define seed 131
using namespace std;
typedef pair<int,int> pii;
typedef unsigned long long ULL;
typedef long long LL;
const int maxn=100005;
int n,m;
struct Edge
{
    int from,to,cap,flow;
    Edge(int a,int b,int c,int d):from(a),to(b),cap(c),flow(d){}
};
vector<Edge>edge;
vector<int>g[50];
int a[50];
int p[50];
int s0=0,t0=45;
int dr[25];
int dc[25];
void init()
{
    for(int i=0;i<=t0;i++)
        g[i].clear();
    edge.clear();
}
void addEdge(int a,int b,int c)
{
    edge.push_back(Edge(a,b,c,0));
    edge.push_back(Edge(b,a,0,0));
    int e=edge.size();
    g[a].push_back(e-2);
    g[b].push_back(e-1);
}
void EK()
{
    for(;;)
    {
        memset(a,0,sizeof(a));
        a[s0]=inf;
        queue<int>que;
        que.push(s0);
        while(!que.empty())
        {
            int u=que.front();
            que.pop();
            int len=g[u].size();
            for(int i=0;i<len;i++)
            {
                Edge& e=edge[g[u][i]];
                if(a[e.to]==0&&e.cap>e.flow)
                {
                    a[e.to]=min(a[u],e.cap-e.flow);
                    p[e.to]=g[u][i];
                    que.push(e.to);
                }
            }
            if(a[t0])
                break;
        }
        if(a[t0]==0)
            break;
        for(int i=t0;i!=s0;i=edge[p[i]].from)
        {
            edge[p[i]].flow+=a[t0];
            edge[p[i]^1].flow-=a[t0];
        }
    }
}
int main()
{
    int t;
    scanf("%d",&t);
    int tot;
    for(int cas=1;cas<=t;cas++)
    {
        tot=0;
        scanf("%d%d",&n,&m);
        init();
        dr[0]=dc[0]=0;
        for(int i=1;i<=n;i++)
            scanf("%d",&dr[i]);
        for(int i=1;i<=m;i++)
            scanf("%d",&dc[i]);
        for(int i=n;i>=1;i--)
        {
            dr[i]=dr[i]-dr[i-1]-m;
            tot+=dr[i];
        }
        for(int i=m;i>=1;i--)
        {
            dc[i]=dc[i]-dc[i-1]-n;
            tot+=dc[i];
        }
        int p;
        int q=n+1;
        for(p=1;p<=n;p++)
        {
            addEdge(s0,p,dr[p]);
            for(int i=n+1;i<=n+m;i++)
                addEdge(p,i,19);
        }
        for(int i=n+1;i<=n+m;i++)
            addEdge(i,t0,dc[i-n]);
        EK();
        printf("Matrix %d\n",cas);
        for(int i=1;i<=n;i++)
        {
            int len=g[i].size();
            int o=0;
            for(int j=1;j<len;j++)
            {
                if(o==0)
                    printf("%d",edge[g[i][j]].flow+1);
                else
                    printf(" %d",edge[g[i][j]].flow+1);
                o++;
            }
            printf("\n");
        }
        if(cas!=t)
            printf("\n");
    }
    return 0;
}