Codeforces 1349 A

最新推荐文章于 2022-10-06 00:05:58 发布

幼儿园大哥~

最新推荐文章于 2022-10-06 00:05:58 发布

阅读量112

点赞数

分类专栏：数论 Codeforces

本文链接：https://blog.csdn.net/weixin_45633221/article/details/110872401

版权

Codeforces 同时被 2 个专栏收录

40 篇文章 0 订阅

订阅专栏

数论

35 篇文章 0 订阅

订阅专栏

这篇博客探讨了一个涉及数学和编程的问题，即给定一个序列，计算序列中两两元素最大公约数（GCD）的最小公倍数（LCM）集合的GCD。通过分解质因数并分析质因子的贡献，可以有效地解决这个问题。代码示例使用了C++，实现了质因子分解和筛法，然后根据质因子幂次计算贡献，最终得出答案。

摘要由CSDN通过智能技术生成

传送门

题目大意

给定长度为 $n$ 的序列 $a_1,a_2...a_n$ .
设可重集 $S=\{{lcm(a_i,a_j)|1\leq i<j\leq n}\}$ .
求出 $g c d (S)$
$1\leq n\leq 10^5,1\leq a_i\leq 2*10^5$ .

思路

分解质因数，记录每个质因子出现的幂次然后考虑贡献
单独考虑一个质因子的贡献：假设 $p$ 这个质因子幂次分别为 $k_1,k_2$ ,那么 $l c m (a, b)$ 中 $p$ 这个质因子的幂次为 $max(k_1,k_2)$
若 $n$ 个数都有质因子 $p$ , $p$ 的贡献应该是第二小的幂次
若 $n - 1$ 个数中都有质因子 $p$ , $p$ 的贡献应该是第一小的幂次
将所有的质因子的贡献乘起来就是答案

代码

vector<int>v[maxn*2];

ll fun(ll a,ll b)
{
    ll ans=1;
    while(b)
    {
        if(b&1)
        ans=1ll*ans*a;
        a=1ll*a*a;
        b>>=1;
    }
    return ans;
}

bool vis[maxn];//标记非素数,0是素数
int primer[maxn/10];//存素数
int cnt=0;//记录素数个数，
void find_primer(){
    for(int i=2;i <=maxn;i++){
        if(!vis[i])primer[cnt++]=i;
        for(int j=0;j<cnt&&primer[j]*i<=maxn;j++){
            vis[i*primer[j]]=1;//筛
            if(i%primer[j]==0)break;//关键！！！找到i*primer[j]的最小质因子primer[j]，退出
        }
    }    
}

int main(){
	int n;
	cin>>n;
	int x;
	find_primer();
	for(int j=1;j<=n;j++){
		cin>>x;
		for(int i=0;i<cnt&&primer[i]*primer[i]<=x;i++){
			if(x%primer[i]==0){
				int cnt1=0;
				while(x%primer[i]==0){
					cnt1++;
					x/=primer[i];
				}
				v[primer[i]].push_back(cnt1);
			}
		}
		if(x>1) v[x].push_back(1);
	}
	ll ans=1;
	for(int i=1;i<=200000;i++){
		sort(v[i].begin(),v[i].end());
		if(v[i].size()==n){
			ans*=fun(i,v[i][1]);
		}
		else if(v[i].size()==n-1){
			ans*=fun(i,v[i][0]);
		}
	}
	printf("%lld\n",ans);
}