[Noi2002]$Savage$

---

其实答案远不到1e6

所以可以枚举！

设答案是m

那\(i,j\)的相遇就可以表示成\(P_ix+C_i=P_jx+C_j+ym\)

移向就是\((P_i-P_j)x-ym=C_j-C_i\)

套扩展欧几里得定理

如果\(C_j-C_i\mod gcd\ !=0\)说明不会相遇

否则的话求一下第一次相遇时间，如果大于一个的寿命也不会相遇

---

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define M 1000010
using namespace std;

int n,m,k,c[M],l[M],p[M],maxx,x,y;

int exgcd(int a,int b,int &x,int &y)
{
    if(!b) {x=1, y=0; return a;}
    int tmp=exgcd(b,a%b,y,x);
    y-=a/b*x; return tmp;
}

bool check(int k)
{
    for(int i=1;i<=n;i++)   
        for(int j=i+1;j<=n;j++) 
        {
            int a=p[i]-p[j], b=c[j]-c[i];
            if(a<0) a=-a, b=-b;
            int g=exgcd(a,k,x,y);
            if(b%g) continue;
            x*=b/g; int f=k/g;
            x=(x%f+f)%f; if(!x) x=f;
            if(x<=min(l[i],l[j])) return 0;
        }
    return 1;
}

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

转载于:https://www.cnblogs.com/ZUTTER/p/10369959.html