本文主要是介绍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);elseprintf(" %d",edge[g[i][j]].flow+1);o++;}printf("\n");}if(cas!=t)printf("\n");}return 0;
}
这篇关于uva11082 Matrix Decompressing 最大流的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!