UVa 106 - Fermat vs. Pythagoras

18 篇文章 0 订阅

对于a2+b2=c2,求出所有满足条件的(a,b,c)三元组,输出a,b,c三个互质的组别的数量和从1到N中所有
的三元组中都没有出现过的数字的个数。

互质毕达哥拉斯三元组,可直接得到:

a = m^2-n^2    b = 2mn    c = m^2+n^2

因为三个变量之间两两互质,所以a,b中必为一奇一偶。(详细证明见:http://www.cnblogs.com/devymex/archive/2010/08/07/1799713.html

因此问题的复杂度就可以降为O(m*n),即:最坏1000*1000。

除互质勾股数外,其他所有勾股数三元组都可以由求出的互质勾股数乘以一定的系数得到,由此可以得到三元组没出现的数字的个数(n - 所有勾股数个数)。

代码如下:

#include <cstdio>
#include <cmath>
#include <cstring>
const int MAXN = 1000001;
bool vis[MAXN];
int gcd(int m, int n)
{
    int k = m%n;
    while(k)
    {
        m = n;
        n = k;
        k = m%n;
    }
    return n;
}
int main()
{
#ifdef test
    freopen("sample.txt", "r", stdin);
#endif
    int n;
    while(scanf("%d", &n) != EOF)
    {
        int a, b, c, sancct = 0, notsanct = n;
        int sqrtn = (int)sqrt((double)n);
        memset(vis, false, sizeof(vis[0])*(n+1));
        for(int i = 1; i <= sqrtn; ++i)
        {
            for(int j = i+1; j <= sqrtn; j+=2)
            {
                if(gcd(i,j) == 1)
                {
                    a = j*j-i*i;
                    b = 2*i*j;
                    c = j*j+i*i;
                    if(c>n)
                        break;
                    ++sancct;
                    int nc = n/c;
                    for(int k = 1; k<=nc; ++k)
                    {
                        vis[k*a]=true;
                        vis[k*b]=true;
                        vis[k*c]=true;
                    }
                }

            }
        }
        for(int i = 1; i <= n; ++i)
            notsanct -= vis[i];
        printf("%d %d\n", sancct, notsanct);
    }
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值