bzoj 4726 [POI2017]Sabota?

题意

题解

二分答案x，然后用树形dp check即可
具体一点，d[i]表示i为根的子树下最多有多少个叛徒
放个代码

void dp(int u){
    if(sz[u]==1){
        d[u]=1;
        return ;
    }
    d[u]=1;
    for(int i=head[u];i;i=edge[i].nxt){
        int v=edge[i].v;
        dp(v);
        if(1.0*d[v]/(sz[u]-1)>X){
            d[u]=sz[u];
            break;
        }
        else
            d[u]=max(d[u],d[v]);
    }
}

就酱。
但 是 这 道 题 卡 常
所以各种地方减常数…比如这个递归版本的dp，我们可以改成迭代版本的，减少回溯时间。迭代版本也很好改，因为本身输入就是拓扑序…（注意观察输入格式）所以倒着扫一遍就相当于拓扑序倒序了
此外什么读入优化啊也要上起
二分次数也要限制
毒瘤啊

#include<cstdio>
#include<cstring>
#include<algorithm>
#define EPS 1e-7
using namespace std;
typedef long long ll;
const int N=500005;
int n,k;
int sz[N];
double X;
int fa[N],d[N];
bool check(){
    for(int i=n;i>=1;i--){
        d[i]=0;
    }
    for(int i=n;i>=1;i--){
        d[i]=max(d[i],1);
        if(d[i]!=sz[i]&&1.0*d[i]/(sz[i]-1)>X)
            d[i]=sz[i];
        d[fa[i]]=max(d[fa[i]],d[i]);
    }
    return d[1]<=k;
}
void Read(int &x){
    x=0;
    char c=getchar();
    while(c<'0'||c>'9')
        c=getchar();
    while(c>='0'&&c<='9')
        x=x*10+c-'0',c=getchar();
}
int main()
{
    Read(n),Read(k);
    //scanf("%d%d",&n,&k);
    for(int i=2;i<=n;i++){
        Read(fa[i]);
        //scanf("%d",&x);
    }
    for(int i=n;i>=1;i--){
        sz[i]++;
        sz[fa[i]]+=sz[i];
    }
    double l=0,r=1;
    for(int i=1;i<=25;i++){
        X=(l+r)/2;
        if(check())
            r=X;
        else
            l=X;
    }
    printf("%.7f\n",X);
}