BZOJ4726: [POI2017]Sabota?

Description

某个公司有n个人, 上下级关系构成了一个有根树。其中有个人是叛徒(这个人不知道是谁)。对于一个人, 如果他

下属(直接或者间接, 不包括他自己)中叛徒占的比例超过x，那么这个人也会变成叛徒，并且他的所有下属都会变

成叛徒。你要求出一个最小的x，使得最坏情况下，叛徒的个数不会超过k。

Input

第一行包含两个正整数n,k(1<=k<=n<=500000)。

接下来n-1行，第i行包含一个正整数p[i+1]，表示i+1的父亲是p[i+1](1<=p[i+1]<=i)。

Output

输出一行一个实数x，误差在10^-6以内都被认为是正确的。

Sample Input

9 3
1
1
2
2
2
3
7
3

Sample Output

0.6666666667

HINT

答案中的x实际上是一个无限趋近于2/3但是小于2/3的数

因为当x取2/3时，最坏情况下3，7，8，9都是叛徒，超过了k=3。

Source

一开始以为是二分答案 然后

1716029	wxh010910	4726	Wrong_Answer	17156 kb	4 ms	C++/Edit	938 B	2016-11-26 17:48:46
1716022	wxh010910	4726	Time_Limit_Exceed	27096 kb	21832 ms	C++/Edit	1676 B	2016-11-26 17:34:16
1716020	wxh010910	4726	Wrong_Answer	19408 kb	12 ms	C++/Edit	1825 B	2016-11-26 17:33:11
1716006	wxh010910	4726	Time_Limit_Exceed	31004 kb	21192 ms	C++/Edit	1715 B	2016-11-26 17:22:24
1716000	wxh010910	4726	Time_Limit_Exceed	27092 kb	21508 ms	C++/Edit	1173 B	2016-11-26 17:10:53
1715998	wxh010910	4726	Time_Limit_Exceed	27092 kb	21108 ms	C++/Edit	1156 B	2016-11-26 17:07:33
1715997	wxh010910	4726	Wrong_Answer	27088 kb	2688 ms	C++/Edit	1156 B	2016-11-26 17:07:10
1715993	wxh010910	4726	Time_Limit_Exceed	27092 kb	21692 ms	C++/Edit	1060 B	2016-11-26 17:04:10

然后发现树形dp可做

dp[i]表示把i和子树全变成叛徒最小的不能满足的k min max注意一下就好了

No.	RunID	User	Memory	Time	Language	Code_Length	Submit_Time
1	1716032	wxh010910	36372 KB	7816 MS	C++	936 B	2016-11-26 17:55:26
2	1715312	ccz181078	29068 KB	11736 MS	C++	964 B	2016-11-25 19:19:02
3	1715543	Claris	6680 KB	18464 MS	C++	690 B	2016-11-26 00:43:06

[TOP]    [STATUS]

#include<bits/stdc++.h>
 
using namespace std;
 
const int maxn=500050;
 
struct edge
{
    int to,nxt;
}e[maxn<<1];
 
int head[maxn],cnt,n,K,size[maxn];
 
double dp[maxn];
 
inline void addedge(int x,int y)
{
    e[++cnt].to=y;
    e[cnt].nxt=head[x];
    head[x]=cnt;
}
 
void dfs(int x)
{
    size[x]++;
    bool flag=false;
    for(int i=head[x];i;i=e[i].nxt) dfs(e[i].to),size[x]+=size[e[i].to],flag=true;
    if(!flag) { dp[x]=1; return ; }
    for(int i=head[x];i;i=e[i].nxt)
    {
        int y=e[i].to;
        dp[x]=max(dp[x],min(dp[y],(double)size[y]/(double)(size[x]-1)));
    }
//  printf("%d %.7lf\n",x,dp[x]);
}
 
int main()
{
    scanf("%d%d",&n,&K);
    int x;
    for(int i=2;i<=n;i++)
        scanf("%d",&x),addedge(x,i);
    dfs(1);
    double ans=0;
    for(int i=1;i<=n;i++)
        if(size[i]>K&&ans<dp[i])
            ans=dp[i];
    printf("%.7lf",ans);
}