poj-2191

最新推荐文章于 2017-04-27 19:59:59 发布

fyfcauc

最新推荐文章于 2017-04-27 19:59:59 发布

阅读量465

点赞数

分类专栏： POJ 数论文章标签： poj

本文链接：https://blog.csdn.net/fyfcauc/article/details/39298417

版权

POJ 同时被 2 个专栏收录

233 篇文章 0 订阅

订阅专栏

数论

26 篇文章 0 订阅

订阅专栏

//404K  0MS G++
#include <cstdio>
#include <cstring>
#include <cmath>

using namespace std;

// bool isPrime(long long N) {
//     for (long long i = 2; i <= N/2+1; i++) {
//         if (N%i == 0) {
//             return 0;
//         }
//     }
//     return 1;
// }

// void preprocess() {
//     long long val = 4;
//     for (int i = 2; i <= 63; i++) {
//         if (isPrime(i)) {        
//             if (!isPrime(val - 1)) {
//                 long long curVal = val - 1;
//                 printf("(2^%d) - 1\n", i);
//                 long long leftPart = curVal;
//                 long long max = (long long)sqrt(leftPart);
//                 for (long long factor = 2; factor <= max; factor++) {
//                         <span style="color:#ff0000;">while (leftPart % factor == 0)
//                         {
//                             printf("%lld * ", factor);
//                             leftPart /= factor;
//                         }</span>
//                 }

//                 printf("%lld = %lld\n", leftPart, curVal);
//             }
//         }
//         val <<=1;
//     }
// }

// int main(int argc, char const *argv[])
// {
//     preprocess();
//     return 0;
// }

int main()
{
int f[10] =
    { 11, 23, 29, 37, 41, 43, 47, 53, 59 };
    char st[10][100] =
    { "23 * 89 = 2047 = ( 2 ^ 11 ) - 1",
            "47 * 178481 = 8388607 = ( 2 ^ 23 ) - 1",
            "233 * 1103 * 2089 = 536870911 = ( 2 ^ 29 ) - 1",
            "223 * 616318177 = 137438953471 = ( 2 ^ 37 ) - 1",
            "13367 * 164511353 = 2199023255551 = ( 2 ^ 41 ) - 1",
            "431 * 9719 * 2099863 = 8796093022207 = ( 2 ^ 43 ) - 1",
            "2351 * 4513 * 13264529 = 140737488355327 = ( 2 ^ 47 ) - 1",
            "6361 * 69431 * 20394401 = 9007199254740991 = ( 2 ^ 53 ) - 1",
            "179951 * 3203431780337 = 576460752303423487 = ( 2 ^ 59 ) - 1" };
    int n;
    scanf("%d", &n);
    for (int i =0; i <9; i++)
        if (f[i] <= n)
            printf("%s\n", st[i]);

    return 0;
}

梅森素数，直接现离线求出所有满足条件的2k次方-1，然后打表过了。

有在线算法，不过太深了，不看了.