[BZOJ 4563]放棋子

[BZOJ 4563]放棋子

题目

给你一个N*N的矩阵，每行有一个障碍，数据保证任意两个障碍不在同一行，任意两个障碍不在同一列，要求你在这个矩阵上放N枚棋子（障碍的位置不能放棋子），要求你放N个棋子也满足每行只有一枚棋子，每列只有一枚棋子的限制，求有多少种方案。

INPUT

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

OUTPUT

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

SAMPLE

INPUT

2
0 1
1 0

OUTPUT

1

解题报告

一道读入都不用打完的题= =

只用读入$n$就好了QAQ

全错位排列往上一打，加个高精就好了QAQ

 #include<iostream>
 #include<cstring>
 #include<cstdio>
 using namespace std;
 struct bign{
     int a[],len;
     bign():len(){memset(a,,sizeof(a));}
     bign(const char *s){*this=s;}
     bign(const int &x){*this=x;}
     inline void operator=(const char *s){
         len=strlen(s);
         for(int i=;i<len;++i)
             a[i]=s[len-i-]-'';
     }
     inline void operator=(const int &x){
         char s[];
         memset(s,,sizeof(s));
         sprintf(s,"%d",x);
         *this=s;
     }
     inline void clear(){while(len>&&!a[len-])--len;}
     inline bign operator+(const bign &x){
         bign ret;
         ret.len=max(len,x.len)+;
         int jw();
         for(int i=;i<len;++i){
             ret.a[i]=a[i]+x.a[i]+jw;
             jw=ret.a[i]/;
             ret.a[i]%=;
         }
         if(jw)
             ret.a[len++]=jw;
         ret.clear();
         return ret;
     }
     inline void operator+=(const bign &x){*this=*this+x;}
     inline bign operator-(const bign &x){
         bign ret;
         ret=*this;
         for(int i=;i<len;++i){
             if(ret.a[i]<x.a[i])
                 ret.a[i]+=,--ret.a[i+];
             ret.a[i]=ret.a[i]-x.a[i];
         }
         ret.clear();
         return ret;
     }
     inline void operator-=(const bign &x){*this=*this-x;}
     inline bign operator-(const int &x){
         bign tmp(x);
         return *this-tmp;
     }
     inline void operator-=(const int &x){*this=*this-x;}
     inline bign operator*(const bign &x){
         bign ret;
         ret.len=len+x.len+;
         int jw;
         for(int i=;i<len;++i){
             jw=;
             for(int j=;j<x.len;++j){
                 ret.a[i+j]+=a[i]*x.a[j]+jw;
                 jw=ret.a[i+j]/;
                 ret.a[i+j]%=;
             }
             ret.a[i+x.len]+=jw;
         }
         ret.clear();
         return ret;
     }
     inline void operator*=(const bign &x){*this=*this*x;}
     inline bign operator*(const int &x){
         bign tmp(x);
         return *this*tmp;
     }
     inline void operator*=(const int &x){*this=*this*x;}
     inline void print(){
         for(int i=len-;i>=;--i)
             printf("%d",a[i]);
     }
 }f[];
 int n;
 int main(){
     scanf("%d",&n);
     f[]=;
     f[]=;
     for(int i=;i<=n;++i)
         f[i]=(f[i-]+f[i-])*(i-);
     f[n].print();
 }