洛谷 P2153 [SDOI2009]晨跑（费用流）

题目描述

Elaxia最近迷恋上了空手道，他为自己设定了一套健身计划，比如俯卧撑、仰卧起坐等 等，不过到目前为止，他坚持下来的只有晨跑。 现在给出一张学校附近的地图，这张地图中包含N个十字路口和M条街道，Elaxia只能从 一个十字路口跑向另外一个十字路口，街道之间只在十字路口处相交。Elaxia每天从寝室出发 跑到学校，保证寝室编号为1，学校编号为N。 Elaxia的晨跑计划是按周期（包含若干天）进行的，由于他不喜欢走重复的路线，所以 在一个周期内，每天的晨跑路线都不会相交（在十字路口处），寝室和学校不算十字路 口。Elaxia耐力不太好，他希望在一个周期内跑的路程尽量短，但是又希望训练周期包含的天 数尽量长。 除了练空手道，Elaxia其他时间都花在了学习和找MM上面，所有他想请你帮忙为他设计 一套满足他要求的晨跑计划。

存在1\rightarrow n1→n的边存在。这种情况下，这条边只能走一次。

输入输出格式

输入格式：

 

第一行：两个数N,M。表示十字路口数和街道数。 接下来M行，每行3个数a,b,c，表示路口a和路口b之间有条长度为c的街道（单向）。

 

输出格式：

 

两个数，第一个数为最长周期的天数，第二个数为满足最长天数的条件下最短的路程长 度。

 

输入输出样例

输入样例#1： 

7 10
1 2 1
1 3 1
2 4 1
3 4 1
4 5 1
4 6 1
2 5 5
3 6 6
5 7 1
6 7 1

输出样例#1： 

2 11

说明

对于30%的数据，N ≤ 20，M ≤ 120。

对于100%的数据，N ≤ 200，M ≤ 20000。

解题思路

以1为源点n为会汇点，将2~n-1的点拆分为i1，i2，每个i1到i2建一条flow为1，dis为0的边，然后读入a,b,c。建边(a2, b1)，flow为1，dis为c（如果a为1，则连（1，b1），若b为n，则连（a2, n））。反向弧的flow为0，dis为-c。然后求最大流和，最小费用最大流。

代码如下

#include <iostream>
#include <cstring>
#include <queue>
#include <vector>
#include <cmath>
#define INF 0x3f3f3f3f
using namespace std;
vector<int> g[405];
struct Line{
	int r, flow, dis;
	Line(int r, int f, int d): r(r), flow(f), dis(d){	}
};
vector<Line> line;
void add_line(int x, int y, int f, int d)
{
	line.push_back(Line(y, f, d));
	g[x].push_back(line.size() - 1);
	line.push_back(Line(x, 0, -d));
	g[y].push_back(line.size() - 1);
}
int pre[405];
int dis[405];
int flow[405];
int edge[405];
int s, t;
bool spfa()
{
	bool vis[405] = {0};
	memset(pre, -1, sizeof(pre));
	memset(edge, -1, sizeof(edge)); 
	memset(dis, 0x7f, sizeof(dis));
	memset(flow, 0x7f, sizeof(flow));
	queue<int> que;
	que.push(s);
	vis[s] = true;
	dis[s] = 0;
	while(!que.empty()){
		int top = que.front();
		que.pop();
		vis[top] = false;
		//cout << top << endl;
		for(int i = 0; i < g[top].size(); i ++){
			int z = g[top][i];
			int r = line[z].r;
			//cout << r << endl;
			if(line[z].flow && dis[top] + line[z].dis < dis[r]){
				pre[r] = top;
				edge[r] = z;
				dis[r] = dis[top] + line[z].dis;
				flow[r] = min(flow[top], line[z].flow);
				if(!vis[r]){
					que.push(r);
					vis[r] = true;
				}
			}
		}
	}
	return pre[t] != -1;
}
void EK()
{
	int max_flow = 0;
	int min_dis = 0;
	while(spfa()){
		//cout << "??" << endl;
		max_flow += flow[t];
		min_dis += flow[t] * dis[t];
		int now = t;
		while(pre[now] != -1){
			line[edge[now]].flow -= flow[t];
			line[edge[now]^1].flow += flow[t];
			now = pre[now];
		}
	}
	cout << max_flow << " " << min_dis << endl;
}
int main()
{
	int n, m;
	while(cin >> n >> m){
		s = 3;
		t = 2*n;
		for(int i = 2; i < n; i ++)
			add_line(2*i, 2*i+1, 1, 0);
		for(int i = 0; i < m; i ++){
			int x, y, d;
			cin >> x >> y >> d;
			add_line(2*x+1, 2*y, 1, d);
		}
		EK();
		for(int i = 1; i <= n; i ++)
			g[i].clear();
		line.clear();
	}
	return 0;
}