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 ∗ t ) − ( C j − P j ∗ t ) ≡ 0 ( M o d R ) (C_i+P_i*t)-(C_j-P_j*t)≡0(Mod \ \ R) (Ci+Pi∗t)−(Cj−Pj∗t)≡0(Mod R)
那么移项可得
( P i − P j ) ∗ t ≡ C j − C i ( M o d R ) (P_i-P_j)*t≡C_j-C_i(Mod \ \ R) (Pi−Pj)∗t≡Cj−Ci(Mod R)
我们可以用扩展欧几里得算出下面式子里的 T T T
( P i − P j ) ∗ T ≡ g c d ( P i − P j , R ) ( M o d R ) (P_i-P_j)*T≡gcd(P_i-P_j,R)(Mod \ \ R) (Pi−Pj)∗T≡gcd(Pi−Pj,R)(Mod R)
那么根据性质, T ∗ ( C j − C i ) / g c d ( P i − P j , R ) T*(C_j-C_i)/gcd(P_i-P_j,R) T∗(Cj−Ci)/gcd(Pi−Pj,R)就是 t t t的一组解。并且 T ∗ ( C j − C i ) / g c d ( P i − P j , R ) + x ∗ R / g c d ( P i − P j , R ) T*(C_j-C_i)/gcd(P_i-P_j,R)+x*R/gcd(P_i-P_j,R) T∗(Cj−Ci)/gcd(Pi−Pj,R)+x∗R/gcd(Pi−Pj,R)就是 t t t的解集、
如果这个解集的最小正解大于 m i n ( l i , l j ) min(l_i,l_j) min(li,lj),说明当圈长为 R R R时,这两个人不会发生矛盾。
当然方程无解也不会发生矛盾,根据裴蜀定理,无解的情况就是 ( C j − C i ) % g c d ( P i − P j , R ) ! = 0 (C_j-C_i)\%gcd(P_i-P_j,R)!=0 (Cj−Ci)%gcd(Pi−Pj,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;
}