人见人爱gcd的解法及简单证明

题目来源：NEFU oj-1221 人见人爱gcd
详细链接已补。

op

首博，请多指教；
&特别感谢ph姐姐的深夜指导；
&显然成立大法好。

大体想法

题目欲求 x²+y² ，同时给出 x+y=a ，依经验得须由 (x+y)²-2xy 得出 x²+y² ，此时求出 xy 即可。
实际上（ph姐姐觉得），xy=gcd(a,b)*b ，故答案为 a²-2b*gcd(a,b) 以下为试证明过程。

简单证明

*不考虑测试数据不合法的情况
*所有字母均代表整数
已知 a=x+y,b=lcm(x,y) ，求证 xy=gcd(a,b)*b 。
设 i=gcd(x,y);
则有 x=p*i,y=q*i;(p,q互质)
则b=lcm(x,y)=pqi,
a=x+y=(p+q)*i;
∴xy=pqi²=i*b;
往证 xy=gcd(a,b)*b ，
即证 gcd(a,b)=i;
gcd(a,b)=gcd(pqi,(p+q)i);
即证 pq 与 p+q 互质；

假设 pq 与p+q 有公因数 j；
则有 pq=pnj,
p+q=mj;
∴p+q-q=(m-n)*j=p;
∴q=n*j;
则 p 与 q 有公因数 j (j>1)；
与前证 p,q 互质矛盾，故假设不成立；

故 pq 与 p+q 互质；
xy=gcd(a,b)*b 得证。

AC代码

#include<bits/stdc++.h>
#define ll long long
using namespace std;

ll gcd(ll a,ll b)
{
    return b?gcd(b,a%b):a;
}

int main()
{
    ll a,b,n;
    while(~scanf("%lld",&n))//cin被卡了
    {
        while(n--)
        {
            scanf("%lld%lld",&a,&b);
            printf("%lld\n",a*a-2*b*gcd(a,b));
        }
    }
    return 0;
}

ed

第一次接触gcd与lcm的非语法题，感觉还蛮有意思。
//虽然昨天晚上被contest暴打；
//也许它就是语法题。