jzoj1158-荒岛野人【扩欧,gcd,同余方程】

最新推荐文章于 2024-09-01 11:01:32 发布

QuantAsk

最新推荐文章于 2024-09-01 11:01:32 发布

阅读量333

点赞数

分类专栏：数论and数学文章标签： Gcd 扩欧 jzoj 同余方程数论

本文链接：https://blog.csdn.net/Mr_wuyongcong/article/details/80970090

版权

数论and数学专栏收录该内容

261 篇文章 5 订阅

订阅专栏

正题

大意

有n个野人，每个野人有一个初始山洞 $C_i$ ，每次向前移动距离 $P_i$ ，寿命 $L_i$ ，如果野人走到了最后一个山洞那么继续就好回到第一个山洞，求至少多少个山洞才可以让野人们不会发生冲突。
这里写图片描述

解题思路

我们可以枚举答案，然后我们要判断是否冲突，我们可以枚举两个野人然后列出式子：

C i + P i x \equiv C j + P j x (m o d l)

$C_i+P_ix\equiv C_j+P_jx(mod\ \ l)$
如果这个式子有解且最小解

x<=Li/Lj x <= L i / L j $x<=L_i/L_j$ 那么就会有冲突。
我们先把式子转换一下

(C i - C j) + (P i - P j) x \equiv 0 (m o d l)

$(C_i-C_j)+(P_i-P_j)x\equiv 0(mod\ \ l)$

(C i - C j) + (P i - P j) x = l y

$(C_i-C_j)+(P_i-P_j)x=ly$

(P i - P j) x - l y = (C i - C j)

$(P_i-P_j)x-ly=(C_i-C_j)$
然后我们进行扩欧

代码

#include<cstdio>
#include<algorithm>
using namespace std;
int x,y,n,c[16],p[16],l[16],maxc,d,A,B,g;
bool flag;
int gcd(int a,int b)//同余
{
    if (b==0)
    {
        x=1;
        y=0;
        return a;
    }
    d=gcd(b,a%b);
    int k=x;
    x=y;
    y=k-a/b*y;
    return d;
}
int main()
{
    scanf("%d",&n);
    for (int i=1;i<=n;i++)
    {
        scanf("%d%d%d",&c[i],&p[i],&l[i]);
        maxc=max(maxc,c[i]);
    }
    while (1)
    {
        flag=false;
        for (int i=1;i<n;i++)
        {
          for (int j=i+1;j<=n;j++)
          {
            A=c[i]-c[j];B=p[j]-p[i];
            if (B<0)
            {
              A=-A;
              B=-B;//取反，因为是   -ly
            }
            d=gcd(B,maxc);
            g=maxc/d;
            g=((x*(A/d))%g+g)%g;//最小解
            if (!(A%d)&&g<=l[i]&&g<=l[j])//判断
              {flag=true;break;}
          }
          if (flag) break;
        }
        if (!flag) break;
        maxc++;//枚举
    }
    printf("%d",maxc);
}