HDU - 5381 - The sum of gcd
左端点相同的区间最多有
log(max(ai))
个
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;
}