loj2061. 「HAOI2016」放棋子
给你一个 的矩阵,每行有一个障碍,数据保证任意两个障碍不在同一行,任意两个障碍不在同一列,要求你在这个矩阵上放 枚棋子(障碍的位置不能放棋子),要求你放 个棋子也满足每行只有一枚棋子,每列只有一枚棋子的限制,求有多少种方案。
Sol
考虑容斥。
前i枚棋子与>=i个障碍冲突的方案数为C(n,i)*(n-i)!=n!*i!
唯一有点麻烦的是要高精。
没错,这题跟障碍在哪里没啥关系

#include<cstdio> #include<iostream> #include<cstdlib> #include<cstring> #include<algorithm> #include<cmath> #define ll __int128 #define maxn 505 using namespace std; int n,s[205],top; struct big{ int l,s[maxn]; }tmp,A,B,ans; big operator +(big a,big b){ int l=max(a.l,b.l); for(int i=1;i<=l;i++){ a.s[i]+=b.s[i]; if(a.s[i]>=10)a.s[i]-=10,a.s[i+1]++,l=max(l,i+1); } a.l=l;return a; } big operator -(big a,big b){ int l=a.l; for(int i=1;i<=l;i++){ a.s[i]-=b.s[i]; if(a.s[i]<0)a.s[i]+=10,a.s[i+1]--; } if(a.s[a.l]==0)a.l--;return a; } big mul(big a,int v){ int l=a.l; for(int i=1;i<=l;i++)a.s[i]=a.s[i]*v; for(int i=1;i<=l;i++){ a.s[i+1]+=a.s[i]/10,a.s[i]%=10; if(a.s[i+1]>0)l=max(l,i+1); } a.l=l; return a; } int main(){ cin>>n; tmp.l=tmp.s[1]=1; for(int i=n;i>=0;i--){ if(i&1)B=B+tmp; else A=A+tmp; tmp=mul(tmp,i); } ans=A-B; for(int i=ans.l;i>=1;i--)cout<<ans.s[i];puts(""); return 0; }
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】博客园社区专享云产品让利特惠,阿里云新客6.5折上折
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· DeepSeek 解答了困扰我五年的技术问题
· 为什么说在企业级应用开发中,后端往往是效率杀手?
· 用 C# 插值字符串处理器写一个 sscanf
· Java 中堆内存和栈内存上的数据分布和特点
· 开发中对象命名的一点思考
· DeepSeek 解答了困扰我五年的技术问题。时代确实变了!
· PPT革命!DeepSeek+Kimi=N小时工作5分钟完成?
· What?废柴, 还在本地部署DeepSeek吗?Are you kidding?
· 赶AI大潮:在VSCode中使用DeepSeek及近百种模型的极简方法
· DeepSeek企业级部署实战指南:从服务器选型到Dify私有化落地
2019-02-21 lct模板
2019-02-21 [ZJOI2013]K大数查询
2019-02-21 叶落归根(hometown)