GCD序列

（一）D. GCD of an Array

题目链接
题意：
给出一个长度为n的序列和q次询问，每次询问把序列中的一个元素扩大x倍，输出整个序列的gcd。
思路：
考虑把每个数质因数分解，用map和multiset实现存储，m[k][x]表示第k个数中含x的公因子的个数，multiset[x]中存储含x的个数序列。如果某个因子个数达到n个，表明每个数字都有该因子，把它算进贡献里面。

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int N=2e5+10,mod=1e9+7;
int n,q;
ll base;
multiset<int> mul[N];
map<int,int> m[N];

ll qpow(ll a,ll b){
    ll res=1;
    while(b){
        if(b&1) res=res*a%mod;
        b>>=1;
        a=a*a%mod;
    }
    return res;
}
void insert(int k,int x,int cnt){
    if(m[k].count(x)){
        if(mul[x].size()==n) base=base*qpow(qpow(x,*mul[x].begin()),mod-2)%mod;
        mul[x].erase(mul[x].find(m[k][x]));//没有find会擦去多个一样的
        m[k][x]+=cnt;
        mul[x].insert(m[k][x]);
        //cout<<x<<","<<mul[x].size()<<endl;
        if(mul[x].size()==n) base=base*qpow(x,*mul[x].begin())%mod;
    }
    else{
        mul[x].insert(cnt);
        m[k][x]=cnt;
        //cout<<x<<"  "<<mul[x].size()<<endl;
        if(mul[x].size()==n) base=base*qpow(x,*mul[x].begin())%mod;
    }
}
void divi(int k,int x){
    for(int i=2;i<=x/i;i++){
        int cnt=0;
        while(x%i==0){
            x/=i;
            cnt++;
        }
        if(cnt) insert(k,i,cnt);
    }
    if(x>1) insert(k,x,1);
}
int main(){
    base=1;
    cin>>n>>q;
    for(int i=1;i<=n;i++){
        int x;cin>>x;
        divi(i,x);
    }
    while(q--){
        int i,x;cin>>i>>x;
        divi(i,x);
        cout<<base<<endl;
    }
return 0;
}

（二）C. Row GCD

题目链接
题意：
给出一个长度为n的序列和q次询问，每次询问把序列中的每个元素加上x，输出整个序列的gcd。
思路：
我们都知道gcd(a,b)=gcd(a,b-a)
对于多个数组而言，gcd(a,b,c,d)=gcd(a,b-a,c-b,d-c)
所以gcd(a+x,b+x,c+x,d+x)=gcd(a+x,b-a,c-b,d-c)
这就是gcd的更相减损性。

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int N=2e5+10,mod=1e9+7;
int n,m;
ll a[N],b[N];

int main(){
    cin>>n>>m;
    for(int i=0;i<n;i++) cin>>a[i];
    for(int i=0;i<m;i++) cin>>b[i];
    if(n==1){
        for(int i=0;i<m;i++) cout<<a[0]+b[i]<<' ';
    }
    else{
        sort(a,a+n);
        ll ans=a[1]-a[0];
        for(int i=1;i<n;i++) ans=__gcd(ans,a[i]-a[i-1]);
        for(int i=0;i<m;i++) cout<<__gcd(ans,a[0]+b[i])<<' ';
    }
return 0;
}