poj 1305 Fermat vs. Pythagoras(毕达哥斯拉三元组)

最新推荐文章于 2019-03-20 20:02:12 发布

原创最新推荐文章于 2019-03-20 20:02:12 发布 · 513 阅读

0 ·

CC 4.0 BY-SA版权

----数论---- 专栏收录该内容

41 篇文章

订阅专栏

本文探讨了毕达哥斯拉三元组，即满足a^2+b^2=c^2的整数数组(a, b, c)。介绍了求解此类三元组的方法，并给出了问题POJ 1305的解决思路，重点在于确定s的取值范围3≤s≤√(2n-1)，以及如何找到最小的t值。建议参考相关资料深入理解。" 115761189,7674386,使用Spring MVC处理FormData上传文件,"['后端开发', 'Spring框架', 'HTTP请求', '文件处理']

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

毕达哥斯拉三元组，即勾股数组。
勾股数组参考资料：http://blog.youkuaiyun.com/magicnumber/article/details/6410043
s的取值范围：3≤s≤√(2n-1)，假设t取最小值，s就是这个范围。具体请查看参考资料。

#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <math.h>
using namespace std;

const int MAXN = 1000100;
bool vis[MAXN];

int gcd(int a, int b)
{
    return b==0?a:gcd(b,a%b);
}

int main()
{
    int n,a,b,c,res,cnt;
    while(~scanf("%d",&n))
    {
        res = cnt = 0;
        memset(vis,0,sizeof(vis));
        int m = sqrt(2*n);
        for(int  t = 1; t <= m; t += 2)
        {
            for(int s = t+2; s <= m; s += 2)
            {
                if(gcd(s,t) != 1) continue;
                a = s*t;
                b = (s*s-t*t)/2;
                c = (s*s+t*t)/2;
                if(c > n) break;
                if(a*a + b*b != c*c)
                    continue;
                ++res;
                for(int k = 1; k*c <= n; ++k)
                    vis[k*a] = vis[k*b] = vis[k*c] = true;
            }
        }
        for(int i = 1; i <= n; ++i)
            if(!vis[i])
                ++cnt;
        printf("%d %d\n",res,cnt);
    }
    return 0;
}