题目大意
请你确定一种点分治策略,使得点分树的深度最小。
n<=1e5
解题思路
根据原题面,定义一个点x的“海胆值”f[x]为max(f[son])+1,son是x在点分树上的儿子。
根据一般点分治,我们可以发现f最大是O(logn)级别的。
还是不好做,再观察一下其他性质。
考虑两个点x,y,其中f[x]=f[y],那么在原树中,x到y的路径必然存在点u,满足f[u]>f[x]。
在原树中,考虑一个点x,他的两颗不同子树如果分别存在两个点a,b,f[a]=f[b],那么f[x]>f[a]。
我们按照这个思路,来设计算法。
设sat[x]表示x这颗子树有哪些f[y]使得y到x(y可以等于x)路径上没有比f[y]大的点。即一个二进制数,由于f的O(logn)级别,这个可以存下来。
那么我们可以贪心地根据x的子树确定f[x]。就用上面那个条件搞一搞就行了。
然后还要维护一下sat[x]。
时间复杂度O(n)
代码
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<set>
#include<map>
using namespace std;
#define fo(i,j,k) for(i=j;i<=k;i++)
#define fd(i,j,k) for(i=j;i>=k;i--)
#define cmax(a,b) (a=(a>b)?a:b)
#define cmin(a,b) (a=(a<b)?a:b)
typedef long long ll;
const int N=1e6+5,M=2e6+5,mo=1e9+7;
int n,ans,i,x,y;
int tt,nxt[N],b[N],fst[N];
void cr(int x,int y)
{
tt++;
b[tt]=y;
nxt[tt]=fst[x];
fst[x]=tt;
}
int f[N],sat[N];
void dfs(int x,int y)
{
int lwb=0;
for(int p=fst[x];p;p=nxt[p])
if (b[p]!=y)
{
dfs(b[p],x);
lwb|=sat[x]&sat[b[p]];
sat[x]|=sat[b[p]];
}
f[x]=0;
while ((1<<f[x])<=lwb||((1<<f[x])&sat[x])) f[x]++;
sat[x]=(sat[x]>>f[x]<<f[x])|(1<<f[x]);
cmax(ans,f[x]);
}
int main()
{
freopen("t11.in","r",stdin);
//freopen("t11.out","w",stdout);
scanf("%d",&n);
fo(i,1,n-1)
{
scanf("%d %d",&x,&y);
cr(x,y);
cr(y,x);
}
dfs(1,0);
printf("%d\n",ans);
}