【线性模方程】BZOJ 1407 [Noi2002]Savage 题解

最新推荐文章于 2021-09-03 08:16:27 发布

「已注销」

最新推荐文章于 2021-09-03 08:16:27 发布

阅读量137

点赞数

分类专栏： BZOJ题解 ======数论====== 模方程文章标签：模方程 BZOJ

本文链接：https://blog.csdn.net/try__jhf/article/details/81369571

版权

BZOJ题解同时被 3 个专栏收录

34 篇文章 0 订阅

订阅专栏

======数论======

8 篇文章 0 订阅

订阅专栏

模方程

2 篇文章 0 订阅

订阅专栏

题目大意

$n$ 个人在 $m$ 个点的环上，每个人初始在点 $c_i$ ，能活 $L_i$ 秒，每秒会向后走 $p_i$ 步，求出最小的 $m$ 使得任意两个人在有生之年都不会在一个点上相遇。

解题分析

转化一下这题，如果两个人 $i,j$ 在第 $x$ 天相遇了，说明什么？

$c_i+x*p_i\equiv c_j+x*p_j(Mod \ m)$

转移

$(p_i-p_j)x \equiv c_i-c_j(Mod \ m)$

就是线性模方程，如果算出x无解（这辈子都碰不上）或 $x>min(L_i,L_j)$ ，则 $i,j$ 一定碰不到，这样可以先找出一个m，再用 $n^2$ 个方程来验证，但是这不满足二分，所以……穷举大法好（ $ans<=10^6$ ）

示例程序

(传送门)

#include<cstdio>
#define LL long long
using namespace std;
int n;
LL a[20],b[20],c[20],s,ans;
inline void readL(LL &x){
    x=0; char ch=getchar();
    while ('0'>ch||ch>'9') ch=getchar();
    while ('0'<=ch&&ch<='9') {x=x*10+ch-'0'; ch=getchar();}
}
void exgcd(LL a,LL b,LL &d,LL &x,LL &y){
    if (!b) {d=a; x=1; y=0;}
    else {exgcd(b,a%b,d,y,x); y-=x*(a/b);}
}
bool _check(LL tem){
    for (int i=1;i<n;i++)
        for (int j=i+1;j<=n;j++){
            LL C=((a[i]-a[j])%tem+tem)%tem,A=((b[j]-b[i])%tem+tem)%tem,g,x,y;
            exgcd(A,tem,g,x,y);
            if (C%g) continue;
            int m=tem/g,k=(x*C/g%m+m)%m;
            if (k<=c[i]&&k<=c[j]) return false;
        }
    return true;
}
int main()
{
    freopen("savage.in","r",stdin);
    freopen("savage.out","w",stdout);
    scanf("%d",&n); s=0;
    for (int i=1;i<=n;i++){
        readL(a[i]); readL(b[i]); readL(c[i]); s=(s<a[i])?a[i]:s;
    }
    for (ans=s;;ans++) if (_check(ans)) {printf("%d",ans); break;}
    return 0;
}