1. 程式人生 > >x mod a=r(N對a,r)

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

解決 for blog printf () style vector 問題: gcd

 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)