HIT2543 Stone IV（一定费用内的最大流）

题目大概说，有n个从0到n-1的城市，要从城市0运送石头到城市1，运送石头的单价是p。城市间的有m条双向路相连，路都有能运送石头的限额c1，如果超过限额运送石头的单价就要提高c2。问在总花费c以内能运送最多多少石头。

不大一样的题。建容量网络很容易：

源点向0连容量INF费用p的边，1向汇点连容量INF费用0的边
每条双向路拆成两条，一条容量c1费用0，另一条容量INF费用c2

而要求在一定费用内的最大流，一个直观的想法就是源点一单位一单位地流，直到费用不够减了，这样流了几单位就是这个费用下的最大流了——

于是可以这样解：在跑最小费用最大流的连续最短路过程中每次找到当前单位费用最少的可行流cost以及该可行流的最大流量flow，通过这个与剩余总价钱c比较来更新即可：

if((long long)cost*flow<=c){
    c-=cost*flow;
    res+=flow;
}else{
    res+=c/cost;
    break;
}

具体见代码：

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<queue>
 4 #include<algorithm>
 5 using namespace std;
 6 #define INF (1<<30)
 7 #define MAXN 1111
 8 #define MAXM 11111*16
 9 struct Edge{
10     int u,v,cap,cost,next;
11 }edge[MAXM];
12 int head[MAXN];
13 int NV,NE,vs,vt;
14 
15 void addEdge(int u,int v,int cap,int cost){
16     edge[NE].u=u; edge[NE].v=v; edge[NE].cap=cap; edge[NE].cost=cost;
17     edge[NE].next=head[u]; head[u]=NE++;
18     edge[NE].u=v; edge[NE].v=u; edge[NE].cap=0; edge[NE].cost=-cost;
19     edge[NE].next=head[v]; head[v]=NE++;
20 }
21 bool vis[MAXN];
22 int d[MAXN],pre[MAXN];
23 bool SPFA(){
24     for(int i=0;i<NV;++i){
25         vis[i]=0;
26         d[i]=INF;
27     }
28     vis[vs]=1;
29     d[vs]=0;
30     queue<int> que;
31     que.push(vs);
32     while(!que.empty()){
33         int u=que.front(); que.pop();
34         for(int i=head[u]; i!=-1; i=edge[i].next){
35             int v=edge[i].v;
36             if(edge[i].cap && d[v]>d[u]+edge[i].cost){
37                 d[v]=d[u]+edge[i].cost;
38                 pre[v]=i;
39                 if(!vis[v]){
40                     vis[v]=1;
41                     que.push(v);
42                 }
43             }
44         }
45         vis[u]=0;
46     }
47     return d[vt]!=INF;
48 }
49 int c;
50 int MCMF(){
51     int res=0;
52     while(SPFA()){
53         int flow=INF,cost=0;
54         for(int u=vt; u!=vs; u=edge[pre[u]].u){
55             flow=min(flow,edge[pre[u]].cap);
56         }
57         for(int u=vt; u!=vs; u=edge[pre[u]].u){
58             edge[pre[u]].cap-=flow;
59             edge[pre[u]^1].cap+=flow;
60             cost+=edge[pre[u]].cost;
61         }
62         if((long long)cost*flow<=c){
63             c-=cost*flow;
64             res+=flow;
65         }else{
66             res+=c/cost;
67             break;
68         }
69     }
70     return res;
71 }
72 int main(){
73     int t,n,m,p,u,v,c1,c2;
74     scanf("%d",&t);
75     while(t--){
76         scanf("%d%d%d%d",&n,&m,&c,&p);
77         vs=n; vt=vs+1; NV=vt+1; NE=0;
78         memset(head,-1,sizeof(head));
79         while(m--){
80             scanf("%d%d%d%d",&u,&v,&c1,&c2);
81             addEdge(u,v,c1,0);
82             addEdge(v,u,c1,0);
83             addEdge(u,v,INF,c2);
84             addEdge(v,u,INF,c2);
85         }
86         addEdge(vs,0,INF,p);
87         addEdge(1,vt,INF,0);
88         printf("%d
",MCMF());
89     }
90     return 0;
91 }