判断素数:
bool isPrime(int a){
if(a<=1)return false;
int sqr=(int)sqrt(1.0*a);
for(int i=2;i<=sqr;i++){
if(a%i==0)return false;
}
return true;
}
//构造素数表
const int MAXN=100010;
int prime[MAXN],pNum=0;//素数表
bool p[MAXN]={false};//p[x]=true,说明x为素数
void Find_Prime(){
for(int i=1;i<MAXN;i++){
if(isPrime(i)){
p[i]=true;
prime[pNum++]=i;
}
}
}
质因数分解:
//质因数
struct factor{
int x;//素数
int cnt;//个数
}fac[10];
//获取质因数
int num=0;
void getFactor(int a){
int sqr=(int)sqrt(1.0*a);
for(int i=0;i<pNum&&prime[i]<=sqr;i++){
if(a%prime[i]==0){
fac[num].x=prime[i];
fac[num].cnt=0;
while(a%prime[i]==0){
fac[num].cnt++;
a/=prime[i];
}
num++;
}
if(a==1)break;
}
if(a!=1){
fac[num].x=a;
fac[num].cnt=1;
num++;
}
}