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

复制代码
 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 }
复制代码

 

posted on   cltt  阅读(246)  评论(0编辑  收藏  举报

编辑推荐:
· go语言实现终端里的倒计时
· 如何编写易于单元测试的代码
· 10年+ .NET Coder 心语,封装的思维:从隐藏、稳定开始理解其本质意义
· .NET Core 中如何实现缓存的预热?
· 从 HTTP 原因短语缺失研究 HTTP/2 和 HTTP/3 的设计差异
阅读排行:
· 分享一个免费、快速、无限量使用的满血 DeepSeek R1 模型,支持深度思考和联网搜索!
· 基于 Docker 搭建 FRP 内网穿透开源项目(很简单哒)
· ollama系列01:轻松3步本地部署deepseek,普通电脑可用
· 25岁的心里话
· 按钮权限的设计及实现
< 2025年3月 >
23 24 25 26 27 28 1
2 3 4 5 6 7 8
9 10 11 12 13 14 15
16 17 18 19 20 21 22
23 24 25 26 27 28 29
30 31 1 2 3 4 5

导航

统计

点击右上角即可分享
微信分享提示