码迷,mamicode.com
首页 > 其他好文 > 详细

HDU 2669 Romantic (扩展欧几里得定理)

时间:2018-04-11 00:04:20      阅读:182      评论:0      收藏:0      [点我收藏+]

标签:nbsp   bsp   a*   main   anti   style   names   div   mil   

题目大意:给两个数a和b,找出一组x,y使得a*x + b*y = 1,如果找不出输出sorry
 
题解:显然是用扩展欧几里得定理求解。
又扩展欧几里得定理有,如果a*x+b*y = d   要使得方程有解必有gcd(a,b)为d的约数。
而此题的d = 1  所以若gcd(a,b)!=1,则应该输出sorry

 

 

#include <bits/stdc++.h>

using namespace std;


long long e_gcd(long long a,long long b,long long &x,long long &y)
{
    if(b==0)
    {
        x=1;y=0;
        return a;
    }
    long long ans=e_gcd(b,a%b,x,y);
    long long temp=x;
    x=y;
    y=temp-a/b*y;
    return ans;
}

long long cal(long long a,long long b,long long c)
{
    long long x,y;
    long long gcd=e_gcd(a,b,x,y);
    if(c%gcd) return -1;
    long long k=c/gcd;
    x*=k,y*=k;//相当于x=x0+b/gcd*c
    if(b<0) b=-b;
    long long ans=x%b;//最小的x
    if(ans<=0) ans+=b;
    return ans;
}

int main()
{
    long long a,b;
    while(cin>>a>>b)
    {
        long long ans=cal(a,b,1);
        if(ans==-1) printf("sorry\n");
        else printf("%lld %lld\n",ans,(1-ans*a)/b);
    }
    return 0;
}

 

HDU 2669 Romantic (扩展欧几里得定理)

标签:nbsp   bsp   a*   main   anti   style   names   div   mil   

原文地址:https://www.cnblogs.com/Fy1999/p/8783489.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!