逆元 快速幂 (b/a)%p=(b%p)*(a^-1%p)=(b%p)*(a^p-2%p)

题目描述
多多喜欢行走，有一天老师问他一个问题：在一个方格点阵中，左上角点的坐标为(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)

#include<iostream>
#include<cstdio>
using namespace std;
typedef long long ll;
const int maxn=2000100;
const int mod=1000000007;

ll fact[maxn],inv[maxn];
int qpow(ll a,ll b)
{
	ll ans=1;
	while(b!=0)
	{
		if(b&1) ans=(ans*a)%mod;
		a=(a*a)%mod;
		b>>=1;
	}
	return ans;
}
void init()
{
	fact[0]=1;
	for(int i=1;i<maxn;i++){
		fact[i]=(fact[i-1]*i)%mod;
	}
	inv[maxn-1]=qpow(fact[maxn-1],mod-2);
	for(int i=maxn-2;i>=0;i--){
		inv[i]=(inv[i+1]*(i+1))%mod;
	}
}
ll C(int n,int m)
{
	return fact[n]*inv[n-m]%mod*inv[m]%mod;
}
int main()
{
	init();
	int n,m;
	int t;
	scanf("%d",&t);
	while(t--)
	{
		scanf("%d%d",&n,&m);
		printf("%lld\n",C(n+m-2,n-1));
	}
    return 0;
}