NOIP2009提高组 最优贸易

传送门

题目

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之间的双向道路。
输出格式

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

过程

during the exam

这是昨晚考试的最优一道题。当时读到题的时候我是懵逼的，因为我以前刷题的时候将这道题读过n遍，然而并没有什么思路。但是既然考试中有这道题，那能得几分是几分呗，反正得分比不得分强。当时用的是floyed的变形，即

        for(k=1;k<=n;k++)
          for(i=1;i<=n;i++)
            for(j=1;j<=n;j++)
              if(b[i][k]&&b[k][j]) b[i][j]=1;
        for(i=1;i<=n;i++)
          for(j=1;j<=n;j++)
            if(b[1][i]&&b[i][j]&&b[j][n]) ans=max(ans,w[j]-w[i]);
        printf("%d",ans);

我本来是奔着30分去的，然而由于手残+样例的恶毒，我在做差的时候写成了w[i]-w[j]，而且样例过了！结果我就不多说了，自己脑补······
今天上网找了n篇题解，终于找到了一篇适合自己读的。

正解思路：

两遍SPFA。第一遍求出在路上经过该点时能得到的最小买入价格，起点为1；第二遍将所有边反过来，求出在路上经过该点时能得到的最大卖出价格，起点为n。最后O(n)计算差值比较即可。
由于我以前用的不是这种邻接表存储，所以在调试程序的过程中又耗费了很长时间······
代码：

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
int n,m,i,x,y,z,num1,num2,h,t,ans;
int w[100001];
struct aa
{
    int from,to,next;
};
aa a[500001];
aa b[500001];
int heada[100001];
int headb[100001];
int maxn[100001];
int minn[100001];
int f[2000001];
bool bo[100001];
void add(aa a[],int head[],int from,int to,int &num)
{
    a[num].from=from;
    a[num].to=to;
    a[num].next=head[from];
    head[from]=num++;
}
int main()
{
    scanf("%d%d",&n,&m);
    for(i=1;i<=n;i++)
      scanf("%d",&w[i]);
    for(i=1;i<=m;i++)
    {
        scanf("%d%d%d",&x,&y,&z);
        if(z==1) {add(a,heada,x,y,num1);add(b,headb,y,x,num2);}
        else {add(a,heada,x,y,num1);add(a,heada,y,x,num1);add(b,headb,x,y,num2);add(b,headb,y,x,num2);}
    }
    memset(minn,100,sizeof(minn));
    memset(maxn,-1,sizeof(maxn));
    h=t=1;
    f[1]=1;
    minn[1]=w[1];
    bo[1]=1;
    while(h<=t)//最小值 
    {
        int tmp=f[h++];
        bo[tmp]=0;
        for(i=heada[tmp];i;i=a[i].next)
          if(minn[tmp]<minn[a[i].to]||w[a[i].to]<minn[a[i].to])
          {
            minn[a[i].to]=min(w[a[i].to],minn[tmp]);
            if(!bo[a[i].to])
            {
                f[++t]=a[i].to;
                bo[a[i].to]=1;
            }
          }
    }
    h=t=1;
    f[1]=n;
    maxn[n]=w[n];
    bo[n]=1;
    while(h<=t)//最大值 
    {
        int tmp=f[h++];
        bo[tmp]=0;
        for(i=headb[tmp];i;i=b[i].next)
          if(maxn[tmp]>maxn[b[i].to]||w[b[i].to]>maxn[b[i].to])
          {
            maxn[b[i].to]=max(maxn[tmp],w[b[i].to]);
            if(!bo[b[i].to])
            {
                f[++t]=b[i].to;
                bo[b[i].to]=1;
            }
          }
    }
    for(i=1;i<=n;i++)
      ans=max(ans,maxn[i]-minn[i]);
    printf("%d",ans);
}