好难好难,将行列当成X和Y,源汇点连接各自的X,Y集,容量为行列的和,相当于从源点流向每一行,然后分配流量给每一列,最后流入汇点,这样执意要推断最后是否满流,就知道有没有解,而解就是每一行流向每一列多少流量。
关键在于怎么推断多解的情况。我想不到啊T_T
题讲解,找到一个长度大于2的环。
想了一想,也就是找到还有剩余流量的环,假设找到了,我就能够把当中一条边的流量转移,由于是一个环,所以它又会达到平衡,不会破坏最大流,可是这样转移后,解就多了一种,所以仅仅要推断是否有一个长度大于2的环就够了。
这里长度为什么要大于2,由于建图时的反向弧会导致 A->B而且B立马->A,这种话,也是一个环,但转移这条环的流量却会破坏最大流。
所以我们在用dfs找环的时候,要注意不能立马走反向弧。
dinic用了当前弧优化。
#include<cstdio> #include<cstring> #include<cmath> #include<iostream> #include<algorithm> #include<map> #include<queue> #include<vector> #include<string> #define eps 1e-12 #define INF 0x7fffffff #define maxn 1005 using namespace std; int n,m,k; int en; int st,ed; int dis[maxn],cur[maxn]; int que[999999]; struct edge { int to,c,next; }; edge e[999999]; int head[maxn]; void add(int a,int b,int c) { e[en].to=b; e[en].c=c; e[en].next=head[a]; head[a]=en++; e[en].to=a; e[en].c=0; e[en].next=head[b]; head[b]=en++; } int bfs() { memset(dis,-1,sizeof(dis)); dis[st]=0; int front=0,rear=0; que[rear++]=st; while(front<rear) { int j=que[front++]; for(int k=head[j];k!=-1;k=e[k].next) { int i=e[k].to; if(dis[i]==-1&&e[k].c) { dis[i] = dis[j]+ 1 ; que[rear++]=i; if(i==ed) return true; } } } return false; } int dfs(int x,int mx) { if(x==ed || mx==0) return mx; int f,flow=0; for(int& i=cur[x];i!=-1;i=e[i].next) { if(dis[x]+1==dis[e[i].to] && (f=dfs(e[i].to,min(mx,e[i].c)))) { e[i].c-=f; e[i^1].c+=f; flow+=f; mx-=f; if(!mx)break; } } return flow; } void init() { en=0; st=0; //源 ed=n+m+1; //汇 memset(head,-1,sizeof(head)); } inline int ReadInt() { char ch = getchar(); int data = 0; while (ch < '0' || ch > '9') { ch = getchar(); } do { data = data*10 + ch-'0'; ch = getchar(); }while (ch >= '0' && ch <= '9'); return data; } int s1,s2; void build() { int x,y,z; for(int i=1;i<=n;i++) { x=ReadInt(); s1+=x; add(st,i,x); } for(int i=1;i<=m;i++) { x=ReadInt(); s2+=x; add(i+n,ed,x); } for(int i=1;i<=n;i++) { for(int j=1;j<=m;j++) { add(i,j+n,k); } } } int dinic() { int tmp=0; int maxflow=0; while(bfs()) { for(int i=st;i<=ed;i++) cur[i]=head[i]; while(tmp=dfs(st,INF)) maxflow+=tmp; } return maxflow; } int TIME; struct node2 { int to,next; }e2[999999]; int en2,head2[maxn]; void add2(int a,int b) { e2[en2].to=b; e2[en2].next=head2[a]; head2[a]=en2++; } bool vis[maxn]; bool loop(int now,int fa) { for(int i=head2[now];~i;i=e2[i].next) { int to=e2[i].to; if(to!=fa) { if(vis[to]) return true; vis[to]=1; if(loop(to,now)) return true; vis[to]=0; } } return false; } bool cal() { en2=0; memset(head2,-1,sizeof(head2)); for(int i=1;i<=n+m;i++) { for(int j=head[i];~j;j=e[j].next) { if(e[j].c&&e[j].to>=1&&e[j].to<=n+m) { add2(i,e[j].to); } } } memset(vis,0,sizeof(vis)); for(int i=1;i<=n;i++) { if(loop(i,-1)) return true; } return false; } int ans[505]; int main() { int ca=1; while(scanf("%d%d%d",&n,&m,&k)!=EOF) { s1=s2=0; init(); build(); int maxflow=dinic(); if(maxflow!=s1||s1!=s2) puts("Impossible"); else { if(cal()) puts("Not Unique"); else { puts("Unique"); for(int i=1;i<=n;i++) { for(int j=head[i];~j;j=e[j].next) { if(e[j].to>n&&e[j].to<=n+m) { ans[e[j].to-n]=k-e[j].c; } } for(int j=1;j<=m;j++) { if(j==1) printf("%d",ans[j]); else printf(" %d",ans[j]); } puts(""); } } } } }