题目传送门:【BZOJ 4726】
本题为权限题,下面有详细题意
题目大意:公司有 n 个人,(除老板以外)每个人有且只有一个直接上司。公司里有一个叛徒(未给明)。若一个人的下属(直接或者间接,不包括他自己)中叛徒占的比例超过 x,他和他的所有下属都会变为叛徒。求出叛徒的个数不会超过 k 的情况下最小的 x。(本题为 SPJ)
输入的第一行包含两个正整数 n , k(1 ≤ k ≤ n ≤ 5*10
5
)。
接下来 n-1 行,第 i 行包含一个正整数 p
i+1
,表示 i+1 的父亲是 p
i+1
(1 ≤ p
i+1
≤ i)。
输出一行,为一个实数 x,x 的误差在 10 −6 以内都被认为是正确的。
9 3
1 1 2 2 2 3 7 3
样例输入 | 样例输出 |
---|---|
9 3 | 0.6666666667 |
1 1 2 2 2 3 7 3 |
HINT:答案中的 x 实际上是一个无限趋近于 2/3 但是小于 2/3 的数;
因为当 x 取 2/3 时,最坏情况下 3,7,8,9 都是叛徒,超过了k=3。
题目分析:
分析整个题意,我们可以发现,题目要求为最坏情况,那么我们可以贪心的思想证明叛徒肯定是某个叶子节点。
同时,每个节点的状态很可能和它的子树的状态有关。于是我们考虑 DP。
我们观察到每个节点的状态有两类:以该节点为根的子树没有全变为叛徒 / 已经全变为叛徒。对于前者,它由两种状态组成:
【1】某一个子树全是叛徒,但是在容忍度之内
【2】某个子树中有叛徒,但是它没有超过容忍度,因此该子树没有全变为叛徒
对于后者,它仅由一种状态组成:
自己的某一个子树中有叛徒,且叛徒比例大于容许比例
于是我们可以使用树形 DP 的思路秒掉这道题。具体的转移方程为:
dp
i
表示无法使i的子树中所有人都变成叛徒的最小x,可以列出方程
dp
i
= max ( dp
i
, min ( dp
j
,
sizej(sizei−1)
) )
(其中 j 是 i 的儿子,size 表示这个节点的子树大小,包括自己)
下面附上代码:
- #include<cstdio>
- #include<algorithm>
- using namespace std;
- const int MX=500005;
- struct Edge{
- int to,next;
- }edge[MX*2];
- int n,k,now=0,head[MX],siz[MX];
- bool leaf[MX];
- double dp[MX],min_d=0;
- //leaf[i]: 判断 i 号节点是否为叶子节点
- //dp[i]: 无法使 i 的子树全变成叛徒的最小 x 值
- inline void adde(int u,int v){
- edge[++now].to=v;
- edge[now].next=head[u];
- head[u]=now;
- }
- void dfs(int u){
- bool son=false;
- siz[u]=1;
- for (int i=head[u];i;i=edge[i].next){
- int v=edge[i].to;
- son=true;
- dfs(v);
- siz[u]+=siz[v];
- }
- if (son)
- for (int i=head[u];i;i=edge[i].next){
- int v=edge[i].to;
- dp[u]=max(dp[u],min(dp[v],1.0*siz[v]/(siz[u]-1)));
- }
- else
- dp[u]=1.0,leaf[u]=true;
- }
- int main(){
- int a;
- scanf(”%d%d”,&n,&k);
- for (int i=2;i<=n;i++){
- scanf(”%d”,&a);
- adde(a,i);
- }
- dfs(1);
- for (int i=1;i<=n;i++)
- if (siz[i]>k)
- min_d=max(dp[i],min_d);
- printf(”%0.10lf”,min_d);
- return 0;
- }