好難好難，將行列當成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("");
}
}
}
}
}