题目链接:http://codeforces.com/gym/102798/problem/D
题意:给出c,要求找出a+b=c,其中要求abc的素数因子相乘得到结果小于c
解题思路:首先排除c为素数的情况,c如果为素数那么结果一定会大于等于c。
当c为合数的时候,c可以分解为素数相乘的形式,如果素数中有素数幂的出现,那么一定使得rad(abc)<c
判断素数幂的方法:
因为c的大小为1e18,可以先预处理1e6内的所有素数
①然后判断c%(prim[i] * prim[i])是否为0,为0则直接输出yes即可
②如果c%prim[i]=0,那么就先将c除去prim[i]
③处理完所有1e6内的素数之后,如果没找到素数幂的出现,那么对c进行处理。剩下的c如果等于1,那么直接为no。如果剩下的c大于1,那么说明还有大于1e6的素数存在,因为c的大小最大为1e18,所以剩下的c中最多有两个素数相乘,如果想要判断素数幂,那么直接求c的根即可,如果c为一个数的平方,那么直接输出yes,否则no
#define _CRT_SECURE_NO_WARNINGS
#include<iostream>
#include<cstdio>
#include<string>
#include<map>
#include<vector>
#include<stdlib.h>
#include<math.h>
using namespace std;
#define ll long long
const int MAXN = 1e6 + 10;
ll c;
int t;
ll prime[MAXN + 1];
void getPrime() {
memset(prime, 0, sizeof(prime));
for (int i = 2; i <= MAXN; i++) {
if (!prime[i])prime[++prime[0]] = i;
for (int j = 1; j <= prime[0] && prime[j] <= MAXN / i; j++) {
prime[prime[j] * i] = 1;
if (i % prime[j] == 0) break;
}
}
}
int main() {
getPrime();
cin >> t;
while (t--) {
cin >> c;
int sc = 0;
for (int i = 1; i <= prime[0] && prime[i] <= c / 2; i++) {
if (c % (prime[i]*prime[i]) == 0) {
sc = 1;
break;
}
else if (c % prime[i] == 0) {
c /= prime[i];
}
}
if (sc == 1)
cout << "yes" << endl;
else if (c > 1) { //如果c不为1,说明c内还会有大于1e6的素数
ll tmp = sqrt(1.0 * c);
if (tmp * tmp == c)
cout << "yes" << endl;
else
cout << "no" << endl;
}
else
cout << "no" << endl;
}
return 0;
}