HDU - 5381 - The sum of gcd（莫队算法）

HDU - 5381 - The sum of gcd
左端点相同的区间最多有 $log(max(a_i))$ 个 $gcd$ ，且随着右端点的变大而变小。右端点相同的区间同理。
于是考虑 $O(nlogn)$ 预处理出来。
然后用莫队维护区间 $gcd$ 和就行了。

#include<bits/stdc++.h>
#define fi first
#define se second
#define mp make_pair
using namespace std;
typedef long long ll;
const int N=1e4+7;
typedef pair<int,int> pii;
vector<pii> vl[N],vr[N];
int a[N],pos[N],L,R,n;
ll ans[N],Ans;
struct Query
{
    int l,r,id;
    bool operator < (const Query & t) const
    {
        return pos[l]<pos[t.l] || (pos[l]==pos[t.l]&&r<t.r);
    }
}Q[N];
int gcd(int a,int b){ return a%b?gcd(b,a%b):b; }
//fi 值,se 下标
void init()
{
    for(int i=0;i<=n+1;++i) vr[i].clear(),vl[i].clear();
    for(int i=1;i<=n;++i)
    {
        int p=a[i];
        vr[i].push_back(mp(p,i));
        for(int j=0;j<vr[i-1].size();++j)
        {
            p=gcd(p,vr[i-1][j].fi);
            if(p!=vr[i].back().fi) vr[i].push_back(mp(p,vr[i-1][j].se));
            else vr[i].back().se=vr[i-1][j].se;
        }
    }
    for(int i=n;i>=1;--i)
    {
        int p=a[i];
        vl[i].push_back(mp(p,i));
        for(int j=0;j<vl[i+1].size();++j)
        {
            p=gcd(p,vl[i+1][j].fi);
            if(p!=vl[i].back().fi) vl[i].push_back(mp(p,vl[i+1][j].se));
            else vl[i].back().se=vl[i+1][j].se;
        }
    }
}
void addL(int L)
{
    ll res=0,cur=L;
    for(int i=0;i<vl[L].size();++i)
    {
        res+=(ll)vl[L][i].fi*(min(vl[L][i].se,R)-cur+1);
        cur=vl[L][i].se+1;
        if(cur>R) break;
    }
    Ans+=res;
}
void addR(int R)
{
    ll res=0,cur=R;
    for(int i=0;i<vr[R].size();++i)
    {
        res+=(ll)vr[R][i].fi*(cur-max(vr[R][i].se,L)+1);
        cur=vr[R][i].se-1;
        if(cur<L) break;
    }
    Ans+=res;
}
void delL(int L)
{
    ll res=0,cur=L;
    for(int i=0;i<vl[L].size();++i)
    {
        res+=(ll)vl[L][i].fi*(min(vl[L][i].se,R)-cur+1);
        cur=vl[L][i].se+1;
        if(cur>R) break;
    }
    Ans-=res;
}
void delR(int R)
{
    ll res=0,cur=R;
    for(int i=0;i<vr[R].size();++i)
    {
        res+=(ll)vr[R][i].fi*(cur-max(vr[R][i].se,L)+1);
        cur=vr[R][i].se-1;
        if(cur<L) break;
    }
    Ans-=res;
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        int q;
        scanf("%d",&n);
        for(int i=1;i<=n;++i) scanf("%d",&a[i]);
        init();
        scanf("%d",&q);
        int sz=sqrt(n);
        for(int i=1;i<=n;++i)
            pos[i]=(i-1)/sz;
        for(int i=0;i<q;++i)
            scanf("%d%d",&Q[i].l,&Q[i].r),Q[i].id=i;
        sort(Q,Q+q);
        L=1;R=0;Ans=0;
        for(int i=0;i<q;++i)
        {
            while(L>Q[i].l) addL(--L);
            while(R<Q[i].r) addR(++R);
            while(L<Q[i].l) delL(L++);
            while(R>Q[i].r) delR(R--);
            ans[Q[i].id]=Ans;
        }
        for(int i=0;i<q;++i) printf("%I64d\n",ans[i]);
    }
    return 0;
}