HDU 2068 排列组合+错排

贴一下原题地址:http://acm.hdu.edu.cn/showproblem.php?pid=2068

今年暑假杭电ACM集训队第一次组成女生队，其中有一队叫RPG，但做为集训队成员之一的野骆驼竟然不知道RPG三个人具体是谁谁.RPG给他机会让他猜猜，第一次猜：R t是公主，P是草儿，G是月野兔;第二次猜：R t是草儿，P是月野兔，G是公主;第三次猜：R t是草儿，P是公主，G是月野兔; ……可怜的野骆驼第六次终于把RPG分清楚了由于RPG的带动，做ACM的女生越来越多，我们的野骆驼想都知道她们，可现在有ñ多人，他要猜的次数可就多了，为了不为难野骆驼，女生们只要求他答对一半或以上就算过关，请问有多少组答案能使他顺利过关。

关于错排公式,高中就有提及.
这里重新放一下推导过程

当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.

由于只要答对一半以上,用C(N,i)(0<=i<=n/2) 奇数额外+1
就能AC

#include <stdio.h>
typedef long long int64;
int64 Com ( int n, int m )  // 取 Cn m 的 组合数
{
      if ( m == 0 )
      {
           return 1; 
      }
      int64 up = 1,down = 1;
      for ( int i = 0; i < m; ++ i )
      {
            up *= n - i;
            down *= i + 1; 
      } 
      return up / down;
}
int64 F[14] = { 0, 0, 1, 2 };
void setNum ()
{
     for ( int i = 4; i < 14; ++ i )
     {
           F[i] = ( i - 1 ) * ( F[ i - 1 ] + F[ i - 2 ] ); 
     } 
}
int main ()
{
    int N;
    int64 sum;
    setNum ();
    while ( scanf("%d", &N), N )
    {
          sum = 1;
          for ( int i = N / 2; i >= 0; -- i )
          {
               sum += Com ( N , i ) * F[i];    // 从N中取i 个人 的错排数
          }
          printf("%I64d\n", sum );
    }
}