wannalfy 挑战赛7 E 珂朵莉与GCD (离线+线段树/树状数组)

链接：https://www.nowcoder.com/acm/contest/56/E

时间限制：C/C++ 5秒，其他语言10秒

空间限制：C/C++ 716800K，其他语言1433600K
64bit IO Format: %lld

题目描述

给你一个长为n的序列a

m次查询

每次查询一个区间的所有子区间的gcd的和mod1e9+7的结果

输入描述:

第一行两个数n,m
之后一行n个数表示a
之后m行每行两个数l,r表示查询的区间

输出描述:

对于每个询问，输出一行一个数表示答案

示例1

输入

5 7
30 60 20 20 20
1 1
1 5
2 4
3 4
3 5
2 5
2 3

输出

30
330
160
60
120
240
100

说明

[1,1]的子区间只有[1,1]，其gcd为30
[1,5]的子区间有：
[1,1]=30,[1,2]=30,[1,3]=10,[1,4]=10,[1,5]=10
[2,2]=60,[2,3]=20,[2,4]=20,[2,5]=20
[3,3]=20,[3,4]=20,[3,5]=20
[4,4]=20,[4,5]=20
[5,5]=20
总共330
[2,4]的子区间有：
[2,2]=60,[2,3]=20,[2,4]=20
[3,3]=20,[3,4]=20
[4,4]=20
总共160
[3,4]的子区间有：
[3,3]=20,[3,4]=20
[4,4]=20
总共60
[3,5]的子区间有：
[3,3]=20,[3,4]=20,[3,5]=20
[4,4]=20,[4,5]=20
[5,5]=20
总共120
[2,5]的子区间有：
[2,2]=60,[2,3]=20,[2,4]=20,[2,5]=20
[3,3]=20,[3,4]=20,[3,5]=20
[4,4]=20,[4,5]=20
[5,5]=20
总共240
[2,3]的子区间有：
[2,2]=60,[2,3]=20
[3,3]=20
总共100

备注:

对于100%的数据，有1 <= n , m , a_i <= 100000

///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////

//////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////

根据一个数的唯一分解，一个r往左边的gcd最多个分成个区间，相同区间内的点到r的gcd是一样的

对于询问（l , r) ，枚举子区间的右端点从l到r，如果子区间右端点为 l到r-1 都处理好了就只需要坐享前面更新的成果，然后更新r为右端点的区间就好了

把询问离线，，因此需要区间更新和区间求和，线段树的就不写了

新学树状数组的区间更新和区间求和，用差分的思想

贴了别人的blog，我也是看了别人的（侵删）http://blog.csdn.net/fsahfgsadhsakndas/article/details/52650026

#include <bits/stdc++.h>
#define mst(a,b)    memset((a),(b), sizeof a)
#define lowbit(a)   ((a)&(-a))
#define IOS         ios::sync_with_stdio(0);cin.tie(0);
#define MP          make_pair
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
const int mod=1e9+7;
const int maxn=1e5+10;
int n,m;
int a[maxn],nx[maxn];
int ans[maxn];
vector<pii>uu[maxn];
int c1[maxn],c2[maxn];
void update(int *bits,int pos,int val){
    for(;pos<=n;pos+=lowbit(pos))
        bits[pos]=(bits[pos]+val)%mod;
}
int get(int *bits,int pos){
    int ret=0;
    for(;pos;pos-=lowbit(pos))ret=(ret+bits[pos])%mod;
    return (ret+mod)%mod;
}
void add(int l,int r,int val){
    update(c1,l,val),update(c2,l,1LL*(l-1)*val%mod);

    update(c1,r+1,-val),update(c2,r+1,-1LL*r*val%mod);
}
int query(int pos){
    return (1LL*pos*get(c1,pos)-get(c2,pos)+mod)%mod;
}
int query(int l,int r){
    return (query(r)-query(l-1)+mod)%mod;
}
int main(){
#ifdef local
freopen("in.txt","r",stdin);
//freopen("out.txt","w",stdout);
#endif
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;++i)scanf("%d",&a[i]),nx[i]=i-1;
    for(int i=1;i<=m;++i){
        int l,r;scanf("%d%d",&l,&r);
        uu[r].push_back(MP(l,i));
    }
    for(int i=1;i<=n;++i){
        for(int j=i;j;j=nx[j])a[j]=__gcd(a[j],a[i]);

        for(int j=i;nx[j];j=nx[j])
            while(nx[j]&&a[j]==a[nx[j]])nx[j]=nx[nx[j]];

        for(int j=i;j;j=nx[j])add(nx[j]+1,j,a[j]);

        for(int j=0;j<uu[i].size();++j){
            pii&k = uu[i][j];
            ans[k.second]=query(k.first,i);
        }
    }
    for(int i=1;i<=m;++i)printf("%d\n",ans[i]);
    return 0;
}