瞬间移动
Time Limit: 4000/2000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)Total Submission(s): 502 Accepted Submission(s): 281
Problem Description
有一个无限大的矩形,初始时你在左上角(即第一行第一列),每次你都可以选择一个右下方格子,并瞬移过去(如从下图中的红色格子能直接瞬移到蓝色格子),求到第
n
行第
m
列的格子有几种方案,答案对
1000000007
取模。
![http://acm.hdu.edu.cn/data/images/C702-1003-1.jpg](http://acm.hdu.edu.cn/data/images/C702-1003-1.jpg)
![http://acm.hdu.edu.cn/data/images/C702-1003-1.jpg](http://acm.hdu.edu.cn/data/images/C702-1003-1.jpg)
Input
多组测试数据。
两个整数 n,m(2≤n,m≤100000)
两个整数 n,m(2≤n,m≤100000)
Output
一个整数表示答案
Sample Input
4 5
Sample Output
10
无非就是个规律。C(m+n-4,m-2);套个模板即可。
#include <iostream>
#include<stdio.h>
using namespace std;
typedef long long ll;
void extendgcd(ll a,ll b,ll &d,ll &x,ll &y)
{
if(b==0){d=a;x=1;y=0;return;}
extendgcd(b,a%b,d,y,x);
y-=x*(a/b);
}
ll GetNi(ll A,ll M)
{
ll rex=0,rey=0;
ll td=0;
extendgcd(A,M,td,rex,rey);
return (rex%M+M)%M;
}
ll C(ll n,ll m,ll p)
{
if(m>n) return 0;
ll up=1,dn=1;
for(int i=0;i<m;i++)
{
up = up*(n-i)%p;
dn = dn*(i+1)%p;
}
return up*GetNi(dn, p)%p;
}
ll lucas(ll n,ll m,ll p)
{
if(m==0) return 1;
return C(n%p,m%p,p)*lucas(n/p,m/p,p) % p;
}
int main()
{
int n,m;
while(scanf("%d%d",&n,&m)!=EOF)
{
ll sum;
ll t1=m+n-4;
ll t2=m-2;
sum=lucas(t1,t2,1000000007);
printf("%d\n",sum);
}
return 0;
}