NOIP2009T3 最优贸易(SPFA)

描述

C 国有 n 个大城市和 m 条道路,每条道路连接这 n 个城市中的某两个城市。任意两个

城市之间最多只有一条道路直接相连。这 m 条道路中有一部分为单向通行的道路,一部分

为双向通行的道路,双向通行的道路在统计条数时也计为 1 条。

C 国幅员辽阔,各地的资源分布情况各不相同,这就导致了同一种商品在不同城市的价

格不一定相同。但是,同一种商品在同一个城市的买入价和卖出价始终是相同的。

商人阿龙来到 C 国旅游。当他得知同一种商品在不同城市的价格可能会不同这一信息

之后,便决定在旅游的同时,利用商品在不同城市中的差价赚回一点旅费。设 C 国 n 个城

市的标号从 1~ n,阿龙决定从 1 号城市出发,并最终在 n 号城市结束自己的旅行。在旅游的

过程中,任何城市可以重复经过多次,但不要求经过所有 n 个城市。阿龙通过这样的贸易方

式赚取旅费:他会选择一个经过的城市买入他最喜欢的商品――水晶球,并在之后经过的另

一个城市卖出这个水晶球,用赚取的差价当做旅费。由于阿龙主要是来 C 国旅游,他决定

这个贸易只进行最多一次,当然,在赚不到差价的情况下他就无需进行贸易。

假设 C 国有 5 个大城市,城市的编号和道路连接情况如下图,单向箭头表示这条道路

为单向通行,双向箭头表示这条道路为双向通行。

20170621110109_23905

假设 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 个整数,表示最多能赚取的旅费。如果没有进行贸易,

样例输入[复制]
5 5 
4 3 5 6 1 
1 2 1 
1 4 1 
2 3 2 
3 5 1 
4 5 2
样例输出[复制]
5
提示

【数据范围】

输入数据保证 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。

 
 
 
这个题简直是难以想到spfa怎么写
先说思路
建正反图,对正图跑一个求最大值的spfa,对反图跑一个最小值的spfa
注意这里的最小值指的是从1号点到这里经过点的最小值
由于是只进行一次,那么又要到达n点
首先条件就是这个点要与终点联通
由于我来的时候肯定要购买最小的,所以肯定对于从1出发要寻找最小值
如果我们反向建图达到了两个目的,一个是证明一个点可以到n同时可以从1到达,第二同时找最大值
如果你直接最大最小值在正图上面跑就会有一个问题,如果这个点不能到达n就gg,所以必须要建立反图
然后跑完之后直接遍历寻找最大值就可以了
那么在写spfa的时候还会遇到一个巨大的问题,如果你只用dis来更新值就会出现双向边回不来的情况,这个时候优先与自己比较在与v比较,可以看代码理解
脑洞啊QAQ
代码:
 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<queue>
 5 #define N 1000006 
 6 using namespace std;
 7 struct node{
 8     int u,v;
 9 }e[N],d[N];
10 int first[N],nxt[N],ans,First[N],Nxt[N],cnt,tot,vis[N],n,m,dis[N],dis1[N],v[N];
11 void add(int u,int v){
12     e[++cnt].u=u;
13     e[cnt].v=v;
14     nxt[cnt]=first[u];
15     first[u]=cnt;
16 }
17 void add1(int u,int v){
18     d[++tot].u=u;
19     d[tot].v=v;
20     Nxt[tot]=First[u];
21     First[u]=tot;
22 }
23 void spfa1(){
24     queue<int>q;
25     memset(vis,0,sizeof vis);
26     for(int i=1;i<=n;i++)dis[i]=0x3f3f3f3f;
27     dis[1]=v[1];
28     vis[1]=1;
29     q.push(1);
30     while(!q.empty()){
31         int u=q.front();
32         q.pop();
33         vis[u]=0;
34         for(int i=first[u];i;i=nxt[i]){
35             int v1=e[i].v;
36             if(dis[v1]>dis[u]){
37                 dis[v1]=min(dis[u],v[v1]);
38                 if(!vis[v1]){
39                     vis[v1]=1;
40                     q.push(v1);
41                 }
42             }
43         }
44     }
45 }
46 void spfa2(){
47     queue<int>q;
48     memset(vis,0,sizeof(vis));
49     for(int i=1;i<=n;i++)dis1[i]=0;
50     dis1[n]=v[n];
51     q.push(n);
52     vis[n]=1;
53     while(!q.empty()){
54         int u=q.front();
55         q.pop();
56         vis[u]=0;
57         for(int i=First[u];i;i=Nxt[i]){
58             int v1=d[i].v;
59             if(dis1[v1]<dis1[u]){
60                 dis1[v1]=max(dis1[u],v[v1]);
61                 ans=max(ans,dis1[v1]-dis[v1]);
62                 if(!vis[v1]){
63                     vis[v1]=1;
64                     q.push(v1);
65                 }
66             }
67         }
68     }
69 }
70 int main(){
71     cin>>n>>m;
72     for(int i=1;i<=n;i++)cin>>v[i];
73     for(int i=1;i<=m;i++){
74         int x,y,z;
75         cin>>x>>y>>z;
76         if(z==1){
77             add1(y,x);
78             add(x,y);
79         }    
80         else {
81             add(x,y),add(y,x);
82             add1(y,x),add1(x,y);
83         }
84     }
85     spfa1();
86     spfa2();
87     cout<<ans;
88 } 

over qwq

原文地址:https://www.cnblogs.com/saionjisekai/p/9569680.html