准备复赛

题目描述
今年的NOIP初赛真是简单，小可可不用吹灰之力就考进了复赛，但是复赛可没有那么简单了，小可可想要好好准备复赛，争取复赛拿个省一。今天小可可在复习树和图的最大匹配时就碰到这样的一个难题：n个节点满足以下性质的不同的树有多少种。

1、树是有标号的，每个节点被标上1到n之间的整数；

2、每个节点最多和其他3个节点相连，但是1号节点最多和其他2个节点相连；

3、这棵树的最大匹配（把树看成二分图后的最大匹配）数为k。

两棵树被认为不同当且仅当存在两个点u、v，在一棵树中u、v之间有边，另一棵树中u、v之间没边。

由于答案可能很大，所以小可可让你输出答案模1000000007 (109 + 7)。

输入
第一行包含两个正整数n，k。

输出
包含一行，为方案数。

样例输入
4 2
样例输出
12
提示
【数据范围】

对于30%的数据，2≤n≤5；

对于60%的数据，2≤n≤20；

对于100%的数据，2≤n≤50。

题解
设f[i][j][0,1]表示用了i个数，然后有j个匹配，根节点是否被匹配。

答案为(f[n][k][0]+f[n][k][1])/n 因为只有1是根，而最后答案多乘了n，所以要除。

代码

#include<bits/stdc++.h>
#define mod 1000000007
#define inv 500000004 
#define ll long long
using namespace std;
ll f[55][55][2],c[55][55];
int n,m;
ll gpow(ll x,ll y){
    ll z=1;
    for(;y;y/=2,x=x*x%mod)if(y&1)z=z*x%mod;
    return z;
}
inline void pre()
{
    c[0][0]=1;
    for (int i=1;i<=n;i++)
    {
        c[i][0]=1;
        for (int j=1;j<=i;j++)
            c[i][j]=(c[i-1][j-1]+c[i-1][j])%mod;
    }
}
int main()
{
    scanf("%d%d",&n,&m);pre();
    f[1][0][0]=1;
    for (int i=2;i<=n;i++)
    {
        for (int j=0;j<=min(i/2,m);j++)
        {
            f[i][j][0]=i*f[i-1][j][1]%mod;
            if (j) f[i][j][1]=i*f[i-1][j-1][0]%mod;
            for (int x=1;x<i;x++)
                for (int k=0;k<=min(x/2,j);k++)
                {
                    int y=i-x-1,l=j-k;
                    if (y/2>=l) f[i][j][0]=(f[i][j][0]+f[x][k][1]*f[y][l][1]%mod*i%mod*inv%mod*c[i-1][x]%mod)%mod;
                    if (y/2<l-1) continue;
                    if (l>=1) f[i][j][1]=(f[i][j][1]+f[x][k][1]*f[y][l-1][0]%mod*i%mod*c[i-1][x]%mod+f[x][k][0]*f[y][l-1][0]%mod*i%mod*inv%mod*c[i-1][x]%mod)%mod;
                }
        }
    }
    printf("%d",(f[n][m][1]+f[n][m][0])%mod*gpow(n,mod-2)%mod);
    return 0;
}