题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5698
题意说我们可以从一个格子跳到右下角的任意一个格子中,那么反过来,跳到某一个格子的方法数显然就是其左上角的所有格子的方法数的和,那么问题又来了,我们怎么才能求出一个格子左上角的所有的格子的方法数和呢,树状数组?线段树?答案都不是,或许可以,不过我们并不需要那么做,观察仔细的同学应该发现每一个格子的方法数就等于其左边和上面的格子的方法数的和,这又是为什么呢?举个例子,比如我们要求的是(x,y)的方法数,那么其左边的格子的方法数和右边格子的方法数分别是以(1,1)和(x-2,y)的矩形和以(1,1)和(x,y-2)的矩形内的所有点的方法数和,不难发现这两个矩形重叠的部分就是以(1,1)到(x-2,y-2)的矩形,刚刚好相当于点(x-2,y-2)的方法数,所以这两个矩形的内的所有点的方法数和就等于(x,y)的方法数。
推出了每一个格子的方法数就等于其左边和上面的格子的方法数的和这个结论后,再仔细思考,我们能够直接通过这个结论得到答案吗,当然不行!!! 1e5的数据范围二维数组连开都开不出来,必须得用其他方法。
我们再来看这个序列,因为有每一个格子的方法数就等于其左边和上面的格子的方法数的和这个性质,所以我们很容易可以把这个和杨辉三角扯上关系,这个题就相当于一个横着的杨辉三角。
那么我们来看看杨辉三角
1
1 1
1 2 1
1 3 3 1
1 4 6 4 1
1 5 10 5 1
……
这个杨辉三角的每一行相当于(x+1)^(row-1)的每一项的系数,那么n行m列的值就等于C(n-1, m-1)。
再来看这个图:
1 0 0 0 0 0
0 1 1 1 1 1
0 1 2 3 4 5
0 1 3 6 10 15
0 1 4 10 20 35
0 1 5 15 35 70
显然从(2,2)右下方开始就是一个横着的杨辉三角,把现在的行和列换成杨辉三角的行和列就可以得到,(x,y)的方法数就等于C(n+m-4, m-2)。那么就等于(n+m-4)! / ((n-2)! * (m-2)!),答案到这里就很接近了,最后因为是要取余,所以别忘了求逆元。
#include <cstdio>#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define LL __int64
const int MOD = 1000000007;
void exGcd(LL a, LL b, LL &x, LL &y)
{
if(b == 0)
{
x = 1;
y = 0;
return ;
}
exGcd(b, a%b, y, x);
y -= (a/b)*x;
}
int main()
{
LL n,m;
while(scanf("%I64d%I64d",&n,&m)!=EOF)
{
LL ans = 1;
for(int i=1; i<=n+m-4; i++)
ans = (ans*i) % MOD;
for(int i=1; i<=n-2; i++)
{
LL x, y;
exGcd(i, MOD, x, y);
x = (x%MOD + MOD) % MOD;
ans = (ans*x) % MOD;
}
for(int i=1; i<=m-2; i++)
{
LL x, y;
exGcd(i, MOD, x, y);
x = (x%MOD + MOD) % MOD;
ans = (ans*x) % MOD;
}
printf("%I64d\n",ans);
}
}