【csp模拟赛1】铁路网络 (network.cpp)
【题目描述】
在暑假来临之际,小 Z 的地理老师布置了一个暑假作业,让同学们暑假期间 了解一下 C 国的铁路发展史。小 Z 在多番查证资料后发现,C 国在铁路发展初期, 铁路网络有着一个严密规整的结构:C 国的 N 个城市按层级分为首都、省会、省 辖市……由于在铁路发展初期,建造铁路花费巨大,为了避免不必要的浪费,铁 路网络的建设保证了任意两个城市相互连通,且仅有一条连通路径。 随后,小 Z 还查阅到历年来铁道部向全国上下发布的大小消息。由于多年来 C 国货币系统保持健康的 3%通胀率带来了稳定的物价上涨,铁道部发布的消息中 不乏对某段连接城市 u 与城市 v 铁路的价格上涨 w 元的通知。于是小 Z 不禁好奇, 在某些特定时间,在某个城市 p 的管辖区域内,任意两个城市间的铁路通勤的费 用和是多少。
【输入格式】
第一行输入两个正整数 N,Q,分别表示城市的数目和操作的数目。 接下来有 N–1 行,第 i 行是两个正整数 p[i], c[i],表示城市 p[i]是城市 i 的父亲结点,且连接 p[i]和 i 的铁路的初始收费为 c[i](1≤c[i]≤1000)。 再接下来有 Q 行,每行是如下两种类型之一: INC u v w (u, v, w 都是整数,且 1≤u, v≤N, 0≤w≤1000) ASK p (p 是整数,且 0≤p≤N) 意义如题目所述。
【输出格式】
对每个 ASK 类型的操作,输出所求的答案。请你输出答案对 2019 取模后的 结果。
【样例输入】
5 5
1 1
2 5
1 2
2 1
INC 2 4 2
INC 3 4 1
ASK 2
INC 2 5 3
ASK 1
【样例输出】
14 84
思路:
题目大意:给你一棵有根树,每条边有边权。实现两种操作:
1.给某一条路径上所 有边的权值加上一个数;
2.询问某棵子树内所有点对的距离和。 本题设置了几个部分分梯度,得部分分的做法有不少, 复杂度也各有千秋。下面给出比较容易理解的一种满分解 法。
从查询的角度思考,如何计算以 i 为根的子树内所有点 对的距离和?
不妨设点 x 与它的父亲 fa[x]相连的边的权值为 p[x],考虑 p[x]会对那些点对产生贡 献?显然是经过 x——fa[x]这条边的那些点对。
记 siz[x]为以 x 为根的子树的大小。则经过 x——fa[x]的点对有 siz[x]×(siz[i] – siz[x])对(如图 5),
于是,子树 i 内的点 x 的贡献就是 p[x]×siz[x]×siz[i] – p[x]×(siz[x])²。黄色的部分与 i 无关。
如果要计算,我们只需把子树 i 内所有点(不包括 i)的自己的黄色部分加起来,然后就可以求出答案。
这是因为 ∑(p[x]×siz[x]×siz[i] – p[x]×(siz[x])²) = siz[i]×∑(p[x]×siz[x]) – ∑(p[x]×(siz[x])²) 而修改操作,就是给一条路径上所有点(准确地说,不含 LCA)的 p 加上一个数。
我们要高效地维护上面黄色部分标出的两种“和”。这里有两种解法:
解法一
:
树链剖分 + DFS 序。 修改操作在树的链上进行,可以考虑使用树链剖分;而查询操作是查询子树的和,可 以考虑用 DFS 序。 树链剖分,把重链的结点在线段树中连续存放; 而 DFS 序,把子树结点连续存放。
两者可以兼得吗? 事实上,进行树链剖分的时候,我们 记 录了每个结点的“重儿子”,求树的 DFS 序 时, 我们对每个结点都先搜索它的重儿子再搜 索 其他儿子。
这样,就可以保证,DFS 序中, 每 条重链是连续存放的,每棵子树也是连续 存 放的,可以使用同一棵线段树存放。(不 妨 参考如图 6 的例子)
在线段树中,很容易 处 理让一段的值都加上一个数和查询一段的 和 这样的操作。 这样做是 O(n (log2n) 2)的。
代码:
#include<bits/stdc++.h>
typedef long long LL;
const int N = 100000;
const int M = 5000;
const int inf = 2147483647;
const int p = 2019;
char buf[1 << 23], *p1 = buf, *p2 = buf, obuf[1 << 23], *O = obuf;
inline void _(char ch){*O++ = ch;}
void print(int x) {if(x < 0) _('-'), x = -x; if(x>9) print(x/10);*O++ = x % 10 ^ 48;}
inline int read(){
int res = 0; char ch = getchar(); bool bo = false;
while(ch < '0' || ch > '9') bo = (ch == '-'), ch = getchar();
while(ch >= '0' && ch <= '9') res = (res << 1) + (res << 3) + (ch ^ 48), ch = getchar();
return bo ? -res : res;
}
int getc() {
char ch = getchar();
while(ch != 'A' && ch != 'I') ch = getchar();
return ch == 'I';
}
struct Edge {
int to;
Edge *nxt;
Edge (int to, Edge *nxt) : to(to), nxt(nxt) {}
}*head[N];
int dep[N], a[N], size[N], siz[N], son[N], Top[N], fa[N], dfn[N], nfd[N], tot;
struct Segment {
struct node {
int sum, tag, siz;
int l, r;
node *L, *R; node(int l, int r) : sum(0), tag(0), siz(0), l(l), r(r), L(NULL), R(NULL) {} void up() {
siz = (L->siz + R->siz) % p;
sum = (L->sum + R->sum) % p;
} void down() {
L->sum = (L->sum + L->siz * tag) % p;
R->sum = (R->sum + R->siz * tag) % p;
L->tag = (L->tag + tag) % p; R->tag = (R->tag + tag) % p;
tag = 0;
} int mid() {return (l + r) >> 1;} }*root; void bulid(node *&o, int l, int r, bool bo) {
o = new node(l, r);
if(l == r) {
int Siz = siz[nfd[l]] % p;
if(bo) Siz = ((LL)Siz * siz[nfd[l]]) % p;
o->siz = Siz;
o->sum = (Siz * a[nfd[l]]) % p;
return;
}
bulid(o->L, l, o->mid(), bo);
bulid(o->R, o->mid() + 1, r, bo);
o->up();
} void chenge(node *o, int l, int r, int tag) {
if(o->l >= l && o->r <= r) {
o->sum = (o->sum + o->siz * tag) % p;
o->tag = (o->tag + tag) % p;
return;
}
if(o->tag) o->down();
if(o->mid() >= l) chenge(o->L, l, r, tag);
if(o->mid() < r) chenge(o->R, l, r, tag);
o->up();
} int ask(node *o, int l, int r) {
if(o->l >= l && o->r <= r) return o->sum;
if(o->tag) o->down();
int ans = 0;
if(o->mid() >= l) ans += ask(o->L, l, r);
if(o->mid() < r) ans += ask(o->R, l, r);
return ans % p;
} }wfx, cdy; void dfs(int x, int f) {
dep[x] = dep[fa[x] = f] + (siz[x] = 1);
for(Edge *i = head[x]; i; i = i->nxt) {
if(i->to == f) continue;
dfs(i->to, x);
siz[x] += siz[i->to];
if(siz[i->to] > siz[son[x]]) son[x] = i->to;
}
} void dfs(int x, int f, int top) {
Top[nfd[dfn[x] = ++tot] = x] = top;
if(son[x]) dfs(son[x], x, top);
for(Edge *i = head[x]; i; i = i->nxt) {
if(i->to == son[x] || i->to == f) continue;
dfs(i->to, x, i->to);
}
} void chenge(int x, int y, int d) {
if(x == y) return;
while(Top[x] != Top[y]) {
if(dep[Top[x]] < dep[Top[y]]) std::swap(x, y);
wfx.chenge(wfx.root, dfn[Top[x]], dfn[x], d);
cdy.chenge(cdy.root, dfn[Top[x]], dfn[x], d);
x = fa[Top[x]];
}
if(x == y) return;
if(dep[x] > dep[y]) std::swap(x, y);
wfx.chenge(wfx.root, dfn[x] + 1, dfn[y], d), cdy.chenge(cdy.root, dfn[x] + 1, dfn[y], d);
} int ask(int k) {
if(siz[k] <= 1) return 0;
int a = wfx.ask(wfx.root, dfn[k] + 1, dfn[k] + siz[k] - 1), b = cdy.ask(cdy.root, dfn[k] + 1, dfn[k] + siz[k] - 1);
return ((siz[k] * a - b) % p + p) % p;
} int main()
{
#ifdef yilnr
#else
freopen("network.in","r",stdin);
freopen("network.out","w",stdout);
#endif
int n = read(), m = read();
for(int i = 2; i <= n; ++i) {
int x = read(); a[i] = read();
head[x] = new Edge (i, head[x]);
}
dfs(1, 0);
dfs(1, 0, 1);
wfx.bulid(wfx.root, 1, n, false);
cdy.bulid(cdy.root, 1, n, true);
for(; m --> 0; ) {
int opt = getc(), l, r, d;
if(opt) l = read(), r = read(), d = read(), chenge(l, r, d);
else printf("%d\n", ask(read()));
}
return 0;
}
【csp模拟赛1】铁路网络 (network.cpp)的更多相关文章
- CSP模拟赛游记
时间:2019.10.5 考试时间:100分钟(连正式考试时间的一半还没有到)题目:由于某些原因不能公开. 由于第一次接触NOIinux系统所以连怎么建文件夹,调字体,如何编译都不知道,考试的前半小时 ...
- 【csp模拟赛4】 珠江夜游 (cruise.cpp)-二分,贪心
Problem 1 珠江夜游 (cruise.cpp) [题目描述] 小 Z 放假后难得来一趟广州游玩,当然要吃遍广州各路美食小吃然后再 到珠江新城看看远近闻名的小蛮腰啦!可当小 Z 一路吃吃吃以后, ...
- 【csp模拟赛4】基站建设 (station.cpp)
[题目描述] 小 Z 的爸爸是一位通信工程师,他所在的通信公司最近接到了一个新的通 信工程建设任务,他们需要在 C 城建设一批新的基站. C 城的城市规划做得非常好,整个城市被规整地划分为 8 行 8 ...
- 【csp模拟赛5】限制 (restrict.cpp)--数学
自己看吧: 爆搜代码: //春水初涨-春林初盛-春风十里-不如你 //----hzwer // 这是啥子题,读不懂-- //题意有问题 -- #include<iostream> #inc ...
- 【csp模拟赛5】购物(shopping.cpp)--常规
多项式,因为每次的x相同,所以把a和b相加就行了,然后找对称轴,找离对称轴最近的整数点,然而我却写了个暴力,没看x #include <iostream> #include <cst ...
- 【csp模拟赛4】旅行计划 (travelling.cpp)--欧拉回路
[题目描述] 小 Z 打算趁着暑假,开启他的旅行计划.但与其他同学不同的是,小 Z 旅 行时并不关心到达了哪个网红景点打了哪些卡.小 Z 更关注沿路的风光,而且 小 Z 觉得,尽管多次到达同一个地方, ...
- 【csp模拟赛1】不服来战 (challenge.cpp)
[题目描述] 最近小 Z 和他的朋友都迷上了一款手机游戏:不服来战. 游戏的设定十分简单,在游戏开始时,会给出一排共 N 个灯,有的灯是开着 的有的是关着的,每个灯都有一个分数.而玩家可以进行任意次操 ...
- 【CSP模拟赛】Freda的迷宫(桥)
题目描述 Freda是一个迷宫爱好者,她利用业余时间建造了许多迷宫.每个迷宫都是由若干房间和走廊构成的,每条走廊都连接着两个不同的房间,两个房间之间最多只有一条走廊直接相连,走廊都是双向通过. 黄昏 ...
- HGOI20191114 CSP模拟赛 反思
Problem A 宇宙魔方 有一个$N \times N \times N$的魔方,每一次操作可以整体转动该魔方,也可以对于一层整体+X. 给出最后魔方的最终状态,其中有一个位置为-1.利用其它位置 ...
随机推荐
- Earth Wind and Fire CodeForces - 1148E (构造)
大意: $n$个石子, 第$i$个石子初始位置$s_i$, 每次操作选两个石子$i,j$, 要求$s_i<s_j$, 任取$d$, 满足$0\le 2d\le s_j-s_i$, 将$s_i,s ...
- docker启动mysql 自定义配置文件
命令行如下: docker run --name mysql56 -p : -v /home/mysql56/data:/var/lib/mysql -v /home/mysql56/conf:/et ...
- 按Excel的模板导出数据
思路:把Excel模板上传到服务器:然后读取文档,写入数据,另存为新的文件 然后就是导出的时候处理数据定义的字符串,按数据导出:注意读取的数据流要处理下,不然会报异常 public Stream Ge ...
- TCP面向字节流和UDP面向报文的区别
TCP面向字节流 打个比方比喻TCP,你家里有个蓄水池,你可以里面倒水,蓄水池上有个龙头,你可以通过龙头将水池里的水放出来,然后用各种各样的容器装(杯子.矿泉水瓶.锅碗瓢盆)接水. 上面的例子中,往水 ...
- 【Caffe学习笔记】一 、环境安装 Caffe + cuda + windows10 + VS2015 安装笔记, win7也适用
1. 下载cuda8.0 cudnn5 anaconda https://developer.nvidia.com/cuda-80-ga2-download-archive https://de ...
- vue-cli3 本地数据模拟后台接口
vue-cli3 本地数据模拟后台接口 原理: 将本地的json数据在前端模拟为后台接口,然后调用接口,完成前端操作.在后台接通后可以直接在api配置文件中修改路径,完成前后台对接. 配置: 1.文件 ...
- Javascript简单教程汇总
什么是函数 一段定义好的代码,并可以反复使用的代码块 函数的作用 提升代码的可复用性,将一段代码进行预定义,需要使用的时候才触发 代码块 形成了一个相对独立的作用域 语法: function 函数名 ...
- ZPL语言完成条形码的打印
近期因为项目的需求,需要使用到打印机来打印业务相关的条形码和其他信息,由于之前有操作其它打印机的经验,Leader就安排我来做这个了(凑哦,这能说我是懵逼的么).于是就开始了我的探索之旅啦,不对,是踩 ...
- Apache httpclient的execute方法调试
因为工作需要,想研究一下execute执行的逻辑. 在这一行调用execute: response = getHttpClient().execute(get); getHttpClient的实现: ...
- 8.Spring整合Hibernate_2_声明式的事务管理(Annotation的方式)
声明式的事务管理(AOP的主要用途之一) (Annotation的方式) 1.加入annotation.xsd 2.加入txManager bean 3.<tx:annotation-drive ...