[bzoj1089]严格n元树

设f[i]表示深度不超过i的方案数，那么有f[0]=1,$f[i]=f[i-1]^{n}+1$，然后用高精度即可（注意深度恰好为d还要用f[d]-f[d-1]才是答案）

#include<bits/stdc++.h>
using namespace std;
struct ji{
    int a[1005];
}a,s,ans;
int n,d;
void jia(){
    for(int i=1;i<=ans.a[0]+1;i++)
        if (ans.a[i]==9)ans.a[i]=0;
        else{
            ans.a[i]++;
            ans.a[0]=max(ans.a[0],i);
            break;
        }
}
void jian(){
    for(int i=1;i<=ans.a[0];i++){
        if (s.a[i]>ans.a[i]){
            ans.a[i+1]--;
            ans.a[i]+=10;
        }
        ans.a[i]-=s.a[i];
    }
    while ((ans.a[0]>1)&&(!ans.a[ans.a[0]]))ans.a[0]--;
}
void cheng(){
    memset(a.a,0,sizeof(a.a));
    a.a[0]=ans.a[0]+s.a[0]-1;
    for(int i=1;i<=s.a[0];i++)
        for(int j=1;j<=ans.a[0];j++)a.a[i+j-1]+=s.a[i]*ans.a[j];
    ans=a;
    for(int i=2;i<=ans.a[0];i++){
        ans.a[i]+=ans.a[i-1]/10;
        ans.a[i-1]%=10;
    }
    while (ans.a[ans.a[0]]>9){
        ans.a[ans.a[0]+1]=ans.a[ans.a[0]]/10;
        ans.a[ans.a[0]++]%=10;
    }
}
int main(){
    scanf("%d%d",&n,&d);
    ans.a[0]=ans.a[1]=1;
    for(int i=1;i<=d;i++){
        s=ans;
        for(int j=1;j<n;j++)cheng();
        jia();
    }
    jian();
    for(int i=ans.a[0];i;i--)printf("%d",ans.a[i]);
}