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 }

View Code