暴力美学——分块

1

CF1017G
给定一棵树，每次操作一个点，如果这个点是白色就把他染黑，否则递归操作所有儿子。或者把一个子树染白。或者询问某个点的颜色（黑白）。

把操作以 $O(\sqrt{n})$ 划分成若干段，把每段涉及到的点建一棵类似虚树的东西，边权是两点之间的白点个数。修改操作在这棵树上改，复杂度 \sqrt n，每段处理完后 O(n) 修改到原树上。

还有一种巧妙的做法。总体想法是我们想要对修改点打标记。问题是如何设计一个标记，使得一个点是黑点当且仅当满足某个条件。发现如果一开始给每个点打上-1的标记，操作一个点就把这个点+1，那么一个点是黑色当且仅当从他开始向根走的某条路径和大于等于0。因此我们树剖一下维护最大后缀和即可。注意子树染白操作，不能简单地打上-1标记，要注意修改子树根的权值使得它满足我们的规定。

2

CF925E
⼀棵 n 个点的有根黑白树，⼀个数组 a 和 m 次操作，每次将某个点反色。若⼀个节点 x 为黑色且其子树内白色节点数量大于 ax，则该节点是“不高兴”的。每次操作后输出“不高兴”的节点数目。

同样对修改分块，相邻的 B 个询问点建立一棵虚树，容易发现在当前块内，虚树上一条边 (u,v) 在原树中的点是一起被 u 的子树影响的，因此可以一起修改贡献。注意 -1(int)>1(unsigned int)，比较时要注意类型转换。

#include<bits/stdc++.h>
#define ll long long
#define pb push_back
#define fir first
#define sec second
#define ld long double
using namespace std;
const int N=100010;
typedef pair <int,int> P;
struct edge {
	int to,next;
}ed[N<<1];
int head[N],sz,qy[N],b[N],d[N],t[N],A[N],B[N],tim,st[N],cnt,fa[N],faa[N],f[N][18],deep[N],vis[N],tmp[N],ans=0,pq[N],lev[N];
vector <P> c[N];
void add_edge(int from,int to)
{
	ed[++sz].to=to;
	ed[sz].next=head[from];
	head[from]=sz;
}
int read()
{
	int x=0;char c=getchar(),flag='+';
	while(!isdigit(c)) flag=c,c=getchar();
	while(isdigit(c)) x=x*10+c-'0',c=getchar();
	return flag=='-'?-x:x;
}
bool cmp(int a,int b) {return A[a]<A[b];}
int Lca(int x,int y)
{
	if(deep[x]<deep[y]) swap(x,y);
	for(int j=17;j>=0;j--) if(deep[f[x][j]]>=deep[y]) x=f[x][j];
	if(x==y) return x;
	for(int j=17;j>=0;j--) if(f[x][j]^f[y][j]) x=f[x][j],y=f[y][j];
	return f[x][0];
}
void dfs0(int u,int ff)
{
	A[u]=++tim;
	fa[u]=ff;
	f[u][0]=ff;
	deep[u]=deep[ff]+1;
	for(int i=head[u];i;i=ed[i].next)
	{
		int v=ed[i].to;
		if(v==ff) continue;
		dfs0(v,u);
	}
	B[u]=++tim;
}
int dfs1(int u,int ff)
{
	int sum=d[u];
	for(int i=head[u];i;i=ed[i].next)
	{
		int v=ed[i].to;
		if(v==ff) continue;
		sum+=dfs1(v,u);
	}
	t[u]-=sum;
	return sum;
}
void solve()
{
	int m=cnt,top=0;
	cnt=0;
	for(int i=1;i<=m;i++) 
	{
		int x=abs(qy[i]);
		if(!vis[x]) b[++cnt]=x,vis[x]=1;
	}
	sort(b+1,b+cnt+1,cmp);
	for(int i=1;i<cnt;i++) 
	{
		int lca=Lca(b[i],b[i+1]);
		if(!vis[lca]) vis[lca]=2,b[++cnt]=lca;
	}
	sort(b+1,b+cnt+1,cmp);
	
	for(int i=1;i<=cnt;i++)
	{
		while(top&&(A[b[i]]<A[st[top]]||A[b[i]]>B[st[top]])) top--;
		faa[b[i]]=st[top];
		st[++top]=b[i];
	}
	
	for(int i=1;i<=cnt;i++)
	{
		int u=b[i],x=fa[u],tot=0;
		while(x^faa[u]) 
		{
			if(!lev[x]) tmp[++tot]=t[x];
			x=fa[x];
		}
		sort(tmp+1,tmp+tot+1);
		c[u].clear();
		pq[u]=-1;
		for(int j=1,p;j<=tot;j=p)
		{
			int num=0;
			for(p=j;tmp[p]==tmp[j]&&p<=tot;p++) num++;
			c[u].pb(P(tmp[j],num));
			if(tmp[j]<0) pq[u]=c[u].size()-1;
		}
	}
	for(int i=1;i<=m;i++)
	{
		int u=qy[i],k;
		if(u>0) //leave
		{
			k=1;
			lev[u]=1;
			if(t[u]<d[u]) ans--;
			d[u]++;
		}
		else //return
		{
			u=-u,k=-1;
			lev[u]=0;
			d[u]--;
			if(t[u]<d[u]) ans++;
		}
		int x=u;
		while(x) 
		{
			if(x!=u)
			{
				if(!lev[x]&&t[x]<d[x]) ans--;
				d[x]+=k;
				if(!lev[x]&&t[x]<d[x]) ans++;
			}
			if(pq[x]>=0&&pq[x]<c[x].size()&&d[x]<=c[x][pq[x]].fir) ans-=c[x][pq[x]].sec,pq[x]--;
			int tmp=c[x].size()-1;
			if(pq[x]<tmp) if(d[x]>c[x][pq[x]+1].fir) ans+=c[x][pq[x]+1].sec,pq[x]++;
			x=faa[x];
		}
		cout<<ans<<' ';
	}
	for(int i=1;i<=cnt;i++) d[b[i]]=0;
	for(int i=1;i<=m;i++) 
	{
		if(qy[i]>0) d[qy[i]]++;
		else d[-qy[i]]--;
	}
	dfs1(1,0);
	for(int i=1;i<=cnt;i++)
	{
		int x=b[i];
		d[x]=vis[x]=0;
	}
}
int main()
{
	int n=read(),m=read(),B=sqrt(n);
	for(int i=2;i<=n;i++) 
	{
		int x=read();
		add_edge(x,i),add_edge(i,x);
	}
	for(int i=1;i<=n;i++) t[i]=read();
	dfs0(1,0);
	for(int j=1;j<=17;j++)
		for(int i=1;i<=n;i++) f[i][j]=f[f[i][j-1]][j-1];
	for(int i=1;i<=m;i++)
	{
		qy[++cnt]=read();
		if(cnt==B||i==m) solve(),cnt=0;
	}
	return 0;
}
/*by DT_Kang*/