<错排>codevs 1697 ⑨要写信

去题目的传送门
题面很长，说白了，就是要你求错排n个数的方案数
表示看完错排的通项公式的推导过程之后，整个人都不好了。但是这个题不用通项公式，只知道递推式就好了。
递推式：f[n]=(n-1)*(f[n-1]+f[n-2])
来一波递推式的推导过程：

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

——–摘自百度百科
一开始看第二步的第二种情况的时候，死活没看明白。这里再解释一下。
为什么是剩下的n-1个元素的错排的方案数呢？
看一个例子：
1 2 3 4 5
现在固定把1放到5的位置，如果5不在1的位置上，我们把4抠掉，相当等于剩下的5个数：
5 2 3 1进行错排（也就是这样所有的错排的情况中，5肯定不会不放在第一个位置）
对了，这个题主要的代码都用来打高精度，n最大是100，long long肯定会爆

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;

int n;
struct hh
{
    int a[1100];
    int cnt;
}f[110];

hh add(hh x,hh y)
{
    hh ans;
    memset(ans.a,0,sizeof(ans.a));
    int l=max(x.cnt,y.cnt);
    for(int i=1;i<=l;++i)
    {
        ans.a[i]+=x.a[i]+y.a[i];
        if(i==l&&ans.a[i]>=10) ans.cnt=++l;
        else ans.cnt=l;
        ans.a[i+1]+=ans.a[i]/10;
        ans.a[i]=ans.a[i]%10;
    }
    return ans;
}
hh muti(int x,hh y)
{
    hh ans;
    memset(ans.a,0,sizeof(ans.a));
    int l=y.cnt;
    for(int i=1;i<=l;++i) ans.a[i]=y.a[i]*x;
    for(int i=1;i<=l;++i)
    {
        ans.a[i+1]+=ans.a[i]/10;
        ans.a[i]=ans.a[i]%10;
    }
    if(ans.a[l+1]) ans.cnt=++l;
    else ans.cnt=l;
    return ans;
}
int main()
{
    scanf("%d",&n);
    f[1].a[1]=0,f[2].a[1]=1;
    f[1].cnt=f[2].cnt=1;
    if(n<3) printf("%d",f[n].a[1]);
    else 
    {
        for(int i=3;i<=n;++i)
          f[i]=muti(i-1,add(f[i-1],f[i-2]));
        for(int i=f[n].cnt;i>=1;--i)
          printf("%d",f[n].a[i]);
    }

    return 0;
}