HDU 5726 GCD 倍增

题目大意，给一列数，多组询问[L,R]，求AL到AR的GCD，以及[1,n]中有多少个子区间的gcd=AL到AR的GCD。
题解：考虑暴力，用线段树维护第一个询问O(nlgn)，第二个询问怎么做？
一个性质是，一个数字如果变成自己的因数，并且严格小于，那么就会变小一半，涨姿势了。
这类性质还有取模，开根，等等，每次成功操作至少减少一半。
其实想了想树剖点分等本质上也是干这种事情，只不过这种性质是人为构造出来的。
有些这类性质需要观察题目。
因此枚举左端点，考虑随着右端点的递增其gcd只会有O(lgn)种，因此计算每个gcd对应的右端点所在的区间，可以用二分做。一个优化是求gcd的时候可以用ST表维护，这是gcd的另一个性质。（类似于max/min，重叠计算不影响答案，因此可以这么搞）。这样最后用一个map存一下每一种gcd的答案即可。
复杂度O(nlg^3n)。
代码：

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<map>
#define LG 20
#define N 100010
#define lint long long
using namespace std;
int f[N][LG],n,lg[N];
map<int,lint> ans;
inline int gcd(int x,int y)
{
    return x?gcd(y%x,x):y;
}
inline int query(int s,int t)
{
    int k=lg[t-s+1];
    return gcd(f[s][k],f[t-(1<<k)+1][k]);
}
inline int solve(int s,int L,int R,int g)
{
    int mid=(L+R)>>1;
    while(L<=R)
    {
        if(query(s,mid)==g) L=mid+1;
        else R=mid-1;mid=(L+R)>>1;
    }
    return R;
}
int main()
{
    int T,Tc=0;scanf("%d",&T);
    for(int i=1;i<N;i++)
    {
        lg[i]=lg[i-1];
        while((1<<(lg[i]+1))<=i) lg[i]++;
    }
    while(T--)
    {
        scanf("%d",&n),ans.clear();
        for(int i=1;i<=n;i++) scanf("%d",&f[i][0]);
        for(int j=1;(1<<j)<=n;j++)
            for(int i=1;i+(1<<j)-1<=n;i++)
                f[i][j]=gcd(f[i][j-1],f[i+(1<<(j-1))][j-1]);
        for(int i=1,j=1,k,g=f[1][0];i<=n;g=f[j=++i][0])
            while(j<=n) k=solve(i,j,n,g),ans[g]+=k-j+1,g=query(i,j=k+1);
        int m;scanf("%d",&m);printf("Case #%d:\n",++Tc);
        while(m--)
        {
            int s,t,g;scanf("%d%d",&s,&t),g=query(s,t);
            printf("%d %lld\n",g,ans[g]);
        }
    }
    return 0;
}