题面传送门

神仙题。

首先乍一看此题非常棘手,不过注意到有一个条件 \(0\le x_{i,j}\le m\),而整个矩阵恰好有 \(m\) 列,这就启发我们考虑将每个元素的上下界求出来,如果我们第一列全填 \(0\),其余每个数都恰好等于它左边的数加 \(1\),那么 \(x_{i,j}\) 刚好取到下界 \(j-1\);如果我们最后一列全填 \(m\),其余每个数都恰好等于它右边的数减 \(1\),那么 \(x_{i,j}\) 刚好取到上界 \(j\),因此对于任意一个第 \(j\) 列的元素 \(x_{i,j}\),它的取值只有两种:\(j-1\) 和 \(j\)。

我们记一类格为填 \(j-1\) 的 \(x_{i,j}\),二类格为填 \(j\) 的 \(x_{i,j}\),那么不难发现对于一个二类格 \(x_{i,j}\),所有形如 \(x_{i-p,j+q}\),\(p\le q\) 的格子(也就是下图中的蓝色区域)都必须是二类格,而不难发现对于任意一个合法的矩阵,都存在唯一的二类格的集合 \(S\),满足任意两个二类格都不在互相所管辖的范围内。因此题目可以转化为,有多少个二类格的集合 \(S\),满足任意两个格子都不在互相管辖的范围内。

考虑 \(dp\),不难发现任意一行,一类格必定是一段后缀,因此记 \(dp_{i,j}\) 表示从下往上考虑到了第 \(i\) 行,一类格前缀长度为 \(j\) 的方案数,那么显然上一行一类格前缀的长度 \(\le j+1\),因此我们不难得到 \(dp\) 转移方程式 \(dp_{i,j}=\sum\limits_{k=j-1}^{m}dp_{i+1,k}\),上式稍微化简一下可以得到 \(dp_{i,j}=dp_{i,j+1}+dp_{i+1,j-1}\),注意,对于 \(dp_{i,m}\) 而言,上式只有 \(dp_{i+1,m-1}\),而实际上 \(dp_{i+1,m}\) 也能转移到 \(dp_{i,m}\),故 \(dp_{i,m}=dp_{i+1,m-1}+dp_{i+1,m}\)。

这样 \(dp\) 是 \(nm\) 的,无法通过,考虑使用组合意义优化,对于 \(n=3,m=3\) 而言,该 \(dp\) 的值等价于下图(这里借用了张题解区的图)中从左上角走到右下角的方案数(提示:如果将整张图旋转 \(\pi\) 那可能比较好理解,因为上面的 \(dp\) 过程是从下往上推的)

稍微将它变形一下可以得到下图:

不难发现这东西就等价于从 \((0,0)\) 出发到达 \(P(n+m+1,n)\) 的方案数,其中不能碰到直线 \(A:y=x+1\) 和直线 \(B:y=x-(m+2)\)

这个东西怎么求呢?一个比较棘手的地方是它涉及两条直线,如果只涉及一条直线那可以像我们的经典问题——求卡特兰数递推式那样做一个对称然后简单求个组合数。因此这里介绍一种思考问题的方法:前缀容斥。注意到对于所有不合法情况,它经过直线 \(AB\) 的情况必然构成一个序列,比方说 \(AABBAABBBA\),把连续段缩一下可以得到 \(ABABA\),我们不妨就从这个缩好的序列入手计算方案数。显然对于每个不合法的序列,它缩好的序列要么以 \(A\) 开头要么以 \(B\) 开头,因此答案就是总方案数减去以 \(A\) 开头的方案数减去以 \(B\) 开头的方案数。

怎样求以 \(A\) 开头的方案数呢?我们不妨把以 \(A\) 开头的方案罗列一下,可以得到:

  • A
  • AB
  • ABA
  • ABAB
  • ...

看到以 \(A\) 开头我们直观地想到做 \(P(n+m+1,n)\) 关于直线 \(A\) 的对称点 \(P'\),但事实上这是不对的。In fact,仔细分析一下就可以发现,对于所有 \((0,0)\) 到 \(P'\) 的路径,如果我们把最后一个与 \(A\) 的交点 \(X\) 找到,然后把 \(X\to P'\) 的折线翻下来,那么最后一段 \(X\to P'\) 必然与 \(A\) 没有交点,因此这样的路径末尾要么是 \(A\)(最后一段与 \(B\) 没有交点)要么是 \(AB\)(最后一段与 \(AB\) 有交点),因此求得的方案数是以 \(A\) 或 \(AB\) 为结尾的方案数,但这样只是以 \(A\) 或 \(AB\) 结尾的啊,还会多算什么 \(BA,BAB,ABA,ABAB\cdots\),别急,如果我们再作 \(P'\) 关于 \(B\) 的对称点 \(P''\),那么所有 \((0,0)\to P''\) 的路径必然经过 \(B\),再把它翻上来就得到了 \((0,0)\to P'\) 的路径,由已知 \((0,0)\to P'\) 所对应的 \((0,0)\to P\) 的路径要么以 \(A\) 为结尾,要么以 \(AB\) 为结尾,前面再填个 \(B\),故 \((0,0)\to P''\) 所对应的路径要么以 \(BA\) 为结尾,要么以 \(BAB\) 结尾,发现就是上面 \(BA,BAB,ABA,ABAB\cdots\),二者相减就得到了 \(A,AB\) 的方案数。

求 \(ABA,ABAB\),以及后面 \(ABABA,ABABAB\cdots\) 的方案数也同理,只需要再做 \(P''\) 关于 \(A\) 的对称点计算方案数,再关于 \(B\) 对称计算方案数,二者相减即可得到 \(ABA,ABAB\) 的方案数,以此类推即可求出以 \(A\) 开头的方案数。求以 \(B\) 开头的也同理。

时间复杂度 \(\mathcal O(n)\)

const int MAXN=6e6+4;
const int MOD=1e9+7;
int n,m,fac[MAXN+5],ifac[MAXN+5];
void init_fac(int n){
for(int i=(fac[0]=ifac[0]=ifac[1]=1)+1;i<=n;i++) ifac[i]=1ll*ifac[MOD%i]*(MOD-MOD/i)%MOD;
for(int i=1;i<=n;i++) fac[i]=1ll*fac[i-1]*i%MOD,ifac[i]=1ll*ifac[i-1]*ifac[i]%MOD;
}
int binom(int x,int y){return 1ll*fac[x]*ifac[y]%MOD*ifac[x-y]%MOD;}
void flip(int &x,int &y,int a){x^=y^=x^=y;x-=a;y+=a;}
int ways(int x,int y){return (x<0||y<0)?0:binom(x+y,x);}
int main(){
scanf("%d%d",&n,&m);init_fac(MAXN);
int x=n+m+1,y=n,ans=ways(x,y);
while(x>=0&&y>=0){
flip(x,y,1);ans=(ans-ways(x,y)+MOD)%MOD;
flip(x,y,-m-2);ans=(ans+ways(x,y))%MOD;
} x=n+m+1,y=n;
while(x>=0&&y>=0){
flip(x,y,-m-2);ans=(ans-ways(x,y)+MOD)%MOD;
flip(x,y,1);ans=(ans+ways(x,y))%MOD;
} printf("%d\n",ans);
return 0;
}

洛谷 P3266 - [JLOI2015]骗我呢(容斥原理+组合数学)的更多相关文章

  1. 【BZOJ4005】[JLOI2015] 骗我呢(容斥,组合计数)

    [BZOJ4005][JLOI2015] 骗我呢(容斥,组合计数) 题面 BZOJ 洛谷 题解 lalaxu #include<iostream> using namespace std; ...

  2. [JLOI2015]骗我呢

    [JLOI2015]骗我呢 Tags:题解 作业部落 评论地址 TAG:数学,DP 题意 骗你呢 求满足以下条件的\(n*m\)的矩阵的个数对\(10^9+7\)取模 对于矩阵中的第\(i\)行第\( ...

  3. 洛谷 P5400 - [CTS2019]随机立方体(组合数学+二项式反演)

    洛谷题面传送门 二项式反演好题. 首先看到"恰好 \(k\) 个极大值点",我们可以套路地想到二项式反演,具体来说我们记 \(f_i\) 为钦定 \(i\) 个点为极大值点的方案数 ...

  4. 洛谷 P1763 状态压缩dp+容斥原理

    (题目来自洛谷oj) 一天,maze决定对自己的一块n*m的土地进行修建.他希望这块土地共n*m个格子的高度分别是1,2,3,...,n*m-1,n*m.maze又希望能将这一些格子中的某一些拿来建蓄 ...

  5. 洛谷 P3263 [JLOI2015]有意义的字符串

    洛谷 首先,看到\((\frac{(b+\sqrt{d})}{2})^n\),很快能够想到一元二次方程的解\(\frac{-b\pm\sqrt{\Delta}}{2a}\). 所以可以推出,\(\fr ...

  6. 洛谷 P6276 - [USACO20OPEN]Exercise P(组合数学+DP)

    洛谷题面传送门 废了,又不会做/ll orz czx 写的什么神仙题解,根本看不懂(%%%%%%%%% 首先显然一个排列的贡献为其所有置换环的乘积.考虑如何算之. 碰到很多数的 LCM 之积只有两种可 ...

  7. 洛谷 P4708 - 画画(Burnside 引理+组合数学)

    洛谷题面传送门 神仙题 %%%%%%%%%%%%%%%%%%%% 题解搬运人来了 首先看到本质不同(无标号)的图计数咱们可以想到 Burnside 引理,具体来说,我们枚举一个排列 \(p\),并统计 ...

  8. 【BZOJ4005】[JLOI2015]骗我呢

    题意: Alice和Bob在经过了数学的洗礼之后,不再喜欢玩对抗游戏了,他们喜欢玩合作游戏.现在他们有一个n×m的网格,Alice和Bob要在一定规则下往网 格里填数字,Alice和Bob都是聪明绝顶 ...

  9. bzoj4005[JLOI2015]骗我呢

    http://www.lydsy.com/JudgeOnline/problem.php?id=4005 神题~远距离orz 膜拜PoPoQQQ大神 #include<cstdio> #i ...

随机推荐

  1. 【UE4 C++】抛物线路径、发射轨道相关

    基于UGameplayStatics Blueprint_PredictProjectilePath_ByObjectType 根据 Object Type,算出抛物线的点集合和检测结果 static ...

  2. Alpha-技术规格说明书

    项目 内容 这个作业属于哪个课程 2021春季软件工程(罗杰 任健) 这个作业的要求在哪里 团队项目-计划-功能规格说明书 一.架构与技术栈 1.整体架构 本项目的整体架构如上图所示.下面我们将对涉及 ...

  3. 2021.8.16考试总结[NOIP模拟41]

    T1 你相信引力吗 肯定是单调栈维护.但存在重复值,还是个环,不好搞. 发现取区间时不会越过最大值,因此以最大值为断点将环断为序列.在栈里维护当前栈中有多少个与当前元素相等的元素,小分类讨论一下. 最 ...

  4. Python课程笔记(一)

    由于新冠状病毒的爆发,不得不在家里上网课,开课已经两个礼拜了,今天上完Python课后,准备整理一下最近学习Python的笔记. 人生苦短,我用Python 一.Hello World 初学一门新的语 ...

  5. 攻防世界 杂项 6.pure_color

    图片隐写 工具 使用StegSolve一把梭 另一种解法 右击图片编辑,画图工具打开,属性设置黑白.

  6. loadRunner运行场景时,事务数为0或是只显示添加的事务的数

    脚本编辑好后,不要着急到controller去执行,注意查看Run-time Settings(运行是设置)-->General(常规)-->Miscellaneous(其他)中查看Aut ...

  7. centos如何上网

    问题 centos如何上网 解决方法 第一步: 打开VMware,选择菜单栏的Edit->Virtual Network Editor(虚拟网络编辑器).点击Restore Defaults(即 ...

  8. 用 Node.js 实现的最简单的 HTTP 服务器

    用 Node.js 实现的最简单的 HTTP 服务器 //app.js var http = require('http'); http.createServer(function(req, res) ...

  9. celery tasks always in pending

    Result backend doesn't work or tasks are always in PENDING state¶All tasks are PENDING by default, s ...

  10. Part 31 AngularJS page refresh problem

    What is the issue : When you navigate to http://localhost/students, you will see list of students as ...