最优贸易-tarjan+拓扑更新

描述
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之间的双向道路。

https://vijos.org/p/1754

tip:

首先就是可以反复走,那么一个连用分量里的点,一定可以取到最大值和最小值,里面走成什么样子可以不关心,于是,我们可以tarjan跑一发,然后缩点,这样只要知道一个圈里面的最大和最小值就可以了。
然后就要知道,所有点从1辗转反侧到达它的途中最小的交易金额,然后求出从每个点到n的最大金额,这样两个数值相减,取max就是ans
怎么求呢。。。
还是要用拓扑排序,这个时候呀,正向走一遍缩点之后的图,更新每个能从1来的点的最小值(记录能不能从1这个点来,拓扑排序过程中只要遇到过1点,就算可以,标记一下就可以)
第二遍拓扑排序,反向走一遍缩点之后的图,相当于把原图中的边反向,更新每个点从n来的最大值(记录和上面是一样的)
最后遍历所有点,当一个点再第一遍拓扑的时候从1来,第二遍拓扑的时候从n来的话,就说明这个点是1-n的某条路径上的某一个点,就满足了题目条件,于是更新ans。。。最后ans取这些点(max-min)中的最大值。。。。

#include <cstdio>
#include <iostream>
#include <queue>
#include <cstring>
using namespace std;
const int maxn = 5e5+10;
queue<int> qf,qb;
int n,m,tot,u,v,k,tt,chuo,top,totf,totb;
int headf[maxn],headb[maxn],inf[maxn],inb[maxn];
int a[maxn],head[maxn],ma[maxn],mi[maxn],finf[maxn],finb[maxn];
int low[maxn],dfn[maxn],belong[maxn],instk[maxn],stk[maxn];
struct node{
    int v,next;
}edges[maxn];

struct forwd{
    int v,next;
}egf[maxn];

struct bak{
    int v,next;
}egb[maxn];

void add(int u,int v){
    edges[tot].v = v;edges[tot].next = head[u];head[u] = tot++;
}
void addf(int u,int v){
    egf[totf].v = v;egf[totf].next = headf[u];headf[u] = totf++;
}
void addb(int u,int v){
    egb[totb].v = v;egb[totb].next = headb[u];headb[u] = totb++;
}

void init(){
    for(int i = 1; i <= n ; i++){
        scanf("%d",&a[i]);
    }
    for(int i = 1; i <= m ; i++){
        scanf("%d%d%d",&u,&v,&k);
        if(k == 2   )   add(v,u);
        add(u,v);
    }
}

void cler(){
    tot=0;
    for(int i = 1; i <= n ; i++){
        head[i] = headf[i] = headb[i] = -1;
        mi[i] = (1<<30);
        ma[i] = -(1<<30);
    }
}

void tsortf(){
    for(int i = 1 ; i <= tt ;i++){
        if(inf[i] == 0){
            qf.push(i);
        }
    }
    while(!qf.empty()){
        int tmp = qf.front();
        qf.pop();
        for(int k = headf[tmp];k != -1 ; k = egf[k].next){
            int to = egf[k].v;
            inf[to]--;
            if(finf[tmp]){
                mi[to] = min(mi[to],mi[tmp]);
                finf[to] = 1;
                if(inf[to] == 0){
                    qf.push(to);
                }
            }
        }
    }
}

void tsortb(){
    for(int i = 1 ; i <= tt ;i++){
        if(inb[i] == 0){
            qb.push(i);
        }
    }
    while(!qb.empty()){
        int tmp = qb.front();
        qb.pop();
        for(int k = headb[tmp];k != -1 ; k = egb[k].next){
            int to = egb[k].v;
            inb[to]--;
            if(finb[tmp]){
                ma[to] = max(ma[to],ma[tmp]);
                finb[to] = 1;
                if(inb[to] == 0){
                    qb.push(to);
                }
            }
        }
    }
}

void tarjan(int u){
    dfn[u] = low[u] = ++chuo;
    instk[u] = 1;stk[top++] = u;
    for(int k = head[u] ;k != -1; k = edges[k].next){
        if(!dfn[edges[k].v]){
            tarjan(edges[k].v);
            low[u] = min(low[u],low[edges[k].v]);
        }
        else if(instk[edges[k].v]){
            low[u] = min(low[u],dfn[edges[k].v]);
        }
    }
    if(low[u] == dfn[u]){
        ++tt;
        int t;
        do{
            t = stk[--top];
            instk[t] = 0;
            belong[t] = tt;
            mi[tt] = min(mi[tt],a[t]);
            ma[tt] = max(ma[tt],a[t]);
            //printf("t = %d  mi[%d] = %d, ma[%d] = %d\n",t,tt,mi[tt],tt,ma[tt]);
        }while(t != u);
    }
}

void sov(){
    for(int i = 1; i <= n ; i++){
        if(!dfn[i]){
            tarjan(i);
        }
    }
    for(int i = 1; i <= n ; i++){
        for(int k = head[i] ; k != -1; k = edges[k].next){
            if(belong[i]!=belong[edges[k].v]){
                addf(belong[i],belong[edges[k].v]);
                inf[belong[edges[k].v]]++;
                addb(belong[edges[k].v],belong[i]);
                inb[belong[i]]++;
            }
        }
    }
    //printf("belong[1] = %d  inf = %d \n",belong[1],inf[belong[1]]);
    finf[belong[1]] = finb[belong[n]] = 1;
    tsortf();
    tsortb();
    int ans = 0;
    for(int i = 1; i <= n; i++){
        if(!finf[i] || !finb[i])    continue;
        ans = max(ans,ma[i]-mi[i]);
    }
    printf("%d\n",ans);
}

int main(){
    scanf("%d%d",&n,&m);
    cler();
    init();
    sov();
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值