洛谷P4831 Scarlet loves WenHuaKe
这道题告诉我们推式子的时候头要够铁。
题意
问一个\(n\times m\)的棋盘,摆上\(n\times 2\)个中国象棋的炮使其两两不能攻击的方案数,对\(998244353\)取模。
\((n\leq m\leq 2000)或(n\leq m\leq 100000且m-n\leq 10)\)。
题解
怎么两个数据范围搞搞。
显然合法方案等价于每行每列炮的数量不超过\(2\),那么每一行就必定放\(2\)个炮了。
我们记\(f(n,m)\)为答案,考虑如何归约到规模更小的问题。
那么我们枚举最后一行炮的个数,分三类情况:
\(1\)、个数为\(0\),归约到\(f(n,m-1)\)。
\(2\)、个数为\(2\)(个数为\(1\)比较麻烦后面再说),那么先枚举放在这一列的是哪两行(\(\times \frac{n(n-1)}{2}\)),接着分类讨论这两行的另一个是否相同:
如果相同,则枚举这是哪一个\((\times (m-1))\),归约到\(f(n-2,m-2)\)。
如果不同,则这两行可以合并(同一行的唯一要求就是两个列不同),只要根据有序性\(\times 2\)即可,于是归约到\(f(n-1,m-1)\)。
\(3\)、个数为\(1\),那么先枚举占了最后一列的是哪一行\((\times n)\),再枚举这一行的另一个在哪一列\((\times (m-1))\),问题就转化为\(n-1\)行\(m-1\)列,其中有一列炮的个数\(\leq 1\)的方案数。
那么考虑容斥,用总方案数减去这一列放了两个的方案数。前者就是\(f(n-1,m-1)\),对于后者,进行与情况\(2\)相似的讨论,也可以进行计算。
可以发现\(n>m\)时\(f(n,m)=0\),于是复杂度就是\(O((m-n)n)\)。
代码里为了方便我将\(m\)减去了\(n\)。
#include<cstdio>
#include<cstring>
const int mod=998244353,inv2=(mod+1)/2;
inline int add(int a,int b)
{
return (a+=b)>=mod?a-mod:a;
}
inline int sub(int a,int b)
{
return (a-=b)<0?a+mod:a;
}
inline int mul(int a,int b)
{
return (long long)a*b%mod;
}
inline int qpow(int a,int b)
{
int res=1;
for(;b;a=mul(a,a),b>>=1)
if(b&1)
res=mul(res,a);
return res;
}
int n,m;
namespace solver1
{
const int N=2005;
int memo[N][N];
inline void init()
{
memset(memo,-1,sizeof(memo));
memo[1][0]=0;
memo[2][0]=1;
memo[3][0]=6;
return;
}
int f(int n,int m)
{
if(m<0)
return 0;
if(n==0)
return 1;
if(~memo[n][m])
return memo[n][m];
int res=0;
//0
res=add(res,f(n,m-1));
//1
res=add(res,mul(mul(n,n+m-1),f(n-1,m)));
if(n>=3)
res=sub(res,mul(mul(n,n+m-1),mul(mul(mul(n-1,n-2),inv2),add(mul(n+m-2,f(n-3,m)),mul(2,f(n-2,m))))));
//2
if(n>=2)
res=add(res,mul(mul(mul(n,n-1),inv2),add(mul(n+m-1,f(n-2,m)),mul(2,f(n-1,m)))));
return memo[n][m]=res;
}
inline void main()
{
init();
printf("%d\n",f(n,m-n));
return;
}
}
namespace solver2
{
const int N=1e5+5;
int memo[N][15];
inline void init()
{
memset(memo,-1,sizeof(memo));
memo[1][0]=0;
memo[2][0]=1;
memo[3][0]=6;
return;
}
int f(int n,int m)
{
if(m<0)
return 0;
if(n==0)
return 1;
if(~memo[n][m])
return memo[n][m];
int res=0;
//0
res=add(res,f(n,m-1));
//1
res=add(res,mul(mul(n,n+m-1),f(n-1,m)));
if(n>=3)
res=sub(res,mul(mul(n,n+m-1),mul(mul(mul(n-1,n-2),inv2),add(mul(n+m-2,f(n-3,m)),mul(2,f(n-2,m))))));
//2
if(n>=2)
res=add(res,mul(mul(mul(n,n-1),inv2),add(mul(n+m-1,f(n-2,m)),mul(2,f(n-1,m)))));
return memo[n][m]=res;
}
inline void main()
{
init();
printf("%d\n",f(n,m-n));
return;
}
}
signed main()
{
scanf("%d%d",&n,&m);
if(n<=2000&&m<=2000)
solver1::main();
else
solver2::main();
return 0;
}
洛谷P4831 Scarlet loves WenHuaKe的更多相关文章
- 洛谷 P1580 yyy loves Easter_Egg I
洛谷 P1580 yyy loves Easter_Egg I 题解: 队列+字符串 #include <cstdio> #include <string> #include ...
- [洛谷2397]yyy loves Maths VI
题目背景 自动上次redbag用加法好好的刁难过了yyy同学以后,yyy十分愤怒.他还击给了redbag一题,但是这题他惊讶的发现自己居然也不会,所以只好找你 题目描述 他让redbag找众数他还特意 ...
- 洛谷P3602 Koishi Loves Segments(贪心,multiset)
洛谷题目传送门 贪心小水题. 把线段按左端点从小到大排序,限制点也是从小到大排序,然后一起扫一遍. 对于每一个限制点实时维护覆盖它的所有线段,如果超过限制,则贪心地把右端点最大的线段永远删去,不计入答 ...
- [CF327E]Axis Walking([洛谷P2396]yyy loves Maths VII)
题目大意:给一个长度为$n(1\leqslant n\leqslant24)$的序列$S$和$k(0\leqslant k\leqslant2)$个数. 求有多少种$S$的排列方式使得其任何一个前缀和 ...
- [洛谷1580]yyy loves Easter_Egg I
题目背景 Soha的出题效率着实让人大吃一惊.OI,数学,化学的题目都出好了,物理的题还没有一道.于是,Huntfire,absi2011,lanlan对soha进行轮番炸,准备炸到soha出来,不料 ...
- 洛谷P2397 yyy loves Maths VI (mode)
P2397 yyy loves Maths VI (mode) 题目背景 自动上次redbag用加法好好的刁难过了yyy同学以后,yyy十分愤怒.他还击给了redbag一题,但是这题他惊讶的发现自己居 ...
- 洛谷P2396 yyy loves Maths VII
P2396 yyy loves Maths VII 题目背景 yyy对某些数字有着情有独钟的喜爱,他叫他们为幸运数字;然而他作死太多,所以把自己讨厌的数字成为"厄运数字" 题目描述 ...
- 洛谷——P2393 yyy loves Maths II
P2393 yyy loves Maths II 题目背景 上次蒟蒻redbag可把yyy气坏了,yyy说他只是小学生,蒟蒻redbag这次不坑他了. 题目描述 redbag给了yyy很多个数,要yy ...
- 洛谷 P4882 lty loves 96! 解题报告
P4882 lty loves 96! 题目背景 众所周知,\(lty\)非常喜欢\(96\)这两个数字(想歪的现在马上面壁去),更甚于复读(人本复)! 题目描述 由于爱屋及乌,因此,\(lty\)对 ...
随机推荐
- JAVA springmvc 转换器
一.有时候springmvc给咱们提供的数据转换并不能转换所有类型比如说由字符串类型转换Date类型,这时候需要我们自定义转换器,帮助springmvc转换我们需要的类型. 二.1)定义我们的转换器: ...
- HDU 2059 龟兔赛跑(超级经典的线性DP,找合适的j,使得每个i的状态都是最好的)
龟兔赛跑 Time Limit:1000MS Memory Limit:32768KB 64bit IO Format:%I64d & %I64u Submit Status ...
- day72
今日内容: 1 创建多表模型(详情见代码) from django.db import models # Create your models here. class Publish(models.M ...
- CAN总线典型特征
CAN总线典型特征 2016-04-12 20:36:54来源: eefocus 关键字:CAN总线 典型特征 收藏 评论(0) 分享到 微博 QQ 微信 LinkedIn CAN总线有如下基本 ...
- sql语句截取字符串
Postgresql 当中有四种方式获取当前时间. 一:now() 通过now()获取的时间是最完整的时间,包括时区,秒也保留到了6位小数. select now(); ...
- 在 R 中估计 GARCH 参数存在的问题(基于 rugarch 包)
目录 在 R 中估计 GARCH 参数存在的问题(基于 rugarch 包) 导论 rugarch 简介 指定一个 \(\text{GARCH}(1, 1)\) 模型 模拟一个 GARCH 过程 拟合 ...
- WPF save listbox config
UI <Grid x:Class="WzlyTool.ReplyContentUI" xmlns="http://schemas.microsoft.com/win ...
- go语言之行--包与变量
一.包的概念 包是go语言中不可缺少部分,在每个go源码的第一行进行定义,定义方式是:package "包名",并且该名称是作为调用该包时候所使用的名称. 包的概念总结: 每个 G ...
- SQLServer数据库还原:无法在已有的mdf文件上还原文件
如果提示无法在已有的mdf文件上还原文件,请修改如下位置
- windows下如何查看进程、端口占用、杀死进程教程
一. 查看所有进程占用的端口 在开始-运行-cmd,输入:netstat –ano 可以查看所有进程 二.查看占用指定端口的程序 当你在用tomcat发布程序时,经常会遇到端口被占用的情况,我们想知道 ...