无平方因子的数（Eratosthenes筛法）

题目：

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

分析：

思路都是紫书上的，只是实现一下；

代码：

#include
#include
#include
using namespace std;
const int maxn=1000002;
int prim[maxn];
int p[maxn];
int n,m,t;
void solve(int k)
{
    memset(p,0,sizeof(p));
    p[1]=1;
    for(int i=2;i<=sqrt(k+0.5);i++)
    {
        if(!p[i])
            for(int j=i*i;j<=k;j+=i)
                p[j]=1;
    }
    t=0;
    for(int i=2;i<=k;i++)
        if(!p[i])
         prim[t++]=i;
}
int judge(int k)
{
    for(int i=0;i     {
        int z=prim[i]*prim[i];
        if(k%z==0)
            return 0；
    }

return 1；
}
int main()
{
    int ans=0;
    scanf("%d %d",&n,&m);
    int k=sqrt(m+0.5);
    solve(k);
    for(int i=n;i<=m;i++)
        if(judge(i))
            ans++;
    printf("%d\n",ans);
    return 0;
}
没找到oj的地方，也不知道有没有错误，仅供参考，欢迎指正