BZOJ 4563 错排+高精度

思路：

把障碍移到对角线

就发现 这是个错位排列问题

用错排公式即可解

s[i]=(s[i-1]+s[i-2])*i

//By SiriusRen
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
struct Bignum{
    int a[],num;
    void init(){memset(a,,sizeof(a));num=;}
    friend Bignum operator+(Bignum a,Bignum b){
        Bignum c;c.init();
        for(int i=;i<=max(a.num,b.num);i++)
            c.a[i]=a.a[i]+b.a[i];
        for(int i=;i<=a.num+b.num;i++){
            c.a[i+]+=c.a[i]/,c.a[i]%=;
            if(c.a[i])c.num=i;
        }return c;
    }
    friend Bignum operator*(Bignum a,Bignum b){
        Bignum c;c.init();
        for(int i=;i<=a.num;i++)
            for(int j=;j<=b.num;j++)
                c.a[i+j-]+=a.a[i]*b.a[j];
        for(int i=;i<=a.num+b.num;i++){
            c.a[i+]+=c.a[i]/,c.a[i]%=;
            if(c.a[i])c.num=i;
        }return c;
    }
}f[],g;
int main(){
    int n;
    f[].a[]=f[].num=g.num=,g.a[]=;
    for(int i=;i<=;i++)g=g+f[],f[i]=(f[i-]+f[i-])*g;
    scanf("%d",&n);
    for(int j=f[n].num;j;j--)printf("%d",f[n].a[j]);
}