NOIP 2009 最优贸易

nkoj 1692

Description

(trade.pas/c/cpp) 
【问题描述】 
　　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 5 6 1
1 2 1
1 4 1
2 3 2
3 5 1
4 5 2

Sample Output

5

呼～撸了一个小时终于做了出来

数据结构:

edge1[]   last[0][]  next[0][] 存正图 ，edge2[] , last[1][] , next[1][] 存反图

minn 表示最小买入价格， maxx表示最大卖出价格

代码如下：

#include<cstdio>   
#include<queue>
#include<cstring>
#include<iostream>
using namespace std;
const int inf=100000000,maxn=100000+5,maxm=500000+5;
int n,m,p=0,s[maxn],next[2][maxm],last[2][maxm];
int minn[maxn],maxx[maxn];
bool vis[maxn];
queue <int> q;
struct EDGE{
	int from,to;
};
EDGE edge1[maxm],edge2[maxm];
void insert(EDGE e[],int x,int y,int z){ 
     //向e[]数组中加入一个x->y的边 
	e[p].from=x;
	e[p].to=y;
	next[z][p]=last[z][x];
	last[z][x]=p;
}
void _read(int &x){
    char t;
    bool mark=false;
    for(;t=getchar(),t<'0'||t>'9';) if(t=='-') mark=1;
    for(x=t-'0',t=getchar();'0'<=t&&t<='9';x=x*10+t-'0',t=getchar());
    x=mark?-x:x;
}

void input(){
	int i,x,y,z;
	_read(n);_read(m);
	for(i=1;i<=n;i++)_read(s[i]);
	for(i=1;i<=m;i++){
		p++;
		_read(x);_read(y);_read(z);
		insert(edge1,x,y,0);
		insert(edge2,y,x,1);
		if(z==2){
			p++;
			insert(edge2,x,y,1);
			insert(edge1,y,x,0);
		}
	}
	m=p;
}

void _spfa1(){   //在反图中找最大 
	q.push(n);
	vis[n]=true;
	maxx[n]=s[n];
	while(!q.empty()){
		int x=q.front();
		q.pop();vis[x]=false;
		for(int i=last[1][x];i;i=next[1][i]){
			EDGE e=edge2[i] ;
			if(maxx[e.from]>maxx[e.to]){
				maxx[e.to]=maxx[e.from];
				if(!vis[e.to]){
					q.push(e.to);vis[e.to]=true;
				}
			}
			if(s[e.to]>maxx[e.to]){
				maxx[e.to]=s[e.to];
				if(!vis[e.to]){
					q.push(e.to);vis[e.to]=true;
				}
			}
		}
	}
}

void _spfa2(){  //在正图找最小 
	q.push(1);
	vis[1]=true;
	minn[1]=s[1];
	while(!q.empty()){
		int x=q.front();
		q.pop();vis[x]=false;
		for(int i=last[0][x];i;i=next[0][i]){
			EDGE e=edge1[i] ;
			if(minn[e.from]<minn[e.to]){
				minn[e.to]=minn[e.from];
				if(!vis[e.to]){
					q.push(e.to);vis[e.to]=true;
				}
			}
			if(s[e.to]<minn[e.to]){
				minn[e.to]=s[e.to];
				if(!vis[e.to]){
					q.push(e.to);vis[e.to]=true;
				}
			}
		}
	}
}
void solve(){
	int i,ans=0;
	for(i=0;i<=n;i++)minn[i]=inf;
	memset(vis,0,sizeof(vis));
	_spfa1();
	memset(vis,0,sizeof(vis));
	_spfa2();
	for(i=1;i<=n;i++)  //找答案 
		ans=max(ans,maxx[i]-minn[i]);
	printf("%d\n",ans);
}
int main(){
	input();
	solve();
}