【Codeforces Round #368 (Div. 2)】Codeforces 707C Pythagorean Triples

Katya studies in a fifth grade. Recently her class studied right
triangles and the Pythagorean theorem. It appeared, that there are
triples of positive integers such that you can construct a right
triangle with segments of lengths corresponding to triple. Such
triples are called Pythagorean triples.

For example, triples (3, 4, 5), (5, 12, 13) and (6, 8, 10) are
Pythagorean triples.

Here Katya wondered if she can specify the length of some side of
right triangle and find any Pythagorean triple corresponding to such
length? Note that the side which length is specified can be a cathetus
as well as hypotenuse.

Katya had no problems with completing this task. Will you do the same?
Input

The only line of the input contains single integer n (1 ≤ n ≤ 109) —
the length of some side of a right triangle. Output

Print two integers m and k (1 ≤ m, k ≤ 1018), such that n, m and k
form a Pythagorean triple, in the only line.

In case if there is no any Pythagorean triple containing integer n,
print  - 1 in the only line. If there are many answers, print any of
them.

构造解。貌似方法很多。
我的方法很麻烦,是用a=pq,b=(p^2-q^2)/2,c=(p^2+q^2)/2推导的。
一种比较简单的方法是对于一个奇数,直接取n^2/2和n^2/2+1,对于一个偶数,将它不断除以二得到奇数再扩大回原来的倍数。

#include<cstdio>
#include<cstring>
#include<cmath>
int main()
{
    long long i,j,k,m,n,p,q,x,y,z;
    scanf("%I64d",&x);
    for (p=1;p*p<=x;p++)
      if (x%p==0)
      {
        q=x/p;
        i=1;
        j=2;
        if ((q*q+p*p)%2==0&&p!=q)
        {
            printf("%I64d %I64d\n",(q*q-p*p)/2,(q*q+p*p)/2);
            return 0;
        }
      }
    for (y=2;y*y<=2*x;y+=2)
      if (x%y==0)
      {
        p=x/y-y/2;
        q=x/y+y/2;
        if (!(p&&q)) continue;
        printf("%I64d %I64d\n",p*q,(p*p+q*q)/2);
        return 0;
      }
    for (p=1;p*p<=2*x;p++)
    {
        q=sqrt(2*x-p*p)+0.1;
        if (q*q==2*x-p&&p!=q&&p&&q)
        {
            printf("%I64d %I64d\n",p*q,(q*q-p*p)/2);
            return 0;
        }
    }
    printf("-1\n");
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值