蓝桥杯——无平方因子数（埃氏筛法）

给出正整数n和m，区间[n,m]内的“无平方因子”的数有多少个？整数p无平方因子当且仅当不存在k > 1,使得p是k * k的倍数。1 <= n <= m <= 10 ^ 12,m - n <= 10 ^ 7

 

#include<iostream>
#include<cstring>

using namespace std;

typedef long long ll;
const ll maxn=100000000;
ll isprime[maxn];
ll vis[maxn];
int first,end;
int count,m;
//埃氏筛法，筛出素数存入数组vis 
void serch(int t)
{
	memset(isprime,0,sizeof(isprime));
	for(int i=2;i*i<=t;i++)
	{
		if(!isprime[i]) 
		{
		vis[m++]=i;
		for(int j=i*i;j<=t;j=j+i)
			isprime[j]=1;
		}	
	}
	return;		
}

int main()
{
	cin>>first>>end;
	serch(end);
	//筛完之后，记得把isprime数组归零 
	memset(isprime,0,sizeof(isprime));
	//然后将素数的平方数值1 
	for(int i=0;i<m;i++)
	{
		for(int j=vis[i]*vis[i];j<=end;j+=vis[i]*vis[i])
			isprime[j]=1;
	}
	//输出[first,end]这个区间范围内的所有无平方因子的数，并统计个数 
	for(int i=first;i<=end;i++)
	{
		if(!isprime[i])
		{
			cout<<i<<" ";
			count++;
		}
	}
	cout<<endl;
	cout<<count;
	return 0;
 }