BZOJ1407: [Noi2002]Savage

BZOJ1407

直接枚举$m$，然后$n^2$判断。
为了方便处理，将所有$C_i$先-1.
对于任意两个野人$i,j$，令他们第$x$天相遇
则有$(C_i+x*p_i)mod \space M=(C_j+x*p_j)mod \space M$
$(C_i-C_j)=x(p_j-p_i) mod \space M$
可以列出方程$(p_j-p_i)*x+M*y=(C_i-C_j)$
方程无解或者x最小非负整数解 <=min(Li,Lj) <script type="math/tex" id="MathJax-Element-2995"><=min(L_i,L_j)</script>，则满足。
$exgcd$解出的$x,y$是$(p_j-p_i)*x+M*y=Gcd((p_j-p_i),M)$
所以求得的$x0=x*\frac{(C_i-C_j)}{gcd}$
通解$X=x0+\frac{M}{gcd}*k$
所以最小非负整数解为$(x0+\frac{M}{gcd})mod \space \frac{M}{gcd}$
注意开始位置要从$Max(C_i)+1$开始

【代码】

#include <cstdio>
#include <iostream>
#include <algorithm>
#define N 100005
#define Mod 999911659
#define INF 0x7fffffff
using namespace std;
typedef long long ll;

ll read()
{
    ll x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
    return x*f;
}

int n,ans,mn;
int C[N],P[N],L[N];

void Exgcd(int a,int b,int &d,int &x,int &y) {
    if(b) {
        Exgcd(b,a%b,d,y,x);
        y-=(a/b)*x;
    }
    else x=1,y=0,d=a;
}

int Check(int i,int j,int M)
{
    int x,y,d;
    int a=(P[j]-P[i]+M)%M,b=(C[i]-C[j]+M)%M;
    Exgcd(a,M,d,x,y);
    if(b%d) return INF;
    x*=b/d;M/=d;
    return (x%M+M)%M;
}

bool Judge(int M)
{
    for(int i=1;i<=n;i++)
    for(int j=1;j<i;j++) if(Check(i,j,M)<=min(L[i],L[j])) 
        return false;
    return true;
}

int main()
{   
    n=read();
    for(int i=1;i<=n;i++) C[i]=read()-1,P[i]=read(),L[i]=read(),mn=max(mn,C[i]);
    for(int x=mn+1;x<=1e6;x++) if(Judge(x)) 
    {
        printf("%d\n",x);
        return 0;
    }
    return 0;
}