hdu 4289 Control
题意就是 要抓捕恐怖分子 在恐怖分子行驶的路上安排SA 但是安排SA要花钱的啊。。
现在就是在能抓到恐怖的前提下花最少的钱
路是双向的 每个点有各自的消耗值 就是求最小的话费使得恐怖分子的起点S与终点T之间不联通
所以就是一个最大流的应用。。。
思路就是拆点网络流 把每个点拆作两个点i和i' 连接i->i'的费用为这个点的花费 如果原图中一个连接a->b
则连接a->b' 求这个图的最大流即可!!!!
贴一个模版代码。。。
#include <cstdio>
#include <iostream>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <string.h>
#include <string>
#include <vector>
#include <queue>
#define eps 1e-8
#define op operator
#define MOD 10009
#define MAXN 10010
#define INF 0x7fffffff
#define FOR(i,a,b) for(int i=a;i<=b;i++)
#define FOV(i,a,b) for(int i=a;i>=b;i--)
#define REP(i,a,b) for(int i=a;i<b;i++)
#define REV(i,a,b) for(int i=a-1;i>=b;i--)
#define MEM(a,x) memset(a,x,sizeof a)
#define ll __int64
using namespace std;
struct Dinic
{
struct Edge
{
int from,to,cap,flow;
Edge() {}
Edge(int from,int to,int cap,int flow):from(from),to(to),cap(cap),flow(flow) {}
};
vector<Edge> edges;
vector<int> G[MAXN];
bool vis[MAXN];//BFS使用
int d[MAXN]; //从起点到i的距离
int cur[MAXN]; //当前弧下标
int n,m,s,t,maxflow; //节点数 边数(包括反向弧) 源点编号和弧点编号
void init(int n)
{
this->n=n;
for(int i=1;i<=n;i++)
G[i].clear();
edges.clear();
}
void addedge(int from,int to,int cap)
{
edges.push_back(Edge(from,to,cap,0));
edges.push_back(Edge(to,from,0,0));
m=edges.size();
G[from].push_back(m-2);
G[to].push_back(m-1);
}
bool bfs()
{
MEM(vis,0);
MEM(d,-1);
queue<int> q;
q.push(s);
d[s]=maxflow=0;
vis[s]=1;
while(!q.empty())
{
int u=q.front(); q.pop();
int sz=G[u].size();
for(int i=0;i<sz;i++)
{
Edge e=edges[G[u][i]];
if(!vis[e.to]&&e.cap>e.flow)
{
d[e.to]=d[u]+1;
vis[e.to]=1;
q.push(e.to);
}
}
}
return vis[t];
}
int dfs(int u,int a)
{
if(u==t||a==0) return a;
int sz=G[u].size();
int flow=0,f;
for(int &i=cur[u];i<sz;i++)
{
Edge &e=edges[G[u][i]];
if(d[u]+1==d[e.to]&&(f=dfs(e.to,min(a,e.cap-e.flow)))>0)
{
e.flow+=f;
edges[G[u][i]^1].flow-=f;
flow+=f;
a-=f;
if(a==0) break;
}
}
return flow;
}
int Maxflow(int s,int t)
{
this->s=s; this->t=t;
int flow=0;
while(bfs())
{
MEM(cur,0);
flow+=dfs(s,INF);
}
return flow;
}
}Dic;
int main()
{
//freopen("ceshi.txt","r",stdin);
int n,m,s,t;
while(scanf("%d%d%d%d",&n,&m,&s,&t)!=EOF)
{
Dic.init(n*2);
for(int i=1;i<=n;i++)
{
int d;
scanf("%d",&d);
Dic.addedge(i,i+n,d);
}
for(int i=1;i<=m;i++)
{
int u,v;
scanf("%d%d",&u,&v);
Dic.addedge(u+n,v,INF);
Dic.addedge(v+n,u,INF);
}
printf("%d\n",Dic.Maxflow(s,t+n));
}
return 0;
}