[NOIP2009 提高组] 最优贸易(C++,tarjan,topo,DP)
创始人
2024-06-01 16:29:44
0

题目描述

$C 国有国有国有 n 个大城市和个大城市和个大城市和 m$ 条道路,每条道路连接这 nnn个城市中的某两个城市。任意两个城市之间最多只有一条道路直接相连。这 mmm 条道路中有一部分为单向通行的道路,一部分为双向通行的道路,双向通行的道路在统计条数时也计为 $1 $条。

$C $国幅员辽阔,各地的资源分布情况各不相同,这就导致了同一种商品在不同城市的价格不一定相同。但是,同一种商品在同一个城市的买入价和卖出价始终是相同的。

商人阿龙来到 CCC 国旅游。当他得知同一种商品在不同城市的价格可能会不同这一信息之后,便决定在旅游的同时,利用商品在不同城市中的差价赚回一点旅费。设 CCC 国 n 个城市的标号从 1∼n1\sim n1∼n,阿龙决定从 $1 $号城市出发,并最终在 nnn 号城市结束自己的旅行。在旅游的过程中,任何城市可以重复经过多次,但不要求经过所有 nnn 个城市。阿龙通过这样的贸易方式赚取旅费:他会选择一个经过的城市买入他最喜欢的商品――水晶球,并在之后经过的另一个城市卖出这个水晶球,用赚取的差价当做旅费。由于阿龙主要是来 CCC 国旅游,他决定这个贸易只进行最多一次,当然,在赚不到差价的情况下他就无需进行贸易。

假设 $C $国有 555个大城市,城市的编号和道路连接情况如下图,单向箭头表示这条道路为单向通行,双向箭头表示这条道路为双向通行。

假设 1n1~n1 n 号城市的水晶球价格分别为 4,3,5,6,14,3,5,6,14,3,5,6,1。

阿龙可以选择如下一条线路:111->222->333->555,并在 $2 号城市以号城市以号城市以 3$ 的价格买入水晶球,在 333号城市以$ 5 $的价格卖出水晶球,赚取的旅费数为 2。

阿龙也可以选择如下一条线路$ 1$->444->555->444->555,并在第$1 次到达次到达次到达 5$ 号城市时以 $1 $的价格买入水晶球,在第 222 次到达$ 4$ 号城市时以$ 6$ 的价格卖出水晶球,赚取的旅费数为$ 5$。

现在给出 $n 个城市的水晶球价格,个城市的水晶球价格,个城市的水晶球价格,m$ 条道路的信息(每条道路所连接的两个城市的编号以及该条道路的通行情况)。请你告诉阿龙,他最多能赚取多少旅费。

输入格式

第一行包含 222 个正整数$ n $和 mmm,中间用一个空格隔开,分别表示城市的数目和道路的数目。

第二行 n 个正整数,每两个整数之间用一个空格隔开,按标号顺序分别表示这 n 个城市的商品价格。

接下来 mmm 行,每行有$ 3 个正整数个正整数个正整数x,y,z$,每两个整数之间用一个空格隔开。如果 z=1z=1z=1,表示这条道路是城市$ x 到城市到城市到城市 y 之间的单向道路;如果之间的单向道路;如果之间的单向道路;如果 z=2$,表示这条道路为城市 $x 和城市和城市和城市y $之间的双向道路。

输出格式

一 个整数,表示最多能赚取的旅费。如果没有进行贸易,则输出 000。

样例 #1

样例输入 #1

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

样例输出 #1

5

提示

【数据范围】

输入数据保证 111 号城市可以到达$ n $号城市。

对于 10%的数据,1≤n≤61≤n≤61≤n≤6。

对于 30%的数据,1≤n≤1001≤n≤1001≤n≤100。

对于 50%的数据,不存在一条旅游路线,可以从一个城市出发,再回到这个城市。

对于 100%的数据,1≤n≤1000001≤n≤1000001≤n≤100000,1≤m≤5000001≤m≤5000001≤m≤500000,1≤x1≤x1≤x,y≤ny≤ny≤n,1≤z≤21≤z≤21≤z≤2,1≤1≤1≤各城市

水晶球价格≤100≤100≤100。

NOIP 2009 提高组 第三题

解题思路:

不需要过多考虑,本题唯一的解法是遍历每一种情况,找出最大差价

关于图的搜索,如果一个点能被重复经过会使问题复杂化,所以我们首先缩点,简化问题

然后考虑如何得出答案

缩点需要保留的信息是环中的最大出手价格和最小入手价格

DFS从城市1开始暴力搜索记录当前行过点中最小的买入价格,每到一个点,计算差价

然后与当前的最大差价比较(差价最小为0,代表不需要贸易)

需要注意的是,我们的目标除了找出最大差价,还需要到达城市n,否则答案没有意义

思路非常简洁,当然判题的回复也非常简洁

TLE(QAQ)

那么就需要考虑如何优化

超时的原因一定是因为暴力搜索重复搜索了太多次,所以关键在于去重

所以想到动态规划

在这里插入图片描述

考虑以上的图,动态规划关键在于从子情况的解推导出下一个情况的解

而我们从111开始DFS会导致直接到达555,我们能够在444处“等待”,但这显然不可能用DFS实现

因为对于444号节点来说,入边是不可见的,不可能判断出是否需要“等待”

这时就需要用到我们的工具,拓扑排序

拓扑排序能够实现在到达一个节点之间已经遍历过之前的所有节点

最后需要提示的一点

在这里插入图片描述

试着考虑图中矩形节点(矩形节点到节点1只有单向边)

这两个节点会影响答案的正确性,因为在拓扑排序 + DP的时候,我们会根据它们的信息更新节点1

解决方法就是在建新图的时候注意将其去除即可(具体如何去除见代码)

Code:

#include 
#include 
#include 
using namespace std;
const int max_n = 1e5;
const int max_m = 5e5;
const int max_value = 100;
const int NaN = 0x3F3F3F3F;int n, m, x, y, z;
struct edge { int v, next; }edges[max_m * 2];
int tot = -1;
int head[max_n + 1];
int value[max_n + 1];
//tarjan
int timeclock = 0, dfn[max_n + 1], low[max_n + 1];
int rsp = -1, stack[max_n + 1], instack[max_n + 1];
int cnt = 0, belong[max_n + 1];
//re_init
edge new_edges[max_m * 2];
int new_head[max_n + 1];
int new_tot = -1;
int invalue[max_n + 1];
int outvalue[max_n + 1];
//topo
int in[max_n + 1];
queuetopoq;
//dp
int ans[max_n + 1];void add_edge(int u, int v) {edges[++tot] = { v,head[u] }; head[u] = tot;
}void add_new_edge(int u, int v) {new_edges[++new_tot] = { v,new_head[u] }; new_head[u] = new_tot;
}void tarjan(int x) {dfn[x] = low[x] = ++timeclock;stack[++rsp] = x;instack[x] = 1;for (int i = head[x]; i != -1; i = edges[i].next) {int v = edges[i].v;if (!dfn[v]) {tarjan(v);low[x] = min(low[x], low[v]);}else if (instack[v]) {low[x] = min(low[x], low[v]);}}if (dfn[x] == low[x]) {cnt++;while (stack[rsp + 1] != x) {int node = stack[rsp--];belong[node] = cnt;invalue[cnt] = min(invalue[cnt], value[node]);outvalue[cnt] = max(outvalue[cnt], value[node]);instack[node] = 0;}}
}void re_init() {for (int i = 1; i <= n; i++) {for (int j = head[i]; j != -1; j = edges[j].next) {int v = edges[j].v;if (belong[i] && belong[v] && belong[i] != belong[v]) {//由于只从1号开始tarjan,从1号不可到达的节点没有归属任何新节点,其belong[x] == 0add_new_edge(belong[i], belong[v]);in[belong[v]]++;}}}
}void topo() {for (int i = 1; i <= cnt; i++) {if (in[i] == 0) {topoq.push(i);}}while (!topoq.empty()) {int node = topoq.front();topoq.pop();ans[node] = max(ans[node], outvalue[node] - invalue[node]);for (int i = new_head[node]; i != -1; i = new_edges[i].next) {int v = new_edges[i].v;ans[v] = max(ans[v], ans[node]);//DPinvalue[v] = min(invalue[v], invalue[node]);//DPin[v]--;if (!in[v]) topoq.push(v);}}
}int main() {memset(head + 1, -1, sizeof(int) * max_n);memset(new_head + 1, -1, sizeof(int) * max_n);memset(invalue + 1, 0x3F, sizeof(int) * max_n);cin >> n >> m;for (int i = 1; i <= n; i++) cin >> value[i];for (int i = 0; i < m; i++) {cin >> x >> y >> z;add_edge(x, y);if (z == 2) add_edge(y, x);}tarjan(1);//只从1号开始tarjanre_init();topo();cout << ans[belong[n]];return 0;
}

相关内容

热门资讯

喜欢穿一身黑的男生性格(喜欢穿... 今天百科达人给各位分享喜欢穿一身黑的男生性格的知识,其中也会对喜欢穿一身黑衣服的男人人好相处吗进行解...
发春是什么意思(思春和发春是什... 本篇文章极速百科给大家谈谈发春是什么意思,以及思春和发春是什么意思对应的知识点,希望对各位有所帮助,...
网络用语zl是什么意思(zl是... 今天给各位分享网络用语zl是什么意思的知识,其中也会对zl是啥意思是什么网络用语进行解释,如果能碰巧...
为什么酷狗音乐自己唱的歌不能下... 本篇文章极速百科小编给大家谈谈为什么酷狗音乐自己唱的歌不能下载到本地?,以及为什么酷狗下载的歌曲不是...
华为下载未安装的文件去哪找(华... 今天百科达人给各位分享华为下载未安装的文件去哪找的知识,其中也会对华为下载未安装的文件去哪找到进行解...
怎么往应用助手里添加应用(应用... 今天百科达人给各位分享怎么往应用助手里添加应用的知识,其中也会对应用助手怎么添加微信进行解释,如果能...
家里可以做假山养金鱼吗(假山能... 今天百科达人给各位分享家里可以做假山养金鱼吗的知识,其中也会对假山能放鱼缸里吗进行解释,如果能碰巧解...
一帆风顺二龙腾飞三阳开泰祝福语... 本篇文章极速百科给大家谈谈一帆风顺二龙腾飞三阳开泰祝福语,以及一帆风顺二龙腾飞三阳开泰祝福语结婚对应...
四分五裂是什么生肖什么动物(四... 本篇文章极速百科小编给大家谈谈四分五裂是什么生肖什么动物,以及四分五裂打一生肖是什么对应的知识点,希...
美团联名卡审核成功待激活(美团... 今天百科达人给各位分享美团联名卡审核成功待激活的知识,其中也会对美团联名卡审核未通过进行解释,如果能...