BZOJ 4563 [Haoi2016]放棋子 错排公式+高精度

Description

给你一个N*N的矩阵，每行有一个障碍，数据保证任意两个障碍不在同一行，任意两个障碍不在同一列，要求你在

这个矩阵上放N枚棋子（障碍的位置不能放棋子），要求你放N个棋子也满足每行只有一枚棋子，每列只有一枚棋子

的限制，求有多少种方案。

Input

第一行一个N，接下来一个N*N的矩阵。N<=200,0表示没有障碍，1表示有障碍，输入格式参考样例

Output

一个整数，即合法的方案数。

Sample Input

2
0 1
1 0

Sample Output

1

HINT

传送门

题目一开始没看清楚……= =

数据是保证每行都有障碍，而且每列每行都只有一个障碍的= =

那么可以看见障碍对答案是毫无影响的，

通过交换一些行始终可以把障碍移动到主对角线的位置，

接着放位置的条件变成了：第x行不能放在x位置，而且所有放的位置是一个排列。

……这不就是错排问题么？

错排公式F(n)=(n-1)*(F(n-1)+F(n-2))

网上有很多说明方法，不赘述了……

然后加个高精度即可。。

因为高精度位数有6000-位，所以不要开数组递推了，用3个变量迭代就好了。

高精度终于会打重载啦哈哈哈= v =

#include<bits/stdc++.h>
using namespace std;
const int
    N=2005,
    maxn=5750;
struct BigNum{
    int len,a[maxn];
    BigNum(){len=0;memset(a,0,sizeof(a));}
    void read(){
        char t[maxn];
        scanf("%s",t+1);
        len=strlen(t+1);
        for (int i=1;i<=len;i++) a[i]=t[len-i+1]-'0';
    }
    void print(){
        for (int i=len;i;i--) printf("%d",a[i]);
        putchar('\n');
    }
    BigNum operator =(int x){
        char t[maxn];
        sprintf(t+1,"%d",x);
        len=strlen(t+1);
        for (int i=1;i<=len;i++) a[i]=t[len-i+1]-'0';
        return (*this);
    }
    BigNum operator =(BigNum x){
        len=x.len;
        for (int i=1;i<=len;i++) a[i]=x.a[i];
        return (*this);
    }
    BigNum operator +(BigNum x){
        BigNum y;
        y.len=max(len,x.len);
        for (int i=1;i<=y.len;i++){
            y.a[i]+=x.a[i]+a[i];
            y.a[i+1]=y.a[i]/10;
            y.a[i]%=10;
        }
        if (y.a[y.len+1]) y.len++;
        return y;
    }
    BigNum operator *(BigNum x){
        BigNum y;
        y.len=x.len+len-1;
        for (int i=1;i<=len;i++)
            for (int j=1;j<=x.len;j++){
                y.a[i+j-1]+=x.a[j]*a[i];
                y.a[i+j]+=y.a[i+j-1]/10;
                y.a[i+j-1]%=10;
            }
        if (y.a[y.len+1]) y.len++;
        return y;
    }
};
int main(){
    int n;scanf("%d",&n);
    if (n==1){puts("0");return 0;}
    BigNum x,y,z,t;
    x=0,y=1,z=1;
    for (int i=3;i<=n;i++)
        t=i-1,z=(x+y)*t,x=y,y=z;
    z.print();
    return 0;
}