标签：障碍物 include num HAOI2016 棋子 lll 错排

[HAOI2016]放棋子

有一\(N\times N\)的棋盘，有N个障碍物放入其中，已经给出其位置，并且保证每一行每一列有且仅有一个障碍物，现在再给你N个棋子，放入其中，保证每一行每一列有且仅有一颗棋子，询问其方案数。

解

根据错排自由定理，我们完全可以把障碍物按照行递增的方式的排序起来，即把障碍物全部摆到对角线上特殊化问题，于是对于每一列摆棋子即变成不能摆对应列号的行，转换模型即一个1~N全排列中，保证每个数都不摆到其对应编号的位置上，此时就变成了一个裸的错排模型，设\(f[n]\)表示长n的错排问题的方案数，与\(f[n-1]\)关系在于多了一个数n放入其中，它显然不能放到第n个位置，于是放在前面有n-1种选择，而被放的位置的数可以选择放到第n个位置，这样就变成\(f[n-2]\)的错排问题，而如果不放到最后一个位置，不难得知等价于\(f[n-1]\)的错排问题，于是我们有

\[f[n]=(n-1)\times(f[n-1]+f[n-2])\]

边界：\(f[1]=0,f[2]=1\)

答案：\(f[n]\)

打个高精即可。

参考代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#define il inline
#define ri register
using namespace std;
struct lll{
    int num[1000];
    lll(){memset(num,0,sizeof(num)),num[0]|=true;}
    il void operator=(string s){
        num[0]=s.size();
        for(int i(1);i<=num[0];++i)
            num[i]=s[num[0]-i]-48;
    }
    il void print(){
        for(ri int i(num[0]);i;--i)
            putchar(num[i]+48);
    }
    il lll operator+(lll x){
        lll y;ri int i;
        for(i=1;i<=num[0]||i<=x.num[0];++i){
            y.num[i]+=num[i]+x.num[i];
            if(y.num[i]>9)y.num[i]-=10,++y.num[i+1];
        }while(!y.num[i]&&i)--i;y.num[0]=i;return y;
    }
    il lll operator*(lll x){
        lll y;ri int i,j,k;
        for(i=1;i<=num[0];++i){
            k=0;
            for(j=1;j<=x.num[0];++j)
                y.num[i+j-1]+=k+num[i]*x.num[j],
                    k=y.num[i+j-1]/10,y.num[i+j-1]%=10;
            y.num[i+x.num[0]]+=k;
        }y.num[0]=num[0]+x.num[0];
        while(y.num[0]&&!y.num[y.num[0]])--y.num[0];
        return y;
    }
}dp[201],lsy;
int main(){
    int n,i;scanf("%d",&n);
    dp[2]="1",lsy="2";
    for(i=3;i<=n;++i,lsy=lsy+dp[2])
        dp[i]=(dp[i-1]+dp[i-2])*lsy;
    dp[n].print();
    return 0;
}

标签：障碍物,include,num,HAOI2016,棋子,lll,错排
来源： https://www.cnblogs.com/a1b3c7d9/p/10801567.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。