51node-1627 瞬间移动

原题链接

1627 瞬间移动

基准时间限制：1 秒 空间限制：131072 KB 分值: 80  难度：5级算法题

 收藏

 关注

有一个无限大的矩形，初始时你在左上角（即第一行第一列），每次你都可以选择一个右下方格子，并瞬移过去（如从下图中的红色格子能直接瞬移到蓝色格子），求到第n行第m列的格子有几种方案，答案对1000000007取模。

Input

单组测试数据。
两个整数n,m（2<=n,m<=100000）

Output

一个整数表示答案。

Input示例

4 5

Output示例

10

枚举到(n,m)的步数k。
k步走到(n,m)的方案=C（n-2，k-1）*C（m-2，k-1）（横坐标增加k次，总和为n-1的方案 * 纵坐标增加k次，总和为m-1的方案）

预处理1！~n！的值计算组合数。所以只需要枚举步数即可。

#include <iostream>
#include <algorithm>
#include <cstring>
#include <cstdio>
#include <vector>
#include <map>
#include <cmath>
#include <queue>
#include <ctime> 
#define maxn 100005
#define INF 1e18
#define MOD 1000000007
using namespace std;
typedef long long ll;

ll d[maxn];
ll solve(ll p){
	
	ll m = MOD - 2, ans = 1;
	while(m){
		if(m&1)(ans *= p) %= MOD;
		(p *= p) %= MOD;
		m >>= 1;
	}
	return ans;
}
int main(){
	
//	freopen("in.txt", "r", stdin);
	d[0] = 1;
	for(int i = 1; i <= 100000; i++)
	  d[i] = d[i-1] * i % MOD;
	int n, m;
	
	scanf("%d%d", &n, &m);
	ll ans = 0;
	int k = min(n - 1, m - 1);
	int p1 = n - 2, p2 = m - 2;
	for(int i = 1; i <= k; i++){
		
		ans += d[p1] * solve(d[i-1]  * d[p1-i+1] % MOD) % MOD 
		        * d[p2] %MOD * solve(d[i-1] * d[p2-i+1] % MOD) % MOD;
		ans %= MOD;
	}
	printf("%I64d\n", ans);
	return 0;
}