NOIP2009 最优贸易
3. 最优贸易
(trade.pas/c/cpp)
【问题描述】
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 条道路的信息(每条道路所连接的两个城市的编号以及该条道路的通行情况)。请你告诉阿龙,他 多能赚取多少旅费。
【输入】
输入文件名为 trade.in。
第一行包含 2 个正整数 n 和 m,中间用一个空格隔开,分别表示城市的数目和道路的数目。
第二行 n 个正整数,每两个整数之间用一个空格隔开,按标号顺序分别表示这 n 个城市的商品价格。
接下来 m 行,每行有 3 个正整数,x,y,z,每两个整数之间用一个空格隔开。如果 z=1,表示这条道路是城市 x 到城市 y 之间的单向道路;如果 z=2,表示这条道路为城市 x 和城市 y 之间的双向道路。
【输出】
输出文件 trade.out 共 1 行,包含 1 个整数,表示 多能赚取的旅费。如果没有进行贸易,则输出 0。
【输入输出样例】
trade.in |
trade.out |
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。
【思路】
本题可以概括为求1-n的一条路,使得路上的max-min最大,但max必须在min之后。
刚开始用dfs解,忽略了max与min的先后关系。而又因为本题可以出现环,所以不能用dfs。
本题可以用两边SPFA完成(较dfs而言SPFA是更新),第一次计算每个节点到1路径上的min,第二次计算每个节点到n路径上的max。ans=max{max[i]-min[i]}
然而还有更优的算法,只进行一次SPFA,维护_min[v]代表到v包含v的路上的最小值,维护f[v]代表到v包含v的路上的p-min所得的最大值。注意更新条件。
【代码1】
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<fstream>
#include<vector>
#include<queue>
#define FOR(a,b,c) for(int a=(b);a<=(c);a++)
using namespace std; const int maxn = +;
const int INF=<<;
vector<int> Gto[maxn],Gbac[maxn];
int n,m,p[maxn],_min[maxn],_max[maxn]; void SPFA_min(int s,int* d) {
int inq[maxn]; fill(inq,inq+n+,);
queue<int> Q;
fill(d,d+n+,INF);
d[s]=p[s]; inq[s]=;
Q.push(s);
while(!Q.empty()) {
int u=Q.front();Q.pop(); inq[u]=;
for(int i=;i<Gto[u].size();i++) {
int v=Gto[u][i];
if(min(d[u],p[v])<d[v]) {
d[v]=min(d[u],p[v]);
if(!inq[v]) {
inq[v]=; Q.push(v);
}
}
}
}
}
void SPFA_max(int s,int* d){
int inq[maxn]; fill(inq,inq+n+,);
queue<int> Q;
d[s]=p[s]; inq[s]=;
Q.push(s);
while(!Q.empty()) {
int u=Q.front();Q.pop(); inq[u]=;
for(int i=;i<Gbac[u].size();i++) {
int v=Gbac[u][i];
if(max(d[u],p[v])>d[v]) {
d[v]=max(d[u],p[v]);
if(!inq[v]) {
inq[v]=; Q.push(v);
}
}
}
}
}
inline void AddEdge(int u,int v) {
Gto[u].push_back(v);
Gbac[v].push_back(u);
}
int main() {
scanf("%d%d",&n,&m);
FOR(i,,n) scanf("%d",&p[i]);
FOR(i,,m) {
int u,v,z;
scanf("%d%d%d",&u,&v,&z);
AddEdge(u,v);
if(z==) AddEdge(v,u);
}
SPFA_min(,_min);
SPFA_max(n,_max);
int ans=;
FOR(i,,n) if(_max[i]&&_min[i]<INF) ans=max(ans,_max[i]-_min[i]);
cout<<ans;
return ;
}
【代码2】
#include<cstdio>
#include<vector>
#include<queue>
#define FOR(a,b,c) for(int a=(b);a<=(c);a++)
using namespace std; const int maxn = +;
const int INF=<<;
vector<int> G[maxn];
int n,m,p[maxn],f[maxn],_min[maxn]; void SPFA(int s) {
int inq[maxn]; fill(inq,inq+n+,);
queue<int> Q;
fill(_min,_min+n+,INF);
_min[s]=p[s]; inq[s]=;
Q.push(s);
while(!Q.empty()) {
int u=Q.front();Q.pop(); inq[u]=;
for(int i=;i<G[u].size();i++) {
int v=G[u][i];
if(min(_min[u],p[v])<_min[v] ||f[v]<p[v]-_min[v]|| f[v]<f[u] ) {
//注意判断条件 只要可以更新v结点
_min[v]=min(_min[v],min(_min[u],p[v]));
f[v]=max(f[v],p[v]-_min[v]);
f[v]=max(f[v],f[u]);
if(!inq[v]) {
inq[v]=; Q.push(v);
}
}
}
}
} int main() {
freopen("trade.in","r",stdin);
freopen("trade.out","w",stdout);
scanf("%d%d",&n,&m);
FOR(i,,n) scanf("%d",&p[i]);
FOR(i,,m) {
int u,v,z;
scanf("%d%d%d",&u,&v,&z);
G[u].push_back(v);
if(z==) G[v].push_back(u);
}
SPFA();
printf("%d",f[n]);
return ;
}
最新文章
- [转] How to debug a ARM Cortex-M hard fault exception
- 深入了解 Oracle Flex ASM 及其优点
- 深入Mysql 导入导出
- session进程和服务
- [windows phone] 教你如何使地图动画缩放
- 第48篇 说说.Net与Java中的字符串
- php计算多个集合的笛卡尔积实例详解
- lufylegend:Lbuttion等UI组件
- STL源码分析-deque
- SQL Server更改排序规则的实现过程
- ntp时钟服务器
- [技巧篇]16.MyEclipse2014安装SVN插件,在线安装
- 【洛谷 P4168】[Violet]蒲公英(分块)
- 线程协作-Semaphore并发限制
- java文本文件读写
- luogu1556 幸福的路
- UESTC 1437
- Operating System-Thread(1)What and Why Thread &;&;进程和线程的对比
- layui 数据table隐藏表头
- 条款16:成对使用new和delete时要采取相同形式