★ HDU 4888 网络流判断解是否唯一

题意：给你一个矩阵，然后给你矩阵每行每列的和，每个格子可以填0~k之间的任何数，让你复原这个矩阵。

分析：这题直接网络流建图不难，关键在于判断是否唯一解的地方。在求完最大流的残余网络中，如果有长度大于2的环，则说明有多个解。因为我们发现，在这个环里面走一遍的话，并不影响最大流。事实上我们可以在这个环里面改变某些边的流量，同时也不影响最大流，因此有多解。对于建图的话，S点向每一行建边，流量为该行的和。每一行向每一列建边，流量为k。每一列向T建边，流量为该列的和，然后跑一边最大流，判环即可，好题。

代码：

//Isap算法，复杂度O(n^2m)
#pragma comment(linker,"/STACK:102400000,102400000")
#include <iostream>
#include <string.h>
#include <stdio.h>
#include <algorithm>
#include <vector>
#include <string>
#include <math.h>
#include <queue>
#include <stack>
#include <map>
#include <set>
using namespace std;
typedef long long ll;   //记得必要的时候改成无符号
const int maxn=1005;
const int maxm=1000005;
const int INF=1000000000;
struct EdgeNode
{
    int from;
    int to;
    int cost;
    int cap;
    int next;
}edge[maxm];
int head[maxn],cnt;
void add(int x,int y,int z)
{
    edge[cnt].from=x;edge[cnt].to=y;edge[cnt].cap=z;edge[cnt].cost=z;edge[cnt].next=head[x];head[x]=cnt++;
    edge[cnt].from=y;edge[cnt].to=x;edge[cnt].cap=0;edge[cnt].cost=0;edge[cnt].next=head[y];head[y]=cnt++;
    //printf("%d    %d    %d\n",x,y,z);
}

void init()
{
    cnt=0;
    memset(head,-1,sizeof(head));
}

int S,T,n,m;
int d[maxn],gap[maxn],curedge[maxn],pre[maxn];
//curedge[]为当前弧数组，pre为前驱数组

int sap(int S,int T,int n)  //n为点数
{
    int cur_flow,flow_ans=0,u,tmp,neck,i;
    memset(d,0,sizeof(d));
    memset(gap,0,sizeof(gap));
    memset(pre,-1,sizeof(pre));
    for(i=0;i<=n;i++)curedge[i]=head[i]; //初始化当前弧为第一条邻接表
    gap[0]=n;
    u=S;
    while(d[S]<n)             //当d[S]>=n时，网络中肯定出现了断层
    {
        if(u==T)
        {
            cur_flow=INF;
            for(i=S;i!=T;i=edge[curedge[i]].to)
            {                           //增广成功，寻找瓶颈边
                if(cur_flow>edge[curedge[i]].cost)
                {
                    neck=i;
                    cur_flow=edge[curedge[i]].cost;
                }
            }
            for(i=S;i!=T;i=edge[curedge[i]].to)
            {                             //修改路径上的边容量
                tmp=curedge[i];
                edge[tmp].cost-=cur_flow;
                edge[tmp^1].cost+=cur_flow;
            }
            flow_ans+=cur_flow;
            u=neck;                     //下次增广从瓶颈边开始
        }
        for(i=curedge[u];i!=-1;i=edge[i].next)
            if(edge[i].cost&&d[u]==d[edge[i].to]+1)
               break;
        if(i!=-1)
        {
            curedge[u]=i;
            pre[edge[i].to]=u;
            u=edge[i].to;
        }
        else
        {
            if(0==--gap[d[u]])break;    //gap优化
            curedge[u]=head[u];
            for(tmp=n,i=head[u];i!=-1;i=edge[i].next)
                if(edge[i].cost)
                   tmp=min(tmp,d[edge[i].to]);
            d[u]=tmp+1;
            ++gap[d[u]];
            if(u!=S)u=pre[u];           //重标号并且从当前点前驱重新增广
        }
    }
    return flow_ans;
}

int a[maxn],b[maxn],chu[maxn][maxn],flag;
int deep[maxn];

void dfs(int u,int dp)
{
    if(deep[u]==-1)deep[u]=dp;
    for(int i=head[u];i!=-1;i=edge[i].next)
    {
        int v=edge[i].to;
        if(edge[i].cost){
            if(deep[v]==-1)dfs(v,dp+1);
            else if(dp+1-deep[v]>2)flag=1;
        }
    }
}

void solve()
{
    for(int i=1;i<=n;i++){
        memset(deep,-1,sizeof(deep));
        dfs(i,0);
        if(flag==1)break;
    }
}

int main()
{
    int k,i,j,suma,sumb,x,y;
    while(~scanf("%d%d%d",&n,&m,&k))
    {
        init();
        suma=sumb=0;
        for(i=1;i<=n;i++){
            scanf("%d",&a[i]);
            suma+=a[i];
        }
        for(i=1;i<=m;i++){
            scanf("%d",&b[i]);
            sumb+=b[i];
        }
        if(suma!=sumb){
            printf("Impossible\n");
            continue;
        }
        S=0; T=n+m+1;
        for(i=1;i<=n;i++){
            add(S,i,a[i]);
            for(j=1;j<=m;j++)add(i,n+j,k);
        }
        for(i=1;i<=m;i++)add(n+i,T,b[i]);
        if(suma!=sap(S,T,T+1))printf("Impossible\n");
        else{
            flag=0;
            solve();
            if(flag==1){
                printf("Not Unique\n");
                continue;
            }
            printf("Unique\n");
            for(i=0;i<cnt;i++){
                if(edge[i].cap){
                    x=edge[i].from; y=edge[i].to;
                    if(x>=1&&x<=n&&y>n&&y<=n+m)chu[x][y-n]=edge[i].cap-edge[i].cost;
                }
            }
            for(i=1;i<=n;i++){
                for(j=1;j<=m;j++){
                    if(j==1)printf("%d",chu[i][j]);
                    else printf(" %d",chu[i][j]);
                }
                printf("\n");
            }
        }
    }
    return 0;
}