一,错排公式
Dn=[n!/e+0.5]只是近似值!4次WA的代价啊!
应该写成这样
当n个编号元素放在n个编号位置,元素编号与位置编号各不对应的方法数用M(n)表示,那么M(n-1)就表示n-1个编号元素放在n-1个编号位置,各不对应的方法数,其它类推.
第一步,把第n个元素放在一个位置,比如位置k,一共有n-1种方法;
第二步,放编号为k的元素,这时有两种情况.1,把它放到位置n,那么,对于剩下的n-2个元素,就有M(n-2)种方法;2,不把它放到位置n,这时,对于这n-2个元素,有M(n-1)种方法;
综上得到
M(n)=(n-1)[M(n-2)+M(n-1)]
特殊地,M(1)=0,M(2)=1
- for (i=3;i<=20;i++)
- s[i] = (i-1)*(s[i-1] + s[i-2]);
n的错排数即为s[n]
例题:HDU 2094
二、组合数
以前是先算出上下两项然后相除,发现会局部越界。
比较科学的写法是乘完一个后马上除掉一个~
typedef long long ll;
using namespace std;
ll c(ll a, ll b)
{
ll ans=1, i;
for(i=1;i<=b;i++)
{
ans*=a-i+1;
ans /=i;
}
return ans;
}