1. 程式人生 > >【NOIP2009】最優貿易

【NOIP2009】最優貿易

!= 連通 最大 sub 旅行 ret 空格 十分 100%

Description

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條道路的信息(每條道路所連接的兩個城市的編號以及該條道路的通行情況) 。請你告訴阿龍,他最多能賺取多少旅費。

Input

第一行包含 2 個正整數 n 和 m,中間用一個空格隔開,分別表示城市的數目和道路的數目。
第二行 n 個正整數,每兩個整數之間用一個空格隔開,按標號順序分別表示這 n 個城市的商品價格。
接下來 m行, 每行有 3 個正整數, x, y, z, 每兩個整數之間用一個空格隔開。 如果 z=1,表示這條道路是城市 x到城市 y之間的單向道路;如果 z=2,表示這條道路為城市 x 和城市y之間的雙向道路。

Output

共1 行, 包含 1 個整數, 表示最多能賺取的旅費。 如果沒有進行貿易,則輸出 0。

Sample Input

5 5
4 3 5 6 1
1 2 1
1 4 1
2 3 2
3 5 1
4 5 2

Sample Output

5

Hint

數據範圍:
輸入數據保證 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。

題解:

  還是一道十分有趣的題目,首先,我們看一下百分之50的數據,保證圖是一個DAG,那麽我們就可以考慮用動態規劃來解決,設dp[i]表示從i點出發,到n之間最大的點權是多大(保證只有能到n的點才能更新答案),那麽顯然dp[i]=max(dp[to],dp[i]),如果i節點能到n節點,dp[i]=max(dp[i],quan[i]),那麽答案就是枚舉一個買的點,那麽答案就是max(dp[i]-quan[i])。

  為什麽要先知道50分的呢?(因為我只想出了50 ^_^),看了題解,50分是提醒我們想到正解的,我們考慮,因為如果我們強聯通分量縮點,那麽圖又變成了一個DAG,我們只要維護每個縮點的最大點權和最小點權,又可以用上述dp方法了。

  正確性是因為顯然如果在一個聯通分量裏面,那麽兩點之間一定都可以到達,並且從一個強連通分量返回另一個強連通分量是不可能的。

代碼:

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <iostream>
#include <queue>
#include <stack>
#define MAXN 500100
#define MAXN2 200100
using namespace std;
struct edge{
    int first;
    int next;
    int to;
}a[MAXN*2];
int dfn[MAXN],low[MAXN],in[MAXN],fa[MAXN],maxx[MAXN],minn[MAXN];
int f[MAXN2],quan[MAXN2],x[MAXN],y[MAXN],z[MAXN];
int num1=0,num2=0,num3=0,ans=0,n,m;
stack<int> s;
void addedge(int from,int to){
    a[++num1].to=to;
    a[num1].next=a[from].first;
    a[from].first=num1;
}

void tarjian(int now){
    dfn[now]=low[now]=num2++;s.push(now);in[now]=1;
    for(int i=a[now].first;i;i=a[i].next){
        int to=a[i].to;
        if(!dfn[to]){
            tarjian(to);
            low[now]=min(low[now],low[to]);
        }
        else if(in[to]) low[now]=min(low[now],dfn[to]);
    }
    if(low[now]==dfn[now]){
        num3++;
        int u=-1;
        while(u!=now){
            u=s.top();s.pop();
            maxx[num3]=max(maxx[num3],quan[u]);
            minn[num3]=min(minn[num3],quan[u]);
            fa[u]=num3;
            in[u]=0;
        }
    }
}

void DP(int now){
    dfn[now]=1;
    if(now==fa[n]) f[now]=maxx[fa[n]];
    for(int i=a[now].first;i;i=a[i].next){
        int to=a[i].to;
        if(!dfn[to]) DP(to);
        f[now]=max(f[to],f[now]);
    }
    if(f[now]){
        f[now]=max(f[now],maxx[now]);
    }
    ans=max(ans,f[now]-minn[now]);
}

void make(){
    for(int i=1;i<=m;i++){
        if(fa[x[i]]!=fa[y[i]]) addedge(fa[x[i]],fa[y[i]]);
    }
}

int main()
{
    scanf("%d%d",&n,&m);
    memset(minn,127,sizeof(minn));
    for(int i=1;i<=n;i++) scanf("%d",&quan[i]);
    for(int i=1;i<=m;i++){
        scanf("%d%d%d",&x[i],&y[i],&z[i]);
        if(z[i]==2)addedge(x[i],y[i]),addedge(y[i],x[i]);
        else addedge(x[i],y[i]);
    }
    while(!s.empty()) s.pop();
    tarjian(1);
    memset(a,0,sizeof(a));num1=0;
    memset(dfn,0,sizeof(dfn));
    make();
    DP(fa[1]);
    printf("%d\n",ans);
    return 0;
}

【NOIP2009】最優貿易