错排问题,又称更列问题,是组合数学中的问题之一。对于它的研究最早可以追溯到十八世纪,当时他被数学家尼古拉·伯努利和欧拉研究,因此在历史上也被称为伯努利--欧拉的错装信封问题。这个问题有许多具体的版本,比如在写信时讲n封信装到n个不同的信封里,有多少种全部装错信封的情况?再比如n个人各写一张贺卡相互赠送,有多少种赠送方法?这些经典的题目都是典型的错排问题。
相信看过上面对于错排问题的简单的介绍,大家也都对它有了一些初步的了解,归结起来,就是考虑一个有n个元素的排列,若一个排列中所有的元素都不在自己原来的位置上,那么这样的排列就称为原排列的一个错排,n个元素的错排数记为D(n)。那么对于这样的排列D(n)有多少种呢?我们一步一步进行分析:
首先,对于D(n),有1~n这样n个元素错排,所以对于第一个元素①,它现在可能的位置有(n-1)个,倘若它在第k个元素的位置上,对于第k个元素而言,它所在的位置就有两种可能—第一种,它处在非第一个元素①位置上,所以对于接下来的排列就相当于是n-1个元素的错排,即D(n-1);第二种,它处在第一个元素①的位置上,所以在排列D(n)中有两个元素找到了位置,那么接下来的队列就相当于是n-2个元素的错排。因此,对于D(n)都有D(n)=(n-1)*(D(n-1)+D(n-2))【特殊的,D(1)=0,D(2)=1】。
得到这个递推式之后,我们进一步进行推导:为了运算的方便,我们设D(n)=n!*N(n),则有:
n!*N(n)=(n-1)*(n-2)!*N(n-2)+(n-1)*(n-1)!*N(n-1) ,两边同时除以(n-1)!,可得:n*N(n)=N(n-2)+(n-1)*N(n-1) ,移项:
N(n)-N(n-1)=(N(n-2)-N(n-1))/n = -(1/n)(N(n-1)-N(n-2)) ,所以N(n-1)-N(n-2)= -(1/(n-1))(N(n-2)-N(n-3)),··· ··· ,
直到N2-N(1)=1/2 ,将每一个式子相加易得:N(n)-N(1)=(1/2!-1/3!+1/4!- ··· ··· +((-1)^(n-1))/(n-1)!+((-1)^n)/n! )
由于N(1)=0,所以N(n)=(1/2!-1/3!+1/4!- ··· ··· +((-1)^(n-1))/(n-1)!+((-1)^n)/n! ) ,于是可以得到:
错排公式为D(n)=n!*(1/2!-1/3!+1/4!- 1/5!+ ··· ··· +((-1)^(n-1))/(n-1)!+((-1)^n)/n! )。
这样,我们就通过简单的推导得到了两个关于错排问题的公式!
现在,我们来具体问题具体分析,了解错排公式如何转化为代码来解决考试中实际遇到的问题,我们这里以HDU Online Judge上的一道题考新郎为例,题目是这样的:
阅读题目后我们不难发现,这道题的本质就是求解排列组合C(n,m)与错排m个元素D(m)的乘积,因此这道题的代码也十分简单,以下提供两种AC程序:
#方法1:递推公式--D(n)=(n-1)*(D(n-1)+D(n-2)) [D(1)=0,D(2)=1]
#include<iostream> #include<cstring> #include<cstdio> using namespace std; long long _cuopai[50]; long long jiecheng[22]={1,1,2,6,24,120,720,5040,40320, 362880,3628800,39916800,479001600, 6227020800,87178291200,1307674368000, 20922789888000,355687428096000,6402373705728000, 121645100408832000,2432902008176640000}; long long cuopai(int x) { if(_cuopai[x]) return _cuopai[x]; if(x==1) return 0; if(x==2) return 1; return _cuopai[x]=(x-1)*(cuopai(x-1)+cuopai(x-2)); } long long c(int y,int z) { return jiecheng[y]/(jiecheng[z]*jiecheng[y-z]); } int main() { memset(_cuopai,0,sizeof(_cuopai)); int a; cin>>a; for(int i=1;i<=a;++i) { int m,n; cin>>m>>n; cout<<c(m,n)*cuopai(n)<<endl; } return 0; }#方法2:通项公式--D(n)=n!*(1/2!-1/3!+1/4!- 1/5!+ ··· ··· +((-1)^(n-1))/(n-1)!+((-1)^n)/n! )
这里可以根据题目做一下变形:
F(n,m)=C(n,m)*D(m)=n!*(1/2!-1/3!+1/4!- 1/5!+ ··· ··· +((-1)^(n-1))/(n-1)!+((-1)^n)/n! )/(n-m)!
#include<iostream> #include<cstdio> using namespace std; int m,n; long long jiecheng[22]={1,1,2,6,24,120,720,5040,40320, 362880,3628800,39916800,479001600, 6227020800,87178291200,1307674368000, 20922789888000,355687428096000,6402373705728000, 121645100408832000,2432902008176640000}; long long cuopai_() { long long sum=0,a=jiecheng[n],b=jiecheng[n-m]; for(int i=2;i<=m;++i) { a/=i; if(i%2==0) sum+=a; else sum-=a; //cout<<"sum["<<i<<"]"<<sum<<endl; } return sum/b; } int main() { int x; cin>>x; for(int i=1;i<=x;++i) { cin>>n>>m; cout<<cuopai_()<<endl; } return 0; }这样就完成了关于错排公式的推论和简单应用~!