题解【洛谷P5959】[POI2018]Plan metra
一道比较神仙的构造题。
首先确定 \(1\) 到 \(n\) 的路径长度,不妨设其长为 \(m\) 。
通过观察发现,\(m\) 就是 \(\min_{1<i<n}\{dist_{1,i} + dist_{i,n}\}\)。
如果所有的 \(| dist_{1, i} - dist_{i, n} |\) 都相等,那么可以特判一下:
- 首先将 \(1\) 与 \(n\) 连接起来,路径长度为 \(| dist_{1, i} - dist_{i, n} |\) ;
- 对于每一个 \(\forall 1<i<n\) :
- 如果\(dist_{1,i}>dist_{i,n}\),那么将 \(i\) 与 \(1\) 连接,长度为 \(dist_{1,i}\) ;
- 如果\(dist_{1,i} \leq dist_{i,n}\) ,那么将 \(i\) 与 \(n\)连接,长度为 \(dist_{i, n}\) 。
若点 \(i\) 满足 \(dist_{1,i} + dist_{i,n} = m\),则 \(i\) 在 \(1\) 到 \(n\) 的路径上。
否则可以求出 \(i\) 到 \((1, n)\) 路径的距离和 \(i\) 在路径上的“投影”的
位置,查找位置可以使用二分查找。
将 \(i\) 挂在路径上对应的点上即可。
如果找不到可以挂起的位置或者 \(1\) 到 \(n\) 的路径上有几个点重复就无解。
注意特判 \(n=2\) 的情况。
代码如下(刚好 \(100\) 行…):
#include <bits/stdc++.h>
#define DEBUG fprintf(stderr, "Passing [%s] line %d\n", __FUNCTION__, __LINE__)
#define itn int
#define gI gi
#define mk make_pair
using namespace std;
inline int gi()
{
int f = 1, x = 0; 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 f * x;
}
const int maxn = 500003;
int sz, ddd, len, n, m = 1000000007, d1[maxn], dn[maxn], dis[maxn], jlx[maxn];
struct Node
{
int id, dis;
} a[maxn];
struct Ans
{
int u, v, w;
} ans[maxn];
inline bool cmp(Node x, Node y) {return x.dis < y.dis;}
inline void sub1_getans()
{
puts("TAK");
printf("1 %d %d\n", n, len);
for (int i = 2; i < n; i+=1)
{
if (d1[i] < dn[i]) printf("1 %d %d\n", i, d1[i]);
else printf("%d %d %d\n", n, i, dn[i]);
}
return;
}
inline int Binary_Search(int x)
{
int l = 1, r = sz;
while (l <= r)
{
int mid = (l + r) >> 1;
if (a[mid].dis == x) return mid;
else if (a[mid].dis < x) l = mid + 1;
else r = mid - 1;
}
return -1;
}
int main()
{
n = gi();
if (n == 2) //特判
{
puts("TAK\n1 2 1");
return 0;
}
d1[1] = dn[n] = 0;
for (int i = 2; i < n; i+=1) d1[i] = gi();
for (int i = 2; i < n; i+=1) dn[i] = gi(), jlx[i] = abs(d1[i] - dn[i]);
bool fl = true;
for (int i = 3; i < n; i+=1) if (jlx[i] != jlx[i - 1]) {fl = false; break;}
if (fl) {len = jlx[2]; sub1_getans(); return 0;} //特判
for (int i = 2; i < n; i+=1) m = min(m, d1[i] + dn[i]);
a[++sz] = (Node){1, 0};
d1[n] = dn[1] = m;
for (int i = 2; i < n; i+=1)
if (d1[i] + dn[i] == m) a[++sz] = (Node){i, d1[i]};
a[++sz] = (Node){n, d1[n]};
sort(a + 1, a + 1 + sz, cmp);
for (int i = 2; i < n; i+=1)
{
dis[i] = d1[i] - (d1[i] + dn[i] - m) / 2;
}
for (int i = 1; i < sz; i+=1)
{
if (a[i].dis == a[i + 1].dis) {puts("NIE"); return 0;}
ans[++ddd] = (Ans){a[i].id, a[i + 1].id, a[i + 1].dis - a[i].dis};
}
for (int i = 2; i < n; i+=1)
{
int nowlen = d1[i] + dn[i];
if (nowlen != m)
{
int find_jl = Binary_Search(dis[i]); //二分挂起的位置
if (((d1[i] + dn[i] - m) & 1) || find_jl == -1) {puts("NIE"); return 0;}
ans[++ddd] = (Ans){a[find_jl].id, i, (d1[i] + dn[i] - m) / 2};
}
}
puts("TAK");
for (int i = 1; i <= ddd; i+=1) printf("%d %d %d\n", ans[i].u, ans[i].v, ans[i].w);
return 0;
}
另:这份代码在洛谷上 AC 了, 但是在 BZOJ 上会 WA, 不知道为什么……
题解【洛谷P5959】[POI2018]Plan metra的更多相关文章
- 【BZOJ5100】[POI2018]Plan metra 构造
[BZOJ5100][POI2018]Plan metra Description 有一棵n个点的无根树,每条边有一个正整数权值,表示长度,定义两点距离为在树上的最短路径的长度. 已知2到n-1每个点 ...
- bzoj5100 [POI2018]Plan metra 构造
5100: [POI2018]Plan metra Time Limit: 40 Sec Memory Limit: 128 MBSec Special JudgeSubmit: 189 Sol ...
- bzoj千题计划249:bzoj5100: [POI2018]Plan metra
http://www.lydsy.com/JudgeOnline/problem.php?id=5100 1.找到d1[i]+dn[i] 最小的点,作为1到n链上的点 2.令链长为D,若abs(d1[ ...
- 题解 洛谷P5018【对称二叉树】(noip2018T4)
\(noip2018\) \(T4\)题解 其实呢,我是觉得这题比\(T3\)水到不知道哪里去了 毕竟我比较菜,不大会\(dp\) 好了开始讲正事 这题其实考察的其实就是选手对D(大)F(法)S(师) ...
- 题解 洛谷 P3396 【哈希冲突】(根号分治)
根号分治 前言 本题是一道讲解根号分治思想的论文题(然鹅我并没有找到论文),正 如论文中所说,根号算法--不仅是分块,根号分治利用的思想和分块像 似却又不同,某一篇洛谷日报中说过,分块算法实质上是一种 ...
- 题解-洛谷P5410 【模板】扩展 KMP(Z 函数)
题面 洛谷P5410 [模板]扩展 KMP(Z 函数) 给定两个字符串 \(a,b\),要求出两个数组:\(b\) 的 \(z\) 函数数组 \(z\).\(b\) 与 \(a\) 的每一个后缀的 L ...
- 题解-洛谷P4229 某位歌姬的故事
题面 洛谷P4229 某位歌姬的故事 \(T\) 组测试数据.有 \(n\) 个音节,每个音节 \(h_i\in[1,A]\),还有 \(m\) 个限制 \((l_i,r_i,g_i)\) 表示 \( ...
- 题解-洛谷P4724 【模板】三维凸包
洛谷P4724 [模板]三维凸包 给出空间中 \(n\) 个点 \(p_i\),求凸包表面积. 数据范围:\(1\le n\le 2000\). 这篇题解因为是世界上最逊的人写的,所以也会有求凸包体积 ...
- 题解-洛谷P4859 已经没有什么好害怕的了
洛谷P4859 已经没有什么好害怕的了 给定 \(n\) 和 \(k\),\(n\) 个糖果能量 \(a_i\) 和 \(n\) 个药片能量 \(b_i\),每个 \(a_i\) 和 \(b_i\) ...
随机推荐
- linux 网络接口,ip地址,路由设定
本文是基于centos 配置DNS条目: 配置文件:/etc/resolv.conf 修改主机名称: 命令:hostname NAME.重启后失效 配置文件:/etc/sysconfig/networ ...
- java-局部变量,成员变量区别
1. 内存中的位置 成员变量: 堆内存 局部变量: 栈内存 2. 生命周期 成员变量:随着对象的创建而存在,随着对象的消失而消失 局部变量:随着方法的调用而存在,随着方法调用完毕而消失 3. 注意事项 ...
- Page Visibility API
在code review时看见同事使用visibilitychange 事件来监听页面的隐藏与显示,之前没有了解过这块,学习一下. document.visibilityState 主要有以下3个状态 ...
- 论文阅读笔记(十五)【CVPR2016】:Top-push Video-based Person Re-identification
Approach 特征由两部分组成:space-time特征和外貌特征.space-time特征由HOG3D[传送门]提取,其包含了空间梯度和时间动态信息:外貌特征采用颜色直方图[传送门]和LBP[传 ...
- 【Debian学徒记事】记一次解决Debian开机1min30s
记一次解决Debian开机1min30s 打开我亲爱的Debian 欸,好像有点慢 [* * * * ] A start job is running for....(*/1min30s) [TIME ...
- Spark学习之路 (八)SparkCore的调优之开发调优[转]
前言 在大数据计算领域,Spark已经成为了越来越流行.越来越受欢迎的计算平台之一.Spark的功能涵盖了大数据领域的离线批处理.SQL类处理.流式/实时计算.机器学习.图计算等各种不同类型的计算操作 ...
- 关于f(x)
有时 z = x + y 有时 0 = x + y 有时单独用f(x) 有时 z = f(x) 很容易分不清. 从集合角度,将f(x)看成映射 即从A集合到B集合的对应关系 这样f(x)可以单独使用, ...
- python的for循环的神奇之处
python的for循环太神奇了: 你可以编写这样的语句: for i in range(10) : j= i**2 print(j) 你也可以编写这样的语句: with open('/path/to ...
- Java集合之Collections 剖析
Collections工具类位于 java.util 包下,是一个比较常用的工具类,关于这个工具类,主要介绍其在使用过程中遇到的大坑!!! [事故现场] 在实际项目开发过程中,在前人代码的基础上,对于 ...
- 500kuai
https://www.bilibili.com/bangumi/media/md11653495/?spm_id_from=666.10.b_62616e67756d695f64657461696c ...