Question

根据youtube link 奇数的素因子可以按如下方式计算：

a = sqrt（N + b ^ 2）

我写下面的程序来做到这一点，但我没有得到2345678917的素数因素。我知道这是素数但是对于其他素数，程序确实返回1和数字本身但是对于这个数字它没有发生。为什么呢？

#include <stdio.h>
#include <math.h>

void foo(unsigned long long x)
{
    int i;
    for (i=1;i<x;i++)
        if (fmod(sqrt(x + i*i), 1) == 0) {
            printf("%f %f\n", (sqrt(x + i*i) - i), (sqrt(x + i*i)+i));
            return;
        }
}

int main(void) {
    foo((unsigned long long)2345678917);
    return 0;
}

Answer 1

首先，i（int）的类型与x（unsigned long long）不匹配。将i的类型更改为unsigned long long即可开始使用。

一旦你这样做，你的程序将高兴地宣布14 * 167548494 = 2345678917.当然，这不是真的，因为两个偶数的乘积不能是奇数。这里的问题是精度损失，所以你需要为整数实现一个平方测试函数，而不是测试浮点平方根是否是整数。

#include <stdio.h>
#include <math.h>

unsigned long long find_sqrt(unsigned long long x)
{
    unsigned long long lo = 1;
    while (4 * lo * lo <= x) lo *= 2;
    unsigned long long hi = 2 * lo;
    while (lo + 1 < hi) {
        unsigned long long mid = lo + (hi - lo) / 2;
        if (mid * mid <= x) lo = mid;
        else hi = mid;
    }
    return lo * lo == x ? lo : 0;
}

void foo(unsigned long long x)
{
    unsigned long long i;
    for (i=1;i<x;i++) {
        unsigned long long sqrt_x_ii = find_sqrt(x + i*i);
        if (sqrt_x_ii) {
            printf("%llu = %llu * %llu\n",
                   x, sqrt_x_ii - i, sqrt_x_ii + i);
            return;
        }
    }
}

int main(void) {
    foo((unsigned long long) 2345678917);
    return 0;
}

用于素因子计算的费马算法

1 个答案: