[NOI2002]荒岛野人【洛谷P2421】

题目描述
克里特岛以野人群居而著称。岛上有排列成环行的M个山洞。这些山洞顺时针编号为1,2,…,M。岛上住着N个野人，一开始依次住在山洞C1,C2,…,CN中，以后每年，第i个野人会沿顺时针向前走Pi个洞住下来。

每个野人i有一个寿命值Li，即生存的年数。

下面四幅图描述了一个有6个山洞，住有三个野人的岛上前四年的情况。三个野人初始的洞穴编号依次为1，2，3；每年要走过的洞穴数依次为3，7，2；寿命值依次为4，3，1。

奇怪的是，虽然野人有很多，但没有任何两个野人在有生之年处在同一个山洞中，使得小岛一直保持和平与宁静，这让科学家们很是惊奇。他们想知道，至少有多少个山洞，才能维持岛上的和平呢？

输入格式
第1行为一个整数N(1<=N<=15)，即野人的数目。

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

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

输入输出样例
输入 #1 复制
3
1 3 4
2 7 3
3 2 1
输出 #1 复制
6
说明/提示
对于50% 的数据：N 的范围是[1…1,000]。

对于另外50% 的数据：N 的范围是[1…100,000]。

对于100% 的数据：C 的范围是[1…1,000,000,000]，N 个整数中每个数的范围是：[0…1,000,000,000]。

题解

不能用二分

1. c[i]+p[i]*T$\equiv$c[j]+ p[j]*T(mod m)
2. T*(p[i]-p[j])$\equiv$c[j]-c[i](mod m)
3. 当对于任意i,j 满足T无解或T>min(p[i],p[j])时，m符合要求
4. 一个小小的优化：从C[i]的最大值开始枚举

#include<cstdio>
#include<iostream>
#include<algorithm>
using namespace std;
const int N=20;
int n,C[N],p[N],l[N];
int x,y,d;
int exgcd(int a,int b,int &x,int &y) {
    if(!b) 
	{
		x=1;
		y=0;
		return a;
	}
	else
	{
		d=exgcd(b,a%b,y,x);
    	y-=a/b*x;
    	return d;
	}
}
bool check(int m) 
{
	int i,j;
    for(i=1;i<=n;++i) 
		for(j=i+1;j<=n;j++) 
		{
        	int a=p[i]-p[j];
			int b=m;
			int c=C[j]-C[i];
			if(a<0)
			{
				a=-a;
				c=-c;
			}
        	d=exgcd(a,b,x,y);
        	if(c%d) 
				continue;
        	int t=b/d;
        	if(t<0)t=-t;
        	x=(x*(c/d)%t+t)%t;
        	if(x<=l[i]&&x<=l[j]) 
				return false;
    	}
    return true;
}
int maxx(int x,int y)
{
	return x>y?x:y;
}
int main() 
{
	//freopen("1757.in","r",stdin);
	//freopen("1757.out","w",stdout);
    scanf("%d",&n);
    int mx=0,i;
    for(i=1;i<=n;i++) 
	{
		scanf("%d%d%d",&C[i],&p[i],&l[i]);
		mx=maxx(mx,C[i]);
	}	
    for(i=mx;i<=1000000;i++)
		if(check(i)) 
		{
			printf("%d\n",i)；
		 	return 0;
		 	}
    return 0;
}