HDU5869 Different GCD Subarray Query(离线线段树)

题意：

找l~r区间中各个字串构成的gcd ，有多少种

思路：

可以分离处区间中有多少不同数字，这可以用离线的线段树方式处理。也就是按r端点排序，从左向右，使得值尽可能的往右靠。

处理gcd的方法可以从左向右推，区间gcd的收敛程度很高，暴力即可。

还有一个问题是在于保存gcd的时候，一定要使其（是从哪个l 推出来的， 这个l尽可能的靠近r）只有这样才能保证在离线处理时，不出错!

#include <iostream>
#include <stdio.h>
#include <vector>
#include <algorithm>
#include <cstring>
using namespace std;
const int maxn=100005;
struct node
{
    int num;
}tree[maxn*4];
int n,m;
struct q
{
    int l,r,id;
    bool operator < (const struct q &rhs) const
    {
        if(r!=rhs.r)
            return r<rhs.r;
        return l<rhs.l;
    }
}ask[maxn];
int has[maxn*10];
vector<pair<int,int> > vec[maxn];
void build(int i,int l,int r)
{
    tree[i].num=0;
    if(l==r)
        return ;
    int mid=(l+r)>>1;
    build(i<<1,l,mid);
    build(i<<1|1,mid+1,r);
}
int gcd(int x,int y)
{
    if(y==0)
        return x;
    else
        return gcd(y,x%y);
}
void update(int i,int l,int r,int pos,int w)
{
    tree[i].num+=w;
    if(l==r)
    {
        return ;
    }
    int mid=(l+r)>>1;
    if(pos<=mid)
        update(i<<1,l,mid,pos,w);
    else
        update(i<<1|1,mid+1,r,pos,w);
}
int query(int i,int l,int r,int L,int R)
{
    if(L<=l&&r<=R)
        return tree[i].num;
    int mid=(l+r)>>1;
    int ans=0;
    if(L<=mid)
        ans+=query(i<<1,l,mid,L,R);
    if(R>mid)
        ans+=query(i<<1|1,mid+1,r,L,R);
    return ans;
}
int ans[maxn];
int main()
{
    while(~scanf("%d%d",&n,&m))
    {
        for(int i=0;i<=n;i++)
            vec[i].clear();
        for(int i=1;i<=n;i++)
        {
            int x;
            scanf("%d",&x);
            int y=i;
            int t;
            for(int j=0;j<vec[i-1].size();j++)
            {
                t = gcd(vec[i-1][j].first,x);
                if(t!=x)
                {
                    vec[i].push_back(make_pair(x,y));
                    x=t; y=vec[i-1][j].second;
                }
            }
            vec[i].push_back(make_pair(x,y));
        }
/*
        for(int i=1;i<=n;i++)
        {
            for(int j=0;j<vec[i].size();j++)
                printf(" --%d %d-- ",vec[i][j].first,vec[i][j].second);
            printf("\n");
        }*/
        build(1,1,n);
        for(int i=1;i<=m;i++)
        {
            scanf("%d%d",&ask[i].l,&ask[i].r);
            ask[i].id=i;
        }
        sort(ask+1,ask+1+m);
        int cur=1;
        memset(has,0,sizeof(has));
        for(int i=1;i<=n;i++)
        {
            for(int j=0;j<vec[i].size();j++)
            {   
                int u=vec[i][j].first,v=vec[i][j].second;
                if(has[vec[i][j].first]<vec[i][j].second)
                {
                    if(has[vec[i][j].first])
                        update(1,1,n,has[vec[i][j].first],-1);
                    update(1,1,n,vec[i][j].second,1);
                    has[vec[i][j].first]=vec[i][j].second;
                }
            } 
            for(int j=cur;ask[j].r<=i&&j<=m;j++)
            {
                ans[ask[j].id]=query(1,1,n,ask[j].l,ask[j].r);
                cur=j+1;
            }
        }
        for(int i=1;i<=m;i++)
        {
            printf("%d\n",ans[i]);
        }
    }
    return 0;
}
/*
5 3
1 3 4 6 9
1 1
2 4
3 3
*/