bzoj4271 烷烃计数 dp

Description

---

求有n个碳原子的烷烃同分异构体个数，例如丁烷有两种同分异构体
$n\le 500$

Solution

---

好像还有生成函数的做法，太神仙了。。

本质上在求n个点无标号无根树的方案数。假设我们求出了n个点无标号有根树的答案，那么我们钦定重心作为树根，无根树就变成有根树了

对于求n个点无标号有根树的方案可以dp，设f[i,j]表示子树size为i，根的度数为j的方案数。

我们可以钦定儿子按照size从小到大排序，枚举最大子树的size为mx，枚举恰好k个儿子的size恰好为mx。由于我们要保证重心为根，因此mx不超过$\frac{n-1}{2}$

记S为$\sum _{j=0}^{m-1}f[mx][j]$，即儿子可选择的总的方案数，那么这一部分的贡献实际上就是将k个儿子各自染上S种颜色中的一种的方案数，等价于将k个无差别球放入S个无差别的盒子里允许空的方案数，即$\left(\genfrac{}{}{0ex}{}{S+k-1}{S-1}\right)=\left(\genfrac{}{}{0ex}{}{S+k-1}{k}\right)$

那么n为奇数时的答案就是${\sum }_{i=0}^{m}f[n][i]$，而当n为偶数时会出现两个重心的情况。记$S^{\prime }=\sum _{j=0}^{m-1}f[\frac{n}{2}][j]$，那么我们少算的方案数等价于将S‘种方案分配给2个子树，也就是$\left(\genfrac{}{}{0ex}{}{S^{\prime }+1}{2}\right)$

当然n为偶数的时候个人还有一个比较感性的理解，我们可以把连接两重心的边拆成一个点，令这个新节点作为重心，那么多出来的方案数就和上面的dp转移一样了

喜闻乐见要写高精度

Code

---

#include <stdio.h>
#include <string.h>
#include <algorithm>
#define rep(i,st,ed) for (int i=st;i<=ed;++i)
#define drp(i,st,ed) for (int i=st;i>=ed;--i)
#define fill(x,t) memset(x,t,sizeof(x))

typedef long long LL;
const int MOD=1000;
const int N=1015;
const int L=205;

char st[L];

struct num {
	int s[L],len;
	inline bool operator ==(num b) {
		num a=*this;
		if (a.len!=b.len) return 0;
		drp(i,a.len,1) if (a.s[i]!=b.s[i]) return 0;
		return 1;
	}
	inline bool operator !=(num b) {
		return !(*this==b);
	}
	inline bool operator <(num b) {
		num a=*this;
		if (a.len<b.len) return 1;
		if (a.len>b.len) return 0;
		drp(i,a.len,1) {
			if (a.s[i]<b.s[i]) return 1;
			else if (a.s[i]>b.s[i]) return 0;
		}
		return 0;
	}
	inline bool operator <=(num b) {
		num a=*this;
		if (a<b||a==b) return 1;
		return 0;
	}
	inline bool operator >(num b) {
		num a=*this;
		if (a.len>b.len) return 1;
		if (a.len<b.len) return 0;
		drp(i,a.len,1) {
			if (a.s[i]>b.s[i]) return 1;
			else if (a.s[i]<b.s[i]) return 0;
		}
		return 0;
	}
	inline bool operator >=(num b) {
		num a=*this;
		if (a>b||a==b) return 1;
		return 0;
	}
	inline num operator +(num b) {
		num a=*this,c=(num) { {0},std:: max(a.len,b.len)};
		int v=0;
		rep(i,1,c.len) {
			c.s[i]=(a.s[i]+b.s[i]+v)%MOD;
			v=(a.s[i]+b.s[i]+v)/MOD;
		}
		if (v) {
			c.len+=1;
			c.s[c.len]=v;
		}
		return c;
	}
	inline num operator -(num b) {
		num a=*this,c=(num) { {0},std:: max(a.len,b.len)};
		rep(i,1,c.len) {
			c.s[i]=a.s[i]-b.s[i];
			if (c.s[i]<0) {
				c.s[i]+=MOD;
				a.s[i+1]-=1;
			}
		}
		while (!c.s[c.len]&&c.len>1) c.len-=1;
		return c;
	}
	inline num operator *(num b) {
		num a=*this,c=(num) { {0},a.len+b.len};
		rep(i,1,a.len) {
			rep(j,1,b.len) {
				c.s[i+j-1]+=a.s[i]*b.s[j];
			}
		}
		rep(i,1,a.len+b.len) {
			c.s[i+1]+=c.s[i]/MOD;
			c.s[i]%=MOD;
		}
		while (!c.s[c.len]&&c.len>1) c.len-=1;
		return c;
	}
	inline num operator /(int b) {
		num a=*this,c=(num) { {0},a.len};
		int v=0;
		drp(i,len,1) {
			int t=v*MOD+a.s[i];
			c.s[i]=t/b;
			v=t%b;
		}
		while (!c.s[c.len]&&c.len>1) c.len-=1;
		return c;
	}
	inline int operator %(int b) {
		LL res=0;
		drp(i,len,1) {
			res=(res*MOD%b+s[i])%b;
		}
		return res;
	}
	inline void read() {
		fill(s,0); len=0;
		scanf("%s",st);
		int v=0,i;
		for (i=strlen(st)-1; i >= 3; i-=3) {
			rep(j,i-2,i) v=v*10+st[j]-'0';
			s[++len]=v;
			v=0;
		}
		rep(j,0,i) v=v*10+st[j]-'0';
		s[++len]=v;
	}
	inline void read1(int x) {
		fill(s,0); len=0;
		do {s[++len]=x%MOD;} while (x/=MOD);
	}
	inline void output() {
		num tmp=*this;
		int i=tmp.len;
		while (!tmp.s[i]&&i>1) i--;
		printf("%d",tmp.s[i]);
		drp(j,i-1,1) {
			int v=tmp.s[j];
			int f[5]= {0};
			rep(k,1,3) {
				f[k]=v%10;
				v/=10;
			}
			drp(k,3,1) printf("%d",f[k]);
		} puts("");
	}
} f[N][N],ans,S,K,ze,one,wjp;

num C(num n,int m) {
	num res=n;
	for (int i=m;i>1;--i) {
		n=n-one;
		res=res*n;
	}
	for (int i=m;i>1;--i) {
		res=res/i;
	}
	return res;
}

void solve(int n,int m) {
	f[1][0]=one;
	rep(size,1,(n-1)/2) {
		S=ze;
		rep(j,0,m-1) S=S+f[size][j];
		drp(i,n,size+1) {
			rep(j,1,m) {
				rep(k,1,j) if (k*size<i) {
					K.read1(k);
					f[i][j]=f[i][j]+f[i-k*size][j-k]*C(S+K-one,k);
				}
			}
		}
	}
}

int main(void) { one.read1(1); ze.len=1;
	int n; scanf("%d",&n);
	solve(n,4);
	rep(i,0,4) ans=ans+f[n][i];
	if (n%2==0) {
		S=ze;
		rep(j,0,3) S=S+f[n/2][j];
		ans=ans+((S+one)*S)/2;
	}
	ans.output();
	return 0;
}