hdu5698瞬间移动-（杨辉三角+组合数+乘法逆元）

瞬间移动

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

Problem Description

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

行第m

列的格子有几种方案，答案对1000000007

取模。

 

Input

多组测试数据。

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

 

Output

一个整数表示答案

 

Sample Input

4 5

 

Sample Output

10

 

 

解题过程：

先看一下杨辉三角的图：

矩阵从a【1】【1】开始，先枚举题目的少数答案：

0  0  0  0  0  0  0
0        1   
0      3      6
0    3    15 21
0  1    20 35 56
0  1  5  15 35 70 126
0    6  21 56 126 252 显然斜着看是一个杨辉三角。

从左上到右下看作一行一行，从左下到右上数该行第几个。
C（x,y） 斜着看，第x行y个
杨辉三角有效部分可用组合数表示为
(0,0) (1,1) (2,2) (3,3) (4,4) (5,5)
(1,0) (2,1) (3,2) (4,3) (5,4)
(2,0) (3,1) (4,2) (5,3)
(3,0) (4,1) (5,2)
(4,0) (5,1)
(5,0)

输入n，m表示n行m列。对应到组合数里。选择杨辉三角部分。
a[n][m]
(2,2) (2,3) (2,4) (2,5) (2,6)
(3,2) (3,3) (3,4) (3,5) (3,6)
(4,2) (4,3) (4,4) (4,5) (4,6)
(5,2) (5,3) (5,4) (5,5) (5,6)
(6,2) (6,3) (6,4) (6,5) (6,6)

 对于从左下到右上这一斜线，n+m相等。

 

将矩形数组a[n][m]和组合数C(x,y)联系起来，我们是利用组合数C(x,y)来算答案,故用x和y表示输入的n和m

对于纵坐标，永远都是y=m-2

对于横坐标，捉摸不定，找出第一列的规律，x=n-2

利用矩阵中n+m相等的条件，推出C(x,y)=C(n+m-4,m-2)。

 

由于n和m巨大，并且题目中的模数p为素数。

求组合数，C(n,m)%p = n! / ( m!*(n-m)! ) %p

用乘法逆元配合快速幂可以解决，本题还不需要用到卢卡斯定理。

 

 AC代码：

 #include<iostream>
 #include<stdio.h>
 #include<algorithm>
 #include<algorithm>
 #include<cstring>
 #define ll long long
 using namespace std;
 const ll p=1e9+;

 ll fact[];

 void init()
 {
     memset(fact,,sizeof(fact));
     fact[]=fact[]=;
     for(ll i=;i<;i++)
         fact[i]=fact[i-]*i%p;

 }

 ll power(ll a,ll b,ll p)
 {
     ll res=;
     while(b)
     {
         if(b%)
             res=res*a%p;
         b=b/;
         a=a*a%p;
     }
     return res%p;
 }

 ll C(ll n, ll m ,ll p)
 {   ///C(n,m) = n! / ( m!*(n-m)! )
     ///数据太大肯定爆，p又是素数。换成求m!*(n-m)!的逆元，又不能一起求，会爆数据，分开求，看做n!/m! * 1/(n-m)!
     ///由于是对p求模，n，m范围在阶乘表范围里
     if(m>n)
         return ;
     return fact[n] * power(fact[m], p-, p)%p * power(fact[n-m], p-, p) % p;
     ///fact * power * power 可能爆long long，第二次就要取模
 }

 int main()
 {
     init();
     ll n,m;
     while(scanf("%lld%lld",&n,&m)!=EOF)
     {
         ll ans=C(n+m-,m-,p);
         printf("%lld\n",ans);
     }
     return ;
 }