CF1493D GCD of an Array (分解定理+思维)

最新推荐文章于 2024-07-24 11:16:44 发布

_九磅十五便士_

最新推荐文章于 2024-07-24 11:16:44 发布

阅读量102

点赞数

分类专栏： codeforces训练 # 思维 # 数论

本文链接：https://blog.csdn.net/m0_51841071/article/details/119461162

版权

codeforces训练同时被 3 个专栏收录

78 篇文章 2 订阅

订阅专栏

思维

54 篇文章 0 订阅

订阅专栏

数论

42 篇文章 1 订阅

订阅专栏

这篇博客介绍了一道编程题目，涉及在序列中进行元素乘法操作并求最大公约数（GCD）的问题。博主利用Multiset数据结构来存储序列中的因子，并通过迭代找到所有数的公共因子，从而计算GCD。文章详细阐述了思路，包括预处理素数、维护Multiset以及处理修改操作的过程。最后，博主给出了完整的C++代码实现。

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

题目传送门

这道题有一个非常牛逼的数据结构就是 $m u l t i s e t$

题意

先给你一个序列，然后给你几组询问，询问的第一个数是在第x个位置上，第二个数p代表，在第x个位置上乘p。请找到这个序列的gcd，并输出。

思路

这道题肯定是让你找到所有数的公共的因子，那么就可以以因子为下标，存到 $m u l t i s e t$ 里面，同时开一个数组标记一下，这个数作为了几个数的因子。如果能找到这个因子，并这个数已经作为全部的因子，答案乘这个数。

#include<iostream>
#include<set>
#define huzanjujutaiqiangle signed main()
//#define int long long 
#define MAXN 250005
//#define MOD 1e9+7
using namespace std;
typedef long long ll;
multiset<int>s[MAXN];
const ll MOD=1e9+7;
int p[MAXN<<1];
bool st[MAXN<<1];
int primes[MAXN<<1];
int cnt;
ll gcd=1;
int _hash[MAXN<<1];
int n,m;

void get_minmu(){
	p[1]=1;
	for(int i=2;i<=500000;i++){
		/*if(i==500000){
			cout<<"ok"<<endl;
		}*/
		//cout<<i<<endl;
		if(!st[i]){
			p[i]=i;
			primes[++cnt]=i;	
		}
		for(int j=1;j<=cnt&&i*primes[j]<=MAXN;j++){
			st[primes[j]*i]=true;
			p[primes[j]*i]=primes[j];
			if(i%primes[j]==0)break;
		}
	}
	
}
int modify(int x,int t){
	while(t!=1){
		if(s[p[t]].find(x)==s[p[t]].end()){
			++_hash[p[t]];
		}
		s[p[t]].insert(x);
		if(_hash[p[t]]==n){
			gcd*=(ll)p[t];
			gcd%=MOD;
			_hash[p[t]]=0;
			for(int i=1;i<=n;i++){
				s[p[t]].erase(s[p[t]].find(i));
				if(s[p[t]].find(i)!=s[p[t]].end()){
					++_hash[p[t]];
				}
			}
		}
		t/=p[t];
	}
}
huzanjujutaiqiangle{
	get_minmu();
	cin>>n>>m;
	for(int i=1;i<=n;i++){
		int t;
		scanf("%d",&t);
		//int t;
		//cin>>t;
		modify(i,t);
	}
	while(m--){
		int x,t;
		scanf("%d %d",&x,&t);
		//cin>>x>>t;
		modify(x,t);
		printf("%lld\n",gcd);
		//cout<<gcd<<endl;
	}
	//cout<<gcd<<endl;
}