NOIP2009 T3最优贸易

题目描述
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。

这题写了两种方法，调了挺久的。。第一种正反图双SPFA，注意我们只应用了SPFA的思想，改变了就入队的思想，正序来求出每个点最小买入价。反向来求出最大卖出价。然后差值最大就是答案。第二种方法是用tarjan缩点把图缩成DAG后进行拓扑排序，最后DP求解。这题的思想很妙，两份代码。

//正反图SPFA
#include<bits/stdc++.h>
using namespace std;

const int MAXN=5e5+2;
const int MINN=1e5+2;

int head[MINN],cnt=0,ru[MINN],chu[MINN],headf[MINN],cntf=0,w[MINN];
bool vis[MINN];
struct edge{
    int to,next;
}e[MAXN<<1],ef[MAXN<<1];

inline void add(int u,int v){e[++cnt]=(edge){v,head[u]},head[u]=cnt;}
inline void addf(int u,int v){ef[++cntf]=(edge){v,headf[u]},headf[u]=cntf;}

queue<int>q;

void spfaz(int x){//最小买入价格 
    memset(ru,0x7f,sizeof(ru));
    memset(vis,0,sizeof(vis));
    ru[x]=w[x];
    vis[x]=1;
    q.push(x);
    while(!q.empty()){
        int u=q.front();
        q.pop();vis[u]=0;
        if(w[u]<ru[u])ru[u]=w[u];
        for(int i=head[u];i;i=e[i].next){
            int v=e[i].to;
            if(ru[u]<ru[v]){
                ru[v]=ru[u];
                if(!vis[v]){
                    q.push(v);
                    vis[v]=1;
                }
            }
        }
    }
}

void spfaf(int x){
    memset(chu,0,sizeof(chu));
    memset(vis,0,sizeof(vis));
    chu[x]=w[x];
    vis[x]=1;
    q.push(x);
    while(!q.empty()){
        int u=q.front();
        q.pop();vis[u]=0;
        if(w[u]>chu[u])chu[u]=w[u];
        for(int i=headf[u];i;i=ef[i].next){
            int v=ef[i].to;
            if(chu[u]>chu[v]){
                chu[v]=chu[u];
                if(!vis[v]){
                    q.push(v);
                    vis[v]=1;
                }
            }
        }
    }
}

int n,m,x,y,z;

int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)scanf("%d",&w[i]);
    for(int i=1;i<=m;i++){
        scanf("%d%d%d",&x,&y,&z);
        if(z==1){
            add(x,y);
            addf(y,x);
        }
        else{
            add(x,y);add(y,x);
            addf(x,y);addf(y,x);
        }
    }
    spfaz(1);
    spfaf(n);
    int ans=0;
    for(int i=1;i<=n;i++){
        ans=max(ans,chu[i]-ru[i]);
    //  printf("%d %d ",chu[i],ru[i]);
    }
    printf("%d\n",ans);
    return 0;
}

//tarjan缩点DP
#include<bits/stdc++.h>
using namespace std;

const int MAXN=5e5+2;
const int MINN=1e5+2;

int n,m,w[MINN],x,y,z,dfn[MINN],low[MINN],head[MINN],cnt=0,cntcol=0,headcol[MAXN],rk[MINN];
int sta[MINN],belong1[MINN],num=0,weizhi=0,times=0,fa[MINN],inde[MINN],xv[MAXN],dq=0;
int max_belong[MINN],low_belong[MINN],dp[MINN],pre[MINN];
bool insta[MINN];

struct edge{
    int to,next;
}e[MAXN<<1],col[MAXN<<1];

inline void add(int u,int v){e[++cnt]=(edge){v,head[u]},head[u]=cnt;}
inline void addcol(int u,int v){col[++cntcol]=(edge){v,headcol[u]},headcol[u]=cntcol;}
//dp[n]=max(dp[n-1],dp[n],w[n]-w[1..n-1])    
void tarjan(int u){
    insta[u]=1;
    sta[++weizhi]=u;
    dfn[u]=low[u]=++times;
    for(int i=head[u];i;i=e[i].next){
        int v=e[i].to;
        if(!dfn[v]){
            tarjan(v);
            if(low[v]<low[u])low[u]=low[v];
        }
        else if(insta[v]&&dfn[v]<low[u])low[u]=dfn[v];
    }
    if(dfn[u]==low[u]){
        int tem=0;num++;
        while(tem!=u){
            tem=sta[weizhi--];
            insta[tem]=0;
            belong1[tem]=num;
            max_belong[num]=max(max_belong[num],w[tem]);
            low_belong[num]=min(low_belong[num],w[tem]);
        }
    }
}

void makeset(int n){for(int i=1;i<=n;i++)fa[i]=i;}

int find(int x){
    return x==fa[x]?x:fa[x]=find(fa[x]);
}

void unset(int x,int y){
    if((x=find(x))==(y=find(y)))return;
    if(rk[x]>rk[y])fa[y]=x;
    else {
        fa[x]=y;
        if(rk[x]==rk[y])rk[y]++;
    }
}//无视掉这个谜之并查集吧 我也不知道怎么想的

queue<int>q;

int main(){
    memset(low_belong,0x7f,sizeof(low_belong));
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)scanf("%d",&w[i]);
    for(int i=1;i<=m;i++){
        scanf("%d%d%d",&x,&y,&z);
        if(z==1){
            add(x,y);
        }
        else{
            add(x,y);add(y,x);
        }
    }
    for(int i=1;i<=n;i++){
        if(!dfn[i])tarjan(i);
    }
    for(int u=1;u<=n;u++){
        for(int i=head[u];i;i=e[i].next){
            int v=e[i].to;
            if(belong1[v]!=belong1[u]){
                addcol(belong1[u],belong1[v]);
        //      printf("b:%d %d\n",belong1[u],belong1[v]);
    //          printf("%d ",belong1[v]);
                inde[belong1[v]]++;
            }
        }
    }
//  for(int i=1;i<=n;i++ )printf("%d \n",belong1[i]);
    for(int i=1;i<=num;i++){
        if(!inde[i]){
            q.push(i);//printf("%d",i);
        }   
    }
    while(!q.empty()){
        int u=q.front();
        q.pop();
        xv[++dq]=u;
        for(int i=headcol[u];i;i=col[i].next){
            int v=col[i].to;
        //  printf("%d",v);
            inde[v]--;
            if(!inde[v]){
                q.push(v);
            }
        }
    }
//  for(int i=1;i<=dq;i++)printf("%d ",xv[i]);
    //dp[n]=max(dp[n-1],dp[n],w[n]-w[1..n-1])
//  printf("%d ",xv[1]);

    for(int i=1;i<=dq;i++){
        for(int j=headcol[xv[i]];j;j=col[j].next){//xv[i]是s 
        //  printf("%d ",low_belong[1]);
            int v=col[j].to;
        //  printf("hd:%d v:%d \n",xv[i],v);
            int minnn=0x7ffff;
            int tem2=0;
            //minnn=min(low_belong[v],minnn);
            low_belong[v]=min(low_belong[xv[i]],low_belong[v]);
            tem2=max(dp[xv[i]],max_belong[v]-low_belong[v]);
            dp[v]=max(tem2,dp[v]);
        }
    }
    printf("%d\n",dp[belong1[n]]);
    return 0;
}