当前位置: 首页 > news >正文

题解:P14038 [PAIO 2025] Adventure Plan

第一次给官方比赛投题,来写发题解 /se

子任务 \(1,3,4\) 随便乱做就好了,没啥技术含量。

\(x_u\) 表示从 \(0\)\(u\) 路径的长度,那么一条有向边 \(u\stackrel{[l,r]}{\to}v\) 等价于一个限制条件 \(l\le x_v-x_u\le r\)。将题目转化为一个差分约束系统,根据三角不等式,只需在差分约束系统建图时连边 \(u\stackrel{r}{\to}v\)\(v\stackrel{-l}{\to}u\) 即可。

容易想到,对于每一次加边操作,可以使用 Bellman-Ford 算法判断负环,没有负环就说明可以加边。最后跑一遍 Bellman-Ford 算法求出所有 \(x_u\),对于一条边 \(u\to v\),将其权值设置为 \(x_v-x_u\) 即可。

时间复杂度 \(O(nm+qn(m+q))\),可以通过子任务 \(2,5\),结合上面的算法可以获得 \(75\) 分。

每次加边都要重新跑一遍 \(O(nm)\) 的 Bellman-Ford 算法,时间开销过大。有没有其他最短路算法支持判断负环、求带有负权边的图的最短路呢?答案是有的。

初始时使用 Floyd-Warshall 算法求出每对点之间的最短路,每次加边时暴力使用这条边更新每对点,若存在 \(\operatorname{dis}(u,u)<0\) 则说明有负环。

时间复杂度 \(O(m+n^3+qn^2)\),可以 AC 本题。

//By: OIer rui_er
#include <bits/stdc++.h>
#define rep(x, y, z) for(int x = (y); x <= (z); ++x)
#define per(x, y, z) for(int x = (y); x >= (z); --x)
#define endl '\n'
using namespace std;
typedef long long ll;template<typename T> void chkmin(T& x, T y) {if(y < x) x = y;}
template<typename T> void chkmax(T& x, T y) {if(x < y) x = y;}const ll inf = 0x3f3f3f3f3f3f3f3fll;int _N;
vector<ll> _dis;
vector<tuple<int, int>> _e;vector<bool> add_roads(int N, int M, int Q,vector<int> U, vector<int> V,vector<int> L, vector<int> R,vector<int> U2, vector<int> V2,vector<int> L2, vector<int> R2) {_N = N;_dis.resize(N);vector<vector<ll>> e(N, vector<ll>(N));rep(i, 0, N - 1)rep(j, 0, N - 1)e[i][j] = i == j ? 0LL : +inf;rep(i, 0, M - 1) {int u = U[i], v = V[i];ll l = L[i], r = R[i];_e.emplace_back(u, v);chkmin(e[u][v], +r);chkmin(e[v][u], -l);}rep(k, 0, N - 1)rep(i, 0, N - 1)rep(j, 0, N - 1)if(e[i][k] < +inf && e[k][j] < +inf)chkmin(e[i][j], e[i][k] + e[k][j]);vector<bool> ans(Q);rep(t, 0, Q - 1) {int u = U2[t], v = V2[t];ll l = L2[t], r = R2[t];vector<vector<ll>> f = e;rep(i, 0, N - 1)rep(j, 0, N - 1)if(f[i][u] < +inf && f[v][j] < +inf)chkmin(f[i][j], f[i][u] + r + f[v][j]);rep(i, 0, N - 1)rep(j, 0, N - 1)if(f[i][v] < +inf && f[u][j] < +inf)chkmin(f[i][j], f[i][v] - l + f[u][j]);bool ok = true;rep(i, 0, N - 1) if(f[i][i] < 0) ok = false;if(ok) {_e.emplace_back(u, v);e = f;}ans[t] = ok;}rep(i, 0, N - 1) _dis[i] = e[0][i];return ans;
}vector<int> assign_times() {int M = _e.size();vector<int> ans(M);rep(i, 0, M - 1) {auto [u, v] = _e[i];ans[i] = _dis[v] - _dis[u];}return ans;
}#ifdef LOCALint main() {ios::sync_with_stdio(false);cin.tie(0); cout.tie(0);int N, M, Q;cin >> N >> M >> Q;vector<int> U(M), V(M), L(M), R(M);rep(i, 0, M - 1) cin >> U[i];rep(i, 0, M - 1) cin >> V[i];rep(i, 0, M - 1) cin >> L[i];rep(i, 0, M - 1) cin >> R[i];vector<int> U2(Q), V2(Q), L2(Q), R2(Q);rep(i, 0, Q - 1) cin >> U2[i];rep(i, 0, Q - 1) cin >> V2[i];rep(i, 0, Q - 1) cin >> L2[i];rep(i, 0, Q - 1) cin >> R2[i];vector<bool> ans = add_roads(N, M, Q,U, V, L, R, U2, V2, L2, R2);rep(i, 0, Q - 1) cout << ans[i] << " \n"[i == Q - 1];vector<int> dis = assign_times();int sz = dis.size();rep(i, 0, sz - 1) cout << dis[i] << " \n"[i == sz - 1];return 0;
}#endif
http://www.hskmm.com/?act=detail&tid=21246

相关文章:

  • 20231414_王仕琪_密码技术密码杂凑算法学习笔记
  • web3D、webGL、webGPU、webGIS、webXR、webCodecs的概念和对比 - 实践
  • Claude code的 thinking on/off差别有多少
  • Ubuntu 25的网络配置
  • 2025.9.26 测试
  • 贝叶斯学习笔记 - 详解
  • Ubuntu 24和25配置apt国内源
  • 实用指南:AWS实战:轻松创建弹性IP,实现固定公网IP地址
  • 完整教程:自然语言处理项目之情感分析(下)
  • 还在为安装PS发愁?这款网页版工具,打开浏览器就能用!
  • 委托相关
  • 清除“请允许观看视频”通知页面的完整指南
  • 千亿芯片公司被股东“抛弃” ,AI芯片第一股前景几何?
  • Java 与智慧港口:航运调度与物流枢纽数字化
  • DeepSeek-V3.2-Exp 发布,训练推理提效,API 同步降价
  • 图片任意切割工具(Python 3.8 实现)
  • 从零构建能自我优化的AI Agent:Reflection和Reflexion机制对比详解与实现
  • 超精简的小型C编译器
  • Day1 Linux 入门:9 个核心命令(whoami/id/pwd 等)
  • 9.29 闲话
  • MMU的作用
  • 大二学计算机系统基础
  • 20250929 之所思 - 人生如梦
  • 9/29
  • 9.29总结
  • lc1040-移动石子直到连续II
  • 2025年9月29日
  • c++算法学习笔记
  • test5
  • 最高人民法院新劳动争议司法解释一 理解与适用