[nowcoder 2020] 牛半仙的妹子Tree

最新推荐文章于 2021-04-01 22:01:31 发布

C202044zxy

最新推荐文章于 2021-04-01 22:01:31 发布

阅读量235

点赞数 1

分类专栏：树链剖分

本文链接：https://blog.csdn.net/C202044zxy/article/details/109249442

版权

树链剖分专栏收录该内容

15 篇文章 0 订阅

订阅专栏

一、题目

点此看题

二、解法

这题是个套路题，两种方法都是套路（我都见过但是不会应用，该反思！）

0x01 树链剖分

这其实是一个套路：如果有关键点和另一堆需要询问某些东西的点，那么可以用树链剖分把信息放到这个点到根的路径上维护
$d[x]+d[u]-2d[lca]\le t_u-t_x$ $d[x]+t_x-2d[lca]\le t_u-d[u]$ 我们实际上要找左式的最小值，如果最小值满足，那么 $u$ 就被污染了。

设 $f$ 为 $x, u$ 的公共祖先，那么可以发现，当 $f = l c a$ 时，左边的 $d[x]+t_x-2d[f]$ 取最小值，所以在线段树上我们维护上式的最小值即可，修改的时候区间修改 $d[x]+t_x$ 为一个统一的值，那么树链剖分就很容易解决了。

0x02 分块

也是套路：如果我们对一道题有两个不同暴力算法（主要是复杂度的计算方式不同），对我们的询问分类，看哪种方法计算他更优就用哪种方法。不需要高级算法，门槛低。

第一种暴力就是 $b f s$ ，如果一个点没有被访问过就在对应的时间扩展他，时间复杂度 $O (n)$ ，可以参考下暴力代码， $\tt because\space words\space failed\space me\space when\space I\space try\space to\space express\space myself$

while(!q.empty() && q.front().t<i)
{
	node t=q.front();q.pop();
	for(int i=f[t.u];i;i=e[i].next)
	{
		int v=e[i].v;
		if(vis[v]) continue;
		vis[v]=1;
		q.push(node{v,t.t+1});
	}
}

第二种暴力是 $s t$ 表暴力 $l c a$ ，设初始讨厌他的妹子数是 $t$ ，那么复杂度是 $O (t)$

对于两个 $2$ 操作之间，如果长度大于 $\sqrt m$ ，那么用第一种算法 $O(n\sqrt m)$ ，否则用第二种算法，每个询问只需要算 $\sqrt m$ 次，那么复杂度 $O(m\sqrt m)$ ，我们用最基础的暴力通过此题。

给出方法 $2$ 的代码（方法 $1$ 有人实践，放心写吧）

#include <cstdio>
#include <cstring>
#include <queue>
#include <map>
using namespace std;
const int M = 200005;
const int sqr = 300;
int read()
{
    int num=0,flag=1;char c;
    while((c=getchar())<'0'||c>'9')if(c=='-')flag=-1;
    while(c>='0'&&c<='9')num=(num<<3)+(num<<1)+(c^48),c=getchar();
    return num*flag;
}
int n,m,tot,cnt,f[M],dfn[M],log[M],dep[M],st[M][20];
int o[M],x[M],vis[M],last; 
struct edge
{
	int v,next;
}e[M];
struct node
{
	int u,t;
};
void dfs(int u,int fa)
{
	dfn[u]=++cnt;dep[u]=dep[fa]+1;st[cnt][0]=u;
	for(int i=f[u];i;i=e[i].next)
	{
		int v=e[i].v;
		if(v==fa) continue;
		dfs(v,u);
		cnt++;st[cnt][0]=u;
	}
}
int Min(int x,int y)
{
	if(dep[x]<dep[y]) return x;
	return y;
}
void init()
{
	for(int i=2;i<=cnt;i++)
		log[i]=log[i>>1]+1;
	for(int j=1;j<20;j++)
		for(int i=1;i+(1<<j)-1<=cnt;i++)
			st[i][j]=Min(st[i][j-1],st[i+(1<<j-1)][j-1]); 
}
int lca(int u,int v)
{
	int l=dfn[u],r=dfn[v];
	if(l>r) swap(l,r);
	int k=log[r-l+1];
	return Min(st[l][k],st[r-(1<<k)+1][k]);
}
void work1(int l,int r)
{
	memset(vis,0,sizeof vis);
	queue<node> q;
	for(int i=l;i<=r;i++)
	{
		while(!q.empty() && q.front().t<i)
		{
			node t=q.front();q.pop();
			for(int i=f[t.u];i;i=e[i].next)
			{
				int v=e[i].v;
				if(vis[v]) continue;
				vis[v]=1;
				q.push(node{v,t.t+1});
			}
		}
		if(o[i]==1 && !vis[x[i]]) vis[x[i]]=1,q.push(node{x[i],i});
		if(o[i]==3) puts(vis[x[i]]?"wrxcsd":"orzFsYo");
	}
}
void work2(int l,int r)
{
	for(int i=l;i<=r;i++)
	{
		if(o[i]==3)
		{
			for(int j=l;j<i;j++)
				if(o[j]==1)
				{
					int dis=dep[x[i]]+dep[x[j]]-2*dep[lca(x[i],x[j])];
					if(dis+j<=i)
					{
						puts("wrxcsd");
						goto In;
					}
				}
			puts("orzFsYo");
			In:;
		}
	}
}
signed main()
{
	//freopen("tgT3.in","r",stdin);
	//freopen("zxy.out","w",stdout);
	n=read();m=read();
	for(int i=1;i<n;i++)
	{
		int u=read(),v=read();
		e[++tot]=edge{v,f[u]},f[u]=tot;
		e[++tot]=edge{u,f[v]},f[v]=tot;
	}
	dfs(1,0);
	init();//每调用调了好久 nmsl 
	for(int i=1;i<=m;i++)
	{
		o[i]=read();x[i]=read();
		if(o[i]==2 || i==m)
		{
			if(i-last-1>sqr) work1(last+1,i);//bfs
			else work2(last+1,i);//st表
			last=i;
		}
	}
}