51nod 1119 机器人走方格 V2（组合数）

1119 机器人走方格 V2 

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

M * N的方格，一个机器人从左上走到右下，只能向右或向下走。有多少种不同的走法？由于方法数量可能很大，只需要输出Mod 10^9 + 7的结果。

Input

第1行，2个数M,N，中间用空格隔开。（2 <= m,n <= 1000000)

Output

输出走法的数量 Mod 10^9 + 7。

Input示例

2 3

Output示例

3

思路：当n,m很小的时候，这是一道很简单的dp，但是现在（2 <= m,n <= 1000000)，无法用dp做。可以用组合数做，从起点到终点一定是走n+m-2步,其中向下n-1步，向左走m-1步，总的组合数就是C(m+n−2,n−1)或C(m+n−2,m−1).

那么对于带模的出发a/b=X(mod M),只要 M 是一个素数，而且 b 不是 M 的倍数，就可以用一个逆元整数 b’，通过 a / b = a * b' (mod M)，来以乘换除。

根据费马小定理说，对于素数 M 任意不是 M 的倍数的 b，都有：b ^ (M-1) = 1 (mod M)
于是可以拆成：b * b ^ (M-2) = 1 (mod M)
于是：a / b = a / b * (b * b ^ (M-2)) = a * (b ^ (M-2)) (mod M)
b的逆元其实就是 b ^ (M-2) (mod M) ——（用快速幂可以求出）

这道题就是根据这么算的。

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
#define ll long long
const ll  mod=1e9+7;
ll n,m;
ll qpow(ll x,ll y){//快速幂 
	long  long sum=1;
	x=x%mod;
	while(y>0){
		if(y&1)
			sum=sum*x%mod;
		x=x*x%mod;
		y>>=1;
	}
	return sum;
}
int main()
{
	scanf("%lld%lld",&m,&n);
	if(n>m) swap(n,m);
	n-=1,m-=1;
	ll a=1,b=1,c=n,d=m+n;
	while(c) {//分子 
		a=(a*d)%mod;
		d--;c--;
	}
	while(n) {//分母 
		b=(b*n)%mod;
		n--;
	}
	ll b1=qpow(b,mod-2);//分母的逆元 
	printf("%lld\n",a*b1%mod);//以乘换除 
	return 0;
}