【图-最短路】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 条道路的信息（每条道路所连接的两个城市的编号以及该条道路的通行情况）。请你告诉阿龙，他最多能赚取多少旅费。

【问题分析】
很容易判断是spfa，问题是在spfa的基础上应该维护一个什么域呢？

最开始考试的时候写的是维护一个max域表示经过此处的最大商品卖出值，min表示最小的买入值。结果样例对了，自己造了一个也对了，就交上去了，结果……………………
那么，我们到底应该维护什么呢？
在对原来的程序调试时(拿到数据之后)发现错误的原因是对于一个城市来说，最大的商品贸易差应该是之前的商品买入最小值与当前卖出值的差。所以想到了用min维护一个到当前节点之前的最小的商品买入价格，用f表示到当前这个节点所能获得的最大利润。
dp方程：f[e[t].to]=max(f[i],w[e[t].to]-mmin[i])

相信大家都能看懂，下面上代码！(很丑陋。。)

#include <iostream>
#include <cstdio>
#include <queue>
#include <climits>
using namespace std;
const int N=100001;
const int M=1000001;
const int INF=INT_MIN;
struct zk { int from,to,next; };  zk e[M];
int n,m,tail;  int w[N],first[N];
queue<int> q;  int f[N],mmin[N];  bool used[N];
void add(int a,int b)
{
    tail++;
    e[tail].from=a; e[tail].to=b;
    e[tail].next=first[a]; first[a]=tail;
    return;
}
void read()
{
    int i,a,b,c;
    scanf("%d%d",&n,&m);
    for (i=1;i<=n;i++)
        scanf("%d",&w[i]);
    for (i=1;i<=n;i++)
    { mmin[i]=w[i]; f[i]=INF; }
    for (i=1;i<=m;i++)
    {
        scanf("%d%d%d",&a,&b,&c);
        add(a,b);   if (c==2)  add(b,a);
    }
    return;
}
void spfa()
{
    int i,t;   bool flag;
    while(!q.empty()) q.pop();
    q.push(1); used[1]=1; f[1]=0;
    while(!q.empty())
    {
        i=q.front(); t=first[i];
        while(t!=0)
        {
            flag=0;
            if (f[e[t].to]<max(f[i],w[e[t].to]-mmin[i]))
            {
                f[e[t].to]=max(f[i],w[e[t].to]-mmin[i]);
                flag=1;
            }
            if (mmin[e[t].to]>mmin[i])
            {
                mmin[e[t].to]=mmin[i];
                flag=1;
            }
            if (flag&&!used[e[t].to])
            {
                q.push(e[t].to);
                used[e[t].to]=1;
            }
            t=e[t].next;
        }
        q.pop();  used[i]=0;
    }
    printf("%d",f[n]);
    return;
}
int main()
{
    freopen("trade.in","r",stdin);
    freopen("trade.out","w",stdout); 
    read();
    spfa();
    return 0;
}