军哥讲的数论

因式分解

 

tn=n;

for(i=2;i*i<=n;i++)//试除2~sqrt(n)

     if(tn%i==0){//如果能被i整除

 p[++cnt]=i;//保存底数

 e[cnt]=0;//保存指数

while(tn%i==0){//计算指数

e[cnt]++;

tn/=i;

} }}

if(tn>1){//存在大于sqrt(n)的素因子

 p[++cnt]=tn;

 e[cnt]=1;

}