uva11752(合数)

最新推荐文章于 2020-03-03 23:47:51 发布

二分查找

最新推荐文章于 2020-03-03 23:47:51 发布

阅读量1k

点赞数

CC 4.0 BY-SA版权

分类专栏：数论文章标签： uva

本文链接：https://blog.csdn.net/yeyeyeguoguo/article/details/47426247

数论专栏收录该内容

62 篇文章

订阅专栏

题意：

如果一个数可以表示成两个或以上的数的幂，就是超级幂；

例如64 = 8^2 = 4^3;

求1到2^64-1所有的的超级幂；题目没有输入；

思路：

所有数的合数（不是素数就是合数）次幂，都是超级幂；

所以求所有每个数的合数次幂；

要注意溢出的边界，和判重；

#include <cstdio>
#include <cstring>
#include <cmath>
#include <set>
using namespace std;
#define MAX 65
typedef unsigned long long ll;

bool prime[MAX];
set<ll> myset;

ll power(int a, int n) {
    if(n == 0) return 1LL;
    ll ans = power(a, n/2);
    ans *= ans;
    if(n%2 == 1) ans *= a;
    return ans;
}

void init() {
    memset(prime, false, sizeof(prime));
    for(int i = 2; i < MAX; ++i) 
		if(!prime[i]) {
        for(int j = 2*i; j < MAX; j += i)
            prime[j] = true;
    }
}

int main(){
    init();
    myset.clear();
    myset.insert(1);
   
    double lnmax = log( pow(2.0, 64.0)-1 );
    for(int i = 2; i < (1<<16); ++i) {
        int maxe = (int)ceil(lnmax/log(i));
        for(int j = 4; j < maxe; ++j) 
			if( prime[j] ) {
            ll temp = power(i, j);
            if(!myset.count(temp))
                myset.insert(temp);
        }
    }
   
    for(set<ll>::iterator it = myset.begin(); it != myset.end(); ++it)
		printf("%llu\n", *it); 
    return 0;
}