UVA-11082 Matrix Decompressing(有上下界的最大流)

题目链接：

Matrix Decompressing

题意：

给一个矩阵的每行和每列的和,(给的是前i行或者列的和);

矩阵中每个元素的值在1到20之间,找出这样的一个矩阵；

思路：

把它转化成一个二分图,每行和每列之间连一条弧，然后设置一个源点和一个汇点，源点与行点相连，汇点与列点相连，求一个最大值,当然这是一个有下界的最大流,需要做-1的处理，同时与源汇相连的边也是要处理的;最后求得的反向边+1就是答案了;

AC代码：

#include <bits/stdc++.h>
/*#include <iostream>
#include <queue>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <cstdio>
*/
using namespace std;
#define Riep(n) for(int i=1;i<=n;i++)
#define Riop(n) for(int i=0;i<n;i++)
#define Rjep(n) for(int j=1;j<=n;j++)
#define Rjop(n) for(int j=0;j<n;j++)
#define mst(ss,b) memset(ss,b,sizeof(ss));
typedef long long LL;
const LL mod=1e9+;
const double PI=acos(-1.0);
const int inf=0x3f3f3f3f;
const int N=1e5+;
 
int n,m,a[],b[],cap[][],path[],flow[];
queue<int>qu;
int bfs()
{
    mst(path,-);
    flow[]=inf;
    path[]=;
    qu.push();
    while(!qu.empty())
    {
        int fr=qu.front();
        qu.pop();
        for(int i=;i<=n+m+;i++)
        {
            if(path[i]==-&&cap[fr][i])
            {
                flow[i]=min(cap[fr][i],flow[fr]);
                path[i]=fr;
                qu.push(i);
            }
        }
    }
    if(path[n+m+]==-)return -;
    return flow[n+m+];
}
void maxflow()
{
    int now,pre;
    while()
    {
        int temp=bfs();
        if(temp==-)break;
        now=n+m+;
        while(now!=)
        {
            pre=path[now];
            cap[pre][now]-=temp;
            cap[now][pre]+=temp;
            now=pre;
        }
    }
}
int main()
{
    int t,Case=;
    scanf("%d",&t);
    while(Case<=t)
    {
        mst(cap,);
        scanf("%d%d",&n,&m);
        Riep(n)scanf("%d",&a[i]);
        Riep(n)
        {
            cap[][i]=a[i]-a[i-]-m;
        }
        Riep(m)scanf("%d",&b[i]);
        Riep(m)
        {
           cap[i+n][n+m+]=b[i]-b[i-]-n;
        }
        for(int i=;i<=n;i++)
            for(int j=;j<=m;j++)
                cap[i][j+n]=;
        maxflow();
        printf("Matrix %d\n",Case++);
        Riep(n)
        {
            Rjep(m-)
                printf("%d ",cap[j+n][i]+);
            printf("%d\n",cap[m+n][i]+);
        }
        printf("\n");
    }
    return ;
}