51nod 1119 组合数

这道题因为数据特别大 递归明显不适用 不过若数据在一百左右可以采用以下的递归模板
51nod 1119

 - void  dfs(int x,int y)
 - {    
 -       if(x=1&&j==1)    return 1;
 -       return  dfs(x+1,y)+dfs(x,y-1);
 - }       

不过呢 我们再看这道题的输入数据 （2 <= m,n <= 1000000) 显然使用递归会超时的
我们在看题意 机器人只能够向下走或者向右走 那么向下就有N-1种 向右就有M-1种 从起点到终点总共有
C(N+M-2,N-1）或者C(N+M-2,M-1)都是一样的
规律是找到了 那么接下来怎么计算这个组合呢？C(N+M-2,N-1）% P

我们知道取余法则里面并没有除法 这时我们就需要用到乘法取逆元 （A/B)%P=(A*-A)%P
这里我们引进一个费马定理：任何一个质数P与不是它的倍数B存在这样的关系：B^P≡1（modP)
那么（A/B)%P≡A/BB(B*B^M-2)≡1(MOD P)
那么B^M-2即为A的乘法逆元
代码

#include<bits/stdc++.h>

using namespace std;
const int p=1e9+7;

long long quick(long long a,long long n)
{
    long long res=1;
    while(n)
    {
        if(n&1)
            res=(res*a)%p;
        a=(a*a)%p;
        n>>=1;
    }
    return res;
}
int main()
{
    int n,m;
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        if(n>m)
            swap(n,m);
        long long a=1,b=1,c=n-1,d=n+m-2;//计算排列组合数C(N+M-2,N-1)
        while(c)//
        {
            a=(a*d)%p;
            d--;
            c--;
        }
        n-=1;
        while(n)
        {
            b=(b*n)%p;
            n--;
        }
        long long x=quick(b,p-2);
        printf("%lld\n",((a%p)*(x%p))%p);
    }
    return 0;
}