雅礼培训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\) 或者 \( ...
随机推荐
- [Rational Rose 2007]解决启动报”解决无法启动此程序因为丢失suite objects.dll“的问题
问题根源1:不是丢失suite objects.dll文件,而是环境变量配置错误或无配置 假如安装目录如:C:\Program Files\Rational 需要配置环境变量的路径为:C:\Progr ...
- CSS3制作的一款按钮特效
/*.btn { width:230px; height:70px; font-size:70px; font-weight:bold; overflow:hidden; font: "He ...
- ios push Payload
https://developer.apple.com/library/content/documentation/NetworkingInternet/Conceptual/RemoteNotifi ...
- elasticsearch最全详细使用教程:入门、索引管理、映射详解、索引别名、分词器、文档管理、路由、搜索详解
一.快速入门1. 查看集群的健康状况http://localhost:9200/_cat http://localhost:9200/_cat/health?v 说明:v是用来要求在结果中返回表头 状 ...
- 解决因为手机设置字体大小导致h5页面在webview中变形的BUG
首先,我们做了一个H5页面,在各种手机浏览器中打开都没问题.我们采用了rem单位进行布局,通过JS来动态计算网页的视窗宽度,动态设置html的font-size,一切都比较完美. 这时候,你自信满满的 ...
- currentStyle和getComputedStyle来获取外部样式
currentStyle和getComputedStyle来获取外部样式 通过document.getElementById(id).style.XXX就可以获取到XXX的值,但意外的是,这样做只能取 ...
- 一个batch如何通过一个网络
一个batch下所有的图片一起经过整个网络,不是说一张图片经过网络后再让下一张进入网络,这样一个batch一起通过网络计算速度比一张一张这样快
- Hopfield 网络(上)
讲的什么 这部分主要对 Hopfield 网络作一大概的介绍.写了其模型结构.能量函数和网络的动作方式.主要参考了网上搜到的一些相关 PPT. 概述 早在 1982 年,Hopfield 发表的文 ...
- Ukulele 常用和弦
- 【离线做法 树状数组】luoguP1972 [SDOI2009]HH的项链
与bzoj3585: mex的线段树做法有着异曲同工之妙 题目描述 HH 有一串由各种漂亮的贝壳组成的项链.HH 相信不同的贝壳会带来好运,所以每次散步完后,他都会随意取出一段贝壳,思考它们所表达的含 ...