hdu 2485(最小费用最大流)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2485
思路:题目的意思是删除最少的点使1,n的最短路大于k。将点转化为边,容量为1,费用为0,然后就是对于那些有道路的城市之间连边,若(u,v)有边,则连边(u+n)->v,容量为inf,费用为花费的时间1,然后就是跑最小费了,若dist[vt]>k,则返回false,最后输出的flow就代表要删除的点的个数。
1 #include<iostream> 2 #include<cstdio> 3 #include<cstring> 4 #include<algorithm> 5 #include<queue> 6 using namespace std; 7 #define MAXN 222 8 #define MAXM 4444444 9 #define inf 1<<30 10 11 struct Edge{ 12 int v,cap,cost,next; 13 }edge[MAXM]; 14 15 int n,m,k,vs,vt,NE; 16 int head[MAXN]; 17 18 void Insert(int u,int v,int cap,int cost) 19 { 20 edge[NE].v=v; 21 edge[NE].cap=cap; 22 edge[NE].cost=cost; 23 edge[NE].next=head[u]; 24 head[u]=NE++; 25 26 edge[NE].v=u; 27 edge[NE].cap=0; 28 edge[NE].cost=-cost; 29 edge[NE].next=head[v]; 30 head[v]=NE++; 31 } 32 33 int dist[MAXN],pre[MAXN],cur[MAXN]; 34 bool mark[MAXN]; 35 bool spfa(int vs,int vt) 36 { 37 memset(mark,false,sizeof(mark)); 38 fill(dist,dist+2*n+1,inf); 39 dist[vs]=0; 40 queue<int>que; 41 que.push(vs); 42 while(!que.empty()){ 43 int u=que.front(); 44 que.pop(); 45 mark[u]=false; 46 for(int i=head[u];i!=-1;i=edge[i].next){ 47 int v=edge[i].v,cost=edge[i].cost; 48 if(edge[i].cap>0&&dist[u]+cost<dist[v]){ 49 dist[v]=dist[u]+cost; 50 pre[v]=u; 51 cur[v]=i; 52 if(!mark[v]){ 53 mark[v]=true; 54 que.push(v); 55 } 56 } 57 } 58 } 59 if(dist[vt]>k)return false; 60 return dist[vt]!=inf; 61 } 62 63 int MinCostFlow(int vs,int vt) 64 { 65 int flow=0,cost=0; 66 while(spfa(vs,vt)){ 67 int aug=inf; 68 for(int u=vt;u!=vs;u=pre[u]){ 69 aug=min(aug,edge[cur[u]].cap); 70 } 71 flow+=aug,cost+=dist[vt]*aug; 72 for(int u=vt;u!=vs;u=pre[u]){ 73 edge[cur[u]].cap-=aug; 74 edge[cur[u]^1].cap+=aug; 75 } 76 } 77 return flow; 78 } 79 80 int main() 81 { 82 int u,v; 83 while(~scanf("%d%d%d",&n,&m,&k)){ 84 if(n==0&&m==0&&k==0)break; 85 NE=0; 86 vs=1,vt=n; 87 memset(head,-1,sizeof(head)); 88 for(int i=2;i<=n-1;i++)Insert(i,i+n,1,0); 89 while(m--){ 90 scanf("%d%d",&u,&v); 91 if(u==vs)Insert(vs,v,inf,1); 92 else Insert(u+n,v,inf,1); 93 } 94 printf("%d\n",MinCostFlow(vs,vt)); 95 } 96 return 0; 97 }