[JLOI2015]骗我呢

Tags:题解

作业部落

评论地址


TAG:数学,DP

题意

骗你呢

求满足以下条件的\(n*m\)的矩阵的个数对\(10^9+7\)取模

对于矩阵中的第\(i\)行第\(j\)列的元素\(x_{i,j}\)都有

  • \(x_{i,j}<x_{i,j+1}\)
  • \(x_{i,j}<x_{i-1,j+1}\)
  • \(0\le x_{i,j}\le m\)

题解

Part 0 前言

不会做啊!(杠了四五个小时!)

谢两位dalao:blog1blog2

以下图片均来自于此篇文章:http://www.cnblogs.com/coco-night/p/9552677.html,如有冒犯请与我联系,谢谢!

Part 1 朴素DP

首先发现一个很好的性质:

每行是递增的并且一行\(m\)个元素,取值只能在\([0,m]\)中选

那么必然该行至多有一个位置与后一个位置相差2,其余的都只相差1

由此可以列出一个简单的\(DP\):

\(dp[i][j]\)表示第\(i\)行没有出现过的数是\(j\)的方案数

\(dp[i][j]=\sum_{k=0}^{j+1}dp[i-1][k]\)

至于上界为什么是\(j+1\)可以手动模拟一下,假设这行\(j\)没有出现过,上一行试一试\(j-1\)、\(j\)、\(j+1\)、\(j+2\),发现大于\(j+1\)的就不合法了

略微优化一下就变成了\(dp[i][j]=dp[i-1][j+1]+dp[i][j-1]\)

Part 2 转化为图形

发现这个\(DP\)像极了组合数公式,把它套用在坐标系里就是这个样子



自上而下第\(i\)行,从左往右第\(j\)列的点就表示\(dp[i][j]\),其指向的点就表示可以转移

这样仍然不太好处理,我们继续转化:



还是不好看,给它对称一下:

Part 3 挖掘组合意义

这么一看,不就是从原点出发,只能向右或向上走,不接触直线A,B,到达点(n+m+1,n)的路径条数吗!

直线\(A:y=x+1\),直线\(B:y=x-(m+2)\)

Part 4 计算

这种格路数计算(如两双手)都可以考虑采用容斥计数

不考虑其他限制,原点到\(x,y\)的方案数是\(C_{x+y}^x\)

考虑不合法方案是什么:如依次经过\(AABBAAAABB\)

把它缩一下:\(ABAB\)

可以发现不合法方案要么以\(A\)开头要么以\(B\)开头

表示为首次跨越的直线是\(A\)还是\(B\)

所以:答案=总方案数 - A开头的方案数 - B开头的方案数

\(x=n+m+1,y=n\),把\((x,y)\)沿\(A\)对称得到\((x',y')=(y-1,x+1)\)

每条从\((0,0)\)到\((x',y')\)的路径都依次对应一条以A结尾或者以AB结尾的路径!

如图:(这个图是我自己画的!)



上面是一条以\(A\)结尾的路径



上面是一条以\(AB\)结尾的路径

所以,总共的不合法方案是

  • A
  • B
  • AB
  • BA
  • ABA
  • BAB
  • ABAB
  • BABA
  • ...

为了减去以\(A\)开头的方案,需要减去以A,AB结尾的方案,加上以BA,BAB结尾的方案,减去....

那么实现方式是:把(x,y)沿A翻折,减去答案;将翻折过的点沿B翻着,加上答案;再沿A翻折...

同理计算以\(B\)开头的方案,就是先沿\(B\)折就好了

具体细节的话沿着\(A\)折是\((x,y)->(y-1,x+1)\),沿着\(B\)折是\((x,y)->(y+(m+2),x-(m+2))\)

复杂度是O(n)的,复杂度瓶颈为预处理阶乘,至于计算每次是像跳棋一样,复杂度为\(O(logn)\)

完美解决本题!

代码

#include<iostream>
using namespace std;
const int P=1e9+7,N=3e6+10;
int n,m,up,inv[N],jc[N],inj[N];
int Calc(int x,int y) {return (x<0||y<0)?0:1ll*jc[x+y]*inj[x]%P*inj[y]%P;}
void flip1(int &x,int &y) {swap(x,y);x--;y++;}
void flip2(int &x,int &y) {swap(x,y);x+=m+2;y-=m+2;}
void add(int &x,int y) {x+=y;if(x>=P) x-=P;}
int main()
{
cin>>n>>m;inv[0]=inv[1]=jc[0]=inj[0]=1;up=max(n,m)*3+1;
for(int i=2;i<=up;i++) inv[i]=(P-1ll*P/i*inv[P%i]%P)%P;
for(int i=1;i<=up;i++) jc[i]=1ll*jc[i-1]*i%P,inj[i]=1ll*inj[i-1]*inv[i]%P;
int x=n+m+1,y=n,ans=Calc(x,y);
while(x>=0&&y>=0)
{
flip1(x,y);add(ans,P-Calc(x,y));
flip2(x,y);add(ans,Calc(x,y));
}
x=n+m+1,y=n;
while(x>=0&&y>=0)
{
flip2(x,y);add(ans,P-Calc(x,y));
flip1(x,y);add(ans,Calc(x,y));
}
return cout<<ans<<endl,0;
}

[JLOI2015]骗我呢的更多相关文章

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

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

  2. 【BZOJ4005】[JLOI2015]骗我呢

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

  3. bzoj4005[JLOI2015]骗我呢

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

  4. [BZOJ4005][JLOI2015]骗我呢-[dp+容斥]

    Description 传送门 Solution 如果单独考虑一行i,则左边位置的数严格比右边位置的数小.而一行有m个位置,它们可以填[0,m]这m+1个数,则必然有一个数不存在. 定义第i行的第j位 ...

  5. 洛谷 P3266 - [JLOI2015]骗我呢(容斥原理+组合数学)

    题面传送门 神仙题. 首先乍一看此题非常棘手,不过注意到有一个条件 \(0\le x_{i,j}\le m\),而整个矩阵恰好有 \(m\) 列,这就启发我们考虑将每个元素的上下界求出来,如果我们第一 ...

  6. 「JLOI2015」骗我呢 解题报告?

    「JLOI2015」骗我呢 这什么神仙题 \[\color{purple}{Link}\] 可以学到的东西 对越过直线的东西翻折进行容斥 之类的..吧? Code: #include <cstd ...

  7. 【LOJ】#2109. 「JLOI2015」骗我呢

    题解 我深思熟虑许久才算是明白个大概的计数问题吧 先是转化成一个矩形,列一条直线y = x,y = x - (m + 1)我们从(0,0)走到(n + m + 1,m + 1)就是答案 因为我们起始相 ...

  8. BZOJ 4005 [JLOI 2015] 骗我呢

    首先,我们可以得到:每一行的数都是互不相同的,所以每一行都会有且仅有一个在 $[0, m]$ 的数没有出现. 我们可以考虑设 $Dp[i][j]$ 为处理完倒数 $i$ 行,倒数第 $i$ 行缺的数字 ...

  9. [JLOI 2015]骗我呢

    传送门 Description 求给\(n*m\)的矩阵填数的方案数 满足: \[ 1\leq x_{i,j}\leq m \] \[ x_{i,j}<x_{i,j+1} \] \[ x_{i, ...

随机推荐

  1. 【AutoFac】依赖注入和控制反转的使用

    在开始之前首先解释一下我认为的依赖注入和控制反转的意思.(新手理解,哪里说得不正确还请指正和见谅) 控制反转:我们向IOC容器发出获取一个对象实例的一个请求,IOC容器便把这个对象实例“注入”到我们的 ...

  2. C# 创建EXCEL图表并保存为图片

    数据表格能够清晰的呈现数据信息,但是我们对于一些繁杂多变的数据想要很直观的看到数据变化走势或者数据的占比时,数据图表会更具代表性,并且在呈现数据信息上也更形象,也能获取更多纯数字信息所不能直接展现的信 ...

  3. Java开发笔记(十六)非此即彼的条件分支

    前面花了大量篇幅介绍布尔类型及相应的关系运算和逻辑运算,那可不仅仅是为了求真值或假值,更是为了通过布尔值控制流程的走向.在现实生活中,常常需要在岔路口抉择走去何方,往南还是往北,向东还是向西?在Jav ...

  4. 【JVM】7、深入理解Java G1垃圾收集器

    本文首先简单介绍了垃圾收集的常见方式,然后再分析了G1收集器的收集原理,相比其他垃圾收集器的优势,最后给出了一些调优实践. 一,什么是垃圾回收 首先,在了解G1之前,我们需要清楚的知道,垃圾回收是什么 ...

  5. vue2.x 在引用插件的时候,npm run dev跑正常 ,npm run build 报错vue-cli Unexpected token: punc (() [

    这是因为,引用的插件在node_modules里,并不在vue-cli的es6编译范围内,所以语法报错,修改方法:

  6. Mysql基本操作指令集锦

    一.MySQL服务的启动.停止与卸载 在 Windows 命令提示符下运行: 启动: net start MySQL 停止: net stop MySQL 卸载: sc delete MySQL 二. ...

  7. 谷歌浏览器扩展程序manifest.json参数详解

    { // Required "manifest_version": 2, // manifest编写规范版本,目前主流2 "name": "My Ex ...

  8. 美团技术沙龙01 - 58到家服务的订单调度&数据分析技术

    1. 2015.4.15 到家服务的订单调度&数据分析技术 58到家· 黄海斌 @xemoaya 2.agenda • 58到家介绍 • 订单管理系统介绍 • 数据分析技术的应用 3.2015 ...

  9. JButton 按钮,JRadioJButton单选按钮,JChectBox复选框

    一. [按钮JButton] //导入Java类 import javax.swing.*;import java.awt.*;import java.awt.event.ActionEvent;im ...

  10. 2019Java查漏补缺(三)

    1.为什么这个public的类的类名必须和文件名相同    是为了方便虚拟机在相应的路径中找到相应的类所对应的字节码文件    2.java8 的一些新特性:     3: 数据库隔离级别 隔离级别 ...