n<=10^18是否能被一个数的平方整除

问n<=10^18是否能被一个数的平方整除
    将n分解为：p1p2pk  其中pi为素数，且pi<=pi+1（这里可以等于）
    若n能被一个数的平方整除，它肯定能被一个素因子的平方p^2整除，我们找到最小的这个p即可  ~~~
    1.如果p不是最后的素因子，在p之后的素数pi>=p，这时就有n>=p^3了，
       即p<=n^(1/3)，所以枚举n^(1/3)内的素数看是否p^2|n即可，而n^(1/3)<=10^6
    2.如果p是最后的素因子，所以n的形式就是p1p2p'p^2，而p'<=p，所以P'^3<=n，p'<=n^(1/3)
       所以可以先将n的这些素因子p1,p2,,p'给去掉，p'<=n^(1/3)，所以用n^(1/3)内的素数去除即可
       最后只剩下n'=p^2，通过判断一个数是否为完全平方数即可


    算法就是预处理10^6内的素数，用这些素数去除n，发现有p^2能整除的return即可

    否则判断剩下的这个n'是否为完全平方数，我用直接开方比较的。。。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<iostream>
#define ll long long
using namespace std;
bool isp[1000010];
int pr[300000], pn;


void init()
{
    pn=0;
    for (int p = 2; p <= 1000000; p++) {
        if (!isp[p]) {
            pr[pn++] = p;
            for (long long j = (long long)p*p; j <= 1000000; j+= p) {
                isp[j] = 1;
            }
        }
    }
}
bool work(ll n)
{
    ll i,j;
    ll a=n;
    for(i=0;i<pn&&pr[i]*pr[i]*pr[i]<=n;i++)
    {
        int coun=0;
        while(n%pr[i]==0)
        {
            n/=pr[i];
            coun++;
            //printf("pr=%d count=%d\n",pr[i],coun);
            if(coun>=2)
            return 0;
        }
    }
    i=(ll)sqrt(n);
   // printf("i=%d\n",i);
    if(n==i*i)
    return 0;
    return 1;
}
int main()
{
   int t;
   ll n,m,cas=1;
   init();
   scanf("%d",&t);
   while(t--)
   {
       scanf("%I64d",&n);
       if(work(n))
       printf("Case %d: Yes\n",cas++);
       else
       printf("Case %d: No\n",cas++);
   }
}