bzoj4003城池攻占 左偏树

4003: [JLOI2015]城池攻占

Time Limit: 20 Sec  Memory Limit: 128 MB
Submit: 1312  Solved: 491
[Submit][Status][Discuss]

Description

小铭铭最近获得了一副新的桌游，游戏中需要用 m 个骑士攻占 n 个城池。

这 n 个城池用 1 到 n 的整数表示。除 1 号城池外，城池 i 会受到另一座城池 fi 的管辖，

其中 fi <i。也就是说，所有城池构成了一棵有根树。这 m 个骑士用 1 到 m 的整数表示，其

中第 i 个骑士的初始战斗力为 si，第一个攻击的城池为 ci。

每个城池有一个防御值 hi，如果一个骑士的战斗力大于等于城池的生命值，那么骑士就可

以占领这座城池；否则占领失败，骑士将在这座城池牺牲。占领一个城池以后，骑士的战斗力

将发生变化，然后继续攻击管辖这座城池的城池，直到占领 1 号城池，或牺牲为止。

除 1 号城池外，每个城池 i 会给出一个战斗力变化参数 ai;vi。若 ai =0，攻占城池 i 以后骑士战斗力会增加 vi；若 ai =1，攻占城池 i 以后，战斗力会乘以 vi。注意每个骑士是单独计算的。也就是说一个骑士攻击一座城池，不管结果如何，均不会影响其他骑士攻击这座城池的结果。

现在的问题是，对于每个城池，输出有多少个骑士在这里牺牲；对于每个骑士，输出他攻占的城池数量。

Input

第 1 行包含两个正整数 n;m，表示城池的数量和骑士的数量。

第 2 行包含 n 个整数，其中第 i 个数为 hi，表示城池 i 的防御值。

第 3 到 n +1 行，每行包含三个整数。其中第 i +1 行的三个数为 fi;ai;vi，分别表示管辖

这座城池的城池编号和两个战斗力变化参数。

第 n +2 到 n + m +1 行，每行包含两个整数。其中第 n + i 行的两个数为 si;ci，分别表

示初始战斗力和第一个攻击的城池。

Output

 输出 n + m 行，每行包含一个非负整数。其中前 n 行分别表示在城池 1 到 n 牺牲的骑士

数量，后 m 行分别表示骑士 1 到 m 攻占的城池数量。

Sample Input

5 5
50 20 10 10 30
1 1 2
2 0 5
2 0 -10
1 0 10
20 2
10 3
40 4
20 4
35 5

Sample Output

2
2
0
0
0
1
1
3
1
1

HINT

 对于 100% 的数据，1 <= n;m <= 300000; 1 <= fi<i; 1 <= ci <= n; -10^18 <= hi,vi,si <= 10^18;ai等于1或者2；当 ai =1 时，vi > 0；保证任何时候骑士战斗力值的绝对值不超过 10^18。

Source

带延迟标记的左偏树
还是自底向上合并一波。
以骑士的战斗力为关键字（小根堆），额外维护延迟标记。弹栈的时候在这个节点把死去骑士的数目++
由于骑士的路径一定是一条链，所以弹出骑士之后直接用当前起始深度减去当前深度即可（sb的我之前还多维护了一个骑士的攻战城池数目，慢得要死）
记得最后要把vectory的骑士弹出去。
维护延迟标记的时候记得要先乘后加，乘的时候加的延迟标记也要乘（和线段树区间修改一样，具体看代码）
常数好大，luogu的时限没过去，科科%>_<%
听说300000点会爆栈，然而我深搜过了。

#include<iostream>#include<cstdlib>#include<cstdio>#include<cstring>#include<algorithm>#include<map>#include<cmath>#define maxn 330000using namespace std;long long read(){    char ch = getchar(); long long x = 0, f = 1;    while(ch < '0' || ch > '9') {if(ch == '-') f = -1; ch = getchar();}    while(ch >= '0' && ch <= '9') {x = x * 10 + ch - '0'; ch = getchar();}    return x * f;}int pre[maxn], top;struct edge {    int to, next;}e[maxn * 2];void add(int u, int v) {    e[top].to = v;    e[top].next = pre[u];    pre[u] = top++;}long long val[maxn], d[maxn], tagm[maxn], taga[maxn], w[maxn];int h[maxn], ls[maxn], rs[maxn], deep[maxn], s[maxn];int sz, n, m, f[maxn], root[maxn], die[maxn], ans[maxn];bool a[maxn];void newnode(long long v, int p) {val[p] = v; h[p] = 1; tagm[p] = 1;}void paint_and(int p, long long v) {val[p] += v; taga[p] += v;}void paint_mult(int p, long long v) {val[p] *= v; tagm[p] *= v; taga[p] *= v;}void paint(int p, bool _, bool __, long long vm, long long va) {if(!p) return;if(_) paint_mult(p, vm);if(__) paint_and(p, va);}void pushdown(int p) {paint(ls[p], 1, 1, tagm[p], taga[p]);paint(rs[p], 1, 1, tagm[p], taga[p]);tagm[p] = 1; taga[p] = 0;}int merge(int u, int v) {if(!u || !v) return u | v;pushdown(u); pushdown(v);if(val[u] > val[v]) swap(u, v);rs[u] = merge(rs[u], v);if(h[rs[u]] > h[ls[u]]) swap(ls[u], rs[u]);h[u] = h[rs[u]] + 1;return u;}void init() {memset(pre, -1, sizeof(pre)); top = 0;n = read(); m = read();for(int i = 1;i <= n; ++i) d[i] = read();for(int i = 2;i <= n; ++i) {f[i] = read(); a[i] = read(); w[i] = read();add(f[i], i);}for(int i = 1;i <= m; ++i) {newnode(read(), i); s[i] = read(); root[s[i]] = merge(root[s[i]], i);}}void work(int p) {deep[p] = deep[f[p]] + 1;for(int i = pre[p]; ~i; i = e[i].next) {work(e[i].to);root[p] = merge(root[p], root[e[i].to]);} while(root[p] && val[root[p]] < d[p]) {++die[p]; ans[root[p]] = deep[s[root[p]]] - deep[p]; pushdown(root[p]);root[p] = merge(ls[root[p]], rs[root[p]]);}paint(root[p], a[p], a[p] ^ 1, w[p], w[p]);}void print() {while(root[1]) {ans[root[1]] = deep[s[root[1]]]; pushdown(root[1]);root[1] = merge(ls[root[1]], rs[root[1]]);}for(int i = 1;i <= n; ++i) printf("%d\n", die[i]);for(int i = 1;i <= m; ++i) printf("%d\n", ans[i]);}int main(){init();work(1);print();return 0;}/*5 550 20 10 10 301 1 22 1 52 1 101 0 1020 210 340 420 435 5*/