#include <cstdio>
#include <set>
#include <cmath>
using namespace std;
typedef unsigned long long ull;
ull maxt = 1 << 16;
int index[100] = { 4, 6, 8, 9, 10, 12, 14, 15, 16, 18, 20, 21, 22, 24, 25, 26,
27, 28, 30, 32, 33, 34, 35, 36, 38, 39, 40, 42, 44, 45, 46,
48, 49, 50, 51, 52, 54, 55, 56, 57, 58, 60, 62, 63, 64
};
set<ull> ans;
int main(int argc, char const *argv[])
{
ans.insert(1);
for (ull i = 2; i < maxt; i++)
{
ull tmp = i * i * i * i;
ull maxi = ceil(64 * (log(2) / log(i))) - 1;
ans.insert(tmp);
for (ull j = 1; index[j] <= maxi; j++)
{
if (index[j] - index[j - 1] == 1) tmp = tmp * i;
else tmp = tmp * i * i;
ans.insert(tmp);
}
}
set<ull>::iterator it = ans.begin();
for (; it != ans.end(); it++) printf("%llu\n", *it);
return 0;
}
当一个数字满足这样的要求的时候。他的指数一定不是素数。因为他的指数可以拆分,所以指数最小是4,所以底数的最大值为2^16-1;
但是这题有个要注意的是,当你判断一个数字是不是越界的时候,可以使用对数来判断,优化的方法。
那么在正式筛选数字前可以通过素数筛选从而来标记出1~64中所有的合数,而对于每一种情况的幂指数就是ceil(log(2^64)/log(i))=ceil(64*(log(2)/log(i)))。