洛谷 P1073 最优贸易

最优贸易 【luogu P1073 [NOIP2009 提高组] 最优贸易】

题目描述：

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

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

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

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

思路：

建一个正图和反图，正图使从1点出发到每一个点的w最小，即购入价最小，反图从n出发,使卖出价（卖出价和购进价是一样的）最大，$ans=max(ans,dp[i]-dis[i])$ $dp[i]$是反图跑出的最大值，$dis[i]$是正图跑出的最小值
代码使用了$dijkstra$+堆优化 （垃圾SPFA毁我青春,天天TLE）

代码：

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<queue>
#include<vector>
#define r register
#define rep(i,x,y) for(r int i=x;i<=y;++i)
#define per(i,x,y) for(r int i=x;i>=y;--i)
using namespace std;
typedef long long ll;
const int V=(1e5+10)*5;
ll n,m,a[V>>1],head1[V>>1],top,head2[V>>1];
ll x,y,k,dis[V>>1],dp[V>>1],tot,ans;
bool vis[V>>1],pd[V>>1];
struct node
{
	ll nxt,to;
}e1[V<<1],e2[V<<1];
void add1(ll x,ll y)
{
	e1[++top]=(node){head1[x],y};
	head1[x]=top;
}
void add2(ll x,ll y)
{
	e2[++tot]=(node){head2[x],y};
	head2[x]=tot;
}
void dijkstra1()
{
	memset(dis,0x3f,sizeof(dis));
	priority_queue< pair<ll,ll> > q;
//	vis[1]=true;
	dis[1]=a[1];
	q.push(make_pair(-dis[1],1));
	while(!q.empty() )
	{
		ll x=q.top().second ;
		q.pop();
		if(vis[x]) continue ;
		vis[x]=true;
		for(r ll i=head1[x];i;i=e1[i].nxt)
		{
			ll y=e1[i].to;
			if(dis[y]>min(dis[x],a[y]))
			{
				dis[y]=min(dis[x],a[y]);
				q.push(make_pair(-dis[y],y));
			}
		}
		
	}
}
void dijkstra2()
{
	memset(dp,0,sizeof(dp));
	priority_queue< pair<ll,ll> > q;
//	pd[n]=true;
	dp[n]=a[n];
	q.push(make_pair(dp[n],n));
	while(!q.empty() )
	{
		ll x=q.top().second ;
		q.pop();
		if(pd[x]) continue ;
		pd[x]=true;
		for(r ll i=head2[x];i;i=e2[i].nxt)
		{
			ll y=e2[i].to;
			if(dp[y]<max(dp[x],a[y]))
			{
				dp[y]=max(dp[x],a[y]);
				q.push(make_pair(dp[y],y));
			}
		}
		
	}
}
int main()
{
	scanf("%lld%lld",&n,&m);
	rep(i,1,n) scanf("%lld",&a[i]);
	rep(i,1,m)
	{
		scanf("%lld%lld%lld",&x,&y,&k);
		add1(x,y);
		add2(y,x);
		if(k==2) add1(y,x),add2(x,y);
	}
	dijkstra1();
	dijkstra2();
	rep(i,1,n) ans=max(ans,dp[i]-dis[i]);
	cout<<ans;
 	return 0;
}