NOIP 2009 最优贸易 luogu1073

(http://www.elijahqi.win/2017/07/10/noip-2009-%E6%9C%80%E4%BC%98%E8%B4%B8%E6%98%93-luogu1073/)
题目描述

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 条道路的信息（每条道路所连接的两个城市的编号

以及该条道路的通行情况）。请你告诉阿龙，他最多能赚取多少旅费。

输入输出格式

输入格式：

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

数目。

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

市的商品价格。

接下来 m 行，每行有 3 个正整数，x，y，z，每两个整数之间用一个空格隔开。如果 z=1，

表示这条道路是城市 x 到城市 y 之间的单向道路；如果 z=2，表示这条道路为城市 x 和城市

y 之间的双向道路。

输出格式：

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

则输出 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
说明

【数据范围】

输入数据保证 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。

NOIP 2009 提高组 第三题

思路 ：两次spfa从起点到终点 第一次找出最小的买入点 第二次找出最大的卖出点

#include<cstdio>
#include<cstring>
#include<queue>
#define N 110000
using namespace std;
queue<int>q;
int num,w[N],f[N],f1[N],flag[N],h[N],h2[N],n,m;
int max1,min1;
struct node{
    int x,y,next;
}data[N*10];
int read(){
    int x=0;char ch=getchar();
    while (ch<'0'||ch>'9') ch=getchar();
    while ('0'<=ch&&ch<='9') {x=x*10+ch-'0';ch=getchar(); }
    return x;
}
void insert1(int x,int y){
    data[++num].x=x;data[num].y=y;data[num].next=h[x];h[x]=num;
    data[++num].x=y;data[num].y=x;data[num].next=h2[y];h2[y]=num;
}
inline int min(int x,int y){
    return x<y?x:y;
}
inline int max(int x,int y){
    return x>y?x:y;
}
void spfa1(){
    for (int i=1;i<=n;++i) f[i]=0x7fffffff;memset(flag,false,sizeof(flag));
    f[1]=w[1];flag[1]=true;q.push(1);
    while(!q.empty()){
        int tmp=q.front();q.pop();
        for (int i=h[tmp];i;i=data[i].next){
            int x=data[i].x,y=data[i].y;
            if (min(w[y],f[x])<f[y]){
                f[y]=min(w[y],f[x]);
                if (flag[y]==false){
                    flag[y]=true;q.push(y);
                }
            }
        }
        flag[tmp]=false;
    }
}
void spfa2(){
    for (int i=1;i<=n;++i) f1[i]=-1;memset(flag,false,sizeof(flag));
    f1[n]=w[n];flag[n]=true;q.push(n);
    while(!q.empty()){
        int tmp=q.front();q.pop();
        for (int i=h2[tmp];i;i=data[i].next){
            int x=data[i].x,y=data[i].y;
            if (max(f1[x],w[y])>f1[y]){
                f1[y]=max(f1[x],w[y]);
                if (flag[y]==false){
                    flag[y]=true;q.push(y);
                }
            }
        }
        flag[tmp]=false;
    }
}
int main (){
    freopen("1073.in","r",stdin);
    freopen("1073.out","w",stdout);
    n=read();m=read();num=0;
    memset(h,0,sizeof(h));
    for (int i=1;i<=n;++i) w[i]=read();
    for (int i=1;i<=m;++i){
        int x=read(),y=read(),t=read();
        if (t==1) insert1(x,y);else insert1(x,y),insert1(y,x);
    }
    /*for (int i=1;i<=num;++i){
        printf("%d %d\n",data[i].x,data[i].y);
    }*/
    spfa1();
    //printf("%d ",min1);
    spfa2();
    int ans=0; 
    for (int i=1;i<=n;++i) {
        ans=max(ans,f1[i]-f[i]);
    } 
    printf("%d",ans);
    return 0;
}