C++素數篩選法
阿新 • • 發佈:2019-01-31
#include<bits/stdc++.h>
using namespace std;
const int arrSize=1001;
int prime[arrSize]; //儲存素數
int primeSize; //限定範圍內素數個數
bool mark[arrSize]; //用來做標記
void primeFilter(){
primeSize=0;
for(int i=2;i<arrSize;i++){
if(!mark[i]){
prime[primeSize++]=i;
for (int j=i*i;j<arrSize;j+=i) mark[j]=true; //從i*i開始
}
}
/*另一種寫法*/
/*for(int i=2;i<arrSize;i++){
if(mark[i]) continue;
prime[primeSize++]=i;
if(i>1000) continue;
for(int j=i*i;j<arrSize;j+=i) mark[j]=true;
}//*/
}
long long numPrime[1001];
int valPrime[1000];
int valSize[1000];
int main(){
long long num;
long long val,temp1,temp2;
primeFilter();
while(scanf("%lld",&num)!=EOF){
scanf("%lld",&val);
for(int i=0;i<1001;i++){
numPrime[i]=0;
valPrime[i]=0;
valSize[i]=0;
}
for (int i=2;i<=num;i++){
if(!mark[i]) numPrime[i]++;
else{
temp2=i;
for(int j=0;j<primeSize;j++){
while(temp2%prime[j]==0){
numPrime[prime[j]]++;
temp2/=prime[j];
}
if(temp2==1) break;
}
}
}
long long index=0;
for(int i=0;i<primeSize;i++){
if(val%prime[i]==0) {
while(val%prime[i]==0){
valSize[index]++;//儲存每個分解出的素數的個數
val/=prime[i];
}
valPrime[index++]=prime[i];//將分解出的素數值儲存下來
if(val==1) break;
}
}
int min=10000;
for(int i=0;i<index;i++){
temp1=numPrime[valPrime[i]]/valSize[i];
if(temp1<min) min=temp1;
}
printf("%d\n",min);
}
return 0;
}