CF - 1350C | 数论

最新推荐文章于 2022-08-21 18:39:33 发布

猹猹猹

最新推荐文章于 2022-08-21 18:39:33 发布

阅读量212

点赞数

分类专栏：思维好题数论基础

本文链接：https://blog.csdn.net/weixin_44855285/article/details/106171484

版权

思维好题同时被 2 个专栏收录

9 篇文章 0 订阅

订阅专栏

数论基础

3 篇文章 0 订阅

订阅专栏

题目

求 $gcd(lcm({a_i,a_j}) | i<j)$ ，其中 $2≤n≤100000，a_i≤2e6$ 。

$g c d (l c m (a, b), l c m (a, c), l c m (a, d)) = l c m (a, g c d (b, c, d))$

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<set>
#include<map>
#include<list>
#include<stack>
#include<queue>
#include<vector>
#define cl (k<<1)
#define cr (k<<1|1)
#define Mid ((a[k].l+a[k].r)>>1)
#define mem(a,x) memset(a,x,sizeof(a))
using namespace std;
#define bug printf("???\n")

typedef long long LL;

const int Inf=0x3f3f3f3f;
const double eps=1e-7;
const int maxn=1e5+50;
//gcd(lcm(a,b), lcm(a,c), lcm(a,d))
// = lcm(a, gcd(b,c,d))

LL a[maxn], suf[maxn];

LL lcm(LL a,LL b){
    return a*b/__gcd(a,b);
}
int main()
{
    int n;
    cin>>n;
    for(int i=1; i<=n; i++){
        scanf("%lld",&a[i]);
    }
    suf[n] = a[n];
    for(int i=n-1; i>=1; i--){
        suf[i] = __gcd(a[i],suf[i+1]);
    }
    LL ans=lcm(a[1], suf[2]);
    for(int i=2; i<=n-1; i++){
        ans = __gcd(ans, lcm(a[i], suf[i+1]));
    }
    printf("%lld\n",ans);
}