瞬间移动
Time Limit: 4000/2000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)Total Submission(s): 490 Accepted Submission(s): 275
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
#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <vector>
using namespace std;
typedef long long ll;
const int Mod = 1e9+7;
/*
方案三:质因数分解,C(n,m)=n!/(m!*(n-m)!)
*/
//利用筛法生成素数
const int MAXN = 1e6;
int prime[MAXN+100]={0};
vector<int> P(){
vector<int> pr;
int i,j;
pr.push_back(2);
for(i=3;i*i<=MAXN;i=i+2){
if(!prime[i]){
pr.push_back(i);
for(j=i*i; j<=MAXN; j+=i)
prime[j]=1;
}
}
while(i<=MAXN){
if(!prime[i])
pr.push_back(i);
i+=2;
}
return pr;
}
//计算n!中的素因子p的指数
int Cal(int x,int p){
int ans=0;
ll rec =p;
while(x>=rec){
ans+=x/rec;
rec*=p;
}
return ans;
}
//计算n的k次方对M取模,二分法
int Pow(ll n, int k, int Mod){
ll ans=1;
while(k){
if(k&1)
ans=(ans*n)%Mod;
n = (n*n)%Mod;
k>>=1; //k=k>>1 k=k/2;
}
return ans;
}
//计算C(n,m)
ll C3(int n,int m){
vector<int> pr = P();
ll ans=1;
int num;
for(int i=0;i<pr.size()&&pr[i]<=n;i++){
num = Cal(n,pr[i])-Cal(m,pr[i])-Cal(n-m,pr[i]);
ans = (ans*Pow(pr[i],num,Mod))%Mod;
}
return ans;
}
int main(){
int n,m;
while(~scanf("%d %d",&n,&m)){
printf("%lld\n",C3(n+m-4,m-2));
}
return 0;
}