CF338D GCD Table 题解

合集 - 题解(11)

1.P5192 Zoj3229 Shoot the Bullet|东方文花帖|【模板】有源汇上下界最大流2023-01-132.Phi的反函数2023-01-153.魔力屏障 (magic) 题解2023-06-01

4.CF338D GCD Table 题解2023-06-05

5.[AGC055B] ABC Supremacy 题解2023-06-096.[AGC055A] ABC Identity 题解2023-06-147.QOJ 6504. CCPC Final 2022 D Flower's Land 2题解2023-07-138.[SDOI2010] 代码拍卖会 题解2023-07-199.PERIODNI - Periodni 题解 & 笛卡尔树讲解 & 树状背包讲解2023-07-2010.Burnside 定理2023-07-3111.ZS Shuffles Cards 题解2023-08-16

收起

CF338D GCD Table 题解

题目描述

你有一个长度为 $k$ 的数列 $a$ ，

询问是否存在 $x\in [1,n]y\in [1,m]$ 使得 $\mathrm{\forall }igcd(x,y+i-1)=a_i$。

解析

我们转换一下可以得到：

$$\mathrm{\forall }i\{\begin{array}{c}x\equiv 0(\mathrm{mod}{a}_i)\\ y+i-1\equiv 0(\mathrm{mod}{a}_i)\end{array}$$

前面一个 $x$ 很好解决，直接最大公倍数。

$y$ 可以转化一下：

$$y\equiv 1-i(\mathrm{mod}{a}_i)$$

经典扩展中国剩余定理。

但是我们因为分开考虑的 $x$ 与 $y$ 得到的不一定是充要条件。

我们需要再次验证一下。

温馨提示

1.在运算过程中会爆 long long 需要龟速乘。

2.在运算过程中最大公倍数也就是 $x$ 会爆long long ，我们需要判断一下有没有超过 $n$ 如果超过直接 NO。

3.计算 $y$ 的时候有可能 $y=0$ 在这个情况下我们考虑 $y=0$ 和 $y=x$ 是等价的我们可以直接赋值成 $x$ 。

考虑 $x$ 还表示了 $a$ 数列的最大公倍数

代码

#include<bits/stdc++.h>
#define ll long long 
using namespace std;
const int N = 1e4 + 10;
int k;
ll n, m, a[N], x = 1,y = 0,M;
ll mmul(ll x,ll k,ll p){
    ll ans = 0,d = x;
    while(k){
        if(k & 1) ans = (ans + d) % p;
        k >>= 1;
        d = d * 2 % p;
    }
    return ans;
}
ll gcd(ll a,ll b){
    if(b == 0) return a;
    return gcd(b, a % b);
}
void exgcd(ll a,ll b,ll &x,ll &y){
    if(b == 0){
        x = 1,y = 0;
        return ;
    }
    exgcd(b, a % b, x, y);
    ll z = x;
    x = y;
    y = z - (a / b) * y;
}
ll jfc(ll a,ll b,ll c){
    ll x = 0,y = 0,g;
    g = gcd(a,b);
    a /= g,b /= g;
    if(c % g) return -1;
    exgcd(a, b, x, y);
    return mmul((x % b + b) % b,c / g,b);
}
void input(){
    cin>>n>>m>>k;
    for(int i = 1; i <= k; ++i){
        cin>>a[i];
    }
}
bool op(){
    for(int i = 1;i <= k; ++i){
        ll now = jfc(x,a[i],(((1 - i + a[i]- y) % a[i] + a[i]) % a[i]));
        if(now == -1) return 0;
        y = (mmul(now,x,x * a[i] / gcd(x,a[i])) + y) % (x * a[i] / gcd(x,a[i]));
        x = x * a[i] / gcd(x,a[i]);
        if(x > n) return 0;
    }    
    if(y == 0) y = x;
    if((y + k - 1) > m) return 0;
    return 1;
}
bool pd(){
    for(int i = 1;i <= k; ++i){
        ll now = gcd(x,y + i * 1ll - 1ll);
        if(now != a[i]) return 0 ;
    }
    return 1;
}
int main(){
    input();
    if(op() && pd()){
        cout<<"YES";    
    }else{
        cout<<"NO";
    }
    return 0;
}