组合数的末十位

组合数的末十位

时间限制：1秒        内存限制：128M

题目描述

求出C(n,r)的最后十位，其中0<r≤n≤30000，输出时不足十位数也按十位输出，此时高位用0表示。C(n,r)=n×(n-1)×……×(n-r+1)/(1×2×3×……×r)。

输入数据为两个以空格隔开的自然数n，r。

输入描述

一行两个整数

输出描述

一行，10位数字。

样例

输入

5 2

输出

0000000010

#include<iostream>
using namespace std;
const int N=3e4+10;
const long long M=1e10;
int n,m,cnt,prime[N],power[N];
bool st[N];
long long ans=1;
void ola(int x){
	for(int i=2;i<=x;i++){
		if(st[i]==0){
			prime[cnt++]=i;
		}
		for(int j=0;j<cnt&&prime[j]*i<=x;j++){
			st[prime[j]*i]=1;
			if(i%prime[j]==0){
				break;
			}
		}
	}
}
void only(int n,int flag){
	for(int i=0;i<cnt;i++){
		long long fac=prime[i],sum=0;
		while(fac<=n){
			sum+=n/fac;
			fac*=prime[i];
		}
		if(flag==1) power[i]+=sum;
		else power[i]-=sum;
	}
}
void mul(int p,int n){
	long long fac=1;
	for(int i=1;i<=n;i++){
		fac=fac*p%M;
	} 
	ans=(ans*fac)%M;
}
int main(){
	cin>>n>>m;
	ola(n);
	only(n,1);
	only(m,-1);
	only(n-m,-1);
	for(int i=0;i<cnt;i++){
		if(power[i]!=0){
			mul(prime[i],power[i]);
		}
	}
	printf("%010lld",ans);
	return 0;
}