hdu 4888 Redraw Beautiful Drawings 最大流

好难好难,将行列当成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("");
                }
            }
        }
    }
}


hdu 4888 Redraw Beautiful Drawings 最大流

上一篇:HDUOJ--4888--Redraw Beautiful Drawings【isap】网络流+判环


下一篇:Elasticsearch升级到1.2(1.x)API的变化-二