[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

---

由于移动障碍点并不会答案造成影响，所以我们把障碍点都移动到\((1,1),(2,2),...,(n,n)\)这些格子上，那么对于每一行每一列而言，不能放置的行与列就是第\(i\)行，第\(i\)列，这样就变成了一个经典的错拍问题

记得高精度

/*program from Wolfycz*/
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define inf 0x7f7f7f7f
using namespace std;
typedef long long ll;
typedef unsigned int ui;
typedef unsigned long long ull;
inline char gc(){
    static char buf[1000000],*p1=buf,*p2=buf;
    return p1==p2&&(p2=(p1=buf)+fread(buf,1,1000000,stdin),p1==p2)?EOF:*p1++;
}
inline int frd(){
    int x=0,f=1; char ch=gc();
    for (;ch<'0'||ch>'9';ch=gc())   if (ch=='-')    f=-1;
    for (;ch>='0'&&ch<='9';ch=gc()) x=(x<<3)+(x<<1)+ch-'0';
    return x*f;
}
inline int read(){
    int x=0,f=1; char ch=getchar();
    for (;ch<'0'||ch>'9';ch=getchar())  if (ch=='-')    f=-1;
    for (;ch>='0'&&ch<='9';ch=getchar())    x=(x<<3)+(x<<1)+ch-'0';
    return x*f;
}
inline void print(int x){
    if (x<0)    putchar('-'),x=-x;
    if (x>9)    print(x/10);
    putchar(x%10+'0');
}
const int digit=4;
const int base=1e4;
const int maxn=2e2;
struct Bignum{
    int v[maxn],len;
    Bignum(){memset(v,0,sizeof(v)),len=1;}
    void write(){
        printf("%d",v[len-1]);
        for (int i=len-2;~i;i--)    printf("%0*d",digit,v[i]);
        putchar('\n');
    }
}f[maxn+10];
Bignum operator +(const Bignum &x,const Bignum &y){
    Bignum z;
    z.len=max(x.len,y.len);
    for (int i=0;i<=z.len;i++)  z.v[i]+=x.v[i]+y.v[i],z.v[i+1]+=z.v[i]/base,z.v[i]%=base;
    while (z.v[z.len])  z.v[z.len+1]+=z.v[z.len]/base,z.v[z.len]%=base,z.len++;
    return z;
}
Bignum operator *(const Bignum &x,int y){
    Bignum z; z.len=x.len;
    for (int i=0;i<=z.len;i++)  z.v[i]+=x.v[i]*y,z.v[i+1]+=z.v[i]/base,z.v[i]%=base;
    while (z.v[z.len])  z.v[z.len+1]+=z.v[z.len]/base,z.v[z.len]%=base,z.len++;
    return z;
}
int main(){
    int n=read();
    for (int i=1;i<=n;i++)  for (int j=1;j<=n;j++)  read();
    f[0].v[0]=1,f[1].v[0]=0;
    for (int i=2;i<=n;i++)  f[i]=(f[i-1]+f[i-2])*(i-1);
    f[n].write();
    return 0;
}

转载于:https://www.cnblogs.com/Wolfycz/p/10060649.html