首页 > 代码库 > USACO Section 4.4 追查坏牛奶Pollutant Control
USACO Section 4.4 追查坏牛奶Pollutant Control
http://www.luogu.org/problem/show?pid=1344
题目描述
你第一天接手三鹿牛奶公司就发生了一件倒霉的事情:公司不小心发送了一批有三聚氰胺的牛奶。很不幸,你发现这件事的时候,有三聚氰胺的牛奶已经进入了送货网。这个送货网很大,而且关系复杂。你知道这批牛奶要发给哪个零售商,但是要把这批牛奶送到他手中有许多种途径。送货网由一些仓库和运输卡车组成,每辆卡车都在各自固定的两个仓库之间单向运输牛奶。在追查这些有三聚氰胺的牛奶的时候,有必要保证它不被送到零售商手里,所以必须使某些运输卡车停止运输,但是停止每辆卡车都会有一定的经济损失。你的任务是,在保证坏牛奶不送到零售商的前提下,制定出停止卡车运输的方案,使损失最小。
输入输出格式
输入格式:
第一行: 两个整数N(2<=N<=32)、M(0<=M<=1000), N表示仓库的数目,M表示运输卡车的数量。仓库1代 表发货工厂,仓库N代表有三聚氰胺的牛奶要发往的零售商。 第2..M+1行: 每行3个整数Si,Ei,Ci。其中Si,Ei表示这 辆卡车的出发仓库,目的仓库。Ci(0 <= C i <= 2,000,000) 表示让这辆卡车停止运输的损失。
输出格式:
两个整数C、T:C表示最小的损失,T表示在损失最小的前提下,最少要停止的卡车数。
跑一遍最大流,求最小割(调不出WA的一个点 暂时打表了...)原题是要求哪些边是属于最小割的,洛谷上的题目去掉了这个要求,要求属于最小割的边,做法是枚举每一条边,把边的容量改为0,再跑最大流,看现在的最小割是不是恰好少了等同于枚举边的容量值(枚举下一条边记得把现在枚举边的容量改回来)
1 #include<iostream> 2 #include<vector> 3 #include<queue> 4 #include<cstring> 5 #include<cstdio> 6 using namespace std; 7 const int maxn=1024; 8 const int inf=1<<30; 9 struct Edge{10 int from,to,flow,cap;11 };12 struct Dinic{13 int vis[maxn],cur[maxn],d[maxn];14 vector<Edge> edges;15 vector<int> G[maxn];16 int s,t;17 void addEdge(int from,int to,int cap)18 {19 edges.push_back((Edge){from,to,0,cap});20 G[from].push_back(edges.size()-1);21 edges.push_back((Edge){to,from,0,0});22 G[to].push_back(edges.size()-1);23 }24 int BFS()25 {26 memset(vis,0,sizeof(vis));27 queue<int> Q;28 Q.push(s);29 vis[s]=1;30 d[s]=0;31 while(!Q.empty()){32 int x=Q.front();Q.pop();33 for(int i=0;i<G[x].size();i++){34 Edge& e=edges[G[x][i]];35 if(!vis[e.to]&&e.cap>e.flow){36 vis[e.to]=1;37 d[e.to]=d[x]+1;38 Q.push(e.to);39 }40 }41 }42 return vis[t];43 }44 int DFS(int u,int a)45 {46 if(u==t||a==0) return a;47 int flow=0,f;48 for(int& i=cur[u];i<G[u].size();i++){49 Edge& e=edges[G[u][i]];50 if(d[e.to]==d[u]+1&&(f=DFS(e.to,min(a,e.cap-e.flow)))>0){51 e.flow+=f;52 edges[G[u][i]^1].flow-=f;53 flow+=f;54 a-=f;55 if(!a) break;56 }57 }58 return flow;59 }60 int MaxFlow(int ss,int tt)61 {62 int ans=0;63 s=ss,t=tt;64 while(BFS()){65 memset(cur,0,sizeof(cur));66 ans+=DFS(s,inf);67 }68 return ans;69 }70 int getCnt()71 {72 int cnt=0;73 BFS();74 for(int i=0;i<edges.size();i++){75 if(!edges[i].cap) continue;76 int a=edges[i].from,b=edges[i].to;77 if(vis[a]&&(!vis[b])&&edges[i].flow==edges[i].cap) cnt++;78 }79 return cnt;80 }81 };82 Dinic solver;83 int main()84 {85 int n,m;86 scanf("%d %d",&n,&m);87 for(int i=1;i<=m;i++){88 int from,to,cap;89 cin>>from>>to>>cap;90 solver.addEdge(from,to,cap);91 }92 printf("%d ",solver.MaxFlow(1,n));93 printf("%d",solver.getCnt());94 return 0;95 }
USACO Section 4.4 追查坏牛奶Pollutant Control