UVA10253 - Series-Parallel Networks（树形DP＋(组合) _ 建模很好）（好题）

题目链接：https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=1194

大致题意: 给你n条边，问你恰好用n条边，能构成几种串并联网络。（串联的各个部分可以任意调换，并联在一起的各个部分也可以任意调换，若通过调换可得，则二者视为等效）

思路：一个n条边的整体总能拆成一些子整体的串联或一些子整体的并联。 这些子整体仍然可以这样拆

每个不同的整体必然有且只有一种唯一的拆开方式。

这样就可以转化成在树上操作，最终拆成只有1条边的整体，就是叶子，然后问题就是一颗n个叶子的树有多少种构树方式，每个孩子节点位置等价，由于位置等价，所以可以按子树的叶子节点从小到大活从大到小DP，这样构树就不会重复，现在问题就是一些子树的包含的叶节点个数一样，这些子树在DP的时候要避免重复计数，只能枚举包含相同叶子的子树有m颗，然后转化成这样的可重组合问题：在每次在n个元素中选m个元素，元素可以重复选，答案就是C(n+m-1, m)

所以DP方程是(dp[x][y]表示含有x个叶子的树，它的每棵子树的叶子个数<=y 有多少种构成方式) 

dp[sum][maxn] = sigma_m  ( C(DP(maxn,maxn-1) - 1 + m, m) * DP(sum - m * maxn, maxn - 1)  )

//#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <iostream>
#include <cstring>
#include <cmath>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <string>
#include <vector>
#include <cstdio>
#include <ctime>
#include <bitset>
#include <algorithm>
#define SZ(x) ((int)(x).size())
#define ALL(v) (v).begin(), (v).end()
#define foreach(i, v) for (__typeof((v).begin()) i = (v).begin(); i != (v).end(); ++ i)
#define reveach(i, v) for (__typeof((v).rbegin()) i = (v).rbegin(); i != (v).rend(); ++ i)
#define REP(i,n) for ( int i=1; i<=int(n); i++ )
#define rep(i,n) for ( int i=0; i< int(n); i++ )
using namespace std;
typedef long long ll;
#define X first
#define Y second
#define PB push_back
#define MP make_pair
typedef pair<int,int> pii;

template <class T>
inline bool RD(T &ret) {
        char c; int sgn;
        if (c = getchar(), c == EOF) return 0;
        while (c != '-' && (c<'0' || c>'9')) c = getchar();
        sgn = (c == '-') ? -1 : 1 , ret = (c == '-') ? 0 : (c - '0');
        while (c = getchar(), c >= '0'&&c <= '9') ret = ret * 10 + (c - '0');
        ret *= sgn;
        return 1;
}
template <class T>
inline void PT(T x) {
        if (x < 0) putchar('-') ,x = -x;
        if (x > 9) PT(x / 10);
        putchar(x % 10 + '0');
}

const int N = 33;
ll dp[N][N]; // x is leaf , y is maxn
ll C(ll n,ll m) {
	ll ans = 1;
	for(int i = 1; i <= m; i ++ ) ans = ans * (n - m + i) / i ;
	return ans ;
}
ll DP(int n, int maxn) {
	ll &ans = dp[n][maxn];
	if( maxn == 1 ) return ans = 1;
   	if( ans != -1 ) return ans ;
	ans = 0;	
	for(int c = 0; c * maxn <= n; c ++ ) {
		ans += C(DP(maxn,maxn-1) - 1 + c, c) * DP(n - c * maxn, maxn - 1) ; 
	}
	return ans ;
}
int main() {
	int n ;
	while( RD(n) ) {
		if( n == 0 ) break;
		memset( dp, -1, sizeof(dp) );
		dp[1][1] = 1;
		if( n == 1 ) {
			puts("1");
		   	continue;
		}
		PT(2 * DP( n, n - 1 ));
		puts("");
	}
}