Codeforces Round #641 (Div. 2) C. Orac and LCM（gcd和lcm的重要理解）

Happig丶

已于 2023-03-22 00:22:54 修改

阅读量288

点赞数 1

文章标签： c语言算法 c++

于 2020-05-14 19:30:28 首次发布

本文链接：https://blog.csdn.net/qq_44691917/article/details/106126238

版权

传送门

题目大意

给出一个序列 $a$ ，求出 $gcd(\{lcm(a_i,a_j) | i<j\})$ 。

做法一

实际上不难发现序列中 $a_1$ 对应的 $l c m$ 对为 $n - 1$ 对， $a_2$ 对应的 $n - 2$ 对…以此类推，关键是如何化简如下的表达式：

$gcd \{ lcm(a_1,a_2)，lcm(a_1,a_3)，...，lcm(a_1,a_n) \}$

GCD和LCM在质因数分解下的意义

假设 $a$ 经过质因数分解后为 $P_1^{k_1}*P_2^{k_2},...*P_n^{k_n}$ ， $b$ 经过质因数分解为 $Q_1^{t_1}*Q_2^{t_2},...*Q_n^{t_n}$

假设质因子 $p$ 为 $a$ 和 $b$ 共有的质因子，在 $a$ 中 $p$ 的幂次为 $k_1$ ，在 $b$ 中 $p$ 的幂次为 $k_2$ ，又因为所有的质因子两两互质，互不影响

那么 $g c d (a, b)$ 的质因子 $p$ 实际上的幂次为 $min\{ k_1,k_2 \}$ ， $l c m (a, b)$ 的质因子 $p$ 实际上的幂次为 $max\{ k_1,k_2 \}$

那么现在我们再看上式，对于所有数的一个公共的质因子 $p$ ，假设幂次分别为 $k_1,k_2,...k_n$ ，上式的意义为：

$min \{ max(k_1,k_2)，max(k_1,k_3)，...，max(k_1,k_n) \}$

然后该式等价于：

$max \{ k_1，min \{k_2,k_3,...,k_n) \} \}$

即化为表达式：

$lcm \{ k_1，gcd \{ k_2,k_3,...,k_n \} \}$

因此我们得到等式：

$gcd \{ lcm(a_1,a_2)，lcm(a_1,a_3)，...，lcm(a_1,a_n) \} = lcm \{ k_1，gcd \{ k_2,k_3,...,k_n \} \}$

其它收获

得知上述等式后，我们对左边和右边 $l c m$ 化 $g c d$ 得到：
在这里插入图片描述
然后两边约去 $a_1$ ：

对于左式的分母，等价于 $gcd\{ a_1,gcd\{ a_2,a_3,...,a_n \} \}$ ，那么我们将 $a_1$ 看做 $x$ ，将 $gcd\{ a_2,a_3,...,a_n \}$ 看作 $y$ ，那么左式实际上变成了：

$y / g c d (x, y)$

那么，实际意义就是将 $y$ 和 $x$ 相同的最大公因数从 $y$ 中约去，而这时再看右式，不难发现实际就是将 $y$ 中的每一部分与 $x$ 相同的最大公因数约去，又因为 $g c d$ 运算和乘法的交换性（即 $k * g c d (a, b) = g c d (k * a, k * b)$ ）最后再求 $g c d$ 即约去了整体与 $x$ 的最大公因数

最后对于本题，我们只需预处理后缀的 $g c d$ 即可

PS：本题的公式还可用另外方法证明，见mrcrack的博客

#include <set>
#include <map>
#include <stack>
#include <queue>
#include <math.h>
#include <cstdio>
#include <string>
#include <cstring>
#include <sstream>
#include <iostream>
#include <algorithm>
#include <unordered_map>
using namespace std;
#define lowbit(x) (x&(-x))
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> P;
const double eps=1e-8;
const double pi=acos(-1.0);
const int inf=0x3f3f3f3f;
const ll INF=1e18;
const int Mod=1e9+7;
const int maxn=2e5+10;

int a[maxn],sub[maxn];

ll gcd(ll a,ll b){
    return b==0?a:gcd(b,a%b);
}

ll lcm(ll a,ll b){
    return a/gcd(a,b)*b;
}

int main(){
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
    int n;
    cin>>n;
    for(int i=1;i<=n;i++) cin>>a[i];
    sub[n]=a[n];
    for(int i=n-1;i>=1;i--)
        sub[i]=gcd(sub[i+1],a[i]);
    ll ans=0;
    for(int i=1;i<=n;i++)
        ans=gcd(ans,lcm(a[i],sub[i+1]));
    cout<<ans<<endl;
    return 0;
}

做法二

根据上述在质因数分解意义下的 $g c d$ 和 $l c m$ ，实际上就是需要求出 $n$ 个数中同时出现最少 $n - 1$ 次的质因子中幂次最小和次小的那个，具体做法要出现 $n$ 次和 $n - 1$ 次来定。

//
// Created by Happig on 2021/1/24.
//
#include <bits/stdc++.h>
#include <unordered_map>
#include <unordered_set>

using namespace std;
#define ENDL "\n"
#define lowbit(x) (x & (-x))
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;
typedef pair<double, double> pdd;
const double eps = 1e-8;
const double pi = acos(-1.0);
const int inf = 0x3f3f3f3f;
const double dinf = 1e300;
const ll INF = 1e18;
const int Mod = 1e9 + 7;
const int maxn = 2e5 + 10;

unordered_map<int, int> fi, se;
unordered_set<int> st;
int minp[maxn], prime[maxn], vis[maxn];
int cnt;

void euler() {
    for (int i = 1; i < maxn; i++) minp[i] = i;
    for (int i = 2; i < maxn; i++) {
        if (minp[i] == i) {
            prime[cnt++] = i;
        }
        for (int j = 0; j < cnt && i * prime[j] < maxn; j++) {
            minp[i * prime[j]] = prime[j];
            if (i % prime[j] == 0) break;
        }
    }
}

void divide(int n) {
    map<int, int> mp;
    int fac, p, num;
    while (n > 1) {
        fac = minp[n];
        n /= fac;
        mp[fac]++;
    }
    for (auto i:mp) {
        p = i.first, num = i.second;
        vis[p]++, st.insert(p);
        if (!fi.count(p)) {
            fi[p] = num;
            continue;
        }
        if (!se.count(p)) {
            if (num < fi[p]) {
                se[p] = fi[p];
                fi[p] = num;
            } else se[p] = num;
            continue;
        }
        if (num <= fi[p]) {
            se[p] = fi[p];
            fi[p] = num;
        } else if (num < se[p]) {
            se[p] = num;
        }
    }
}

ll qkp(ll x, ll n) {
    ll ans = 1;
    while (n) {
        if (n & 1) ans = ans * x;
        x *= x;
        n >>= 1;
    }
    return ans;
}

int main() {
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    ios_base::sync_with_stdio(0), cin.tie(0), cout.tie(0);
    int n;
    cin >> n;
    euler();
    for (int i = 1, x; i <= n; i++) {
        cin >> x;
        divide(x);
    }
    ll ans = 1;
    for (auto p:st) {
        //cout << p << ": " << fi[p] << " " << se[p] << ENDL;
        if (vis[p] >= n - 1) {
            if (vis[p] == n - 1) ans *= qkp(p, fi[p]);
            else ans *= qkp(p, se[p]);
        }
    }
    cout << ans << ENDL;
    return 0;
}

Happig丶

关注

1
点赞
踩
1

收藏

觉得还不错? 一键收藏
0
评论
Codeforces Round #641 (Div. 2) C. Orac and LCM（gcd和lcm的重要理解）

传送门实际上不难发现序列中a1对应的lcm对为n-1对，a2对应的n-2对…以此类推，关键是如何化简如下的表达式：gcd { lcm(a1,a2)，lcm(a1,a3)，…，lcm(a1,an) }GCD和LCM在质因数分解下的意义假设a经过质因数分解后为P1k1*P2k2,…*Pnkn，b经过质因数分解为Q1t1*Q2t2,…*Qmtm假设质因子p为a和b共有的质因子，在a中p的幂次为k1，在b中p的幂次为k2，又因为所有的质因子两两互质，互不影响那么gcd(a,b)的质因子p实际上的幂次
复制链接

扫一扫