测试地址:抱歉...我能找到的测试地址只有一个私人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;
}