JZYZOJ1372 [noi2002]荒岛野人 扩展欧几里得

http://172.20.6.3/Problem_Show.asp?id=1372

想法其实很好想，但是我扩展欧几里得还是用得不熟练，几乎是硬套模板，大概因为今天一个下午状态都不大好。
扩展欧几里得算法计算的是 : ab互质时ax+by=1或ab不互质时ax+by=gcd(a,b)（废话）的一个整数解，可以据此推导一个方程是否有解。
然后我理解这个基本概念理解了一个下午，非常智障了。
这道题也是模板，两两对比即可。

代码

#include<iostream>  
#include<cstdio>  
#include<cstring>  
#include<algorithm>  
#include<cmath>
using namespace std;
const int maxn=100010;
long long n,p[20]={},l[20]={},c[20]={};
int exgcd(int a,int b,long long &x,long long &y){
    if(!b){
        x=1;y=0;
        return a;
    }
    int d=exgcd(b,a%b,x,y);
    long long w=x;x=y;
    y=w-y*(a/b);
    return d;
}
int main(){
    scanf("%d",&n);long long a,b,m,d,x,y,now=n;
    for(int i=1;i<=n;i++){
        scanf("%I64d%I64d%I64d",&c[i],&p[i],&l[i]);
        now=max(now,c[i]);
    }
    for(int k=now;;k++){
        int f=0;
        for(int i=1;i<=n;i++){
            for(int j=i+1;j<=n;j++){
                a=p[i]-p[j];b=k;m=c[j]-c[i];
                d=exgcd(a,b,x,y);
                if(m%d)continue;
                b/=d;m/=d;x=x*m;x%=b;
                if(x<0)x+=abs(b);
                if(x<=min(l[i],l[j])){f=1;break;}
            }
            if(f)break;
        }
        if(!f){
            printf("%d
",k);
            break;
        }
    }
    return 0;
}