最优贸易[Tanjan][拓扑排序][DP]

题目描述

C国有n个大城市和m 条道路，每条道路连接这 n个城市中的某两个城市。任意两个城市之间最多只有一条道路直接相连。这 m 条道路中有一部分为单向通行的道路，一部分为双向通行的道路，双向通行的道路在统计条数时也计为 1条。

C国幅员辽阔，各地的资源分布情况各不相同，这就导致了同一种商品在不同城市的价格不一定相同。但是，同一种商品在同一个城市的买入价和卖出价始终是相同的。

商人阿龙来到 C 国旅游。当他得知同一种商品在不同城市的价格可能会不同这一信息之后，便决定在旅游的同时，利用商品在不同城市中的差价赚回一点旅费。设 C 国 n 个城市的标号从 1~ n，阿龙决定从 1号城市出发，并最终在 n 号城市结束自己的旅行。在旅游的过程中，任何城市可以重复经过多次，但不要求经过所有 n 个城市。阿龙通过这样的贸易方式赚取旅费：他会选择一个经过的城市买入他最喜欢的商品――水晶球，并在之后经过的另一个城市卖出这个水晶球，用赚取的差价当做旅费。由于阿龙主要是来 C 国旅游，他决定这个贸易只进行最多一次，当然，在赚不到差价的情况下他就无需进行贸易。

假设 C国有 55个大城市，城市的编号和道路连接情况如下图，单向箭头表示这条道路为单向通行，双向箭头表示这条道路为双向通行.

假设 1~n 号城市的水晶球价格分别为 4,3,5,6,1。

阿龙可以选择如下一条线路：1->2->3->5，并在 2号城市以3 的价格买入水晶球，在 3号城市以5的价格卖出水晶球，赚取的旅费数为 2。

阿龙也可以选择如下一条线路1->4->5->4->5，并在第1次到达5 号城市时以 1的价格买入水晶球，在第 2 次到达4 号城市时以6的价格卖出水晶球，赚取的旅费数为5。

现在给出 n个城市的水晶球价格，m 条道路的信息（每条道路所连接的两个城市的编号以及该条道路的通行情况）。请你告诉阿龙，他最多能赚取多少旅费。

输入格式：

第一行包含 2 个正整数n和 m，中间用一个空格隔开，分别表示城市的数目和道路的数目。

第二行 n 个正整数，每两个整数之间用一个空格隔开，按标号顺序分别表示这 n 个城市的商品价格。

接下来 m 行，每行有3个正整数x,y,z每两个整数之间用一个空格隔开。如果 z=1，表示这条道路是城市x到城市y之间的单向道路；如果z=2，表示这条道路为城市 x和城市y之间的双向道路。

输出格式：

一 个整数，表示最多能赚取的旅费。如果没有进行贸易，则输出 0。

输入样例#1： 

5 5 
4 3 5 6 1 
1 2 1 
1 4 1 
2 3 2 
3 5 1 
4 5 2 

输出样例#1： 

5

n<=100000,m<=500000

分析

最大值问题一般用dp或贪心解决

不知为何,我一下就想到了tarjan

一个强联通分量中的最大收益就是这个分量中的最大值减去最小值

同时,这个图变成了有向无环图

我们就可以dp了

dp[i] 表示到第i个强联通分量的最大收益

dp[i]=max(max(dp[i],dp[j]),Max[i]-Rmin[i])

Max[i]表示第i个强联通分量的最大值,Rmin[i]表示1--i路径上的最小值

Rmin[i]=min(Rmin[j],Min[i])

Min[i]表示第i个强联通分量的最小值

建议用拓扑排序来递推,dfs容易写挂

拓扑排序的优势是,当遍历到i时,i一定是最优的

...哎,想了好久,不好表述啊,反正是很优秀的,感性理解一下

#include<bits/stdc++.h>
#define N 100005
#define M 500005
using namespace std;
int first[N],next[M*2],to[M*2],tot;
int id[N],dfn[N],low[N],sta[N],insta[N],top,sign,sum;
int a[N],n,m,vis[N],du[N],b[N],ret;
int Min[N],Max[N],dp[N],Rmin[N];
int f2[N],n2[N],t2[N],tot2;
int read(){
	int cnt=0;char ch=0;
	while(!isdigit(ch))ch=getchar();
	while(isdigit(ch))cnt=cnt*10+(ch-'0'),ch=getchar();
	return cnt; 
}
void add(int x,int y){
	next[++tot]=first[x],first[x]=tot,to[tot]=y;
}
void add1(int x,int y){
	n2[++tot2]=f2[x],f2[x]=tot2,t2[tot2]=y;
}
void tarjan(int cur){
	dfn[cur]=low[cur]=++sign,sta[++top]=cur,insta[cur]=1;
	for(int i=first[cur];i;i=next[i]){
		int t=to[i];
		if(!dfn[t]) tarjan(t),low[cur]=min(low[cur],low[t]);
		else if(insta[t]&&low[cur]>dfn[t]) low[cur]=dfn[t];
	}
	if(low[cur]==dfn[cur]){
		sum++;do{
		insta[sta[top]]=0,id[sta[top]]=sum;
		Min[sum]=min(Min[sum],a[sta[top]]);
		Max[sum]=max(Max[sum],a[sta[top]]);
		}while(sta[top--]!=cur);
	}
}
void dfs(int cur){
	queue<int> q;
	for(int i=1;i<=n;i++)
		for(int j=first[i];j;j=next[j])
			if(id[to[j]]!=id[i]) add1(id[i],id[to[j]]),du[id[to[j]]]++;
	for(int i=1;i<=sum;i++) if(du[i]==0) q.push(i);
	while(!q.empty()){
		int x=q.front();q.pop();
		b[++ret]=x;
		for(int i=f2[x];i;i=n2[i]){
			if(!(--du[t2[i]])) q.push(t2[i]);
		}
	}
	for(int i=1;i<=ret;i++){
		int cur=b[i];
		for(int j=f2[cur];j;j=n2[j]){
		int t=t2[j];
		Rmin[t]=min(Rmin[cur],Min[t]);
		dp[t]=max(max(dp[t],dp[cur]),Max[t]-Rmin[t]);
		}
	}
}
int main(){
	memset(Min,127,sizeof(Min));
	memset(Max,128,sizeof(Max));
	memset(Rmin,127,sizeof(Rmin));
	n=read(),m=read();
	for(int i=1;i<=n;i++)a[i]=read();
	for(int i=1;i<=m;i++){
		int x=read(),y=read(),z=read();
		add(x,y);if(z==2)add(y,x);
	}
	for(int i=1;i<=n;i++) if(!dfn[i]) tarjan(i);
	dfs(1);cout<<dp[id[n]];
	return 0;
}