1. 程式人生 > >1877. [SDOI2009]晨跑【費用流】

1877. [SDOI2009]晨跑【費用流】

clu use post AC while tdi nod %d 最短

Description

Elaxia最近迷戀上了空手道,他為自己設定了一套健身計劃,比如俯臥撐、仰臥起坐等 等,不過到目前為止,他 堅持下來的只有晨跑。 現在給出一張學校附近的地圖,這張地圖中包含N個十字路口和M條街道,Elaxia只能從 一 個十字路口跑向另外一個十字路口,街道之間只在十字路口處相交。Elaxia每天從寢室出發 跑到學校,保證寢室 編號為1,學校編號為N。 Elaxia的晨跑計劃是按周期(包含若幹天)進行的,由於他不喜歡走重復的路線,所以 在一個周期內,每天的晨跑路線都不會相交(在十字路口處),寢室和學校不算十字路 口。Elaxia耐力不太好, 他希望在一個周期內跑的路程盡量短,但是又希望訓練周期包含的天 數盡量長。 除了練空手道,Elaxia其他時間 都花在了學習和找MM上面,所有他想請你幫忙為他設計 一套滿足他要求的晨跑計劃。

Input

第一行:兩個數N,M。表示十字路口數和街道數。 接下來M行,每行3個數a,b,c,表示路口a和路口b之間有條長度為c的街道(單向)。 N ≤ 200,M ≤ 20000。

Output

兩個數,第一個數為最長周期的天數,第二個數為滿足最長天數的條件下最短的路程長 度。

Sample Input

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

Sample Output

2 11 一道近乎最小費用最大流的模板題
唯一和模板不同的就是我們要限制每個點只能到一次。
這一看就是常規拆點啊
所以
我們將點裂成兩個,然後在兩點間連一個容量為1費用為0的邊,用來限制每個點只能走一次
然後再在x+n和y中的連邊,把容量設為1,費用設為邊長
跑一邊從(n+1)到n的最小費用最大流(因為1節點可以重復走所以要n+1)
最大流和最小費用即為答案。
#include<iostream>
#include<cstring>
#include<cstdio>
#include<queue>
using namespace std;
struct node
{
	int Flow;
	int Cost;
	int to;
	int next;
}edge[100010];
queue<int>q;
int INF,dis[601],pre[601];
int head[601],num_edge;
bool used[601];

void add(int u,int v,int l,int c)
{
	edge[++num_edge].to=v;
	edge[num_edge].Flow=l;
	edge[num_edge].Cost=c;
	edge[num_edge].next=head[u];
	head[u]=num_edge;
}

bool SPFA(int s,int e)
{
	memset(pre,-1,sizeof(pre));
	memset(dis,0x7f,sizeof(dis));
	q.push(s);
	dis[s]=0;
	used[s]=true;
	while (!q.empty())
	{
		int x=q.front();
		q.pop();
		for (int i=head[x];i!=0;i=edge[i].next)
			if (edge[i].Flow>0 && dis[x]+edge[i].Cost<dis[edge[i].to])
			{
				dis[edge[i].to]=dis[x]+edge[i].Cost;
				pre[edge[i].to]=i;
				if (!used[edge[i].to])
				{
					used[edge[i].to]=true;
					q.push(edge[i].to);
				}
			}
		used[x]=false;
	}
	return (dis[e]!=INF);
}

void MCMF(int s,int e)
{
	int d=INF,Ans=0,Fee=0;
	while (SPFA(s,e))
	{
		for (int i=e;i!=s;i=i=edge[((pre[i]-1)^1)+1].to)
			d=min(edge[pre[i]].Flow,d);
		for (int i=e;i!=s;i=i=edge[((pre[i]-1)^1)+1].to)
		{
			edge[pre[i]].Flow-=d;
			edge[((pre[i]-1)^1)+1].Flow+=d;
		}
		Ans+=d;
		Fee+=dis[e]*d;
	}
	printf("%d %d",Ans,Fee);
}

int main()
{
	memset(&INF,0x7f,sizeof(INF));
	int n,m,u,v,l;
	scanf("%d%d",&n,&m);
	for (int i=1;i<=n;++i)
	{
		add(i,i+n,1,0);
		add(i+n,i,0,0);
	}
	for (int i=1;i<=m;++i)
	{
		scanf("%d%d%d",&u,&v,&l);
		add(u+n,v,1,l);
		add(v,u+n,0,-l);
	}
	MCMF(n+1,n);
}

1877. [SDOI2009]晨跑【費用流】