hdu4888 最大流，dinic模板

源点向行节点建边，容量为行的和，列节点向汇点建边，容量为列的和，中间行节点向列节点建边，容量为k，若最大流等于总和，则存在。

关键在于如何判断唯一性，在残流网络中，若存在>2的环（u->v和反向边v->u构成的环不行），就不唯一。

#include<iostream>
#include<string>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<iomanip>
#include<map>
#include<algorithm>
#include<queue>
#include<set>
#include<stack>
#define inf 1000000000
#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;
struct Edge
{
    int to,c,next;
}e[200000*2];
int st,ed,cnt;
int head[1000];
void add(int a,int b,int c)
{
    e[cnt].to=b;
    e[cnt].c=c;
    e[cnt].next=head[a];
    head[a]=cnt++;
    e[cnt].to=a;
    e[cnt].c=0;
    e[cnt].next=head[b];
    head[b]=cnt++;
}
int vis[1000],dis[1000],que[10000];
int cur[1000];
bool 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;
}
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;
}
//以上为dinic模板

int n,m;
int d1[505],d2[505];

bool dfs(int u,int fa)
{
    for(int i=head[u];i!=-1;i=e[i].next)
    {
        Edge& ee=e[i];
        if(ee.to==fa)
            continue;
        if(ee.c)
        {
            if(vis[ee.to])
                return true;
            vis[ee.to]=1;
            if(dfs(ee.to,u))
                return true;
            vis[ee.to]=0;
        }
    }
    return false;
}
int mp[505][505];
int main()
{
    int k;
    int tot1,tot2;
    while(~scanf("%d%d%d",&n,&m,&k))
    {
        tot1=tot2=0;
        memset(head,-1,sizeof(head));
        cnt=0;
        st=0,ed=n+m+1;
        for(int i=1;i<=n;i++)
        {
            scanf("%d",&d1[i]);
            add(st,i,d1[i]);
            tot1+=d1[i];
            for(int j=m;j>=1;j--)
                add(i,n+j,k);
        }
        for(int i=1;i<=m;i++)
        {
            scanf("%d",&d2[i]);
            tot2+=d2[i];
            add(n+i,ed,d2[i]);
        }
        if(tot1!=tot2)
        {
            printf("Impossible\n");
            continue;
        }
        int tot=dinic();
        if(tot==tot1&&tot==tot2)
        {
            bool flag=false;
            memset(vis,0,sizeof(vis));
            for(int i=1;i<=n;i++)
            {
                if(vis[i]==0)
                if(dfs(i,-1))
                {
                    flag=true;
                    break;
                }
            }
            if(flag)
                printf("Not Unique\n");
            else
            {
                printf("Unique\n");
                for(int i=1;i<=n;i++)
                {
                    for(int j=head[i];j!=-1;j=e[j].next)
                    {
                        if(e[j].to!=st&&e[j].to!=ed)
                            mp[i][e[j].to-n]=k-e[j].c;
                    }
                }
                for(int i=1;i<=n;i++)
                {
                    printf("%d",mp[i][1]);
                    for(int j=2;j<=m;j++)
                    {
                        printf(" %d",mp[i][j]);
                    }
                    printf("\n");
                }
            }
        }
        else
            printf("Impossible\n");
    }
    return 0;
}