标签:
最大流+拆点
#include<cstdio> #include<cstring> #include<string> #include<cmath> #include<vector> #include<queue> #include<algorithm> using namespace std; const int maxn=1000+10; const int INF=0x7FFFFFFF; struct Edge { int from,to,cap,flow; }; vector<Edge>edges; vector<int>G[maxn]; bool vis[maxn]; int d[maxn]; int cur[maxn]; int n,m,s,t; //求出层次网络 bool BFS() { memset(vis,0,sizeof(vis)); queue<int>Q; Q.push(s); d[s]=0; vis[s]=1; while(!Q.empty()) { int x=Q.front(); Q.pop(); for(int i=0; i<G[x].size(); i++) { Edge& e=edges[G[x][i]]; if(!vis[e.to]&&e.cap>e.flow) { vis[e.to]=1; d[e.to]=d[x]+1; Q.push(e.to); } } } return vis[t]; } //加边 void AddEdge(int from,int to,int cap) { Edge r; r.from=from; r.to=to; r.cap=cap; r.flow=0; edges.push_back(r); Edge d; d.from=to; d.to=from; d.cap=0; d.flow=0; edges.push_back(d); m=edges.size(); G[from].push_back(m-2); G[to].push_back(m-1); } //每个阶段来一次DFS增广 int DFS(int x,int a) { if(x==t||a==0) return a; int flow=0,f; for(int i=cur[x]; i<G[x].size(); i++) { Edge& e=edges[G[x][i]]; if(d[x]+1==d[e.to]&&(f=DFS(e.to,min(a,e.cap-e.flow)))>0) { e.flow+=f; edges[G[x][i]^1].flow-=f; flow+=f; a-=f; if(a==0) break; } } return flow; } //多个阶段,多次建立层次网络。 int Maxflow(int ss,int tt) { int flow=0; while(BFS()) { memset(cur,0,sizeof(cur)); flow+=DFS(ss,INF); } return flow; } int N,F,D; char S[1000]; int main() { while(~scanf("%d%d%d",&N,&F,&D)) { edges.clear(); for(int i=0; i<maxn; i++) G[i].clear(); s=0; t=F+N+N+D+1; for(int i=1; i<=F; i++) { int x; scanf("%d",&x); AddEdge(s,i,x); } for(int i=1; i<=D; i++) { int x; scanf("%d",&x); AddEdge(F+N+N+i,t,x); } for(int i=1; i<=N; i++) { scanf("%s",S); for(int j=0; S[j]; j++) if(S[j]==‘Y‘) AddEdge(j+1,F+i,1); } for(int i=1;i<=N;i++) AddEdge(F+i,F+N+i,1); for(int i=1; i<=N; i++) { scanf("%s",S); for(int j=0; S[j]; j++) if(S[j]==‘Y‘) AddEdge(F+N+i,F+N+N+j+1,1); } printf("%d\n",Maxflow(s,t)); } return 0; }
标签:
原文地址:http://www.cnblogs.com/zufezzt/p/4741483.html