https://www.luogu.org/problemnew/show/P1073

C国有 n n个大城市和 mm 条道路,每条道路连接这 nn个城市中的某两个城市。任意两个城市之间最多只有一条道路直接相连。这 mm 条道路中有一部分为单向通行的道路,一部分为双向通行的道路,双向通行的道路在统计条数时也计为 1 1条。C C国幅员辽阔,各地的资源分布情况各不相同,这就导致了同一种商品在不同城市的价格不一定相同。但是,同一种商品在同一个城市的买入价和卖出价始终是相同的。商人阿龙来到 CC 国旅游。当他得知同一种商品在不同城市的价格可能会不同这一信息之后,便决定在旅游的同时,利用商品在不同城市中的差价赚回一点旅费。设 CC 国 n 个城市的标号从 1~ n1 n,阿龙决定从 1 1号城市出发,并最终在 nn 号城市结束自己的旅行。在旅游的过程中,任何城市可以重复经过多次,但不要求经过所有 nn 个城市。阿龙通过这样的贸易方式赚取旅费:他会选择一个经过的城市买入他最喜欢的商品――水晶球,并在之后经过的另一个城市卖出这个水晶球,用赚取的差价当做旅费。由于阿龙主要是来 CC 国旅游,他决定这个贸易只进行最多一次,当然,在赚不到差价的情况下他就无需进行贸易。假设 C C国有 55个大城市,城市的编号和道路连接情况如下图,单向箭头表示这条道路为单向通行,双向箭头表示这条道路为双向通行。假设 1~n1 n 号城市的水晶球价格分别为 4,3,5,6,14,3,5,6,1。阿龙可以选择如下一条线路:11->22->33->55,并在 2 2号城市以 33 的价格买入水晶球,在 33号城市以 5 5的价格卖出水晶球,赚取的旅费数为 2。阿龙也可以选择如下一条线路 11->44->55->44->55,并在第1 1次到达 55 号城市时以 1 1的价格买入水晶球,在第 22 次到达 44 号城市时以 66 的价格卖出水晶球,赚取的旅费数为 55。现在给出 n n个城市的水晶球价格,mm 条道路的信息(每条道路所连接的两个城市的编号以及该条道路的通行情况)。请你告诉阿龙,他最多能赚取多少旅费。

题意

第一眼看觉得是先缩点然后DAG图上跑一边拓扑排序,除了敲起来手有点酸之外没有什么难度。就直接过了。

#include <map>
#include <set>
#include <ctime>
#include <cmath>
#include <queue>
#include <stack>
#include <vector>
#include <string>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <sstream>
#include <iostream>
#include <algorithm>
#include <functional>
using namespace std;
#define For(i, x, y) for(int i=x;i<=y;i++)
#define _For(i, x, y) for(int i=x;i>=y;i--)
#define Mem(f, x) memset(f,x,sizeof(f))
#define Sca(x) scanf("%d", &x)
#define Sca2(x,y) scanf("%d%d",&x,&y)
#define Sca3(x,y,z) scanf("%d%d%d",&x,&y,&z)
#define Scl(x) scanf("%lld",&x);
#define Pri(x) printf("%d\n", x)
#define Prl(x) printf("%lld\n",x);
#define CLR(u) for(int i=0;i<=N;i++)u[i].clear();
#define LL long long
#define ULL unsigned long long
#define mp make_pair
#define PII pair<int,int>
#define PIL pair<int,long long>
#define PLL pair<long long,long long>
#define pb push_back
#define fi first
#define se second
typedef vector<int> VI;
int read(){int x = 0,f = 1;char c = getchar();while (c<'0' || c>'9'){if (c == '-') f = -1;c = getchar();}
while (c >= '0'&&c <= '9'){x = x * 10 + c - '0';c = getchar();}return x*f;}
const double eps = 1e-9;
const int maxn = 1e5 + 10;
const int maxm = 5e5 + 10;
const int INF = 0x3f3f3f3f;
const int mod = 1e9 + 7;
int N,M,K;
int val[maxn];
struct Edge{int to,next;
}edge[maxm * 2];
int head[maxn],tot;
void init(){for(int i = 0 ; i <= N + 1; i ++) head[i] = -1;tot = 0;
}
void add(int u,int v){edge[tot].to = v;edge[tot].next = head[u];head[u] = tot++;
}
int Low[maxn],dfn[maxn],Stack[maxn],Belong[maxn],num[maxn];
int Index,top,scc;
bool Instack[maxn];
void Tarjan(int u){int v;Low[u] = dfn[u] = ++Index;Stack[top++] = u;Instack[u] = true;for(int i = head[u]; ~i; i = edge[i].next){v = edge[i].to;if(!dfn[v]){Tarjan(v);if(Low[u] > Low[v]) Low[u] = Low[v];}else if(Instack[v] && Low[u] > dfn[v]){Low[u] = dfn[v];}}if(Low[u] == dfn[u]){scc++;do{v = Stack[--top];Instack[v] = false;Belong[v] = scc;num[scc]++;}while(v != u);}
}
int MAX[maxn],MIN[maxn],ind[maxn],ans[maxn];
vector<int>P[maxn];
int main(){Sca2(N,M); init();for(int i = 1; i <= N ; i ++) Sca(val[i]);for(int i = 1; i <= M ; i ++){int u,v,w; Sca3(u,v,w);add(u,v);if(w == 2) add(v,u);}for(int i = 1; i <= N ; i ++) if(!dfn[i]) Tarjan(i);for(int i = 1; i <= scc; i ++){MAX[i] = -INF; MIN[i] = INF;}for(int i = 1; i <= N ; i ++){int u = Belong[i];MAX[u] = max(MAX[u],val[i]);MIN[u] = min(MIN[u],val[i]);for(int j = head[i]; ~j; j = edge[j].next){int v = Belong[edge[j].to];if(u == v) continue;P[u].push_back(v);ind[v]++;    }}queue<int>Q;for(int i = 1; i <= scc; i ++) if(!ind[i]) Q.push(i);while(!Q.empty()){int u = Q.front(); Q.pop();ans[u] = max(ans[u],MAX[u] - MIN[u]);for(int i = 0; i < P[u].size(); i ++){int v = P[u][i];ind[v]--;MIN[v] = min(MIN[v],MIN[u]);ans[v] = max(ans[v],ans[u]);if(!ind[v]) Q.push(v);}}Pri(ans[Belong[N]]);return 0;
}

缩点 + 拓扑排序

后来发现还有一种更加新奇的解法,构造分层图,将原图变为一个带权图,权的意思是商人走这条路要得到的钱(负数代表花费),一开始的图所有边权都是0,表示这个商人什么也不买,走来走去不亏也不赚。然后我们对于每个点v,构造一条边通向v + N,边权为-val[v],表示商人在这个点买了一个水晶球,1 + N ~ N + N就是出现的第二个图,对于第二图,也像第一条路一样构造出相同的边权为0的边,表示商人买了水晶球之后依然可以不买也不卖的在这个图上走来走去,同理,我们再v + N ~ v + N + N构造一条边权为val[v]的边,表示商人买完之后在这个点又卖出去了,同第二个图一样,构造出第三个图,最后题意要求到达终点,所以可取的点只有N和3 * N,即第一个图和第三个图的终点,跑一边最短长路即可。

由于边权有正有负的,不能用Dijkstra跑,无向图也不可以拓扑排序dp,所以只能SPFA(大概是这个做法的唯一缺点)

#include <map>
#include <set>
#include <ctime>
#include <cmath>
#include <queue>
#include <stack>
#include <vector>
#include <string>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <sstream>
#include <iostream>
#include <algorithm>
#include <functional>
using namespace std;
#define For(i, x, y) for(int i=x;i<=y;i++)
#define _For(i, x, y) for(int i=x;i>=y;i--)
#define Mem(f, x) memset(f,x,sizeof(f))
#define Sca(x) scanf("%d", &x)
#define Sca2(x,y) scanf("%d%d",&x,&y)
#define Sca3(x,y,z) scanf("%d%d%d",&x,&y,&z)
#define Scl(x) scanf("%lld",&x);
#define Pri(x) printf("%d\n", x)
#define Prl(x) printf("%lld\n",x);
#define CLR(u) for(int i=0;i<=N;i++)u[i].clear();
#define LL long long
#define ULL unsigned long long
#define mp make_pair
#define PII pair<int,int>
#define PIL pair<int,long long>
#define PLL pair<long long,long long>
#define pb push_back
#define fi first
#define se second
typedef vector<int> VI;
int read(){int x = 0,f = 1;char c = getchar();while (c<'0' || c>'9'){if (c == '-') f = -1;c = getchar();}
while (c >= '0'&&c <= '9'){x = x * 10 + c - '0';c = getchar();}return x*f;}
const double eps = 1e-9;
const int maxn = 3e5 + 10;
const int maxm = 2e6 + 10;
const int INF = 0x3f3f3f3f;
const int mod = 1e9 + 7;
int N,M,K;
int head[maxn],tot;
struct Edge{int to,next,dis;
}edge[maxm];
void init(){for(int i = 0 ; i <= 3 * N + 1; i ++) head[i] = -1;tot = 0;
}
void add(int u,int v,int w){edge[tot].to = v;edge[tot].next = head[u];edge[tot].dis = w;head[u] = tot++;
}
int val[maxn];
int dis[maxn],vis[maxn];
int SPFA(int s,int t){for(int i = 0 ; i <= 3 * N + 1; i ++) dis[i] = -INF;dis[1] = 0;queue<int>Q;Q.push(s);while(!Q.empty()){int u = Q.front(); Q.pop();vis[u] = 0;for(int i = head[u]; ~i ; i = edge[i].next){int v = edge[i].to;if(dis[v] < dis[u] + edge[i].dis){dis[v] = dis[u] + edge[i].dis;if(!vis[v]){vis[v] = 1;Q.push(v);}}} }if(dis[t] == -INF) dis[t] = 0;return dis[t];
}
int main(){Sca2(N,M); init();for(int i = 1; i <= N ; i ++){Sca(val[i]);add(i,i + N,-val[i]);add(i + N,i + N + N,val[i]);} for(int i = 1; i <= M ; i ++){int u,v,w; Sca3(u,v,w);add(u,v,0);add(u + N,v + N,0);add(u + N + N,v + N + N,0);if(w == 2){add(v,u,0);add(v + N,u + N,0);add(v + N,u + N + N,0);} }int t = 3 * N + 1,s = 1;add(N,t,0); add(3 * N,t,0);Pri(SPFA(s,t));return 0;
}

转载于:https://www.cnblogs.com/Hugh-Locke/p/10330356.html

洛谷P1073 Tarjan + 拓扑排序 // 构造分层图相关推荐

  1. 洛谷 - P4009 汽车加油行驶问题(分层图最短路/最小费用最大流)

    题目链接:点击查看 题目大意:给出一个n*n的矩阵表示道路,途中有一些加油站,现在要从点(1,1)到达点(n,n),问最小花费,其中的一些规则如下: 汽车只能沿着网格边行驶,装满油后可以行驶K条边,出 ...

  2. 洛谷4400 BlueMary的旅行(分层图+最大流)

    qwq 首先,我们观察到题目中提到的每天只能乘坐一次航班的限制,很容易想到建分层图,也就是通过枚举天数,然后每天加入一层新的点. (然而我一开始想的却是erf) 考虑从小到大枚举天数,然后每次新建一层 ...

  3. 图论--拓扑排序--判断一个图能否被拓扑排序

    拓扑排序的实现条件,以及结合应用场景,我们都能得到拓扑排序适用于DAG图(Directed Acyclic Graph简称DAG)有向无环图, 根据关系我们能得到一个线性序列,实现的方式是DFS,具体 ...

  4. 洛谷 P1073 最优贸易 (分层图状态转移+SPFA,求最长路径;另附某dalao的超短代码:暴力+动规)

    题目链接1 题目链接2 另附某dalao的超短代码:暴力+动规 P1073 最优贸易 题目描述 C国有 n 个大城市和 m 条道路,每条道路连接这 n 个城市中的某两个城市. 任意两个城市之间最多只有 ...

  5. 队爷的讲学计划(tarjan +拓扑排序)

    题目描述 队爷为了造福社会,准备到各地去讲学.他的计划中有n个城市,从u到v可能有一条单向道路,通过这条道路所需费用为q.当队爷在u城市讲学完之后,u城市会派出一名使者与他同行,只要使者和他在一起,他 ...

  6. 二分图——洛谷P1155 双栈排序

    https://daniu.luogu.org/problem/show?pid=1155 二分图染色+模拟 1.首先考虑一个简单情况--单栈排序,显然有这样的一个事实: a[i]和a[j] 不能压入 ...

  7. 洛谷 P1073 最优贸易

    题目描述 CC C 国有 n n n 个大城市和 m mm 条道路,每条道路连接这 nnn 个城市中的某两个城市.任意两个城市之间最多只有一条道路直接相连.这 mmm 条道路中有一部分为单向通行的道路 ...

  8. 【算法】图(一)拓扑排序的实现 图的邻接表算法 判断是否图G中存在环

    文章目录 用list来表示图,判断是否存在环 邻接表实现拓扑排序 用DFS(邻接矩阵) 来实现拓扑排序. 判断无向图顶点是否全部连通 判断图G中从顶点u到v是否存在简单路径 输出图G中从顶点u到v的所 ...

  9. Robots on a Grid CodeForces - 1335F(拓扑排序+正反建图+判环)

    There is a rectangular grid of size n×m. Each cell of the grid is colored black ('0') or white ('1') ...

最新文章

  1. 逻辑回归 logistic regression
  2. 动态分辨率是什么意思_b站么么直播最新动态里都有啥 b站什么意思
  3. apache+mod_wsgi+django的环境配置
  4. 如何加快Json 序列化?有哪些方法?
  5. Android N 新特性 + APP开发注意事项
  6. IOS开发基础知识--碎片5
  7. SSL、OPENSSL、SSH、OPENSSH
  8. 监控 MySQL的多种方法
  9. php时间格式转换成时间戳,php把时间格式转换为时间戳的案例
  10. SpringBoot用HttpClient调用快递物流查询API接口
  11. 系统学习机器学习之组合多分类器
  12. 使用logisim设计简易CPU
  13. 跳转微信功能(学习总结)
  14. jquery开发知识总结2
  15. vue 组件第一次不渲染问题_解决vue页面渲染但dom没渲染的操作
  16. Fedora下的应用软件大集合
  17. 华清远见嵌入式培训_第一周回顾与反思
  18. 阿里云视频云推出低代码音视频工厂vPaaS
  19. libvirt 安装
  20. java计算机毕业设计英语课程学习网站源码+数据库+系统+lw文档+mybatis+运行部署

热门文章

  1. keil优化等级设置
  2. BREW做的第一个程序--Hello world!
  3. 如何实现从wgs-84到beijing54的坐标转换
  4. vue render函数
  5. Node --- Promise中的多异步协作
  6. 世界顶级的程序员们告诉你:这些书都是你应该读的
  7. LVM逻辑卷详解及创建
  8. 甲骨文称 Java 序列化的存在是个错误,计划删除
  9. python 进程编程速成
  10. 【第二十七章】 springboot + zipkin(brave-okhttp实现)