几道数论题

第二类Stirling数

n个有类别的球放到m个相同的盒子中，要求不能有空盒，有多少种放法？

这道题就是第二类Stirling数的经典例子，表示为S(n,m)。

它有递推关系S(n,m)=m*S(n-1,m)+S(n-1,m-1) 其中n,m>=1

当盒子数大于球数的时候即m>n有S(n,m)=0

         

盒子和球的问题

这个问题在高中的排列组合中就经常碰到，下面就列个表分析一下它们的性质：

N个球	M个盒子	是否允许空盒	方案数
有区别	有区别	是	M^N
有区别	有区别	否	M!*S(N,M)
有区别	无区别	是	∑­1MS(N,i)
有区别	无区别	否	S(N,M)
无区别	有区别	是	C(N+M-1,N)
无区别	有区别	否	C(N-1,M-1)
无区别	无区别	是	Xn­的系数
无区别	无区别	否	Xn­的系数

1、这个很好理解。

2、这个可以通过盒子没有区别来理解，当盒子有区别时，盒子产生了M!种可能了。

3、这个跟4是一个意思，不过盒子可以为空，就当这个盒子不存在。

4、这个是Stirling数的性质。

5、这个可以用隔板的思想来理解，有N+M-1个位置，其中先用M-1个位置放置隔板，剩下的放置球，隔板隔开的就是盒子里面的球，关键是如何选择隔板的位置，相当于如何选择球的位置。

6、当盒子不能为空时，隔板的放法就有限制了，只能选择球之间的位置，而且每个这样的位置至多只能放置一块这样的隔板。

7、在球和盒子都没有区别的分析中，可能读者还深感疑惑，下面将深入分析：

这个问题可以很直接第地转化为N可以分解为M个数的和，有多少种分法。

设式子(1+x+x²+x³+x⁴…)表示一个数可以分解为多少个1的可选式子，x的指数部分/1表示分解为1的个数；

设式子(1+x²+x⁴+x⁶+x⁸…)表示一个数可以分解为多少个2的可选式子，x的指数部分/2表示分解为2的个数；

…

设式子(1+x^2r+x^4r+x^6r+x^8r…)表示一个数可以分解为多少个r的可选式子，x的指数部分/r表示分解为r的个数；

这样，(1+x+x²+x³+x⁴…)* (1+x²+x⁴+x⁶+x⁸…)*…* (1+x^2r+x^4r+x^6r+x^8r…)展开后xⁿ的系数(假设r>=n)就是n分解为若干数的和的种类数。r的作用就是控制其中分解的最大数字：

例如对6进行分解，如果要求最大的数不能超过2，则只有(1,1,1,1,1,1)，(2,1,1,1,1)，(2,2,1,1)，(2,2,2)这4种分法，对应于上式，则分别是x⁶*1，x⁴*x²，x²*x⁴，1*x⁶，对应于格子图(Ferrer图)，将格子其按红色的线旋转，得到下面的格子图，这就是分解不超过两个数的格子图。也就是说将N分解成不超过M=2个数的分法，等价于将N分解成的数不超过M的分法，等价于将无区别的N个球放入无区别的M个盒子中，盒子可以为空。

 

 

 

8、当盒子不能为空时候，取N个无区别的球放入无区别的M个盒子，盒子可以为空的放法，减去，取N个无区别的球放入无区别的M-1个盒子，盒子可以为空的放法，即得，取N个无区别的球放入无区别的M个盒子，盒子不可以为空。

关于自然数的分解，我这里写了个程序，有兴趣的读者可以编译运行一下，不过运行70以上的data，可能就要等上10s并且准备160m的空间存放输出文件了。

 #include <stdio.h> #define MAXN 100 int a[MAXN],r[MAXN]; void number_divide(int n) { int i,k; k=0;r[0]=n; do{ if(r[k]==0){ printf("%d=%d",a[0],a[1]); for(i=2;i<=k;i++) printf("+%d",a[i]); printf("/n"); while(k>0 && a[k]==1) k--; if(k>0){a[k]--;r[k]++;} } else { a[k+1]=a[k]<r[k]?a[k]:r[k]; r[k+1]=r[k]-a[k+1]; k++; } }while(k>0); } int main() { freopen("out.txt","w",stdout); int data=70; a[0]=data; number_divide(data); return 0; }

  

贝尔数 (bell数--集合划分数)

给定一个n表示{1,2,...,n}这样的集合，比如给定3表示{1,2,3}它可以构造出

{{1}{2}{3}}

{{1,2}{3}}

{{1,3}{2}}

{{2,3}{1}}

{{1,2,3}}

这样5个集合。题目给出n，算出能构造多少个这样的集合。

 

贝尔数就是集合的划分数。它具有以下性质：

1

1 2

2 3 5

5 7 10 15

...

它的第一行第一项是1

第n行第一项是第n-1行的最后一项，其他则为左边的数加上面的数，每一行的最后那个数，也就是下一行的第一个数，就是这道题的答案。这道题也可以用上面讲的Stirling数来解决，盒子分析的第三种情况，也是这道题的答案，读者自己分析一下吧。

 

卡特兰(Catalan)数

令h(0)=1且h(n)=C(2n,n)/(n+1) (n=1,2,3,…)

这种数可以解决括号化问题，出栈次序问题，多边形划分为三角形问题，二叉树种类问题。

1、括号化问题来自矩阵连乘：P=a1×a2×a3×……×an，依据乘法结合律，不改变其顺序，只用括号表示成对的乘积，试问有几种括号化的方案？(h(n-1)种)百度词条里面我自己修改编辑的。

2、出栈次序问题来自栈本身固有的先进先出的性质，这种性质决定了栈不能产生全部的序列。这个问题可以衍生为排队入场问题，即电影票5元一张，售票员在没有观众进场前没有备有零钱，但是排队的2n个人当中有n个是持5元钞票，另外n个是持10元钞票，要令到售票员能顺利让所有人看到电影，有多少种排队方法。h(n)种

3、多边形划分是指将一个n+2条边的凸多边形分成三角形区域的方法数问题。衍生出来的问题有方块走法问题，即n个方块*n个方块 的地图中，从最左上角走到最右下角(下图中1走到2)，路程最短的走法有多少种。h(n)种

1
			2

4、给定顶点组成二叉树的问题。

给定n个顶点，能构成多少种不同的二叉树。h(n)种