UVA 10253 - Series-Parallel Networks(数论+计数问题+递推)

来源:互联网 发布:arcgis js 事件关闭 编辑:程序博客网 时间:2024/05/16 00:47

题目链接:10253 - Series-Parallel Networks

白书的例题。
这题也是需要把问题进行转化,一个并联可以分为几个串联,然后串联可以分成边。
如此一来,最后叶子结点种数会是n,问题转化为去分配叶子结点,使得总和为n。
书上有两种方法,一种直接去递归,利用组合数学的方式去计算答案。
一种是推出递推式:
设dp[i][j]为一共j个叶子结点的树,子树的叶子最多的为i个的情况。然后对于一颗树,枚举恰好包含i个叶子的子树为p棵,那么相当于从f[i]颗树中选出p棵树的方案数,是可重复选择的组合,组合数为:C(f[i] + p - 1, p)种,然后每种子树对应的情况数为dp[i - 1][j - p * i]
所以状态转移方程为dp[i][j] = sum{C(f[i] + p - 1, p) * d(i - 1, j - p * i)},最后答案为dp[n - 1][n],注意边界的设定,详见代码
代码:
#include <stdio.h>#include <string.h>#include <stdlib.h>#include <algorithm>using namespace std;long long C(long long n, long long m) {    double ans = 1;    for (int i = 0; i < m; i++) {        ans *= n - i;    }    for (int i = 0; i < m; i++)    ans /= i + 1;    return (long long)(ans + 0.5);}const int N = 31;int n;long long f[N];long long dp[N][N];int main() {f[1] = 1;    for (int i = 0; i < N; i++) dp[i][0] = 1;    for (int i = 1; i < N; i++) {dp[i][1] = 1; dp[0][i] = 0;}    for (int i = 1; i < N; i++) {        for (int j = 2; j < N; j++) {            dp[i][j] = 0;            for (int p = 0; p * i <= j; p++) {            dp[i][j] += C(f[i] + p - 1, p) * dp[i - 1][j - p * i];            }        }        f[i + 1] = dp[i][i + 1];    }    while (~scanf("%d", &n) && n) {    printf("%lld\n", n == 1 ? 1 : 2 * f[n]);    }    system("pause");    return 0;}


1 0