HOJ 3826 Squarefree number（素数线性筛法）

素数线性筛法
题目意思：
给出一个数 n <= 10^18, 如果这个数的某个约数是 完全平方数，输出 “No”，否则输出 “Yes”
本题要点：
1、n 很多，用 long long 。然后开个 MaxN = 10^6 的数组， 在 MaxN 的范围内，收集所有的素数到
数组 prime, m 记录素数的个数。 这个过程，使用的是线性筛法，时间复杂度 O(n)
2、然后 用 n 去除 MaxN = 10^6 以内的素数， 除完后，如果 n > 1, 说明 n 含有 大于 MaxN = 10^6 的素因子。
此时， 最多含有两个这样的大素因子（假设为 a, b 把）， 如果 a == b ,那么容易判断 sqrt(n) * sqrt(n) == n,
否则，a != b.

#include <cstdio>
#include <cstring>
#include <iostream>
#include <cmath>
using namespace std;
const int MaxN = 1e6 + 10;
typedef long long LL;
int T;
LL n;
LL prime[MaxN + 1];
int m;	// MaxN 范围内，素数个数
LL vis[MaxN];	// vis[k] 存放 k的最小素因子

void init()	//线性筛法
{
	for(int i = 2; i <= MaxN; ++i)
	{
		if(vis[i] == 0)
		{
			vis[i] = i, prime[++m] = i;
		}
		for(int j = 1; j <= m; ++j)
		{
			if(prime[j] > vis[i] || prime[j] * i > MaxN)
				break;
			vis[i * prime[j]] = prime[j];
		}
	}
}

void solve(int t)
{
	printf("Case %d: ", t);
	LL tmp = n;
	for(int i = 1; i <= m; ++i)
	{
		
		if(tmp % prime[i] == 0)
		{
			int cnt = 0;
			while(tmp % prime[i] == 0)
			{
				tmp /= prime[i];
				++cnt;
			}
			if(cnt >= 2)
			{
				printf("No\n");
				return;
			}
		}
	}
	if(tmp > 1)
	{
		LL k = sqrt(tmp);
		if(k * k == tmp)
		{
			printf("No\n");
			return;
		}
	}
	printf("Yes\n");
}

int main()
{
	init();
	scanf("%d", &T);
	for(int t = 1; t <= T; ++t)
	{
		scanf("%lld", &n);
		solve(t);
	}
	return 0;
}

/*
2
30
75
*/

/*
Case 1: Yes
Case 2: No
*/