poj-2191

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, 然后打表过了。 

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

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值