UVa 10253 Series-Parallel Networks

      《训练指南》中的第二种算法,其实本质上就是个背包。d[i][j]表示,在子树的节点数最大为i的情况下,j个节点的解。当之前的i-1,i-2,....0的结果都已知的时候,d[i][j]自然可根据下式求解:

d[i][j]=sum{C(f(i)+p-1,p)*d[i-1][j-p*i] | p*i<=j}

其中f(i)表示恰好有i个节点的子树的数量。而C(f(i)+p-1,p)则表示有p棵i节点子树形成的组合数。然后运用乘法原理,剩余的节点数为j-p*i且包含的子树最多只有i-1个节点。继而从i=2时开始遍历,将所有的d[i][j]求出,而f(i)=d[i-1][i]。最后的解f(n)自然等于d[n-1][n]。

        还要注意的一点是边界的问题。当d[i][j]中的j=0时,则d[i][0]=1,这与求解背包时的边界相似。d[0][1]=1,表示叶子节点。而其他的d[0][i]则都为0,因为有i(i>1)个节点,但没有子树的最大节点数为0,显然是不可能的。最后,在输出结果时,f(1)始终为1,但其他的f(i)要乘以2,因为正如书中所说,根节点为并联节点或串联节点将决定两个不同的序列。


注:关于组合数的求解为什么能这样做,希望有高手能解答!

#include <iostream>
#include <cstdio>
#define MAX 30+2
using namespace std;

long long d[MAX][MAX];
long long f[MAX];
long long C(long long n,long long m)//求解组合数C(n,m)
{
    double ans=1;
    for(int i=0;i<m;++i) ans*=n-i;
    for(int i=1;i<=m;++i) ans/=i;
    return (long long)(ans+0.5);

}

int main()
{
    //freopen("data.txt","r",stdin);
    for(int i=0;i<MAX;++i) d[i][0]=d[1][i]=1;//注意边界
    for(int i=2;i<MAX;++i) d[0][i]=0;
    d[0][1]=1;
    for(int i=2;i<MAX;++i){
        f[i]=d[i-1][i];
        for(int j=1;j<i;++j) d[i][j]=d[i-1][j];//最大子树的节点数大于总节点数,则继承
        for(int j=i;j<MAX;++j)
        for(int p=0;p*i<=j;++p)
            d[i][j]+=C(f[i]+p-1,p)*d[i-1][j-p*i];
    }
    f[1]=1;
    for(int i=2;i<MAX;++i) f[i]*=2;//修改f(i)
    
    int N;
    while(cin>>N&&N){
        cout<<f[N]<<endl;
    }

    return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值