【Luogu】P2351吊灯（脑洞后模拟）

　　题目链接

　　这题要智商qwq。玩不来玩不来。

　　观察到（个P,能观察到的全都是dalao）x是解的充要条件是至少有n/x个节点的size是x的倍数。

　　证明请看这里

　　然后这题就变模拟了呀。

　　

#include<cstdio>
#include<cstring>
#include<cctype>
#include<algorithm>
#include<cstdlib>
#define maxn 2000010
using namespace std;
inline long long read(){
    long long num=,f=;
    char ch=getchar();
    while(!isdigit(ch)){
        if(ch=='-')    f=-;
        ch=getchar();
    }
    while(isdigit(ch)){
        num=num*+ch-'';
        ch=getchar();
    }
    return num*f;
}

int father[maxn];
int size[maxn];
int s[maxn],tot;

struct Edge{
    int next,to;
}edge[maxn];
int head[maxn],num;
inline void add(int from,int to){
    edge[++num]=(Edge){head[from],to};
    head[from]=num;
}

int d[maxn];
void dfs(int x,int fa){
    size[x]=;
    for(int i=head[x];i;i=edge[i].next){
        int to=edge[i].to;
        if(to==fa)    continue;
        dfs(to,x);
        size[x]+=size[to];
    }
    d[size[x]]++;
    return;
}

int main(){
    int n=read();
    for(int i=;i<=n;++i)
        if(n%i==)    s[++tot]=i;
    for(int i=;i<=n;++i){
        father[i]=read();
        add(father[i],i);
    }
    dfs(,);
    printf("Case #1:\n");
    for(int i=;i<=tot;++i){
        int now=;
        for(int j=s[i];j<=n;j+=s[i]){
            now+=d[j];
            if(now>=n/s[i]){
                printf("%d\n",s[i]);
                break;
            }
        }
    }
    for(int i=;i<=;++i){
        memset(d,,sizeof(d));
        for(int j=;j<=n;++j)    father[j]=(father[j]+)%(j-)+;
        printf("Case #%d:\n",i+);
        for(int j=;j<=n;++j)    size[j]=;
        for(int j=n;j>=;--j){
            size[father[j]]+=size[j];
            d[size[j]]++;
        }
        for(int j=;j<=tot;++j){
            int now=;
            for(int k=s[j];k<=n;k+=s[j]){
                now+=d[k];
                if(now>=n/s[j]){
                    printf("%d\n",s[j]);
                    break;
                }
            }
        }
    }
    return ;
}