标签:des style http color os io for ar
5 6 5 3 5 2 3 4 12 1 5 5 4 2 3 2 4 4 3 2 1
3
题目大意:
解题思路:先是n,m表示n个城市,m条边,接下来src,sink,表示一个小偷从起点到终点,接下来n行表示每个城市放置警察需要的花费,接下来m行表示m条无向边。问你阻断小偷需要多少花费?
解题代码:将n个城市拆点构边,根据最小割定理,阻断起点到终点的连同最小割的花费也就是求最大流。
#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
#include <vector>
#include <cmath>
#include <algorithm>
using namespace std;
const int INF=(1<<30);
const int maxn=410,maxm=201000;
struct edge{
int u,v,f,next;
edge(int u0=0,int v0=0,int f0=0){
u=u0;v=v0;f=f0;
}
}e[maxm];
int src,sink,cnt,head[maxn];
void adde(int u,int v,int f){
e[cnt].u=u,e[cnt].v=v,e[cnt].f=f,e[cnt].next=head[u],head[u]=cnt++;
e[cnt].u=v,e[cnt].v=u,e[cnt].f=0,e[cnt].next=head[v],head[v]=cnt++;
}
void init(){
cnt=0;
memset(head,-1,sizeof(head));
}
queue <int> q;
bool visited[maxn];
int dist[maxn];
void bfs(){
memset(dist,0,sizeof(dist));
while(!q.empty()) q.pop();
visited[src]=true;
q.push(src);
while(!q.empty()){
int s=q.front();
q.pop();
for(int i=head[s];i!=-1;i=e[i].next){
int d=e[i].v;
if(e[i].f>0 && !visited[d]){
q.push(d);
dist[d]=dist[s]+1;
visited[d]=true;
}
}
}
}
int dfs(int u,int delta){
if(u==sink) return delta;
else{
int ret=0;
for(int i=head[u];delta && i!=-1;i=e[i].next){
if(e[i].f>0 && dist[e[i].v]==dist[u]+1){
int d=dfs(e[i].v,min(e[i].f,delta));
e[i].f-=d;
e[i^1].f+=d;
delta-=d;
ret+=d;
}
}
if(!ret) dist[u]=-2;
return ret;
}
}
int maxflow(){
int ret=0;
while(true){
memset(visited,false,sizeof(visited));
bfs();
if(!visited[sink]) return ret;
ret+=dfs(src,INF);
}
return ret;
}
int n,m;
void input(){
init();
scanf("%d%d",&src,&sink);
sink+=n;
for(int i=1;i<=n;i++){
int x;
scanf("%d",&x);
adde(i,i+n,x);
}
for(int i=0;i<m;i++){
int u,v;
scanf("%d%d",&u,&v);
adde(u+n,v,INF);
adde(v+n,u,INF);
}
}
void solve(){
printf("%d\n",maxflow());
}
int main(){
while(scanf("%d%d",&n,&m)!=EOF){
input();
solve();
}
return 0;
}
HDU 4289 Control (网络流-最小割),布布扣,bubuko.com
标签:des style http color os io for ar
原文地址:http://blog.csdn.net/a1061747415/article/details/38586835