[cf674E]Bear and Destroying Subtrees

令$f_{i,j}$表示以$i$为根的子树中，深度小于等于$j$的概率，那么$ans_{i}=sum_{j=1}^{dep}(f_{i,j}-f_{i,j-1})j$

大约来估计一下$f_{i,j}$的大小，较坏情况下是$lfloorfrac{n-1}{j} floor$个深度为$j$的节点（若选择有公共部分，虽然会增加节点数但并实际边的数量减少），即可以认为$f_{i,j}ge (1-frac{1}{2^{j}})^{lfloorfrac{n-1}{j} floor}$

其在$jge 60$时，可以认为$f_{i,j}=1$，代入$ans_{i}$的式子即$j>60$的部分不需要计算

因此此时状态数为$o(Dn)$（其中$D=60$），接下来考虑如何维护

如果对其暴力dp，转移为$f_{i,j}=frac{1}{2}f_{i,j}(1+f_{son,j-1})$，那么当插入节点的$k$，显然会修改$f_{fa_{k},0}$、$f_{fa_{fa_{k}},1}$等位置，因此插入复杂度为$o(D)$，查询时对该位置$f$求和同样为$o(D)$

总复杂度为$o(Dn)$，可以通过

#include<bits/stdc++.h>
using namespace std;
#define N 500005
#define D 50
vector<int>v;
int V,n,p,x,fa[N];
double ans,f[N][D+5];
void New(int k){
    fa[++V]=k;
    for(int i=0;i<=D;i++)f[V][i]=1;
}
int main(){
    scanf("%d",&n);
    New(0);
    for(int i=1;i<=n;i++){
        scanf("%d%d",&p,&x);
        if (p==1){
            New(x);
            for(int j=1,k=x;(j<=D)&&(k);j++,k=fa[k])v.push_back(k);
            while (v.size()){
                f[fa[v.back()]][v.size()]/=1+f[v.back()][(int)v.size()-1];
                v.pop_back();
            }
            f[x][0]/=2;
            for(int j=1,k=x;(j<=D)&&(k);j++,k=fa[k])f[fa[k]][j]*=1+f[k][j-1];
        }
        else{
            ans=0;
            for(int j=1;j<=D;j++)ans+=(f[x][j]-f[x][j-1])*j;
            printf("%.9f
",ans);
        }
    }
}