【基础编程】HDOJ2409 组合+错排

【1】题目分析
从N个新郎中选出M个，这是组合，同时M个新郎全部找错了新娘，是错排
错排的公式D(n) = (n-1)*(D(n-1)+D(n-2))
【2】代码

#include<stdio.h>
#define MAX 21
int main(){
    int b,n,m,i,j;
    long long sum1;
    long long a[MAX];
    a[0] = 0;
    a[1] = 0;
    a[2] =1;
    for(i = 3;i < MAX;i++){
        
        a[i] = (i-1) * (a[i-1]+a[i-2]);
    }
    while(EOF!=scanf("%d",&b)){
        while(b--){
            sum1 = 1;
            scanf("%d %d",&n,&m);
            
            for(i=n-m+1;i<=n;i++){
            
                sum1 *= i;
            }
            for(j = m;j >=1;j--){
                
                sum1 /= j;
            } 
            printf("%I64d\n",sum1*a[m]);
        }
    }
    return 0;
} 

【3】关于错排公式推导过程
当n个编号元素放在n个编号位置，元素编号与位置编号各不对应的方法数用D(n)表示，那么D(n-1)就表示n-1个编号元素放在n-1个编号位置，各不对应的方法数，其它类推.
第一步，把第n个元素放在一个位置，比如位置k，一共有n-1种方法；
第二步，放编号为k的元素，这时有两种情况：⑴把它放到位置n，那么，对于剩下的n-1个元素，由于第k个元素放到了位置n，剩下n-2个元素就有D(n-2)种方法；⑵第k个元素不把它放到位置n，这时，对于这n-1个元素，有D(n-1)种方法；
综上得到
D(n) = (n-1) [D(n-2) + D(n-1)]
特殊地，D(1) = 0, D(2) = 1.
下面通过这个递推关系推导通项公式：
为方便起见，设D(k) = k! N(k), k = 1, 2, …, n,
则N(1) = 0, N(2) = 1/2.
n ≥ 3时，n! N(n) = (n-1) (n-1)! N(n-1) + (n-1)! N(n-2)
即 nN(n) = (n-1) N(n-1) + N(n-2)
于是有N(n) - N(n-1) = - [N(n-1) - N(n-2)] / n = (-1/n) [-1/(n-1)] [-1/(n-2)]…(-1/3) [N(2) - N(1)] = (-1)^n / n!.
因此
N(n-1) - N(n-2) = (-1)^(n-1) / (n-1)!,
N(2) - N(1) = (-1)^2 / 2!.
相加，可得
N(n) = (-1)^2/2! + … + (-1)^(n-1) / (n-1)! + (-1)^n/n!
因此
D(n) = n! [(-1)^2/2! + … + (-1)^(n-1)/(n-1)! + (-1)^n/n!].
此即错排公式。
【4】容斥原理
用容斥原理也可以推出错排公式：
正整数1, 2, 3, ……, n的全排列有 n! 种，其中第k位是k的排列有 (n-1)! 种；当k分别取1, 2, 3, ……, n时，共有n*(n-1)!种排列是至少放对了一个的，由于所求的是错排的种数，所以应当减去这些排列；但是此时把同时有两个数不错排的排列多排除了一次，应补上；在补上时，把同时有三个数不错排的排列多补上了一次，应排除；……；继续这一过程，得到错排的排列种数为
D(n) = n! - n!/1! + n!/2! - n!/3! + … + (-1）^n*n!/n! = ∑(k=2~n) (-1)^k * n! / k!,
即D(n) = n! [1/0! - 1/1! + 1/2! - 1/3! + 1/4! + … + (-1)^n/n!].
其中，∑表示连加符号，k=2~n是连加的范围；0! = 1，可以和1!相消。