HDU 5974 A Simple Math Problem 数论 gcd

    今天做山科的题，结果这套题是大连区域赛的题...又想起了当年的黑历史，好吧。这道题当时打现场赛的时候倒是也A了，不过当时我是先跑的暴力TLE了一遍，然后队友剪了剪枝过的。现在想想看来数据出的还是不够卡时间，当年的债终究还是要换的。现在过去半年，我也知道有公式了，虽然也是推了好长时间才推出来= =

    好吧，下面说正题。题意就是说已知x+y=a，lcm(x,y)=b。给出你a和b的值，问能不能求出x和y的值，x和y都是整数。设gcd(x,y)=g，设存在m，n，使mg=x,ng=y。则b=lcm(x,y)=mng，a=mg+ng。因为x和y的公约数为g，则m和n必定互质，所以gcd(a,b)=g=gcd(x,y)。则联立方程x+y=a和xy=b*gcd(x,y)，可以得到方程x^2-ax+b*gcd(a,b)=0。之后要做的就是求出方程的解就可以了，二次方程利用求根公式什么的高中知识也不必多说了。另外注意一点x,y是整数，所以判断一下，如果得出的结果为小数的话，也是要输出No Solution的。（表示一开始没注意这个问题，直接解方程然后WA了...）。

    下面AC代码：

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;

long long gcd(long long m,long long n)
{
    while(m>0)
    {
        long long c = n % m;
        n = m;
        m = c;
    }
    return n;
}

long long x;

int main()
{
    long long p;
    long long a,b,g;
    double d,dd;
    while(scanf("%lld%lld",&a,&b)!=EOF)
    {
        //cout<<sqrt(a)<<endl;
        g=b*gcd(a,b);
        p=a*a-4*1*g;
        //cout<<sqrt(p)<<" and "<<a-sqrt(p)<<endl;
        d=sqrt(p);
        dd=a-sqrt(p);
        if(d-(int)(d)>1e-9||dd-(int)(dd)>1e-9)
        {
            cout<<"No Solution"<<endl;
            continue;
        }
        if(p>=0)
        {
            x=(a-sqrt(p))/2;
            cout<<x<<" ";
            x=a-x;
            cout<<x<<endl;
        }
        else
            cout<<"No Solution"<<endl;
    }
    return 0;
}