唯一分解定理:任何一个大于1的自然数n都可以用某些素数的次方的乘积来表示
一,我们要得到素数表,使用欧拉筛
ll num[maxn];
ll prime[maxn];
void getprime(ll n)
{
ll i,n,j,cou=0;
memset(num,0,sizeof num);//全是素数
for(i=2;i*i<=n;i++)
{
if(!num[i])prime[cou++]=i;
for(j=0;j<cou;j++)//用现有的素数去筛
{
num[prime[j]*i]=1;
if(i%prime[j]==0)break;
}
}
}
二,对n进行分解
void getfactor(int n)
{
cnt=0;;
for (i=1;prime[i]*prime[i]<=n;i++)
{
if(n%prime[i]==0)
{
cnt++;
p[cnt]=prime[i];
a[cnt]=0;//统计上面这个素数被乘的次数
while(n%prime[i]0)
{
a[cnt]++;
n/=prime[i];
}
}
if(n1)break;//除完了
}
if(n>1)
{
cnt++;
p[cnt]=n;
a[cnt]=1;
}
}
三,对n!进行分解
n!=123····n
所以对n!分解的时候的质因子肯定包含1-n的所有质因子,只不过每个质因子被乘的次数不清楚
eg. 8!=12345678
其中 2的倍数2 4 6 8
4的倍数4 8
8的倍数8
所以质因数2被乘的次数即指数部分为4+2+1=7
所以cnt(2)= 8/2+8/2/2+8/2/2/2=4+2+1=7
ll getcnt(ll n,ll p)
{
if(n==0)return 0;
return getcnt(n/p,p)+n/p;//返回本次计算的加上递归下个计算的
}