hdu5968(组合数取模Lucas定理)

瞬间移动

Time Limit: 4000/2000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others)
Total Submission(s): 1460    Accepted Submission(s): 703

Problem Description

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

 

Input

多组测试数据。

两个整数 n,m(2≤n,m≤100000)

 

Output

一个整数表示答案

 

Sample Input

  
  

   
   4 5
  
  

 

Sample Output

  
  

   
   10
  
  

 

思路：这个矩形形成了一个杨辉三角，坐标（n,m）对应杨辉三角中的坐标（n+m-3,m-1），又因为 杨辉三角中第n行的m个数可表示为 C(n-1，m-1)，所以结果就为C(n+m-4,m-2)，再使用Lucas定理求出组合数取模即可。

#include<iostream>
#include<stdio.h>
#include<string.h>
#include<math.h>
using namespace std;
typedef long long LL;
#define MOD  1000000007
LL n,m;
LL p=MOD;
LL quick_mod(LL a, LL b)  
{  
    LL ans = 1;  
    a %= p;  
    while(b)  
    {  
        if(b & 1)  
        {  
            ans = ans * a % p;  
            b--;  
        }  
        b >>= 1;  
        a = a * a % p;  
    }  
    return ans;  
}  
  
LL C(LL n, LL m)  
{  
    if(m > n) return 0;  
    LL ans = 1;  
    for(int i=1; i<=m; i++)  
    {  
        LL a = (n + i - m) % p;  
        LL b = i % p;  
        ans = ans * (a * quick_mod(b, p-2) % p) % p;  
    }  
    return ans;  
}  
  
LL Lucas(LL n, LL m)  
{  
    if(m == 0) return 1;  
    return C(n % p, m % p) * Lucas(n / p, m / p) % p;  
}
int main()
{
	
	while(~scanf("%lld%lld",&n,&m))
	{
		printf("%lld\n",Lucas(n+m-4,m-2));
	}
}