最初是在PKU的acm上刷题,刚接触1061青蛙的约会这道题的时候写了个简单暴力的方法,在一定范围内感觉挺对的,但是提交后WA了,于是在一起K题的战友那里第一次听说了这个算法——欧几里得扩展算法。费了很大的劲之后,终于有些了解这个算法了,在这记录下。
所谓扩展欧几里得算法就是给定一组a,b,求aX+bY=gcd(a,b)的解(gcd表示最大公约数,根据数论中的知识,这个方程一定有解)。下面是求解的代码:
public static int x,y;
public static int extendgcd(int a,int b){
if(b==0){
x=1;
y=0;
return a;
}
int d=extendgcd(b,a%b);
int temp=x;
x=y;
y=temp-a/b*x;
return d;
}
为什么要这样求呢,推算过程如下:aX+bY=gcd(a,b)=gcd(b,a%b)=bX'+(a%b)Y'=bX'+(a-a/b*b)Y'=aY'+b(X'-a/bY')(/为整除),所以,X=Y‘,Y=X'-a/bY'。
经过上面的步骤,求出了aX+bY=gcd(a,b)(0)的一组解(x,y),但在1061这个题中,要求解的是aX+bY=c(1)的最小整数解,到这只是迈出了第一步。
将式子1两边同时除以d=gcd(a,b),如果d不能整除c,说明式1没有整数解,除完后得a'X+b'Y=c'(2),将式子0两边同时除以d,得a'X+b'Y=1,显然(2)式的解就是(c'x,c'y),所有整数解的集合为(c'x+b'k,c'y-a'k)(其中k是任意整数)。因为要求最小整数解,所以还要再处理一下,以c'x+b'k为例,它的最小整数值为:(c'x%b'+b')%b',第一步是为了让它成为绝对值小于b'的整数,然后因为可能会变为负数所以加b',又因为可能原本是正数,加b'后大于b',所以又对b'取模。完整的题目1061的AC代码如下:
import java.util.Scanner;
public class Main {
public static long t,p;
public static void main(String[] args) {
Scanner in=new Scanner(System.in);
long x=in.nextLong();
long y=in.nextLong();
long m=in.nextLong();
long n=in.nextLong();
long L=in.nextLong();
long a=m-n,b=L,c=y-x;
if(a<0){
a=-a;
c=-c;
}
long d=extendgcd(a,b);
if((m==n&&x!=y)||c%d!=0)
System.out.println("Impossible");
else{
b=b/d;
c=c/d;
t=c*t;
System.out.println((t%b+b)%b);
}
}
public static long extendgcd(long a,long b){
if(b==0){
t=1;
p=0;
return a;
}
long d=extendgcd(b,a%b);
long temp=t;
t=p;
p=temp-a/b*t;
return d;
}
}