补题报告————组合数的末十位

错误题目：组合数的末十位。

思路：我们之前学过，类似C(n,r)这种算式我们有三种做法：

1. 杨辉三角；
2. 直接用公式；

但是我试了杨辉三角后发现根本不行，因为杨辉三角所需的是一个二维数组，但是数据量不允许开那么大，杨辉三角如下：

1

1   1

1   2   1

1   3   3   1

1   4   6   4   1

              1   5   10  10  5   1

 ··· ··· ··· ··· ···

所以只能直接用公式了，但是当时我只想到要用高精度除法，却没想到应该是高精度除高精度，结果交上去时间超限······

这却思路应该是这样：

1. 首先需要利用欧拉筛(就是线性筛)筛一遍；
2. 让然后进入唯一分解定理，进行约分(这里是阶乘形式的唯一分解定理)
3. 遍历素数表(这里面是次幂)；
4. 把数乘上素数表里的次幂；
5. 取余；
6. 加前导零输出（要用printf）。

AC代码：

#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;

}