HDUOJ--4888--Redraw Beautiful Drawings【isap】网络流+判环

标签：网络最大流

链接：http://acm.hdu.edu.cn/showproblem.php?pid=4888

题意：一个矩阵，限定每行行和、列和，每个格子数字不超过k，问矩阵是否存在，如存在判断有单解还是多解。

思路：之前多校的题目，那时候还不会网络流，现在A掉了，矩阵的建图模型，判断网络流是否可行只要判断最大流是否等于总行和或总列和即可，判环是看的别人的解题报告，方法是使用dfs查找残余网络中是否有还存在容量的弧形成了环，如果有，说明可以通过这个环改变容量网络内部的增广路方式，而源汇的流量是不会变的，就说明存在多解。如果没有环，就是单一解。

建图：源点向每个行节点连弧，容量为该行行和，每个列节点向汇点连边，容量为每个列和，每个行节点与每个列节点之间连边，容量为k。

细节：我用的我之前AC过题目的isap模板，但是WA了，我debug之后发现错误不会在建图和判环里，然后就和别人的isap代码对照，发现了一个不一样的地方，改之AC，不懂为何，不一样的地方代码中有注释。我的理解nn个点，源点src的层次是0，则最大的层次应该就是nn-1，但是minm初始化为nn-1是WA，初始化为nn就AC了。理论上应该是初始化大于等于nn-1都是可以的，以后还是就初始化为nn吧避免莫名其妙的WA。

#include<cstring>
#include<string>
#include<fstream>
#include<iostream>
#include<iomanip>
#include<cstdio>
#include<cctype>
#include<algorithm>
#include<queue>
#include<map>
#include<set>
#include<vector>
#include<stack>
#include<ctime>
#include<cstdlib>
#include<functional>
#include<cmath>
using namespace std;
#define PI acos(-1.0)
#define MAXN 50100
#define eps 1e-7
#define INF 0x7FFFFFFF
#define LLINF 0x7FFFFFFFFFFFFFFF
#define seed 131
#define MOD 1000000007
#define ll long long
#define ull unsigned ll
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1

struct node{
    int u,v,w,next;
}edge[500000];
int head[820],dist[820],cur[820],fa[820],num[820],vis[820];
int n,m,k,cnt,nn,src,sink;
void add_edge(int a,int b,int c){
    edge[cnt].u = a;
    edge[cnt].v = b;
    edge[cnt].w = c;
    edge[cnt].next = head[a];
    head[a] = cnt++;
}
void bfs()
{
    int x,i,j;
    queue<int> q;
    memset(dist,-1,sizeof(dist));
    memset(num,0,sizeof(num));
    q.push(sink);
    dist[sink] = 0;
    num[0] = 1;
    while(!q.empty()){
        x = q.front();
        q.pop();
        for(i=head[x];i!=-1;i=edge[i].next){
            if(dist[edge[i].v]<0){
                dist[edge[i].v] = dist[x] + 1;
                num[dist[edge[i].v]]++;
                q.push(edge[i].v);
            }
        }
    }
}

int augment()
{
    int x=sink,a=INF;
    while(x!=src){
        a = min(a,edge[fa[x]].w);
        x = edge[fa[x]].u;
    }
    x=sink;
    while(x!=src){
        edge[fa[x]].w -= a;
        edge[fa[x]^1].w += a;
        x = edge[fa[x]].u;
    }
    return a;
}

int isap()
{
    int i,x,ok,minm,flow=0;
    bfs();
    for(i=0;i<=nn+5;i++) cur[i] = head[i];
    x=src;
    while(dist[src]<nn){
        if(x==sink){
            flow += augment();
            x = src;
        }
        ok=0;
        for(i=cur[x];i!=-1;i=edge[i].next){
            if(edge[i].w && dist[x]==dist[edge[i].v]+1){
                ok=1;
                fa[edge[i].v] = i;
                cur[x] = i;
                x = edge[i].v;
                break;
            }
        }
        if(!ok){
            minm = nn;	//minm = nn - 1 就WA
            for(i=head[x];i!=-1;i=edge[i].next)
                if(edge[i].w && dist[edge[i].v]<minm)   minm=dist[edge[i].v];
            if(--num[dist[x]]==0)break;
            num[dist[x]=minm+1]++;
            cur[x]=head[x];
            if(x!=src)  x=edge[fa[x]].u;
        }
    }
    return flow;
}
bool dfs(int u,int pre){
    int i,j;
    if(vis[u])  return true;
    vis[u] = 1;
    for(i=head[u];i!=-1;i=edge[i].next){
        if(edge[i].w>0&&edge[i].v!=pre&&dfs(edge[i].v,u))
            return true;
    }
    vis[u] = 0;
    return false;
}
int row[410],col[410];
int ans[420][420];
int main(){
    int i,j;
    int sumr,sumc;
    while(scanf("%d%d%d",&n,&m,&k)!=EOF){
        memset(head,-1,sizeof(head));
        sumr = sumc = 0;
        cnt = 0;
        src = 0;
        sink = n + m + 1;
        nn = sink + 1;
        for(i=1;i<=n;i++){
            scanf("%d",&row[i]);
            sumr += row[i];
            add_edge(src,i,row[i]);
            add_edge(i,src,0);
            for(j=1;j<=m;j++){
                add_edge(i,j+n,k);
                add_edge(j+n,i,0);
            }
        }
        for(i=1,j=n+1;i<=m;j++,i++){
            scanf("%d",&col[i]);
            sumc += col[i];
            add_edge(j,sink,col[i]);
            add_edge(sink,j,0);
        }
        if(sumr!=sumc){
            puts("Impossible");
            continue;
        }
        int flag = 0;
        int flow = isap();
        if(flow!=sumr){
            puts("Impossible");
            continue;
        }
        memset(vis,0,sizeof(vis));
        for(i=1;i<=n;i++){
            if(dfs(i,-1)){
                flag = 1;
                break;
            }
        }
        if(flag){
            puts("Not Unique");
            continue;
        }
        puts("Unique");
        memset(ans,0,sizeof(ans));
        for(i=1;i<=n;i++){
            for(j=head[i];j!=-1;j=edge[j].next){
                int u = edge[j].v;
                if(u>n&&u<=n+m){
                    ans[i][u-n] = k - edge[j].w;
                }
            }
        }
        for(i=1;i<=n;i++){
            for(j=1;j<=m;j++){
                if(j>1) printf(" ");
                printf("%d",ans[i][j]);
            }
            printf("\n");
        }
    }
    return 0;
}

HDUOJ--4888--Redraw Beautiful Drawings【isap】网络流+判环

标签：网络最大流

原文地址：http://blog.csdn.net/zzzz40/article/details/38765293