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

x mod a=r(N对a,r)

时间:2018-07-14 14:46:36      阅读:139      评论:0      收藏:0      [点我收藏+]

标签:span   typedef   ios   ace   lin   const   return   问题:   can   

 1 //模数不一定互质,互质才可以用孙子定理。
 2 /*
 3 https://blog.csdn.net/zmh964685331/article/details/50527894
 4 uu遇到了一个小问题,可是他不想答。你能替他解决这个问题吗?
 5  问题:给你k对a和r是否存在一个正整数x使每队a和r都满足:x mod a=r,求最小正解x或无解。
 6 */
 7 
 8 #include <iostream>
 9 #include <cstdio>
10 #include <queue>
11 #include <algorithm>
12 #include <vector>
13 #include <set>
14 #include <string>
15 #include <cstring>
16 typedef long long ll;
17 using namespace std;
18 const int N=1e5+9;
19 ll a[N],r[N];
20 int n;
21 ll egcd(ll a,ll b,ll &x,ll &y)
22 {
23     ll d=a;
24     if(!b)
25     {
26         x=1,y=0;
27     }
28     else{
29         d=egcd(b,a%b,y,x);
30         y-=x*(a/b);
31     }
32     return d;
33 }
34 ll  solve()
35 { // M:上个式子的模数,下个式子的 a.  R:每一次的解.
36     ll  M=a[1],R=r[1],x,c,y,d;
37     for(int  i=2;i<=n;i++)
38     {
39     d=egcd(M,a[i],x,y);
40     c=R-r[i];
41     if(c%d) return -1;
42     x=c/d*x%(a[i]/d);//a[i]下个式子的b.范围在【0,b-1].所以取余a[i].
43     R-=x*M;//特解X0,新的余数。
44     M=M*a[i]/d;//lcm(a,b);新的模数。
45     R%=M;
46     }
47     return (R%M+M)%M;
48 }
49 int main()
50 {
51     while(~scanf("%d",&n))
52     {
53         for(int i=1;i<=n;i++)
54         {
55             scanf("%lld%lld",&a[i],&r[i]);
56         }
57         printf("%lld\n",solve());
58     }
59     return 0;
60 }

 

x mod a=r(N对a,r)

标签:span   typedef   ios   ace   lin   const   return   问题:   can   

原文地址:https://www.cnblogs.com/tingtin/p/9309329.html

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