当前位置: 首页 > 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.rkmt.cn/news/14160.html

相关文章:

  • 20231414_王仕琪_密码技术密码杂凑算法学习笔记
  • 堆设置了8G,java进程却占用了12G内存
  • IntelliJ IDEA 中 Shared Build Process Heap Size 的重要性与配置
  • 企业数字化转型战略规划:从愿景到落地的完整路径
  • Java 与智慧港口:航运调度与物流枢纽数字化
  • 千亿芯片公司被股东“抛弃” ,AI芯片第一股前景几何?
  • DeepSeek-V3.2-Exp 发布,训练推理提效,API 同步降价
  • 9.29 闲话
  • US$16 Yanhua BMW F/G Chassis Odometer Wiring Harness
  • US$980 Xhorse VVDI2 BMW OBD + CAS4 +FEM/BDC Functions Full BMW License
  • 2025年9月29日
  • 最高人民法院新劳动争议司法解释一 理解与适用
  • 逆元 组合数问题
  • 一、驱动基础知识速览(迅为RK3568)
  • MonoDETR(2)
  • 记录---window.close()失效 + Chrome浏览器调试线上代码
  • 启发式合并 [USACO22DEC] Making Friends P
  • 加密的病例单
  • 【多线程】什么是原子操作(Atomic Operation)? - 详解
  • 复刻江协旋钮控制模块
  • c语言switch和if语句
  • Qt(制作一个方便的文本编辑器)
  • tldr的安装与利用
  • 题解:P7126 [Ynoi2008] rdCcot
  • 实用指南:汽车地带AutoZone EDI需求分析及对接指南
  • 航司网站url后缀参数FECU分析
  • 优化 if/else 的四种设计模式
  • 多corner综合
  • Day11-C:\Users\Lenovo\Desktop\note\code\JavaSE\Basic\src\com\oop\demo06
  • OpenLayers地图交互 -- 章节十一:拖拽材料交互详解