筛法的几种变形

我们知道筛法一般就是开辟一个大的数组，然后循环可以对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集训学到了一些东西，让我们学会要变通思考，不过关键还是多练多思考，继续努力。