[agc009d]Uninity

题目大意

请你确定一种点分治策略，使得点分树的深度最小。
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);
}