雅礼培训4.3 Problem A 【点分治】
题目简述
一个\(N\)个节点的树,有\(M\)个炸弹分布在一些节点上,有各自的威力,随着其他点距离增大对其他点的伤害呈等差减小,直至为0
问每个点受到的伤害
题解
QAQ考场代码没处理好有些炸弹威力很大这个事实,,数组爆掉。。。
AC算法直接变暴力分,,,
点分治即可
我是每次将子树内所有的炸弹统计到根来,再用一个差分数组求出各个深度受到的伤害,累加入每个节点的答案
但是由于可能会出现伤害来自同一个子树的情况,我们再对每个子树做一遍撤销
常数略大,,经玄学优化强行卡入时限
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<vector>
#include<cmath>
#include<ctime>
#define LL long long int
#define REP(i,n) for (int i = 1; i <= (n); i++)
#define Redge(u) for (int k = h[u],to; k; k = ed[k].nxt)
using namespace std;
const int maxn = 500005,maxm = 600005,INF = 100000000;
inline int read(){
int out = 0; char c = getchar();
while (c < 48 || c > 57) c = getchar();
while (c >= 48 && c <= 57){out = (out << 3) + (out << 1) + c - 48; c = getchar();}
return out;
}
inline void write(LL x){
LL tmp = 0; int cnt = 0;
while (x) tmp = (tmp << 3) + (tmp << 1) + x % 10,x /= 10,cnt++;
while (cnt--) putchar(tmp % 10 + 48),tmp /= 10;
putchar('\n');
}
int h[maxn],ne = 2,n,m;
struct EDGE{int to,nxt;}ed[maxn * 2];
int F[maxn],Siz[maxn],rt,sum,vis[maxn];
int siz[maxn],st[maxm],top,d[maxn],fa[maxn],nd[maxn],ndi;
int md;
LL ans[maxn],D[maxn];
vector<int> power[maxn];
inline void build(int u,int v){
ed[ne] = (EDGE){v,h[u]}; h[u] = ne++;
ed[ne] = (EDGE){u,h[v]}; h[v] = ne++;
}
void getRT(int u){
Siz[u] = 1 + power[u].size();
F[u] = 0;
Redge(u) if (!vis[to = ed[k].to] && to != fa[u]){
fa[to] = u;
getRT(to);
Siz[u] += Siz[to];
F[u] = max(F[u],Siz[to]);
}
F[u] = max(F[u],sum - Siz[u]);
if (F[u] < F[rt]) rt = u;
}
void dfs1(int u){
Siz[u] = 1 + power[u].size(); siz[u] = 1;
nd[++ndi] = u;
md = max(md,d[u]);
Redge(u) if (!vis[to = ed[k].to] && to != fa[u]){
fa[to] = u; d[to] = d[u] + 1;
dfs1(to);
Siz[u] += Siz[to];
siz[u] += siz[to];
}
}
void dfs4(int u){
nd[++ndi] = u;
for (int j = 0; j < power[u].size(); j++){
if (power[u][j] > d[u] + 1) st[++top] = power[u][j] - d[u] - 1;
}
Redge(u) if (!vis[to = ed[k].to] && to != fa[u]){
dfs4(to);
}
}
void dfs5(int u){
ans[u] -= D[d[u]];
Redge(u) if (!vis[to = ed[k].to] && to != fa[u]){
dfs5(to);
}
}
void cal(int u){
top = 0;
for (int j = 0; j < power[u].size(); j++)
if (power[u][j] > 2) st[++top] = power[u][j] - 2;
for (int i = 0; i <= siz[u]; i++) D[i] = 0;
ndi = 0;
Redge(u) if (!vis[to = ed[k].to] && to != fa[u]){
dfs4(to);
}
for (int i = 1; i <= top; i++){
if (st[i] <= md){
D[1] += st[i];
D[2] += -1 - st[i];
D[st[i] + 2] += 1;
}
else {
D[1] += st[i];
D[2] += -1 - st[i];
}
}
for (int i = 1; i <= siz[u]; i++) D[i] += D[i - 1];
for (int i = 1; i <= siz[u]; i++) D[i] += D[i - 1];
ans[u] -= D[1];
for (int i = 1; i <= ndi; i++) ans[nd[i]] -= D[d[nd[i]]];
}
void solve(int u){
vis[u] = true; siz[u] = 1; Siz[u] = 1 + power[u].size();
ndi = 0; md = 0;
Redge(u) if (!vis[to = ed[k].to]){
fa[to] = u; d[to] = 1;
dfs1(to);
siz[u] += siz[to];
Siz[u] += Siz[to];
}
md += 2;
if (Siz[u] == siz[u]) return;
top = 0;
for (int j = 0; j < power[u].size(); j++) st[++top] = power[u][j];
for (int i = 1; i <= ndi; i++){
int v = nd[i];
for (int j = 0; j < power[v].size(); j++){
if (power[v][j] > d[v]) st[++top] = power[v][j] - d[v];
}
}
for (int i = 0; i <= siz[u]; i++) D[i] = 0;
for (int i = 1; i <= top; i++){
if (st[i] <= md){
D[0] += st[i];
D[1] += -1 - st[i];
D[st[i] + 1] += 1;
}
else {
D[0] += st[i];
D[1] += -1 - st[i];
}
}
for (int i = 1; i <= siz[u]; i++) D[i] += D[i - 1];
for (int i = 1; i <= siz[u]; i++) D[i] += D[i - 1];
ans[u] += D[0];
for (int i = 1; i <= ndi; i++) ans[nd[i]] += D[d[nd[i]]];
Redge(u) if (!vis[to = ed[k].to]){
if (siz[to] == Siz[to]) continue;
cal(to);
}
Redge(u) if (!vis[to = ed[k].to]){
if (siz[to] == Siz[to]) continue;
sum = Siz[to]; F[rt = 0] = INF;
getRT(to);
solve(rt);
}
}
int main(){
//double t = clock();
//freopen("1.in","r",stdin);
//freopen("a.out","w",stdout);
n = read(); m = read();
int pos,w;
for (int i = 2; i <= n; i++) build(i,read());
for (int i = 1; i <= m; i++){
pos = read(); w = read();
power[pos].push_back(w);
}
F[rt = 0] = INF; sum = n + m;
getRT(1);
solve(rt);
for (int i = 1; i <= n; i++) write(ans[i]);
//cerr << (clock() - t) / CLOCKS_PER_SEC << 's' << endl;
return 0;
}
雅礼培训4.3 Problem A 【点分治】的更多相关文章
- 雅礼培训 Problem B 【图论 + 贪心】
题意 A和B在树上轮流选点,记A的联通块个数为\(x\),B的联通块个数为\(y\) A使\(x - y\)最大,B使\(x - y\) 二人采取最优策略,求\(x-y\) 题解 树联通块个数 = 点 ...
- 雅礼培训 Problem A 【线段树】
题意 维护一段区间,支持求区间最大值,区间且,区间或 \(n,q<=2*10^5\) 题解 我们用线段树维护区间最大值 对于and和or运算, and实质就是强行把一些位改为0 or实质就是强行 ...
- Loj#6503-「雅礼集训 2018 Day4」Magic【分治NTT】
正题 题目链接:https://loj.ac/p/6503 题目大意 \(n\)张卡\(m\)种,第\(i\)种卡有\(a_i\)张,求所有排列中有\(k\)对相邻且相同的卡牌. \(1\leq n\ ...
- LOJ_6045_「雅礼集训 2017 Day8」价 _最小割
LOJ_6045_「雅礼集训 2017 Day8」价 _最小割 描述: 有$n$种减肥药,$n$种药材,每种减肥药有一些对应的药材和一个收益. 假设选择吃下$K$种减肥药,那么需要这$K$种减肥药包含 ...
- 雅礼集训【Day6-1】字符串
雅礼集训[Day6-1]字符串 假设我们有串\(a\),我们设\(a'\)为\(a\)翻转后按为取反过后的串. 我们只考虑前一半的,长为\(m\)的串.如果前半截匹配了\(a\)或者\(a'\),则\ ...
- 「雅礼集训 2017 Day7」事情的相似度
「雅礼集训 2017 Day7」事情的相似度 题目链接 我们先将字符串建后缀自动机.然后对于两个前缀\([1,i]\),\([1,j]\),他们的最长公共后缀长度就是他们在\(fail\)树上对应节点 ...
- NOIP2018 模拟赛(二十二)雅礼NOI
Preface 这次的题目都是NOI+的题,所以大家的分数都有点惨烈. 依靠T1大力骗分水到Rank2 所以想看正解的话看这里吧 A. 「雅礼NOI2018模拟赛(一) Day1」树 看一眼题目感觉十 ...
- 「雅礼集训 2017 Day2」解题报告
「雅礼集训 2017 Day2」水箱 我怎么知道这种题目都能构造树形结构. 根据高度构造一棵树,在树上倍增找到最大的小于约束条件高度的隔板,开一个 \(vector\) 记录一下,然后对于每个 \(v ...
- 「雅礼集训 2017 Day1」 解题报告
「雅礼集训 2017 Day1」市场 挺神仙的一题.涉及区间加.区间除.区间最小值和区间和.虽然标算就是暴力,但是复杂度是有保证的. 我们知道如果线段树上的一个结点,\(max=min\) 或者 \( ...
随机推荐
- 【NumPy学习指南】day5 改变数组的维度
我们已经学习了怎样使用reshape函数,现在来学习一下怎样将数组展平. (1) ravel 我们可以用ravel函数完成展平的操作: In: b Out: array([[[ 0, 1, 2, 3] ...
- 关于sqlserver帐号被禁用问题
若发现sqlsrver所有帐号不小心被禁用了,这个时候怎么办?用重装吗?不用,仔细看小白是怎么一步一步解开这个谜题的.首先需要Windows帐号设置里重新添加一个新帐号.并将其添加到管理员组里面,然后 ...
- 洛谷 P2292 [HNOI2004]L语言
题目描述 标点符号的出现晚于文字的出现,所以以前的语言都是没有标点的.现在你要处理的就是一段没有标点的文章. 一段文章T是由若干小写字母构成.一个单词W也是由若干小写字母构成.一个字典D是若干个单词的 ...
- 洛谷 P1744 采购特价商品
题目背景 <爱与愁的故事第三弹·shopping>第一章. 题目描述 中山路店山店海,成了购物狂爱与愁大神的“不归之路”.中山路上有n(n<=100)家店,每家店的坐标均在-1000 ...
- 安卓统一推送联盟融云成唯一IM云服务企业
10月16日,安卓统一推送联盟在北京正式成立,来自中国信息通信研究院,华为.小米.OPPO等手机厂商,BAT等互联网巨头公司等75家机构及企业代表参加了联盟成立大会,融云也受邀参会并成为首批成员单位中 ...
- mysql利用binlog恢复数据详细例子
模拟数据恢复的案例 有些时候脑瓜就会短路,难免会出错 场景:在生产环境中,我们搭建了mysql主从,备份操作都是在从备份数据库上 前提:有最近一天或者最近的全备 或者最近一天相关数据库的备份 最重要的 ...
- 浏览器输入一个url到整个页面显示出来经历了哪些过程?
https://cloud.tencent.com/developer/article/1396399 https://www.cnblogs.com/haonanZhang/p/6362233.ht ...
- oracle数据比对工具
上半年的工作重心主要是机房搬迁,免不了要经常的数据比对,保证主备库数据一致,为了节约工作时间,提高工作效率,开发了这个数据比对小工具.用起来还可以.有需要的QQ私聊(1603039990),方便大家, ...
- Java Miniui实现批量上传文件demo 201906221520
可能需要的jar包: 需要miniui(类似easyui). Test2019062201.jsp <%@ page language="java" contentType= ...
- No-12.函数进阶
函数进阶 目标 函数参数和返回值的作用 函数的返回值 进阶 函数的参数 进阶 递归函数 01. 函数参数和返回值的作用 函数根据 有没有参数 以及 有没有返回值,可以 相互组合,一共有 4 种 组合形 ...