[题解][YZOJ50113] 枇杷树
简要题意
\(m\) 个操作,每次操作都会产生一个树的版本 \((\)从 \(0\) 开始\()\).
一次操作把 \(x_i\) 版本的树的点 \(u\) 和 \(y_i\) 版本的树的点 \(v\) 连一条权值是 \(w\) 的边。\((\) \(y_i\) 上的全部点的编号加 \(siz_{x_i}\) \()\)
求每个版本:
\]
\(m\le 300,u,v\leq 10^{18}, siz \leq 2 \times 10 ^{18}\)
解题思路
首先,将答案考虑为每条边乘上两边的子树大小。
然后,考虑产生的新版本 (记为 \(z\) ) 的答案 \(ans_z\) 的组成:
- 原来版本内的路径和:\(ans_x+ans_y\);
- 边 \((u,v,w)\) 的贡献:\(siz_x\times siz_y\times w\);
- 版本间的贡献 ( 记 \(f_x(u)\) 为版本 \(x\) 中,所有点到 \(u\) 的路径和):\(f_x(u)\times siz_y+f_y(v)\times siz_x\).
那么现在的主要问题就是每次如何得到 \(f_x(u)\)。( 记 \(dis_x(u,v)\) 表示版本 \(x\) 中,\(u\) 到 \(v\) 的路径长度)
还是分析当得到一个新的版本 \(z\) 时,\(f_z(a)\) 的组成:
- 若 \(a\) 原来属于版本 \(x\),\(f_x(a)+f_y(v)+(w+dis_x(a, u))\times siz_y\).
- 若 \(a\) 原来属于版本 \(y\),\(f_y(a)+f_x(u)+(w+dis_y(a, v))\times siz_x\).
此时,\(f_x(u)/f_y(v)\) 算作已知,故只需考虑 \(dis_z(a,b)\) 如何计算:
- \(a,b\) 都在版本 \(x/y\),调用 \(dis_{x/y}(a,b)\);
- \(a,b\) 不在同一版本,\(dis_x(a,u)+w+dis_y(b_v)\).
此时整个问题在递归中解决了,现在来分析一下复杂度。
对于 \(f_x(u)\),版本数是 \(O(m)\) 级别,\(u\) 的取值是 \(O(m)\) 级别,总的状态就是 \(O(m^2)\) 级别,转移 \(O(1)\)。
\(dis_x(a, b)\) 的计算平均是 \(O(1)\) 的,所以计算所有的 \(f_x(u)\) 是 \(O(m^2)\) 的。
计算 \(ans\) 时, \(f_x(u)\) 会被调用 \(O(m)\) 次,所以总的时间复杂度是 \(O(m^3)\)。
代码
inline int Dis(int z, int a, int b){
if(a == b) return 0;
if(dis[z].find({a, b}) != dis[z].end()) return dis[z][{a, b}];
int x = X[z], y = Y[z], w = W[z]; LL u = U[z], v = V[z];
if(a <= siz[x] && b <= siz[x]) return dis[z][{a, b}] = Dis(x, a, b);
if(siz[x] < a && siz[x] < b) return dis[z][{a, b}] = Dis(y, a - siz[x], b - siz[x]);
int sum = w; if(a > siz[x]) swap(a, b);
Plus(sum, Dis(x, a, u)), Plus(sum, Dis(y, b - siz[x], v));
return dis[z][{a, b}] = sum;
}
inline int F(int z, LL a){
if(z == 0 && a == 1) return 0;
if(f[z].find(a) != f[z].end()) return f[z][a];
int x = X[z], y = Y[z], w = W[z]; LL u = U[z], v = V[z];
if(a <= siz[x]){
int sum = 0;
Plus(sum, F(x, a)), Plus(sum, F(y, v));
Plus(sum, siz[y] % mod * Mod(w + Dis(x, a, u) - mod) % mod);
return f[z][a] = sum;
}
int sum = 0;
Plus(sum, F(y, a - siz[x])), Plus(sum, F(x, u));
Plus(sum, siz[x] % mod * Mod(w + Dis(y, a - siz[x], v) - mod) % mod);
return f[z][a] = sum;
}
inline void solve(int x, int y, LL u, LL v, int w, int z){
X[z] = x, Y[z] = y, U[z] = u, V[z] = v, W[z] = w;
Plus(Ans[z], Mod(Ans[x] + Ans[y] - mod));
Plus(Ans[z], siz[x] % mod * (siz[y] % mod) % mod * w % mod);
Plus(Ans[z], Mod(siz[y] % mod * F(x, u) % mod + siz[x] % mod * F(y, v) % mod - mod));
siz[z] = siz[x] + siz[y];
}
参考
https://www.cnblogs.com/qjbqjb/p/15905279.html
[题解][YZOJ50113] 枇杷树的更多相关文章
- 2016 华南师大ACM校赛 SCNUCPC 非官方题解
我要举报本次校赛出题人的消极出题!!! 官方题解请戳:http://3.scnuacm2015.sinaapp.com/?p=89(其实就是一堆代码没有题解) A. 树链剖分数据结构板题 题目大意:我 ...
- noip2016十连测题解
以下代码为了阅读方便,省去以下头文件: #include <iostream> #include <stdio.h> #include <math.h> #incl ...
- BZOJ-2561-最小生成树 题解(最小割)
2561: 最小生成树(题解) Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 1628 Solved: 786 传送门:http://www.lyd ...
- Codeforces Round #353 (Div. 2) ABCDE 题解 python
Problems # Name A Infinite Sequence standard input/output 1 s, 256 MB x3509 B Restoring P ...
- 哈尔滨理工大学ACM全国邀请赛(网络同步赛)题解
题目链接 提交连接:http://acm-software.hrbust.edu.cn/problemset.php?page=5 1470-1482 只做出来四道比较水的题目,还需要加强中等题的训练 ...
- 2016ACM青岛区域赛题解
A.Relic Discovery_hdu5982 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Jav ...
- poj1399 hoj1037 Direct Visibility 题解 (宽搜)
http://poj.org/problem?id=1399 http://acm.hit.edu.cn/hoj/problem/view?id=1037 题意: 在一个最多200*200的minec ...
- 网络流n题 题解
学会了网络流,就经常闲的没事儿刷网络流--于是乎来一发题解. 1. COGS2093 花园的守护之神 题意:给定一个带权无向图,问至少删除多少条边才能使得s-t最短路的长度变长. 用Dijkstra或 ...
- CF100965C题解..
求方程 \[ \begin{array}\\ \sum_{i=1}^n x_i & \equiv & a_1 \pmod{p} \\ \sum_{i=1}^n x_i^2 & ...
随机推荐
- java-jsp-learnning
简介:JSP(全称Java Server Pages) 一种动态网页开发技术.它使用JSP标签在HTML网页中插入Java代码.标签通常以<%开头以%>结束. JSP是一种Java ser ...
- 怎么将 byte 转换为 String?
可以使用 String 接收 byte[] 参数的构造器来进行转换,需要注意的点是要使用 的正确的编码,否则会使用平台默认编码,这个编码可能跟原来的编码相同,也 可能不同.
- 初识Spring(为什么要使用Spring?)
Spring,英文翻译是春天的意思,而在Java中,是一个开放源代码的设计层面框架(手动滑稽,程序员的春天),他解决的是业务逻辑层和其他各层的松耦合问题,因此它将面向接口的编程思想贯穿整个系统应用.S ...
- java-jdbc-all
jdbc相关解析 JDBC(Java DataBase Connectivity,Java数据库连接)是一种用于执行SQL语句的Java API,可以为多种关系数据库提供统一访问,它由一组用Java语 ...
- 关于TP框架
TP的特性有哪些? 1.多表查询非常方便,在model中几句代码就可以完成对多表的关联操作. 2.融合了smarty模板,使前后台分离 3.支持多种缓存技术,尤其对memcache技术支持非常好 4. ...
- mybatis-01-基本流程
mybatis执行流程 1. 加载配置文件并初始化(SqlSession) 配置文件来源于两个地方,一个是配置文件(主配置文件conf.xml,mapper文件*.xml), 一个是java代码中 ...
- 微信小程序加密数据(encryptedData)解密中的PHP代码,php7.1报错
问题描述 最近在开发微信小程序涉及到加密数据(encryptedData)的解密,用的是PHP代码,在运行后报错mcrypt_module_ xxx is deprecated,提示方法已过时了 经研 ...
- Eureka Server 的 Instance Status 一直显示主机名问题
Eureka Server 的 Instance Status 一直显示主机名问题 注册中心启动后,访问 http://localhost:8761/ 后: 如何调整为具体所在的服务器 IP 呢? 解 ...
- centos和redhat的区别
CentOS(Community Enterprise Operating System,中文意思是:社区企业操作系统)是Linux发行版之一,它是来自于Red Hat Enterprise Linu ...
- Java学习day14
可变参数用作方法的形参,方法参数的个数就可变 格式:修饰符 返回值类型 方法名(数据类型...变量名){ } 方法内的形参只能有一个,这里的变量是一个数组 public static <T> ...