Project-Euler-009素勾股数

最新推荐文章于 2024-07-21 16:04:23 发布

ivan_zcy

最新推荐文章于 2024-07-21 16:04:23 发布

阅读量243

点赞数

本文链接：https://blog.csdn.net/Ivan_zcy/article/details/83031827

版权

009题：

题意：

毕达哥拉斯三元组是三个自然数a < b < c组成的集合，并满足

a2 + b2 = c2

例如，32 + 42 = 9 + 16 = 25 = 52。

有且只有一个毕达哥拉斯三元组满足 a + b + c = 1000。求这个三元组的乘积abc

思路：

暂时想到了二个：

一.两层嵌套for循环分别从1～1000遍历a和b，只要满足a + b < 1000 && a^2 + b ^ 2 == （1000 - a - b）^2即可，时间复杂度大概1e6

二.利用素勾股数公式：满足gcd(n, m) == 1的互质数，存在a = 2 * n * m; b = m^2 - n^2; c = m^2 + n ^2，这是只需要遍历n和m（取值范围都是1～32），只要满足a+b+c == 1000即可，时间复杂度大概是1e3

代码：

#include <stdio.h>
#include <math.h>
#include <string.h>
#include <inttypes.h>
#include <algorithm>
#define ll int64_t
#define d int32_t
#define f double
#define r return
#define mem(a) memset(a, 0, sizeof(a));
using namespace std;

//求gcd
d gcd(d a, d b) {
    r b == 0?a : gcd(b, a % b);
}

void work() {
    for (d i = 1; i <= 32; i++) {
        for (d j = i + 1; j <= 32 && i * i + j * j < 1000; j++) {
            if (gcd(j, i) == 1) {
                d a = 2 * i * j;
                d b = j * j - i * i;
                d c = j * j + i * i;
                d sum = a + b + c;
                if (1000 % sum == 0) {
                    d t = 1000 / sum;
                    printf("%" PRId64"\n", (ll)t * t * t * a * b * c);
                    r;
                }
            }
        }
    }
    r;
}

d main () {
    work();
    r 0;
}

转载请注明出处！！！

如果有写的不对或者不全面的地方可通过主页的联系方式进行指正，谢谢