NOIP2009 最优贸易

Description

　　C国有n个大城市和m条道路，每条道路连接这n个城市中的某两个城市。任意两个城市之间最多只有一条道路直接相连。这m条道路中有一部分为单向通行的道路，一部分为双向通行的道路，双向通行的道路在统计条数时也计为1条。
　　C国幅员辽阔，各地的资源分布情况各不相同，这就导致了同一种商品在不同城市的价格不一定相同。但是，同一种商品在同一个城市的买入价和卖出价始终是相同的。
　　商人阿龙来到C国旅游。当他得知同一种商品在不同城市的价格可能会不同这一信息之后，便决定在旅游的同时，利用商品在不同城市中的差价赚回一点旅费。设C国n个城市的标号从1-n，阿龙决定从1号城市出发，并最终在n号城市结束自己的旅行。在旅游的过程中，任何城市可以重复经过多次，但不要求经过所有n个城市。阿龙通过这样的贸易方式赚取旅费：他会选择一个经过的城市迈入他最喜欢的商品——水晶球，并在之后经过的另一个城市卖出这个水晶球。用赚取的差价当作旅费。由于阿龙主要是来C国旅游，他决定这个贸易只进行最多一次。当然，在赚不到差价的情况下它就无需进行贸易。
　　假设C国有5个大城市，城市的编号和道路连接情况如下图，单向箭头表示这条道路为单向通行。双向箭头表示这条道路为双向通行。
　　假设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条道路的信息（每条道路所连接的两个城市的编号以及该条道路的通行情况）。请你告诉阿龙，他最多能赚钱多少旅费。

Input

　　输入文件名为trade.in。
　　第一行包含2个正整数n和m，中间用一个空格隔开，分别表示城市的数目和道路的数目。
　　第二行n个正整数，每两个正整数之间用一个空格隔开，按标号顺序分别表示这n个城市的商品价格。
　　接下来m行，每行有3个正整数，x，y，z，每两个整数之间用一个空格隔开。如果z=1，表示这条道路是城市x到城市y之间的单向道路；如果z=2，表示这条道路为城市x和城市y之间的双向道路。

Output

　　输出文件trade.out共1行，包含1个整数，表示最多能赚取的旅费。如果没有进行贸易，则输出0。

Sample Input

5 5

4 3 6 5 1

1 2 1

1 4 1

2 3 2

3 5 1

4 5 2

Sample Output

4

Hint

【数据范围】
　　输入数据保证1号城市可以到达n号城市。
　　对于10%的数据，1≤n≤6。
　　对于30%的数据，1≤n≤100。
　　对于50%的数据，不存在一条旅游路线，可以从一个城市出发，再回到这个城市。
　　对于100%的数据，1≤n≤100000，1≤m≤500000，1≤x，y≤n，1≤z≤2，1≤各城市水晶球价格≤100。

缩点+DAG的dp，我一开始居然sb的用dfs跑，tle飞还不知到原因。。。

顺便复习一下tarjan和DAG拓扑排序。

#include<bits/stdc++.h>
using namespace std;
#define Inc(i,L,r) for(register int i=(L);i<=(r);++i)
const int N = 1e5+10;
int n,m,w[N];
struct Edge{
	int cnt,h[N],to[N<<1],next[N<<1];
	inline void add(int x,int y){
		next[++cnt]=h[x];
		to[cnt]=y;
		h[x]=cnt;
	}
}e;
inline void init(){
	scanf("%d%d",&n,&m);
	Inc(i,1,n)scanf("%d",&w[i]);
	Inc(i,1,m){
		int x,y,z;scanf("%d%d%d",&x,&y,&z);
		e.add(x,y);if(z==2)e.add(y,x);
	}
}
int sign,tj[N],ll[N];
int top,s[N],istk[N];
int scc,bl[N];
int mn[N],mx[N];
vector<int>vec[N];
void Tarjan(int u){
	tj[u]=ll[u]=++sign;
	s[++top]=u;
	istk[u]=1;
	int v;
	for(int p=e.h[u];p;p=e.next[p]){
		v=e.to[p];
		if(!tj[v]){
			Tarjan(v);
			ll[u]=min(ll[u],ll[v]);
		}
		else if(istk[v])ll[u]=min(ll[u],tj[v]);
	}
	if(ll[u]==tj[u]){
		scc++;
		mx[scc]=-(1<<30),mn[scc]=1<<30;
		do{
			v=s[top--];
			istk[v]=0;
			vec[bl[v]=scc].push_back(v);
			mx[scc]=max(mx[scc],w[v]);
			mn[scc]=min(mn[scc],w[v]);
		}while(v!=u);
	}
}
bool inq[N];
int f[N][2];
queue<int>Q;
inline void bfs(){
	Q.push(bl[1]);
	f[bl[1]][0]=mn[bl[1]];
	f[bl[1]][1]=mx[bl[1]]-mn[bl[1]];
	while(!Q.empty()){
		int nc=Q.front();Q.pop();
		inq[nc]=0;
		int size=vec[nc].size();
		Inc(i,0,size-1){
			int x=vec[nc][i];
			for(int p=e.h[x];p;p=e.next[p]){
				int to=e.to[p];
				if(bl[to]==nc||inq[bl[to]])continue;
				inq[bl[to]]=1;
				Q.push(bl[to]);
				f[bl[to]][0]=min(mn[bl[to]],f[nc][0]);
				f[bl[to]][1]=max(mx[bl[to]]-f[bl[to]][0],f[nc][1]);
			}
		}
	}
}
inline void solv(){
	Tarjan(1);
	bfs();
	cout<<f[bl[n]][1]<<"\n";
}
int main(){
	init();
	solv();
	return 0;
}