2016 多校联赛 第一场 HDU 5726 GCD

题意：

n个数给出q 次询问，每个询问有多少个子区间 与给定的（l ,r）区间的GCD相同；

很容易想到类似RMQ 的预处理，或者线段树处理出区间的GCD; 复杂度 O(N*log (n);  每次询问一个区间复杂度 logn,  

接着也很容易想到枚举每个区间的左端点很右端点，这样的复杂度是  N*N*log n   显然这样是过不了的，

这里要知道一个gcd特点，就是以区间起点为起点的gcd 是由大到小，切成阶梯状， 例如；

１８　１２　２４　６　４　２　３　

（１　２）ｇｃｄ＝＝　６；

（１　３）ｇｃｄ＝＝６

（１　４）ｇｃｄ＝＝６

（１　４）ｇｃｄ＝＝２

（１　５）ｇｃｄ＝＝２；

（１　６）ｇｃｄ　＝１；

还有一个特点就是　如果　（ＸＹ）ｇｃｄ＝＝ｗ＆＆（Ｘ，Ｚ）ｇｃｄ＝＝ｗ，（Ｚ＞＝Ｙ）　那么

（Ｘ，Ｙ＋１）（Ｘ，Ｙ＋２），．．．．（Ｘ，Ｚ）　都＝＝ｗ；　

这样的话我们就可以枚举左端点，找到最长的右端点，由于ＧＣＤ递减，可以用二分查找。

找到后，那么左端点到右端点之前的区间ｇｃｄ　＝＝ｗ，继续向后找，左端点不动，右端点后移一个数，

 相当于找每个阶梯的起点和二分终点 ，

　用线段树写了超时，

/*
#include<bits/stdc++.h>
using namespace std;
const int maxn = 100010;
int g[maxn<<2];
int a[maxn];
int ans[maxn];
map<int,long long>m;
int  gcd(int  x,int y)
{
    //  printf("y==%d",y);
    return x%y==0?y:gcd(y,x%y);
}
void pushup(int rt)
{
    g[rt] = gcd(g[rt<<1],g[rt<<1|1]);
}
void build(int l,int r,int rt)
{
    if(l==r)
    {
        g[rt]=a[l];
        return;
    }
    int mid = (l+r)>>1;
    build(l,mid,rt<<1);
    build(mid+1,r,rt<<1|1);
    pushup(rt);
}
int query(int L ,int R, int l,int r,int rt)
{
    if(L<=l&&r<=R)
    {
        return g[rt];
    }
    int mid = (l+r)>>1;
    int x  = 0;
    int  y = 0;
    if(L<=mid)
        x = query(L,R,l,mid,rt<<1);
    if(R>mid)
        y = query(L,R,mid+1,r,rt<<1|1);
    //  printf("%d %d  \n",x,y);
    if(x==0||y==0)return max(x,y);
    return gcd(x,y);
}
int main()
{
    int T;
    scanf("%d",&T);
    int F = 1;
    while(T--)
    {
        int n;
        scanf("%d",&n);
        for(int i = 1; i<=n; i++)
        {
            scanf("%d",&a[i]);
        }
        build(1,n,1);
        int l,r;
        int q;
        scanf("%d",&q);
        for(int i = 1; i<=q; i++)
        {
            scanf("%d%d",&l,&r);
            ans[i] = query(l,r,1,n,1);
            m.insert(make_pair(ans[i],0));
        }
        for(int i = 1; i<=n; i++)
        {
            int beg= i;
            for(int j =beg; j<=n;)
            {
                int  val = query(beg,j,1,n,1);
                l = j;
                r = n;

                while(l<r)
                {// printf("l-----%d r------%d\n",l,r);
                    int mid = (l+r+1)>>1;
                    if(query(beg,mid,1,n,1)<val)
                    {
                        r = mid-1;
                    }
                    else
                    {
                        l = mid;
                    }
                }
               if(m.find(ans[i])!=m.end())
                m[val]+=(l- j+1);
                j = r+1;
            }
        }
        printf("Case #%d:\n",F++);
        for(int i = 1; i<=q; i++)
        {
            printf("%d %lld\n",ans[i],m[ans[i]]) ;
        }
    }
    return 0;
}*/

#include<bits/stdc++.h>
using namespace std;
const int maxn = 100010;
int g[maxn<<2];
int a[maxn];
int ans[maxn];
map<int,long long>m;
int  gcd(int  x,int y)
{
    //  printf("y==%d",y);
    return x%y==0?y:gcd(y,x%y);
}
int maxrmq[100009][20];
/*
void pushup(int rt)
{
    g[rt] = gcd(g[rt<<1],g[rt<<1|1]);
}
void build(int l,int r,int rt)
{
    if(l==r)
    {
        g[rt]=a[l];
        return;
    }
    int mid = (l+r)>>1;
    build(l,mid,rt<<1);
    build(mid+1,r,rt<<1|1);
    pushup(rt);
}
int query(int L ,int R, int l,int r,int rt)
{
    if(L<=l&&r<=R)
    {
        return g[rt];
    }
    int mid = (l+r)>>1;
    int x  = 0;
    int  y = 0;
    if(L<=mid)
        x = query(L,R,l,mid,rt<<1);
    if(R>mid)
        y = query(L,R,mid+1,r,rt<<1|1);
    //  printf("%d %d  \n",x,y);
    if(x==0||y==0)return max(x,y);
    return gcd(x,y);
}
*/

void RMQ(int num)
{
    for(int j = 1; j < 20; ++j)
        for(int i = 1; i <= num; ++i)
            if(i + (1 << j) - 1 <= num)
            {
                maxrmq[i][j] = gcd(maxrmq[i][j - 1], maxrmq[i + (1 << (j - 1))][j - 1]);
            }
}
int query(int start,int end )
{
    int k=log10(end-start+1.0)/log10(2.0);
    int maxdata=gcd(maxrmq[start][k],maxrmq[end-(1<<k)+1][k]);
    return maxdata;
}
int main()
{
    int T;
    scanf("%d",&T);
    int F = 1;
    while(T--)
    {
        m.clear();
        int n;
        scanf("%d",&n);
      //  printf("%d\n",maxrmq[0][0]);
        for(int i = 1; i<=n; i++)
        {
            scanf("%d",&a[i]);
          maxrmq[i][0] = a[i];

        }
        RMQ(n);
        int l,r;
        int q;
        scanf("%d",&q);
        for(int i = 1; i<=q; i++)
        {
            scanf("%d%d",&l,&r);
            ans[i] = query(l,r);
             m.insert(make_pair(ans[i],0));
        }
        for(int i = 1; i<=n; i++)
        {
            int beg= i;
            for(int j =beg; j<=n;)
            {
                int  val = query(beg,j);
                l = j;
                r = n;
                while(l<r)
                {
                    int mid = (l+r+1)>>1;
                    if(query(beg,mid)<val)
                    {
                        r = mid-1;
                    }
                    else
                    {
                        l = mid;
                    }
                } if(m.find(val)!=m.end())
                    m[val]+=(l- j+1);
                j = r+1;
            }
        }
        printf("Case #%d:\n",F++);
        for(int i = 1; i<=q; i++)
        {
            printf("%d %lld\n",ans[i],m[ans[i]]) ;
        }
    }
    return 0;
}