标签:
2 3 3 1 2 2 2 2 3 2 2 3 1 2 2 3 3 1 2 10 2 2 3 2 2 3 1 2 2
Case #1: happy Case #2: unhappy
1
4 5
1 4 10 2
4 1 6 2
2 1 1 1
4 3 6 2
3 2 4 2
时,输出的是happy。可是实际上,选择1、4节点,Y=2,X=8,显然是unhappy。
附上非正解代码:(AC了只能说明题目数据太水了)
#include <cstdio> #include <cstring> #define N 205 #define sf(x) scanf("%d",&x) int x[N],y[N]; int main(){ int t; sf(t); for(int cas=1;cas<=t;cas++){ memset(x,0,sizeof x); memset(y,0,sizeof y); printf("Case #%d: ",cas); int n,m; sf(n); sf(m); for(int i=1;i<=m;i++){ int u,v,d,b; sf(u);sf(v);sf(d);sf(b); x[u]+=d; y[v]+=d+b; } int ok=1; for(int i=1;i<=n;i++) if(x[i]>y[i])ok=0; if(ok)puts("happy"); else puts("unhappy"); } }
正解是无源无汇带上下界判断是否有可行流。
求s到t的最大流,如果源点汇点连接的边全部满流则有可行解。
参考国家集训队论文《一种简易的方法求解流量有上下界的网络中网络流问题》
#include <iostream> #include <algorithm> #include <cstdio> #include <cstring> #include <cmath> #include <queue> #include <vector> #include <map> #define sf(a) scanf("%d",&a) using namespace std; #define N 500 #define M 100001 #define inf 0x3f3f3f3f struct edge{ int to,next,cap,flow; }e[M]; int head[N],cnt; int gap[N],dep[N],cur[N]; void init(){ cnt=0; memset(head, -1, sizeof head); } void add(int u,int v,int w,int rw=0){ e[cnt]=(edge){v,head[u],w,0}; head[u]=cnt++; e[cnt]=(edge){u,head[v],rw,0}; head[v]=cnt++; } int q[N]; void bfs(int st,int ed){ memset(dep,-1,sizeof dep); memset(gap,0,sizeof gap); gap[0]=1; int front=0,rear=0; dep[ed]=0; q[rear++]=ed; while(front!=rear){ int u=q[front++]; for(int i=head[u];~i;i=e[i].next){ int v=e[i].to; if(dep[v]!=-1)continue; q[rear++]=v; dep[v]=dep[u]+1; gap[dep[v]]++; } } } int s[N]; int sap(int st,int ed,int n){ bfs(st,ed); memcpy(cur,head,sizeof head); int top=0; int u=st; int ans=0; while(dep[st]<n){ if(u==ed){ int Min=inf; int inser; for(int i=0;i<top;i++) if(Min>e[s[i]].cap-e[s[i]].flow){ Min=e[s[i]].cap-e[s[i]].flow; inser=i; } for(int i=0;i<top;i++){ e[s[i]].flow+=Min; e[s[i]^1].flow-=Min; } ans+=Min; top=inser; u=e[s[top]^1].to; continue; } bool flag=false; int v; for(int i=cur[u];~i;i=e[i].next){ v=e[i].to; if(e[i].cap-e[i].flow&&dep[v]+1==dep[u]){ flag=true; cur[u]=i; break; } } if(flag){ s[top++]=cur[u]; u=v; continue; } int Min=n; for(int i=head[u];~i;i=e[i].next) if(e[i].cap-e[i].flow &&dep[e[i].to]<Min){ Min=dep[e[i].to]; cur[u]=i; } gap[dep[u]]--; if(!gap[dep[u]])return ans; gap[dep[u]=Min+1]++; if(u!=st)u=e[s[--top]^1].to; } return ans; } int main(){ int t; sf(t); for(int cas=1;cas<=t;cas++){ printf("Case #%d: ",cas); int n,m; sf(n); sf(m); int s=0,t=n+1,in[N]; memset(in,0,sizeof in); init(); for(int i=1;i<=m;i++){ int u,v,d,b; sf(u);sf(v);sf(d);sf(b); add(u,v,b); in[v]+=d; in[u]-=d; } int need=0; for(int i=1;i<=n;i++){ if(in[i]>0){ add(s,i,in[i]); need+=in[i]; } else add(i,t,-in[i]); } int ans=sap(s, t, t+1); if(need==ans)puts("happy"); else puts("unhappy"); } }
【HDU 4940】Destroy Transportation system(数据水/无源无汇带上下界可行流)
标签:
原文地址:http://www.cnblogs.com/flipped/p/5792694.html