HDU 5726

RMQ应用题。

引用别人的解题思路：

#include<cstdio>
#include<cmath>
#include<algorithm>
#include<map>
#define MAXN 100000+5
using namespace std;
map<int,long long> m;
int n,a[MAXN],q;
int dp[MAXN][17];
int gcd(int a,int b){return b==0?a:gcd(b,a%b);}
void RMQ_init()
{
    for(int i=1;i<=n;i++) dp[i][0]=a[i];
    double j_max=(log(n)/log(2.0));
    for(int j=1;j<=j_max;j++)
    {
        for(int i=1;i<=n;i++)
            if(i+(1 << j)-1 <= n) dp[i][j]=gcd(dp[i][j-1],dp[i+(1<<(j-1))][j-1]);
    }
}
int query(int l,int r)
{
    int k=log(r-l+1)/log(2.0);
    return gcd(dp[l][k],dp[r-(1<<k)+1][k]);
}
void pretreat()//这个预处理是需要好好理解的地方
{
    m.clear();
    for(int l=1;l<=n;l++)
    {
        int r=l;
        while(r <= n)
        {
            int d=query(l,r);
            int st=r,ed=n,mid;
            while(ed-st>1)
            {  
                mid=(st+ed)>>1;
                if(query(l,mid) == d) st=mid;
                else ed=mid;
            }
            m[d]+=st-r+1;
            r=st+1;
        }
    }  
}
int main()
{
    int t;
    scanf("%d",&t);
    for(int kase=1;kase<=t;kase++)
    {
        scanf("%d",&n);
        for(int i=1;i<=n;i++) scanf("%d",&a[i]);
        RMQ_init();
        pretreat();
        scanf("%d",&q);
        int l,r;
        printf("Case #%d:
",kase);
        for(int i=1;i<=q;i++)
        {
            scanf("%d%d",&l,&r);
            int ans=query(l,r);
            printf("%d %I64d
",ans,m[ans]);
        }
    }
}