我们知道筛法一般就是开辟一个大的数组,然后循环可以对2,3,5等的倍数剔除,通过昨天的acm集训,我了解了几种变形,很有趣。
1):一般筛法:
bool prime[MAX];
memset(prime,true,sizeof(prime));
for(k=4;k<MAX;k+=2)
prime[k]=false;
for(i=3;i<=SQRTMAX;++i)
{
if(prime[i])
{
for(int j=i*i;j<MAX;j+=2*i)
prime[j]=false;
}
}
j从i*i开始,显然之前比i*i更小的数已经被剔除了,但是此法注意i的大小可能导致i*i越界,这个时候开大范围或者进行判断分为i*i和i+i两种情况先把2倍数剔除,之后可以j+=2*i.
附:HDU1397
#include<iostream>
#include<cmath>
#include<cstring>
using namespace std;
const int MAX=32770;
const int SQRTMAX=181;
const int SUBMAX=10000;
bool prime[MAX];
int subprime[SUBMAX];
void creatprime(void)
{
int i,k;
memset(prime,true,sizeof(prime));
for(k=4;k<MAX;k+=2)
prime[k]=false;
for(i=3;i<=SQRTMAX;++i)
{
if(prime[i])
{
for(int j=i*i;j<MAX;j+=2*i)
prime[j]=false;
}
}
subprime[0]=0;
for(i=2;i<MAX;++i)
if(prime[i])
subprime[++subprime[0]]=i;
}
int main(void)
{
int n,count;
creatprime();
cin>>n;
while(n!=0)
{
count=0;
if(prime[n-2])
count++;
for(int i=1;i<=subprime[0]&&subprime[i]<=n/2;i++)
if(prime[subprime[i]]&&prime[n-subprime[i]])
count++;
cout<<count<<endl;
cin>>n;
}
return 0;
}
2):求解一个数的因子之和:在一般筛的基础上稍变形,prime[i]累加因子即可.
memset(prime,0,sizeof(prime));
for(i=1;i<=MAX/2;++i)
{
for(j=i+i;j<MAX;j+=i)
prime[j]+=i;
}
附:HDU1215
#include<stdio.h>
#include<string.h>
#define MAX 500001
int prime[MAX];
void creatprime(void)
{
int i,j;
memset(prime,0,sizeof(prime));
for(i=1;i<=MAX/2;++i)
{
for(j=i+i;j<MAX;j+=i)
prime[j]+=i;
}
}
int main()
{
int n,ncase;
creatprime();
scanf("%d",&ncase);
while(ncase--)
{
scanf("%d",&n);
printf("%d\n",prime[n]);
}
return 0;
}
3):求解一个数的最大质因子(不超过他本身在质数表的第几个)例如9的最大因子3排名第二个,15的最大因子5排名第三个:
memset(prime,0,sizeof(prime));
for(i=2;i<=MAX;++i)
{
if(prime[i]==0)
{
count++;
for(j=i;j<MAX;j+=i)
prime[j]=count;
}
}
附: HDU2136
#include<stdio.h>
#include<string.h>
#define MAX 1000000
int prime[MAX];
void creatprime(void)
{
int i,j,count=0;
memset(prime,0,sizeof(prime));
for(i=2;i<=MAX;++i)
{
if(prime[i]==0)
{
count++;
for(j=i;j<MAX;j+=i)
prime[j]=count;
}
}
}
int main()
{
int n;
creatprime();
while(scanf("%d",&n)!=EOF)
{
printf("%d\n",prime[n]);
}
return 0;
}
4):求解一个数的最大质因子(如果是质数则是它本身)
memset(prime,0,sizeof(prime));
for(int i=2;i<MAX;++i)
{
if(prime[i]==0)
{
for(int j=i;j<MAX;j+=i)
prime[j]=i;
}
}
昨天的acm集训学到了一些东西,让我们学会要变通思考,不过关键还是多练多思考,继续努力。