D. GCD of an Array(multiset+map)

D. GCD of an Array(multiset+map)

D. GCD of an Array

总结：1.multiset使用的一些细节
2.分解质因数的想法
3.(b/a)%mod=( b * a^(mod-2) ) % mod;

multiset

multiset是多集合，与set区别在于，它是可以重复插入，就这点区别. 有以下的操作方法： insert（传入数值,插入数值） find（传入数值,寻找值返回最先出现的地址） erase（传入数值,删除一数值，如果有多个同样的值，全删除).
(传地址，删除地址对应的值) ps：那我要只删除一个值（含有重复值）怎么做呢. 那就这样 a.erase(a.find(elem));

#include<iostream>
#include<set>
using namespace std;
multiset<int> a;
void print(){
     //默认排序是从小到大.
    for (multiset<int>::iterator i = a.begin(); i != a.end();i++){
        cout << *i << " ";
    }
    cout << endl;
}
int main (){
    a.insert(3);
    a.insert(5);
    a.insert(4);
    a.insert(5);
    a.insert(3);
    a.insert(6);
    a.insert(4);
    print();
    a.erase(3);
    print();
    a.erase(a.find(4));
    print();
    getchar();getchar();
    return 0;
}

题意：长度为n的b数组,m次操作，每次操作一次打印出这一数组的最大公约数.
思路：用multiset创建 一维数组 a[N]，
那么a[k] 这个集合含义：质数为k的，存放每个数（b[i]）的幂次数。
中间用到了逆元：(b/a)%mod=( b * a^(mod-2) ) % mod;取余传送门

代码如下：

#include<iostream>
#include<algorithm>
#include<math.h>
#include<cstring>
#include<vector>
#include<queue>
#include<set>
#include<map>
#define ll long long
ll gcd(ll a,ll b){ return b? gcd(b,a%b):a;}
const int N=2e5+10;
const ll mod=1e9+7;
ll read(){
    ll s = 0, f = 1; char ch = getchar();
    while(!isdigit(ch)){
        if(ch == '-') f = -1;
        ch = getchar();
    }
    while(isdigit(ch)) s = (s << 3) + (s << 1) + (ch ^ 48), ch = getchar();
    return s * f;
}
using namespace std;
int n, m;
ll sum = 1;
multiset<int> a[N]; // 质数为k的，存放每个数（b[i]）的幂次数
map<int, int> mp[N];  //mp[i][p]=cnt  b[i] 质数为p  p的幂次数
ll fastpow(ll a,ll b){
    ll sum = 1;
    while(b){
        if(b&1){
            sum = (sum * a) % mod;
        }
        a = (a*a)%mod;
        b >>= 1;
    }
    return sum;
}
ll niuyuan(ll a){
    return fastpow(a, mod - 2);
}
void insert(int k,int i,int cnt){ //b[k] i质数 cnt幂次数
        if(mp[k][i]){
                //gcd去掉
            if(a[i].size()==n)
                sum = sum * niuyuan(fastpow(i,*a[i].begin())) % mod;
                //去掉原有的
            a[i].erase(a[i].find(mp[k][i]));
            //cout << *a[i].find(mp[k][i]) << " " << mp[k][i] << endl;
            mp[k][i] += cnt;
            a[i].insert(mp[k][i]);
            if(a[i].size()==n)
                sum = sum *fastpow(i,*a[i].begin()) % mod;
        }else{
            a[i].insert(cnt);
            mp[k][i] = cnt;
            if(a[i].size()==n)
                sum = (sum * fastpow(i, *a[i].begin())) % mod;
        }
}
void divide(ll k,ll ans){ //b[k] 
    ll an = sqrt(ans);
    for (int i = 2; i <=an;i++){
        if(ans%i==0){
            int cnt = 0;
            while(ans%i==0){
                cnt++;
                ans /= i;
                //cout << ans << endl;
            }
            insert(k, i, cnt);
        }
    }
    if(ans>1)
        insert(k,ans, 1);
}
void solve(){
    n = read();
    m = read();
    ll ans;
    for (int i = 1; i <= n;i++){
        ans = read();
        divide(i, ans);
    }
    int x;
    for (int i = 1; i <= m;i++){
        x = read();
        ans = read();
        divide(x, ans);
        cout << sum << endl; 
    }
    //cout << sum << endl;
}
int main (){
    solve();
    
    getchar();
    getchar();
}