BZOJ4726: [POI2017]Sabota?

一开始看这题…二分好像很兹瓷啊，虽然n有50w，但是给了20s啊，不虚啊

然后T了…

emmmmmmmmmmmm
所以还是要dp
经过观察可以发现，对于最坏的情况一定是叶子开始叛变，且如果一个人没有被他的下属们策反，他的祖先们一定也不会叛变，所以最坏情况下，叛变的点一定是一棵完整的子树
所以我们令f[i]表示以i为根的子树不叛变的情况下最小的x
得 $f[i]=max[min(f[j],siz[j]/(siz[i]-1))]$
最终要使所有大小>k的子树全部不叛变，所以对这些子树的f[i]取max

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
#define inf 1e9
using namespace std;
 
const int maxn = 510000;
 
int n,k;
struct edge
{
    int y,nex;
    edge(){}
    edge(const int _y,const int _nex){y=_y;nex=_nex;}
}a[maxn<<1]; int len,fir[maxn];
inline void ins(const int x,const int y){a[++len]=edge(y,fir[x]);fir[x]=len;}
 
long double f[maxn];
int siz[maxn];
void dfs(const int x)
{
    siz[x]=0;
    for(int k=fir[x];k;k=a[k].nex)
        dfs(a[k].y),siz[x]+=siz[a[k].y];
    if(!siz[x]) f[x]=1.0;
    else
    {
        f[x]=0.0;
        for(int k=fir[x];k;k=a[k].nex)
            f[x]=max(f[x],min((long double)siz[a[k].y]/siz[x],f[a[k].y]));
    }
    siz[x]++;
    return;
}
 
int main()
{
    scanf("%d%d",&n,&k);
    for(int i=2;i<=n;i++)
    {
        int x; scanf("%d",&x);
        ins(x,i);
    }
    dfs(1);
    long double ans=0;
    for(int i=1;i<=n;i++) if(siz[i]>k)
        if(f[i]>ans) ans=f[i];
    printf("%.10Lf\n",ans);
     
    return 0;
}