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

hdu 1576 A/B

时间:2014-11-19 17:56:09      阅读:157      评论:0      收藏:0      [点我收藏+]

标签:style   blog   http   io   ar   color   os   sp   on   

思路:a/b=x; a=b*x; b*x-9973*y=n;  用扩展欧几里得;

对于等式ax+by=c,abc皆为整数且c是gcd(a, b)的倍数,且(x1, y1)是方程ax+by=gcd(a, b)一组整数解,则(x1*(c/gcd(a, b)), y1*(c/gcd(a, b)))是方程ax+by=c的一组解。

bubuko.com,布布扣
 1 #include <cstdio>
 2 #include <cstring>
 3 #include <algorithm>
 4 #define LL __int64
 5 using namespace std;
 6 
 7 int t;
 8 int n,b;
 9 
10 void extend(LL a,LL b,LL &d,LL &x,LL &y)
11 {
12     if(!b)
13     {
14         x=1;
15         y=0;
16         d=a;
17     }
18     else
19     {
20         extend(b,a%b,d,y,x);
21         y-=x*(a/b);
22     }
23 }
24 
25 int main()
26 {
27     scanf("%d",&t);
28     while(t--)
29     {
30         scanf("%d%d",&n,&b);
31         LL d,x,y;
32         extend((LL)b,9973,d,x,y);
33         x*=n;
34         x=(x%9973+9973)%9973;
35         printf("%I64d\n",x);
36     }
37     return 0;
38 }
View Code

 

hdu 1576 A/B

标签:style   blog   http   io   ar   color   os   sp   on   

原文地址:http://www.cnblogs.com/fanminghui/p/4108654.html

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