牛客-Walk （组合数学+逆元）

链接：https://ac.nowcoder.com/acm/contest/5929/K
来源：牛客网
 

题目描述

多多喜欢行走，有一天老师问他一个问题：在一个方格点阵中，左上角点的坐标为(1, 1)，行坐标从上到下依次递增，列坐标从左到右依次递增，每次行走可以向上、下、左、右移动一格。现在要从(1, 1)点走到(N, M)点，在行走步数最少的情况下，有多少种行走方法？(答案可能过大，请对答案取模1000000007)

输入描述:

第一行输入一个正整数 T，代表询问次数 (1 ≤ T ≤ 100000)

接下来 T 行，每行输入两个正整数 N，M (1 ≤ N ≤ 106，1 ≤ M ≤ 106)

输出描述:

对于每次询问，输出一个正整数，代表在行走步数最少的情况下，从(1, 1)点走到(N, M)点的方法总数 (答案可能过大，请对答案取模1000000007)

 

解题思路：预处理+逆元

从(1, 1)到(n, m)的最短路为向下走n – 1步，向右走m –1步，所以方法总数为C(n + m – 2, n – 1) 

因为是多组输入样例T很大，所以要进行预处理。因为这里1e9+7是一个大质数，所以可以采取下面的方法求解（逆元）：

对 n!%p 进行预处理，即 f[i] 表示 i!%p；

对 (n-m)!^(p-2)%p 进行预处理，即 a[i] 表示 i!^(p-2)%p。

取余运算规则：a ^ b % p = ((a % p)^b) % p 

附上多种方法求解排列组合链接：https://blog.csdn.net/zengaming/article/details/63681754

#include <iostream>
#include <cstdio>//a ^ b % p = ((a % p)^b) % p 
#include <cstring>
#include <algorithm>
#include <cmath>
#include <map>
#include <set>
#include <stack>
#include <queue>
#include <vector>
#include <string>
#define cla(a, sum) memset(a, sum, sizeof(a))
#define rap(i, m, n) for(int i=m; i<=n; i++)
#define rep(i, m, n) for(int i=m; i>=n; i--)
#define bug printf("???\n")
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<ll, ll> P;
const int Inf = 0x3f3f3f3f;
const double eps = 1e-8;
const int maxn = 2e6+5;
const int MOD=1e9+7;
const int N=1e6;
template <typename T> void read(T &x){
    x = 0; int f = 1; char ch = getchar();
    while (!isdigit(ch)) {if (ch == '-') f = -1; ch = getchar();}
    while (isdigit(ch)) {x = x * 10 + ch - '0'; ch = getchar();}
    x *= f;
}

ll t,n,m;
ll f[maxn],a[maxn];

ll quick(ll x,ll y){
	ll ans=1;
	while(y){
		if(y&1){
			ans=ans*x%MOD;
		}
		x=x*x%MOD;
		y/=2;
	}
	return ans%MOD;
}

void init(){//打表预处理 
	f[0]=f[1]=1;
	for(ll i=2;i<=N*2;i++){
		f[i]=(f[i-1]*i)%MOD;
	}
	a[0]=a[1]=1;
	for(ll i=2;i<=N;i++){
		a[i]=quick(f[i],MOD-2);
	}
}

int main()
{
	init();
    cin>>t;
    while(t--){
    	cin>>n>>m;
    	ll sum=f[n+m-2]%MOD*a[n-1]%MOD*a[m-1]%MOD;
    	cout<<sum<<endl;
	}
	return 0;
}