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
//该样例对应于题目描述中的例子。
解法:
分析两个野人的情况,设山洞数目为m,如果i和j在第x年相遇,那么可以得到一个同余方程ci+x*pi=cj+x*pj(mod m),
即(pi-pj)*x=cj-ci(mod m)。所以要使两个野人不能相遇,只需该同余方程无解或者x比其中一个野人的寿命长。推广到
n个野人两两不相遇,只需任意两个野人满足上述条件,因为n很小,所以可以枚举。
我们一开始将m赋值为max{c[i]},每次判断m是否可行,如果可行即为答案,否则就将m加1。在判断m是否可行时要用
到扩展欧几里得算法求。
注意:扩展欧几里得算法求满足ax=gcd(a,b) (mod b)的x的最小正整数解时,要将x对b/gcd(a,b)取模,而不是对b取
模。
形如ax=b (mod c)可以转化为ax+cy=b (mod c),然后a,b,c同时mod gcd(a,c)再求。
//BZOJ 1407
#include <bits/stdc++.h>
using namespace std;
///Ci+x*Pi同余(Cj+x*Pj)%m
const int limitn = 20;
int C[limitn], P[limitn], L[limitn];
int limit, n, a, b, c, x, y;
inline int gcd(int a, int b){
if(!b) return a;
else return gcd(b,a%b);
}
inline void exgcd(int a, int b, int &x, int &y){
if(!b){x=1;y=0;}
else exgcd(b,a%b,y,x), y-=x*(a/b);
}
bool solve(int m)
{
for(int i=1; i<n; i++){
for(int j=i+1; j<=n; j++){
a = P[i]-P[j];
b = m;
c = C[j]-C[i];
int t=gcd(a, b);
if(c%t) continue;
a/=t, b/=t, c/=t;
exgcd(a, b, x, y);
b=abs(b);
x=((c*x)%b+b)%b;
if(!x) x+=b;
if(x>L[i]||x>L[j]) continue;
else return false;
}
}
return true;
}
int main()
{
scanf("%d", &n);
for(int i=1; i<=n; i++){
scanf("%d%d%d", &C[i],&P[i],&L[i]);
if(C[i]>limit) limit = C[i];
}
for(int i=limit; ; i++){
if(solve(i)){
printf("%d\n", i);
return 0;
}
}
return 0;
}