Description
克里特岛以野人群居而著称。岛上有排列成环行的M个山洞。这些山洞顺时针编号为1,2,…,M。岛上住着N个野人,一开始依次住在山洞C1,C2,…,CN中,以后每年,第i个野人会沿顺时针向前走Pi个洞住下来。每个野人i有一个寿命值Li,即生存的年数。下面四幅图描述了一个有6个山洞,住有三个野人的岛上前四年的情况。三个野人初始的洞穴编号依次为1,2,3;每年要走过的洞穴数依次为3,7,2;寿命值依次为4,3,1。
奇怪的是,虽然野人有很多,但没有任何两个野人在有生之年处在同一个山洞中,使得小岛一直保持和平与宁静,这让科学家们很是惊奇。他们想知道,至少有多少个山洞,才能维持岛上的和平呢?
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
Data Constraint
Hint
该样例对应于题目描述中的例子。
本题是拓展欧几里得算法
枚举m,假设i,j两个野人会在第x年在同一个山洞(摔♂跤),就会满足:Ci+Pi*x≡Cj+Pj*x(%m) →Ci+Pi*x-Cj-Pj*x=my →(Pi-Pj)x-my=Cj-Ci
设(Pi-Pj)=a(野人i每次跑的距离可以比野人j的少),-m=b(可能为负数,要取abs),Cj-Ci=c,就得到ax+by=c,只要是在0~Min(Li,Lj)之间有整数解,就会在同一个山洞不可描述。
然后用拓展欧几里得算法来解方程了,之后求出x是否
#include<iostream>
#include<cmath>
#include<cstdio>
#include<cstring>
#define N 16
using namespace std;
int c[N],p[N],l[N];
int n,m,x,y,k;
int gcd(int a,int b) { return b==0?a:gcd(b,a%b); }
void gcd(int a,int b,int &x,int &y)//拓展欧几里得算法
{
if(b==0) { x=1; y=0; return;}
gcd(b,a%b,x,y);
int t=x;
x=y; y=t-a/b*y;
}
bool check(int m)
{
for (int i=1;i<=n;i++)
for (int j=i+1;j<=n;j++)
{
int a=p[i]-p[j],b=m,cc=c[j]-c[i],d=gcd(a,b),x,y;
if (cc%d==0)
{
a/=d; b/=d; cc/=d;
gcd(a,b,x,y);
b=abs(b);
x=((x*cc)%b+b)%b;
while (!x) x+=b;
if (x<=min(l[i],l[j])) 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]);
m=max(c[i],m);
}
while (!check(m)) m++;//一直到不会出现野人在同一山洞为止
printf("%d",m);
return 0;
}