POJ 2891 扩展欧几里得

初始写了个很裸很暴力的中国剩余定理，并且不会判-1的条件，眼一闭无耻交了一个，TLE，不是WA真给面子 — —|||

看了官方题解，咨询了娃娃师父方知扩展欧几里得方为正解。。。

不过迭代求x的时候不加上那个mod a2/GCD(a1,a2) discuss里一组样例测不过，不解。。。

题意是说，给一组a[i],r[i],求一个最小的m，满足 for_each(i=0: n) a[i]%m=r[i].即 a[i]*x+r[i]=m,不存在输出-1。我们首先得到 a1 * x + r1 = m 和 a2 * y + r2 = m，联立得到a1 * x - a2 * y = r2 - r1,利用ex_gcd解之，先解得a1*X-a2*Y=GCD(a1,a2)时的X Y值和a1,a2的gcd，x y存在整数解的条件是 gcd | r2-r1.是否输出-1就是从这判断。有解即为x=(r2-r1)/gcd*X tmp=a2/gcd, x=( ( r2 - r1 ) / gcd * X % tmp + tmp ) % tmp;接着输入下一组a3,r3,问题变为求满足LCM(a1,a2)%m1=m和a3%m1=r3的新的m值：m1LCM(a1,a2) * x + a3 * y = r3 - m这样只要a1 a2 r1 r2 轮流迭代即可,不用开数组。题目也没说k多大，暗示本题不用开数组。每次更新a1为LCM(a1,a2),r1为m即可。不过注意求x时保持通用性，java高精解决import java.io.*;import java.util.*;import java.math.*;import java.util.Map;import java.util.HashMap;import java.util.Map.Entry;import java.text.*;//SimpleDateFormat;public class Main{public static final BigInteger one=BigInteger.ONE,zero=BigInteger.ZERO;public static final int sup=100010;public static BigInteger x,y;public static BigInteger ex_gcd(BigInteger a,BigInteger b){if(0==b.compareTo(zero)){x=one;y=zero;return a;}BigInteger r=ex_gcd(b,a.mod(b));BigInteger tmp=x;x=y;y=tmp.subtract(a.divide(b).multiply(y));//tmp-a/b*yreturn r;}public static void main(String pzjay[]){Scanner pz=new Scanner(new BufferedInputStream(System.in));int k;BigInteger a1,r1,a2,r2;BigInteger Gcd;boolean ok;BigInteger tmp;while(pz.hasNext()){ok=true;k=pz.nextInt();a1=pz.nextBigInteger();r1=pz.nextBigInteger();for(int i=1;i<k;++i){a2=pz.nextBigInteger();r2=pz.nextBigInteger();if(false==ok)continue;Gcd=ex_gcd(a1,a2);//pz.nextBigInteger());//tmp记录ab的公约数if(0!=r2.subtract(r1).mod(Gcd).compareTo(zero))//Gcd|r2-r1{ok=false;continue;}tmp=a2.divide(Gcd);//x=((r2-r1)/Gcd*x%tmp+tmp)%tmp // 先乘x就不对, discuss里一组数据过不去x=r2.subtract(r1).divide(Gcd).multiply(x).mod(tmp).add(tmp).mod(tmp);r1=a1.multiply(x).add(r1);a1=a1.multiply(a2).divide(Gcd);}if(false == ok)System.out.println("-1");elseSystem.out.println(r1);}}