C - Orac and LCM CodeForces - 1350C 数论

https://vjudge.net/contest/374073#problem/C
第一种方法是传统的质因数分解。找到一个质因数，若在n个式子中出现n次，如果次小值的指数为p，则最终结果贡献为n^p；若出现n-1次，最小值指数为p,则最终贡献为n P;若小于n-1次，则贡献为0。
第二种方法是公式法，用到的公式为GCD（lcm(a1,a2),lcm(a1,a3),…lcm(a1,am) ) = lcm(a1,gcd(a2,a3…am))
下面为证明

GCD（lcm(a1,a2),lcm(a1,a3),...lcm(a1,am) ) = lcm(a1,gcd(a2,a3...am))
证明 ：
从GCD(lcm(a1,a2),lcm(a1,a3))=lcm(a1,gcd(a2,a3))开始
左式可表达成gcd(a1*a2/gcd(a1,a2),a1*a3/gcd(a1,a3)
		=a1*gcd(a2/gcd(a1,a2),a3/gcd(a1,a3))
右式可表达成a1*gcd(a2,a3)/gcd(a1,a2,a3)

故目标为证明『gcd(a2/gcd(a1,a2),a3/gcd(a1,a3)) = gcd(a2,a3)/gcd(a1,a2,a3)』

由于每个数可以由若干个质因数的指数乘积表示，因此
设a1 = p1^e1 * p2^e2...pm^em
  a2 = p1^f1 * p2^f2...pm^fm
  a3 = p1^q1 * p2^q2...pm^qm
  A = a2/gcd(a1,a2) = p1^(f1-min(e1,f1))*...pm^(fm-min(em,fm))
  					= p1^(min(0,f1-e1))*...pm^(min(0,fm-em))
  B = a3/gcd(a1,a3) = p1^(min(0,g1-e1))*...pm^(min(0,gm-em))
  
  『gcd(A,B) = p1^min(0,f1-e1,g1-e1)*...pm^min(0,fm-em,gm-em)』①
  
  再求gcd(a2,a3) = p1^min(f1,g1)*...pm^min(fm,gm)
  	  gcd(a1,a2,a3)=p1^min(e1,f1,g1)*...pm^min(em,fm,gm)
  因此『gcd(a2,a3)/gcd(a1,a2,a3) = p1^(min(f1,g1)-min(e1,f1,g1))*...pm^(min(fm,gm)-min(em,fm,gm))②』
  因此，当ei=min(fi,gi,ei)且fi=min(fi,gi)时②式的指数为fi-ei
  		 ei=min(fi,gi,ei)且gi=min(fi,gi)时②式的指数为gi-ei
  		 fi=min(fi,gi,ei)且gi=min(fi,gi)时②式的指数为0
  因此①式和②式相等，既gcd(a2/gcd(a1,a2),a3/gcd(a1,a3)) = gcd(a2,a3)/gcd(a1,a2,a3)
  所以得证GCD(lcm(a1,a2),lcm(a1,a3))=lcm(a1,gcd(a2,a3))
  含有更多的数时证明方法同理。

该题代码如下

//#include<bits/stdc++.h>//poj,hdu can't use it
#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
#include <stack>
#include <map>
#include <cmath>
#include <algorithm>
#include <iomanip>
using namespace std;
/************************head file************************/
#define inf 0x3f3f3f3f
#define SI(x) scanf("%d",&x)
#define set0(x) memset(x,0,sizeof(x))
#define ms(x,n) memset(x,n,sizeof(x));

typedef long long ll;
const int mod = 998244353;
const int maxn = 100010;
using namespace std;
/***********************defination************************/
ll f[maxn],s[maxn],t[maxn],ans;
ll gcd(ll a,ll b){
	return b?gcd(b,a%b):a;
}
ll lcm(ll a,ll b){
	return a*b/gcd(a,b);
}

int main()
{
    int n,i;
	cin>>n;
	for(i = 1;i <= n;i++)
        cin>>f[i];
	s[n] = f[n];
	for(i = n-1;i >= 1;i--)s[i]=gcd(s[i+1],f[i]);
	for(i = 1;i < n;i++)t[i]=lcm(f[i],s[i+1]);
	ans = t[1];
	for(i = 2;i < n;i++)ans=gcd(ans,t[i]);
	cout<<ans<<endl;
	return 0;
}