HDU 3618 Good Plan(费用流)

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=3618

Problem Description

FJ has two same houses for rant. Now he has n (1 ≤ n ≤ 1000) piece of order, the orders are given in the form: 

s t v

means that someone want to rant a house from the day s to t paying v yuan totally (including the day s and t, 0 ≤ s ≤ t ≤ 400, 0 ≤ v ≤ 100,0000). 

A house can be only rant to one person, and FJ should either accept an order totally or reject it. 

 

Input

The first line of input file is a single integer T - The number of test cases. For each test case, the first line is a single integer n then there n lines, each line gives an order. 

 

Output

For each data set, print a single line containing an integer, the maximum total income for the data set. 

 

此题用来测试了一下SPFA费用流的模版。建图：每个s到t连一条容量为1，代价为v的边，每一天与下一天连一条容量为2，代价为0的边，S与第0天连边，T与最后一天连边，容量为2，代价为0，最大费用最大流即为答案。

#include<cstdio>   
#include<cstring>   
#include<algorithm>   
#include<queue>   
#define INF 0x7fffffff   
using namespace std;   
const int MAXN=500,MAXM=4000;   
struct mcflow{   
       int N,M,S,T;   
       int head[MAXN],next[MAXM],to[MAXM],cost[MAXM],f[MAXM],ecnt;   
       int pre[MAXN],d[MAXN];   
       bool vis[MAXN];   
       #define CL(x) memset(x, 0, sizeof(x));   
       void init(){   
            CL(head); CL(next); CL(to); CL(cost);   
            ecnt = 2;   
            T = S = N = M = 0;   
       }   
       void make_edge(int a, int b, int cc, int d){   
            next[ecnt] = head[a];   
            head[a] = ecnt;   
            to[ecnt] = b;   
            f[ecnt] = cc;   
            cost[ecnt++] = d;   
       }   
       void insert(int a, int b, int cc, int d){   
            make_edge(a, b, cc, d);   
            make_edge(b, a, 0, -d);   
       }   
       bool spfa(){   
            CL(vis);   
            for(int i = 0; i <= T; ++i)   
                d[i] = INF;   
            queue<int> Q;   
            Q.push(S);   
            vis[S]=1;   
            d[S] = 0;   
            while(!Q.empty()){   
                int u = Q.front(); Q.pop();   
                vis[u] = 0;   
                for(int p = head[u]; p; p = next[p]){   
                    if(f[p] > 0 && d[to[p]] > d[u] + cost[p]){   
                        d[to[p]] = d[u] + cost[p];   
                        pre[to[p]] = p;   
                        if(!vis[to[p]]){   
                            vis[to[p]]=1;   
                            Q.push(to[p]);   
                        }   
                    }   
                }   
            }   
            return d[T] < INF;   
       }   
       void min_cost_flow(int &Flow, int &fee){   
            Flow = fee = 0;   
            while(spfa()){   
                fee += d[T];   
                int u = T, tmp = INF;   
                while(u != S){   
                    tmp=min(f[pre[u]], tmp);   
                    u=to[pre[u]^1];   
                }   
                u=T;   
                while(u != S){   
                    f[pre[u]] -= tmp;   
                    f[pre[u]^1] += tmp;   
                    u = to[pre[u]^1];   
                }   
                Flow += tmp;   
            }   
       }   
       int mincost(){   
           int ret, tmp;   
           min_cost_flow(tmp, ret);   
           return ret;   
       }   
       int maxflow(){   
           int ret, tmp;   
           min_cost_flow(ret, tmp);   
           return ret;   
       }   
}G;   
       
int main()   
{   
    int T,m,s,t,v,i;   
    scanf("%d",&T);   
    while(T--){   
        G.init();   
        scanf("%d",&m);   
        while(m--){   
            scanf("%d%d%d",&s,&t,&v);   
            G.insert(s,t+1,1,-v);   
        }   
        G.S = 401; G.T = 402;   
        G.insert(G.S,0,2,0); G.insert(400,G.T,2,0);   
        for(i = 0; i < 400; ++i) G.insert(i,i+1,2,0);   
        printf("%d\n",-G.mincost());   
    }   
    return 0;   
}

---

更新（2015年11月7日）：

上面的模板只在每次流量增量是1的时候才是对的。增量不是1就会出错。