【NOI2002T4】荒岛野人-扩展欧几里得

测试地址：抱歉...我能找到的测试地址只有一个私人OJ才有，所以就不贴出来了...

题目大意：一个岛是环状的，环上排列有M个洞穴，顺时针编号为1~M，有N(不超过15)个野人，第i个野人一开始住在洞穴Ci中，每一年要顺时针迁移Pi个洞穴，走Li年后就会死去。求满足在野人有生之年都不存在两个野人同住一个洞穴的情况下，最少的洞穴总数。保证M不超过10^6。

做法：因为答案不满足单调性（可能存在M=6时成立而M=7时却不成立的情况），所以我们只能从小到大枚举洞穴数M，查看其是否能满足条件，如果能满足条件，这个M就一定是最佳答案了。那么，怎么判断一个M满不满足条件呢？我们来讨论不满足条件的充要条件是什么，很显然就是两个野人在有生之年住在了同一个洞穴里，即下列同余方程有解且解小于等于Li和Lj：Ci+Pi*x=Cj+Pj*x(mod M)，即(Pi-Pj)*x+M*y=Cj-Ci，这就是经典的二元不定方程了，用扩展欧几里得算法先求出(Pi-Pj)*x+M*y=gcd(Pi-Pj,M)的解（为了简便，下文称gcd(Pi-Pj,M)为G），然后判断Cj-Ci能否整除G，如果不能，说明原方程没有整数解，如果能，则将x做以下处理：x=x*(Cj-Ci)/G,x=x mod (M/G)，这时候x就是一个最小的非负整数解了，这时候再判断解满不满足x≤Li且x≤Lj即可。对每两个野人都这么检查一遍，最后的复杂度是O(N^2*M)，可以通过全部数据。注意M不能小于任何一个Ci，以及负数取模的问题。

以下是本人代码：

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
#define ll long long
using namespace std;
int n;
ll c[20],p[20],l[20],g,m=0;

ll exgcd(ll a,ll &x,ll b,ll &y)
{
  ll x0,x1,y0,y1,t0,t1,q;
  x0=1,y0=0,x1=0,y1=1;
  q=a/b;
  t0=a,a=b,b=(t0%b+b)%b;
  while(b!=0)
  {
    t0=x0,t1=y0;
	x0=x1,y0=y1,x1=t0-q*x1,y1=t1-q*y1;
	x=x1,y=y1;
	q=a/b;
	t0=a,a=b,b=(t0%b+b)%b;
  }
  return a;
}

bool check()
{
  ll x,y,g;
  for(int i=1;i<=n;i++)
    for(int j=i+1;j<=n;j++)
	{
	  g=exgcd(p[i]-p[j],x,m,y);
	  if (((c[j]-c[i])%g+g)%g==0)
	  {
	    x=x*(c[j]-c[i])/g;
		x=(x%(m/g)+(m/g))%(m/g);
		if (x<=l[i]&&x<=l[j]) return 0;
	  }
	}
  return 1;
}

int main()
{
  scanf("%d",&n);
  for(int i=1;i<=n;i++)
  {
    scanf("%lld%lld%lld",&c[i],&p[i],&l[i]);
    m=max(m,c[i]);
  }
  
  while(!check())
  {
    m++;
  }
  
  printf("%lld",m);
  
  return 0;
}