CF966E-May Holidays【虚树,分块】

正题

题目链接:https://codeforces.ml/contest/966/problem/E


题目大意

nnn个点的一棵树,每个点有一个tit_iti,每次修改一个点是否为关键点,每次修改完后要求有多少个点满足该点不是关键点且子树中关键点数量超过tit_iti


解题思路

对操作分块,对于每个块首我们处理出树的形状,那么显然在这个块中我们需要处理的点的个数级别为n\sqrt nn的。那么一个点修改后会对往上的节点造成影响,除了虚树上的点以外,还有一些隐在虚树边上的点。

对于虚树上每一个非关键点,我们把上面的点按照zi−tiz_i-t_iziti排序(ziz_izi表示点iii子树种关键点数量)。然后我们可以对于每条边维护一个指针,显然每次修改是指针的移动次数是常数级别的。

注意排序的时候用基数排序就可以做到O(nn)O(n\sqrt n)O(nn)的解决该题。


codecodecode

#pragma GCC optimize(2)
%:pragma GCC optimize(3)
%:pragma GCC optimize("Ofast")
%:pragma GCC optimize("inline")
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cctype>
#include<cmath>
#include<vector>
using namespace std;
const int N=1e5+10,T=20;
struct node{int to,next;
}a[N*2];
int n,m,num,cnt,tot,ans,q[N],ls[N],t[N],dep[N];
int f[N][T+1],ft[N],z[N],s[N],dfn[N],w[N],val[N];
int c[N*2],sa[N],l[N],r[N],mid[N],pos[N],rfn[N];
bool ins[N],v[N];
vector<int> p;
int read(){int x=0,f=1;char c=getchar();while(!isdigit(c)){if(c=='-')f=-f;c=getchar();}while(isdigit(c)){x=(x<<1)+(x<<3)+c-'0';c=getchar();}return x*f;
}
void addl(int x,int y){a[++tot].to=y;a[tot].next=ls[x];ls[x]=tot;return;
}
void dfs(int x){dfn[x]=++cnt;dep[x]=dep[f[x][0]]+1;for(int i=ls[x];i;i=a[i].next){int y=a[i].to;f[y][0]=x;dfs(y);}return;
}
void solve(int x){z[x]=v[x];for(int i=ls[x];i;i=a[i].next){int y=a[i].to;solve(y);z[x]+=z[y];}return;
}
bool cmp(int x,int y)
{return dfn[x]<dfn[y];}
int LCA(int x,int y){if(dep[x]>dep[y])swap(x,y);for(int i=T;i>=0;i--)if(dep[f[y][i]]>=dep[x])y=f[y][i];if(x==y)return y;for(int i=T;i>=0;i--)if(f[x][i]!=f[y][i])x=f[x][i],y=f[y][i];return f[x][0];
}
void Insert(int x,int y){//加入一条链 int now=y;l[y]=num;r[y]=num-1;ins[y]=ins[x]=1;w[y]=0;while(f[now][0]!=x){now=f[now][0];++num;if(!v[now])pos[now]=y;}ft[y]=x;return;
}
void Add(int x){if(!cnt){s[++cnt]=x;return;}int lca=LCA(x,s[cnt]);while(cnt>1&&dep[s[cnt-1]]>dep[lca])Insert(s[cnt-1],s[cnt]),cnt--;if(dep[s[cnt]]>dep[lca])Insert(lca,s[cnt]),cnt--;if((!cnt)||(s[cnt]!=lca))s[++cnt]=lca;s[++cnt]=x;return;
}
void Sort(){for(int i=1;i<=n*2;i++)c[i]=0;for(int i=1;i<=n;i++)c[z[i]+n]++;for(int i=1;i<=n*2;i++)c[i]+=c[i-1];for(int i=1;i<=n;i++)sa[c[z[i]+n]--]=i;return;
}
void Move(int x){if(x==1)return;while(mid[x]>l[x]&&rfn[mid[x]-1]+w[x]>0)mid[x]--,ans+=val[mid[x]];while(mid[x]<=r[x]&&rfn[mid[x]]+w[x]<=0)ans-=val[mid[x]],mid[x]++;return;
}
void Change(int x,int W){if(v[x]&&z[x]>1)ans++;if(!v[x]&&z[x]>0)ans--;v[x]=(W==1);w[x]+=W;z[x]+=W;Move(x);x=ft[x];while(x){w[x]+=W;z[x]+=W;Move(x);if(W==1&&!v[x]&&z[x]==1)ans++;if(W==-1&&!v[x]&&z[x]==0)ans--;x=ft[x];}return;
}
int main()
{n=read();m=read();for(int i=2;i<=n;i++)addl(read(),i);dep[1]=1;dfs(1);for(int j=1;j<=T;j++)for(int i=1;i<=n;i++)f[i][j]=f[f[i][j-1]][j-1];for(int i=1;i<=n;i++)t[i]=read();for(int i=1;i<=m;i++)q[i]=read();int Seq=sqrt(m)+1,L=0,R=0;while(1){L=R+1;R+=Seq;R=min(R,m);ans=cnt=num=0;solve(1);p.clear(); for(int i=1;i<=n;i++)z[i]-=t[i],ans+=((z[i]>0)&(!v[i]));for(int i=L;i<=R;i++)p.push_back(abs(q[i]));sort(p.begin(),p.end(),cmp);if(p[0]!=1)s[++cnt]=1;for(int i=0;i<p.size();i++)if((!i)||p[i]!=p[i-1])Add(p[i]);while(cnt>1)Insert(s[cnt-1],s[cnt]),cnt--;Sort();for(int i=1;i<=n;i++){//处理边上节点 if(!pos[sa[i]])continue;int x=sa[i];if(r[pos[x]]<l[pos[x]]||z[x]!=rfn[r[pos[x]]])rfn[++r[pos[x]]]=z[x],val[r[pos[x]]]=1;else val[r[pos[x]]]++;pos[sa[i]]=0;}for(int i=1;i<=n;i++){//处理虚树上节点 if(!ins[i])continue;mid[i]=r[i]+1;for(int j=l[i];j<=r[i];j++)if(rfn[j]>0){mid[i]=j;break;}ins[i]=0;} for(int i=L;i<=R;i++){Change(abs(q[i]),q[i]>0?1:-1);printf("%d ",ans);}if(R==m)break;}return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/320453.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

DFS序讲解

我们经常会遇到树的问题&#xff0c;但树是非线性的结构&#xff0c;操作起来始终还是麻烦&#xff0c;如果我们能把树改造成线性结构&#xff0c;有什么方法&#xff1f;对&#xff0c;就是今天要讲的DSF序&#xff1b; dfs序呢&#xff0c;就是把一棵树区间化&#xff0c;我们…

纪中C组模拟赛总结(2019.8.9)

成绩&#xff1a; 注&#xff1a;rankrankrank是有算其他dalaodalaodalao的 rankrankranknamenamenamescorescorescoreT1T1T1T2T2T2T3T3T3T4T4T4111lyflyflyf250250250100100100100100100202020303030777hkyhkyhky200200200100100100100100100000000424242cyzcyzcyz1201201201…

虚树(Virtual Tree)学习笔记

虚树(Virtual Tree)学习笔记 一道题目(BZOJ-2286消耗战) Description 在一场战争中&#xff0c;战场由n个岛屿和n-1个桥梁组成&#xff0c;保证每两个岛屿间有且仅有一条路径可达。现在&#xff0c;我军已经侦查到敌军的总部在编号为1的岛屿&#xff0c;而且他们已经没有足够多…

利用Asp.Net Core的MiddleWare思想处理复杂业务流程

最近利用Asp.Net Core 的MiddleWare思想对公司的古老代码进行重构&#xff0c;在这里把我的设计思路分享出来&#xff0c;希望对大家处理复杂的流程业务能有所帮助。背景一个流程初始化接口&#xff0c;接口中根据传入的流程类型&#xff0c;需要做一些不同的工作。1.有的工作是…

P4284-[SHOI2014]概率充电器【树形dp,换根法,数学期望】

正题 题目链接:https://www.luogu.com.cn/problem/P4284 题目大意 nnn个点的一棵树&#xff0c;每个点有pip_ipi​概率通电&#xff0c;每个边有一定概率可以导电。求期望有电的节点个数。 解题思路 发现考虑通电概率很难&#xff0c;所以我们可以考虑不通电的概率。 如果只…

【每日一题】4月7日题目精讲 树

树 时间限制&#xff1a;C/C 1秒&#xff0c;其他语言2秒 空间限制&#xff1a;C/C 131072K 其他语言262144K 64bit IO Format:%lld 题目描述 shy有一颗树&#xff0c;树有n个结点。有k种不同颜色的染料给树染色。一个染色方案是合法的&#xff0c;当且仅当对于所有相同颜色的点…

【贪心】Radar Installation(poj 1328)

Radar Installation poj 1328 题目大意&#xff1a; 在平面直角坐标系的一二象限上有n个小岛&#xff0c;现在让你在x坐标上布置雷达&#xff0c;每个雷达可以侦测以它为原心&#xff0c;半径为m的圆内的所有小岛&#xff0c;现在问侦测完这n个小岛最少要多少个雷达 输入样…

HDU5765 Bonds (高维前缀和)

HDU5765 Bonds (高维前缀和) 题意&#xff1a;\(n(n<20)\)个点\(m\)条边无向图&#xff0c;求每条边出现在多少个\(Bond\)里。一个图的\(cut\)指&#xff0c;对于一个图\(G\)的边集的某个子集\(E\)&#xff0c;如果删除\(E\)中的所有边&#xff0c;原图不连通。一个图的\(Bo…

F# 4.5提供Spans、Match!等特性

F# 4.5预览版现已发布&#xff0c;其中提供了一系列新特性&#xff0c;包括对.NET Core 2.1的新原生类型Span<T>的支持、新关键字Match!等。类型Span意在实现底层代码指针操作的安全性和可预测性&#xff0c;这可使得很多情况下不必再分配内存&#xff0c;进而改进了内存…

P4172-[WC2006]水管局长【LCT,最小生成树】

正题 题目链接:https://www.luogu.com.cn/problem/P4172 题目大意 nnn个点的一张图&#xff0c;有两个操作 询问两个点之间的一条路径使得最长边最短删除一条边 解题思路 显然这条边一定是在最小生成树上的&#xff0c;所以我们需要维护支持删边的最小生成树。 显然LCTLCTL…

牛客OI周赛15-普及组

1. 咪咪游戏 bool f1,f2,f; f1用于记录m f2用于记录q f用于判断串的好坏。1为坏&#xff0c;0为好 注意几种情况&#xff1a; mqm mmq 具体看代码&#xff1a; #include<bits/stdc.h> using namespace std; int main() {int n;cin>>n;string s;while(n--){bool f10…

【单调栈】向右看齐(luogu 2947)

向右看齐 luogu 2947 题目大意&#xff1a; 有n头奶牛&#xff0c;它们各有各的高度&#xff0c;现在要求出每一头奶牛右边第一头比他高的奶牛是那一头 输入样例 6 3 2 6 1 1 2 输出样例 3 3 0 6 6 0 输入说明 6 头奶牛的身高分别为 3, 2, 6, 1, 1, 2. 输出…

Codeforces 786B Legacy (线段树优化建图)

Codeforces 786B Legacy (线段树优化建图) 题意&#xff1a;\(n\)个点&#xff0c;有\(3\)种连边操作&#xff1a;1.将\(u\)指向\(v\)&#xff1b;2.将\(v\)指向编号在区间\([l,r]\)的点&#xff1b;3.将\([l,r]\)中的所有点指向\(v\) 做法&#xff1a;线段树优化建图。拓展一些…

Abp + Grpc 如何实现用户会话状态传递

0.背景在实际项目当中&#xff0c;我采用的是 Abp 框架&#xff0c;但是 Abp 框架官方并没有针对 Grpc 进行模块封装。基于此我结合 Abp 与 MagicOnion 封装了一个 Abp.Grpc 模块&#xff0c;它包括服务端和调用端两部分的包。通过这两个包&#xff0c;你可以很方便地在 Abp 框…

P4197-Peaks【Kruskal重构树,主席树】

正题 题目链接:https://www.luogu.com.cn/problem/P4197 题目大意 nnn个点的一张无向图&#xff0c;每个点有一个hih_ihi​&#xff0c;边有权值。 qqq次询问从vvv出发不走权值超过xxx的路径能到达的第kkk大hih_ihi​是多少。 解题思路 构一颗KruskalKruskalKruskal重构树之后…

多项式相关操作学习笔记

多项式相关操作学习笔记 标签&#xff1a; 多项式 说在前边 记录一下相关的多项式操作&#xff0c;顺便存个模板。&#xff08;多点求值之后的部分&#xff0c;有点写不动了。。。留坑留坑 多项式 定义 给定一个环\(R\)(\(R\)通常是交换环&#xff0c;可以是有理数、实数或者复…

恢复数列

题目链接 比赛链接 时间限制&#xff1a;C/C 1秒&#xff0c;其他语言2秒 空间限制&#xff1a;C/C 262144K&#xff0c;其他语言524288K Special Judge,64bit IO Format: %lld 题目描述 小y的数学作业不小心被泼上了墨水。有道题看不清了&#xff0c;现在他想请你帮他恢复这道…

【单调栈】Largest Rectangle in a Histogram(luogu-SP1805/poj 2559)

Largest Rectangle in a Histogram luogu-SP1805 poj 2559 题目大意&#xff1a; 有n个并排的矩阵&#xff0c;高度为aia_iai​&#xff0c;宽度为1&#xff0c;现在让求包含于这些矩阵的并集内部的最大矩阵的面积 输入样例 7 2 1 4 5 1 3 3 4 1000 1000 1000 1000 0输出…

【翻译】asp.net core中使用MediatR

这篇文章来自&#xff1a;https://ardalis.com/using-mediatr-in-aspnet-core-apps本文作为翻译&#xff0c;有一些单词翻译成中文可能会有一些误解&#xff08;对于读者&#xff09;或者错误&#xff08;对于作者&#xff09;的地方&#xff0c;所以在文章中你可以看到一些单词…

P2056-[ZJOI2007]捉迷藏【点分树,堆】

正题 题目链接:https://www.luogu.com.cn/problem/P2056 题目大意 nnn个点的一棵树&#xff0c;开始全是黑点&#xff0c;有操作 取反一个点的颜色求最远的黑点之间的距离 解题思路 根据点分治每个点和分散开来的重心连边&#xff0c;然后每个点往上只会有logloglog层节点。…