错排公式详释，HDU1465，HDU4535，HDU2048，HDU2049

定义：

错排问题，是组合数学中的问题之一。考虑一个有n个元素的排列，若一个排列中所有的元素都不在自己原来的位置上，那么这样的排列就称为原排列的一个错排。n个元素的错排数记为D(n)。

经典的递推：

递归关系式为：

D(n)=(n-1)(D(n-1)+D(n-2))
D(1)=0,D(2)=1
可以得到:
错排公式为 f(n) = n![1-1/1!+1/2!-1/3!+……+(-1)^n*1/n!] 
其中,n!=1*2*3*.....*n,
特别地,有0!=0,1!=1.

解释： 

n 个不同元素的一个错排可由下述两个步骤完成： 
第一步，“错排” 1 号元素（将 1 号元素排在第 2 至第 n 个位置之一），有 n - 1 种方法。 
第二步，“错排”其余 n - 1 个元素，按如下顺序进行。视第一步的结果，若1号元素落在第 k 个位置，第二步就先把 k 号元素“错排”好， k 号元素的不同排法将导致两类不同的情况发生：

1、 k 号元素排在第1个位置，留下的 n - 2 个元素在与它们的编号集相等的位置集上“错排”，有 f(n -2) 种方法；
2、 k 号元素不排第 1 个位置，这时可将第 1 个位置“看成”第 k 个位置(也就是说本来准备放到k位置为元素，可以放到1位置中),于是形成（包括 k 号元素在内的） n - 1 个元素的“错排”，有 f(n - 1) 种方法。据加法原理，完成第二步共有 f(n - 2)+f(n - 1) 种方法。 
根据乘法原理， n 个不同元素的错排种数 
f(n) = (n-1)[f(n-2)+f(n-1)] (n>2) 。 

HDU 含错排公式的题目：

HDU1465 不容易系列之一：

错排公式模板题。

AC代码：

/**
  * 行有余力，则来刷题！
  * 博客链接：http://blog.csdn.net/hurmishine
  *
*/
#include <iostream>
#include <cstdio>
using namespace std;

int main()
{
    //freopen("data.txt","r",stdin);
    int n;
    long long a[25];
    a[1]=0,a[2]=1;
    for(int i=3;i<=20;i++)
        a[i]=(i-1)*(a[i-2]+a[i-1]);
    long long b[25];
    while(cin>>n)
    {
        cout<<a[n]<<endl;
    }
    return 0;
}

HDU 4535 吉哥系列故事

错排+取模

AC代码：

/**
  * 行有余力，则来刷题！
  * 博客链接：http://blog.csdn.net/hurmishine
  *
*/
#include <iostream>
#include <cstdio>
using namespace std;
const int mod=1e9+7;
const int maxn=105;
int main()
{
    //freopen("data.txt","r",stdin);
    int n;
    long long a[maxn];
    a[1]=0,a[2]=1;
    for(int i=3;i<=maxn;i++)
        a[i]=((i-1)*(a[i-2]+a[i-1]))%mod;
    int T;
    cin>>T;
    while(T--)
    {
        cin>>n;
        cout<<a[n]<<endl;
    }
    return 0;
}

HDU 2048 神、上帝以及老天爷 

错排+全排。

#include <cstdio>
#include <iostream>
using namespace std;
typedef long long LL;
LL fac[25];
LL d[25];
void getFac()
{
    fac[0]=fac[1]=1;
    for(int i=2; i<=20; i++)
    {
        fac[i]=fac[i-1]*i;
    }
}
void getD()
{
    d[1]=0,d[2]=1;
    for(int i=3;i<=20;i++)
    {
        d[i]=(i-1)*(d[i-2]+d[i-1]);
    }
}
int main()
{
    int T,n;
    getFac();
    getD();
    cin>>T;
    while(T--)
    {
        cin>>n;
        printf("%.2lf%%\n",d[n]*100.0/fac[n]);
    }
    return 0;
}

HDU 2049 不容易系列之(4)——考新郎

排列组合+错排

#include <cstdio>
#include <iostream>
using namespace std;
typedef long long LL;
LL d[25];
void getD()
{
    d[1]=0,d[2]=1;
    for(int i=3; i<=20; i++)
    {
        d[i]=(i-1)*(d[i-1]+d[i-2]);
    }
}
LL C(int n,int m)
{
    if(m>n/2) m=n-m;
    LL ans=1;
    for(int i=n; i>n-m; i--)
        ans*=i;
    for(int i=m; i>0; i--)
        ans/=i;
    return ans;
}
int main()
{
    int T;
    int n,m;
    cin>>T;
    getD();
    while(T--)
    {
        cin>>n>>m;
        //cout<<C(n,m)<<endl;
        cout<<d[m]*C(n,m)<<endl;
    }
    return 0;
}