BZOJ1407荒岛野人savage——扩展欧几里得

Input
第1行为一个整数N(1<=N<=15)，即野人的数目。
第2行到第N+1每行为三个整数Ci, Pi, Li表示每个野人所住的初始洞穴编号，每年走过的洞穴数及寿命值。
(1<=Ci,Pi<=100, 0<=Li<=10^6 )
Output

仅包含一个数M，即最少可能的山洞数。输入数据保证有解，且M不大于10^6。
Sample Input
3

1 3 4

2 7 3

3 2 1
Sample Output
6

//该样例对应于题目描述中的例子。

---

野人不会在相同的位置，那也就表明不会出现下列情况：（R表示一整圈的长度。）

$(C_i+P_i\ast t)-(C_j-P_j\ast t)\equiv 0（ModR）$

那么移项可得

$(P_i-P_j)\ast t\equiv C_j-C_i（ModR）$

我们可以用扩展欧几里得算出下面式子里的$T$

$(P_i-P_j)\ast T\equiv gcd(P_i-P_j,R)（ModR）$

那么根据性质，$T\ast (C_j-C_i)/gcd(P_i-P_j,R)$就是$t$的一组解。并且$T\ast (C_j-C_i)/gcd(P_i-P_j,R)+x\ast R/gcd(P_i-P_j,R)$就是$t$的解集、

如果这个解集的最小正解大于$min(l_i,l_j)$，说明当圈长为$R$时，这两个人不会发生矛盾。

当然方程无解也不会发生矛盾，根据裴蜀定理，无解的情况就是$(C_j-C_i)\%gcd(P_i-P_j,R)!=0$

至于上面那个解集的证明，详情请看这篇博客

#include<bits/stdc++.h>
using namespace std;
int read(){
    char c;int x;while(c=getchar(),c<'0'||c>'9');x=c-'0';
    while(c=getchar(),c>='0'&&c<='9') x=x*10+c-'0';return x;
}
int n,Max,Min[18][18];
struct node{
    int c,p,t;
}F[18];
void exgcd(int a,int b,int &x,int &y){
    if(!b){
        x=1;y=0;return;
    }
    exgcd(b,a%b,y,x);
    y-=a/b*x;
}
int gcd(int x,int y){
    return y?gcd(y,x%y):x;
}
int check(int x,int y,int Md){
    int xx=F[x].p-F[y].p;
    int a=0,b=0,c=gcd(xx,Md);
    int yy=F[y].c-F[x].c;
    exgcd(xx,Md,a,b);
    if(yy%c) return 0;
    a=a*yy/c;c=abs(c);a=(a%(Md/c)+(Md/c))%(Md/c);
    return (a>Min[x][y])^1;
}
int main()
{
    n=read();int i,j,k;
    for(i=1;i<=n;i++){
        F[i].c=read();F[i].p=read();F[i].t=read();
        Max=max(Max,F[i].c);
    }
    for(i=1;i<=n;i++)
     for(j=i+1;j<=n;j++) Min[i][j]=min(F[i].t,F[j].t);
    for(k=Max;k<=1000000;k++){
        int flag=0;
        for(i=1;i<=n;i++){
         if(flag) break;
         for(j=i+1;j<=n;j++){
          flag|=check(i,j,k);
          if(flag) break;   
          }
        }
        if(!flag){
            printf("%d",k);break;
        }
    }
    return 0;
}