模拟题 疯狂bb的zy prufer序列 DP

题目大意:给定n个点的度数上限,求出由1~n个点能构成的树有多少种。n<=100
这道题题解给的其实是一个非常恶心的生成函数+NTT,时间复杂度n^3logn,但是有另一种n^4的DP并且由于小常成功虐了正解,这里我们给出代码复杂度小的多的第二种题解。
首先我们要引入一个新的知识点叫prufer序列,每一棵树都唯一对应着一个prufer序列,每一个prufer序列也都唯一对应着一棵树,具体操作就是每次找到树上度数为一的点中标号最小的一个,然后将其连接的点加入prufer序列,并在原树中删掉这个点,直到树上只剩两个点时停止操作,得到的长度为n-2的序列即为prufer序列。prufer序列有一个性质就是一个点出现的次数即为这个点在树上的度数-1。
这样的话我们只需搞出一个DP算出来有多少种prufer序列即可。因为我们要算出1~n的所有答案,所以我们要维护一个状态表示选择了几个数,dp[i][j][k]表示前i种点选了j种一共k个点在prufer中的方案总数,组合数转移即可。这里要注意的是我可以选择一个点出现0次,这表示它的度数为1.最后第s个答案即为dp[n][s][s-2]

#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<iomanip>
#include<ctime>
#include<cmath>
#include<cstring>
#include<string>
#include<algorithm>
using namespace std;
const int mod=1004535809;
long long dp[101][101][101];
long long c[101][101];
int a[101];
int main()
{
    for(int i=0;i<=100;i++) c[i][0]=1;
    for(int i=1;i<=100;i++)
        for(int j=1;j<=100;j++)
            c[i][j]=(c[i-1][j-1]+c[i-1][j])%mod;
    int n;
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        scanf("%d",&a[i]);
    dp[0][0][0]=1;
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
            for(int k=0;k<n;k++)
            {
                if(dp[i][j][k]==0) continue;
                dp[i+1][j][k]+=dp[i][j][k],dp[i+1][j][k]%=mod;
                for(int m=0;m<=a[i+1]-1 && k+m<=n;m++)
                    dp[i+1][j+1][k+m]+=(dp[i][j][k]*c[k+m][m]),dp[i+1][j+1][k+m]%=mod;
            }
    for(int i=1;i<=n;i++)
    {
        if(i==1) cout<<n;
        else printf(" %lld",dp[n][i][i-2]);
    }
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值