2021-05-31

珂朵莉与GCD

题目链接

题意：给你n个元素，m次询问，每次询问[l,r]这个区间所有的子区间的gcd的和%p

n<=1e5

固定右端点r，如果[1,r-1]的子区间gcd和都求出来了，那么对于r这个点，将其贡献加入到[1,r-1]这个区间里就是[1,r]的所有子区间gcd的和。循环r的时候不要忘记更新右端点为r的询问。
我们知道对于右端点r来说，其往左看的不同gcd的个数最多为$lo{g}_{2}a[i]$个，因为最多情况是将a[i]看成2的倍数这样分解。所以我们可以对其暴力求gcd的和，用区间加和区间求和来优化。线段树来维护[l,r]这个区间所有的子区间的gcd和，更新只要将加入a[r]后一样的gcd压缩就行，用区间加将其加到区间上。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=4e5+10;
const int p=1e9+7;
ll tree[maxn],lazy[maxn];
void pushup(int x)
{
    tree[x]=(tree[x<<1]+tree[x<<1|1])%p;
}
void pushdown(int l,int r,int x)
{
    int mid=l+r>>1;
    tree[x<<1]=(tree[x<<1]+lazy[x]*(mid-l+1)%p)%p;
    tree[x<<1|1]=(tree[x<<1|1]+lazy[x]*(r-mid)%p)%p;
    lazy[x<<1]=(lazy[x<<1]+lazy[x])%p;lazy[x<<1|1]=(lazy[x<<1|1]+lazy[x])%p;
    lazy[x]=0;
}
void modify(int l,int r,int x,int st,int ed,ll c)
{
    if(st<=l&&r<=ed) {
        tree[x]+=(r-l+1)*c;
        lazy[x]+=c;
        return ;
    }
    int mid=l+r>>1;
    pushdown(l,r,x);
    if(mid>=st) modify(l,mid,x<<1,st,ed,c);
    if(mid<ed) modify(mid+1,r,x<<1|1,st,ed,c);
    pushup(x);
}
ll query(int l,int r,int x,int st,int ed)
{
    if(st<=l&&r<=ed) {
        return tree[x];
    }
    int mid=l+r>>1;
    pushdown(l,r,x);
    ll ans=0;
    if(mid>=st) ans=query(l,mid,x<<1,st,ed);
    if(mid<ed) ans=(ans+query(mid+1,r,x<<1|1,st,ed))%p;  //不%全wa
    return ans;
}
ll gcd(ll a,ll b)
{
    if(!b) return a;
    return gcd(b,a%b);
}
ll a[maxn],ans[maxn],pre[maxn];
pair<int,int>b[maxn];
vector<pair<int,int>>v[maxn];
int main()
{
    ios_base::sync_with_stdio(false);
    cin.tie(nullptr);
    cout.tie(nullptr);
    int n,m;cin>>n>>m;
    for(int i=1;i<=n;i++)
    cin>>a[i],pre[i]=i-1;     //pre记录gcd(pre[i]+1,a[i])!=gcd(pre[i],a[i])
    for(int i=1;i<=m;i++)
    {
        cin>>b[i].first>>b[i].second;
        v[b[i].second].push_back({b[i].first,i});
    }
    for(int r=1;r<=n;r++)
    {
        int i=r;
        while(i) {
            a[i]=gcd(a[i],a[r]);
            i=pre[i];
        }
        i=r;
        while(pre[i]) {
            while(pre[i]&&a[i]==a[pre[i]]) 
            pre[i]=pre[pre[i]];        //将相等区间合并
            i=pre[i];
        }
        i=r;
        while(i) {
            modify(1,n,1,pre[i]+1,i,a[i]);
            i=pre[i];
        }
        for(auto i:v[r]) {
            ans[i.second]=query(1,n,1,i.first,r);   //[l,r]
        }
    }
    for(int i=1;i<=m;i++)
    cout<<ans[i]<<endl;
    return 0;
}