codeforces 1349A（数学）

题意描述

思路

题意要求求出$gcd(lcm(a_i,a_j))|i<j$，由于$lcm(a_i,a_j)=a_i\ast a_j/gcd(a_i,a_j)$，故得到式子$gcd(a_i\ast a_j/gcd(a_i,a_j))$。接下来对式子进行化简。

$gcd(a_i\ast a_j/gcd(a_i,a_j))$
=$a_i\ast gcd(a_j/gcd(a_i,a_j))$
=$a_i\ast gcd(a_j)/gcd(gcd(a_i,a_j))$
=$a_i\ast gcd(a_j)/gcd(a_i,gcd(a_j))$

关于化简步骤的证明

第一步到第二步：
$gcd(a\ast k_i)=a\ast gcd(k_i)$，$1\le i\le n$
显然，无论$i$取何值，$a$都是$a\ast k_i$的公因子，故$gcd(a\ast k_i)=a\ast gcd(k_i)$，$1\le i\le n$
即$gcd(a_i\ast a_j/gcd(a_i,a_j))$=$a_i\ast gcd(a_j/gcd(a_i,a_j))$，$1\le i\le n$

第二步到第三步：
$gcd(a_i/b_i)$=$gcd(a_i)/gcd(b_i)$，$1\le i\le n$
保证$a_i$整除$b_i$
移项得$gcd(a_i/b_i)\ast gcd(b_i)$=$gcd(a_i)$，$1\le i\le n$
因为保证$a_i$整除$b_i$，故两部分$a_i/b_i$和$b_i$最大公约数的成绩一定与$gcd(a_i)$相同。
即$a_i\ast gcd(a_j/gcd(a_i,a_j))$=$a_i\ast gcd(a_j)/gcd(gcd(a_i,a_j))$，$1\le i\le j\le n$。

这样就可以扫一遍数组来记录下来每一个部分的最大公约数。

AC代码

#include<bits/stdc++.h>
#define x first
#define y second
#define PB push_back
#define mst(x,a) memset(x,a,sizeof(x))
#define all(a) begin(a),end(a)
#define rep(x,l,u) for(ll x=l;x<u;x++)
#define rrep(x,l,u) for(ll x=l;x>=u;x--)
#define sz(x) x.size()
#define IOS ios::sync_with_stdio(false);cin.tie(0);
using namespace std;
typedef unsigned long long ull;
typedef pair<int,int> PII;
typedef pair<long,long> PLL;
typedef pair<char,char> PCC;
typedef long long ll;
const int N=2*1e5+10;
const int M=1e6+10;
const int INF=0x3f3f3f3f;
const int MOD=1e9+7;
/*
独立思考
不要看测试样例
找性质
试着证明
写完后不盲目交
*/
ll a[N],ggcd[N];
void solve(){
    int n;cin>>n;
    rep(i,1,n+1) cin>>a[i];
    ll ans=0;
    rep(i,1,n+1){
        ggcd[i]=__gcd(a[i],ggcd[i-1]);
        ans=__gcd(a[i]*ggcd[i-1]/ggcd[i],ans);
    }
    cout<<ans<<endl;
}
int main(){
    IOS;
    //freopen("test.txt", "r", stdin);
    //freopen("test.txt", "w", stdout);
    //int t;cin>>t;
    //while(t--)
        solve();
    return 0;
}