小 D 特别喜欢玩游戏。这一天,他在玩一款填数游戏。
这个填数游戏的棋盘是一个 n×m 的矩形表格。玩家需要在表格的每个格子中填入一个数字(数字 0 或者数字 1 ),填数时需要满足一些限制。
下面我们来具体描述这些限制。
为了方便描述,我们先给出一些定义:
我们用每个格子的行列坐标来表示一个格子,即(行坐标,列坐标)。(注意: 行列坐标均从 0 开始编号)
一条路径是合法的当且仅当:
(1)这条路径从矩形表格的左上角的格子(0,0)出发,到矩形的右下角格子 (n−1,m−1)结束;
(2)在这条路径中,每次只能从当前的格子移动到右边与它相邻的格子,或者 从当前格子移动到下面与它相邻的格子。
对于一条合法的路径 P,我们可以用一个字符串 w(P) 来表示,该字符串的长度为 n+m−2,其中只包含字符“R”或者字符“D”, 第 i 个字符记录了路径 P 中第 i 步的移动 方法,“ R”表示移动到当前格子右边与它相邻的格子,“ D”表示移动到当前格子下面 与它相邻的格子。
同时,将每条合法路径 P 经过的每个格子上填入的数字依次连接后,会得到一个长 度为 n+m−1 的 01 字符串,记为 s(P) 。
游戏要求小 D 找到一种填数字 0、 1 的方法,使得对于两条路径P1,P2,如果w(P1) > w(P2),那么必须s(P1) ≤ s(P2)。
但是仅仅是找一种方法无法满 足小 D 的好奇心,小 D 更想知道这个游戏有多少种玩法,也就是说,有多少种填数字 的方法满足游戏的要求?
小 D 能力有限,希望你帮助他解决这个问题,即有多少种填 0、1 的方法能满足题目要求。由于答案可能很大,你需要输出答案对 10^9 + 7 取模的结果。
输入
输入共一行,包含两个正整数 n,m,由一个空格分隔,表示矩形的大小。其中 n 表示矩形表格的行数,m 表示矩形表格的列数。
输出
输出共一行,包含一个正整数,表示有多少种填 0、1 的方法能满足游戏的要求。 注意:输出答案对 10^9+7 取模的结果。
输入样例#1
2 2
输出样例#1
12
输入样例#2
3 3
输出样例#2
112
输入样例#3
5 5
输出样例#3
7136
数据规模与约定
n<=8,m<=10^6。
傻逼规律题,我竟然不会做。
我太弱了……
题目中的限制等价于:对于所有的点,如图示的红色路径的数字串要 ≥ 蓝色路径的数字串。
我们再把这个限制进行转换:
(1)每一条自右上至左下的对角线的数字总是先出现一堆 0,再出现一堆 1。
(2)如果 (a, b+1) = (a+1, b),则以(a+1, b+1)为左上角,以(n-1, m-1)为右下角的矩形,每一条自右上至左下的对角线的数字全部相同。
可以证明以上的转换都是等价的。
然后就会发现一件很有意思的事情:对于这条对角线 2,肯定是000,001,011,111。这条对角线上一定有 2 个元素相同,也就说限制(2)会限制掉很大一个区域!
于是我们就可以根据这个设计算法了:
先手推 n = 1,n = 2, n = 3 的情况,将它们特判掉。
当 n = 1 时,ans = 2m。
当 n = 2 时,ans = 4*3m-1。
当 n = 3 时,ans = 112*3m-3。
然后:如果交换 n 和 m ,答案不会变。所以我们不妨假设 n <= m。
计算出上图中对角线 1 为 00, 11 的方案数。除此之外对角线 1 一定为 01。
对于对角线 2,先考虑它为 000,111 的情况,然后考虑它为 001,011 的情况。
接下来考虑如何具体计数。
假如对角线 1 为 00, 11:
将剩余的对角线分为三类:两个端点没有被限制(形如 1),对答案的贡献为 4;两个端点被限制了一个(形如 2),对答案的贡献为 3;两个端点都被限制(形如 3),对答案的贡献为 2。将所有对角线的贡献乘起来即可。可以用快速幂。
假如对角线 2 为 000, 111:
图中所标出的对角线是不受限制的。蓝色限制与紫色限制相互组合,可以发现和上面那种情况是差不大多的(除了边缘的行列都被限制)。因此,依然还是按照上文所提到的方法进行计数。
假如对角线 2 为 001,011。因为是对称的,我们只需要弄清 001 的求解方式即可类比推导出 011 的求解方法。
紫色是限制区域。如果绿色矩形中有一个对角线相同,则限制区域会如下图所示。
对比上面那种情况,我们只是将蓝色矩形往下平移了一点。具体的计数方式还是跟上面那种情况是一样的。
还要讨论绿色矩形内对角线都互不相同的方案数。
注意要分类讨论 n = m 和 n < m 的 情况。
我觉得我的题解可能描述的不是很具体……因为把这道题解释清楚是一件对我来说不太容易的事情……但是我觉得关键的地方我都写上了。
还是请大家多多见谅啦,如果有疑问就留言提出来吧,我会尽力解答的。
#include
#include
using namespace std;
const int MOD = int(1E9) + 7;
int pow_mod(int b, int p) {
int ret = 1;
while( p ) {
if( p & 1 ) ret = 1LL*ret*b%MOD;
b = 1LL*b*b%MOD;
p >>= 1;
}
return ret;
}
int main() {
int n, m;
scanf("%d%d", &n, &m);
if( n > m ) swap(n, m);
if( n == 1 ) {
printf("%d\n", pow_mod(2, m));
return 0;
}
else if( n == 2 ) {
printf("%lld\n", 4LL*pow_mod(3, m-1)%MOD);
return 0;
}
else if( n == 3 ) {
printf("%lld\n", 112LL*pow_mod(3, m-n)%MOD);
return 0;
}
else {
int ans = 0;
if( n == m ) {
ans = (ans + 2LL*pow_mod(4, n-2)%MOD*pow_mod(2, n-1)%MOD)%MOD;
ans = (ans + 2LL*5LL*pow_mod(4, n-4)%MOD*pow_mod(2, n-1)%MOD)%MOD;
for(int i=4;i<=m;i++) {
if( i == m ) ans = (ans + 4LL*3LL*pow_mod(2, n-2)%MOD)%MOD;
else ans = (ans + 4LL*5LL*pow_mod(4, n-i-1)%MOD*pow_mod(2, n-1)%MOD)%MOD;
}
ans = (ans + 3LL*pow_mod(2, n-2))%MOD;
for(int i=4;i<=n;i++) {
if( i == n ) ans = (ans + 4LL*3LL*pow_mod(2, m-2)%MOD)%MOD;
else ans = (ans + 4LL*5LL*pow_mod(4, m-i-1)%MOD*pow_mod(2, m-1)%MOD)%MOD;
}
ans = (ans + 3LL*pow_mod(2, n-2))%MOD;
}
else {
ans = (ans + 2LL*pow_mod(4, n-2)%MOD*pow_mod(3, m-n)%MOD*pow_mod(2, n-1)%MOD)%MOD;
ans = (ans + 2LL*5LL*pow_mod(4, n-4)%MOD*pow_mod(3, m-n)%MOD*pow_mod(2, n-1)%MOD)%MOD;
for(int i=4;i<=m;i++) {
if( i == m ) ans = (ans + 3LL*3LL*pow_mod(2, n-2)%MOD)%MOD;
else if( i > n ) ans = (ans + 3LL*4LL*pow_mod(3, m-i-1)%MOD*pow_mod(2, n-1)%MOD)%MOD;
else if( i == n ) ans = (ans + 4LL*4LL*pow_mod(3, m-n-1)%MOD*pow_mod(2, n-1)%MOD)%MOD;
else if( i < n ) ans = (ans + 4LL*5LL*pow_mod(4, n-i-1)%MOD*pow_mod(3, m-n)%MOD*pow_mod(2, n-1)%MOD)%MOD;
}
ans = (ans + 3LL*pow_mod(2, n-2))%MOD;
for(int i=4;i<=n;i++) {
if( i == n ) ans = (ans + 4LL*4LL*pow_mod(3, m-n-1)%MOD*pow_mod(2, n-1)%MOD)%MOD;
else ans = (ans + 4LL*5LL*pow_mod(4, n-i-1)%MOD*pow_mod(3, m-n)%MOD*pow_mod(2, n-1)%MOD)%MOD;
}
ans = (ans + 4LL*pow_mod(3, m-n-1)%MOD*pow_mod(2, n-1)%MOD)%MOD;
}
printf("%d\n", 2LL*ans%MOD);//*2:左上角的填数方式
}
}
就是这样,新的一天里,也请多多关照哦(ノω<。)ノ))☆.。