首页 > 代码库 > HDU 4888

HDU 4888

2024-07-15 13:38:00   222人阅读

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

这题时限给的太短了，以至于卡了一些奇怪的东西，我开始先输入再建图一直T，一边输入一边连边就AC了

题意：n*m的矩阵，每个点可以放0-k的数，给你每行每列的和，问可行否，可行的话解是否唯一，唯一则输出那组解

网络流行列建边，容量为k，这样一条边代表一个点，起点和每行连容量为此行和的边，列和终点连容量为每列和的边，跑最大流，看是否满流

如果能在残量网络找到环，则答案不唯一（按环走一遍不改变结果），注意找环的时候不能马上找其反向边，所以dfs要加一个参数记录当前节点的父亲

输出答案的时候只要看行列每边的流量，流量=k-剩余流量

#include <iostream>#include <cstdio>#include <cstring>using namespace std ;const int INF=0xfffffff ;struct node{    int s,t,cap,nxt ;}e[400005] ;int m,n,cnt,head[100005],level[100005],q[100005] ;void add(int s,int t,int cap){    e[cnt].s=s ;e[cnt].t=t ;e[cnt].cap=cap ;e[cnt].nxt=head[s] ;head[s]=cnt++ ;    e[cnt].s=t ;e[cnt].t=s ;e[cnt].cap=0 ;e[cnt].nxt=head[t] ;head[t]=cnt++ ;}bool build(int s,int t){    int front=0,rear=0 ;    memset(level,-1,sizeof(level)) ;    q[rear++]=s ;    level[s]=1 ;    while(front<rear)    {        int u=q[front++] ;        for(int i=head[u] ;i!=-1 ;i=e[i].nxt)        {            int tt=e[i].t ;            if(level[tt]==-1 && e[i].cap>0)            {                level[tt]=level[u]+1 ;                if(tt==t)return true ;                q[rear++]=tt ;            }        }    }    return false ;}int find(int s,int t,int flow){    if(s==t)return flow ;    int ret=0,a ;    for(int i=head[s] ;i!=-1 ;i=e[i].nxt)    {        int tt=e[i].t ;        if(level[tt]==level[s]+1 && e[i].cap>0)        {            a=find(tt,t,min(e[i].cap,flow-ret)) ;            e[i].cap-=a ;            e[i^1].cap+=a ;            ret+=a ;            if(ret==flow)                return ret ;        }    }    if(!ret)level[s]=-1 ;    return ret ;}int dinic(int s,int t){    int flow,ret=0 ;    while(build(s,t))        while(flow=find(s,t,INF))            ret+=flow ;    return ret ;}int S,T ;int vis[100005] ;int JG(int u,int fa){    if(vis[u])return 1 ;    vis[u]=1 ;    for(int i=head[u] ;i!=-1 ;i=e[i].nxt)    {        int v=e[i].t ;        if(v!=fa && v!=S && v!=T && e[i].cap)            if(JG(v,u))return 1 ;    }    vis[u]=0 ;    return 0 ;}int A[505],B[505] ;int res[505][505] ;int main(){    int k ;    while(~scanf("%d%d%d",&n,&m,&k))    {        cnt=0 ;        memset(head,-1,sizeof(head)) ;        int sum1=0,sum2=0 ;        S=0 ;T=n+m+1 ;        for(int i=1 ;i<=n ;i++)        {            scanf("%d",&A[i]) ;            sum1+=A[i] ;            add(S,i,A[i]) ;            for(int j=1 ;j<=m ;j++)                add(i,j+n,k) ;        }        for(int i=1 ;i<=m ;i++)        {            scanf("%d",&B[i]) ;            sum2+=B[i] ;            add(i+n,T,B[i]) ;        }        int ans=dinic(S,T) ;        if(ans==sum1 && ans==sum2)        {            int flag=0 ;            memset(vis,0,sizeof(vis)) ;            for(int i=1 ;i<=n ;i++)            {                if(JG(i,-1))                {                    flag=1 ;                    break ;                }            }            if(flag)            {                puts("Not Unique") ;            }            else            {                puts("Unique") ;                memset(res,0,sizeof(res)) ;                for(int i=1 ;i<=n ;i++)                {                    for(int j=head[i] ;j!=-1 ;j=e[j].nxt)                    {                        if(e[j].t<=n+m && e[j].t>n)                            res[i][e[j].t-n]=k-e[j].cap ;                    }                }                for(int i=1 ;i<=n ;i++)                {                    for(int j=1 ;j<=m ;j++)                     {                         if(j>1)printf(" ") ;                         printf("%d",res[i][j]) ;                    }                    putchar(‘\n‘) ;                }            }        }        else        {            puts("Impossible") ;        }    }    return 0 ;}
View Code

 


联系
我们