给出正整数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;
}