Codeforces 1493D GCD of an array 数论，数据结构

题意：给一个长度为n的数组，q次询问，每次修改a[i]乘等x，求每次修改完之后数组的gcd
思路
①gcd由各个数的因子最小值构成，其中一个因子的数量如果到达n的话，那么这n个因子里的数量最小值就会对gcd有贡献
②记录一个数有哪些因子可以通过map来记录，同时支持修改，删除和查最小值的数据结构就是multiset。
③每次a[i]*x，可以先对x进行质因数分解，然后对原来的map和multiset进行维护。
维护思路：

• 如果之前已经有n个这个因子了，那么先对gcd除掉这个因子的贡献，然后加入这个数之后再在set中取因子的最小数对gcd进行累加
• 如果没有n个这个因子，那么直接从set中取出来这个数，然后修改map中的值再塞回set中
• 如果之前a[i]没有这个因子，那么直接修改map，然后塞到set里。如果有n个数了，修改gcd

#include <bits/stdc++.h>
using namespace std;
#define LL long long
template<typename T>void Read(T &x){x=0;char ch=getchar();LL f=1;while(!isdigit(ch)){if(ch=='-')f*=-1;ch=getchar();}while(isdigit(ch)){x=x*10+ch-48;ch=getchar();}x*=f;}
inline LL QuickPow(LL base,LL n,LL Mod=0){LL ret(1);while(n){if(n&1){ret*=base;if(Mod)ret%=Mod;}base*=base;if(Mod)base%=Mod;n>>=1;}return Mod?ret%Mod:ret;}
const int kMax = 2e5 + 10;
const int kMod = 1e9 + 7;
LL inv(LL num){
	return QuickPow(num,kMod-2,kMod);
}
multiset<LL>s[kMax];
map<LL,LL>mp[kMax];
LL gcd = 1;
int n,q;
void add(LL x,LL cnt,LL num)
{
	if(mp[x].count(num))
	{
		if(s[num].size() == n)
		{
			LL tmp = QuickPow(num,*s[num].begin(),kMod)%kMod;
			gcd = (gcd%kMod*inv(tmp))%kMod;
		}
		s[num].erase(s[num].find(mp[x][num]));
		mp[x][num] += cnt;
		s[num].insert(mp[x][num]);
		if(s[num].size() == n)
		{	
			gcd = (gcd%kMod*QuickPow(num,*s[num].begin(),kMod))%kMod;
		}
	}
	else
	{
		mp[x][num] = cnt;
		s[num].insert(cnt);
		if(s[num].size() == n)
		{
			gcd = (gcd%kMod*QuickPow(num,*s[num].begin(),kMod))%kMod;
		}
	}
}
void divide(LL x,LL num)
{
	for(int i = 2;i*i <= num;i++)
	{
		if(num%i==0)
		{
			int cnt = 0;
			while(num%i == 0)
			{
				cnt++;
				num /= i;
			}
			add(x,cnt,i);
		}
	}
	if(num > 1)
	add(x,1,num);
}

void solve()
{
	cin>>n>>q;
	for(int i = 1;i <= n;i++)
	{
		LL tmp;
		Read(tmp);
		divide(i,tmp);
	}
	for(int i = 1;i <= q;i++)
	{
		LL x,val;
		scanf("%lld%lld",&x,&val);
		divide(x,val);
		cout<<gcd<<'\n';
	}
}
signed main()
{
	//int _;
	//for(scanf("%d",&_);_;_--)
		solve();
    return 0;
}