HDU - 5726 ​​​​​​​RMQ+二分

GCD

 HDU - 5726 

题目大意：给出了n个数，m次查询，每次查询gcd（l,r）并且输出有多少个区间的gcd==gcd(l,r)；

解题思路：做的时候误以为要求查询的区间在（l,r）内，并且也没有意识到固定左端点移动右端点gcd是单调不增的。

我们可以固定左端点，然后根据上边那个性质，二分右端点。然后统计贡献。

#include<bits/stdc++.h>
using namespace std;
#define LL long long

const int N = 1e5+5;
map<int,LL>mmp;
int a[N];
int f[N][20];
void build(int n)
{
    for(int i=1;i<=n;i++)f[i][0]=a[i];

    for(int i=1;i<=19;i++)
    {
        for(int j=1;j+(1<<i)-1<=n;j++)
        {
            f[j][i]=__gcd(f[j][i-1],f[j+(1<<(i-1))][i-1]);
        }
    }
}

int ask(int l,int r)
{
    int k=log2(r-l+1);
    return __gcd(f[l][k],f[r-(1<<k)+1][k]);
}

void solve(int n)
{
    mmp.clear();
    for(int i=1;i<=n;i++)
    {
        int now=f[i][0];
        int L=i,R=i;
        while(L<=n)
        {
            int l=L,r=n;
            int mid;
            while(l<=r)
            {
                int m=(l+r)>>1;
                if(ask(i,m)==now)l=m+1,mid=m;
                else r=m-1;
            }
            mmp[now]+=mid-L+1;
            L=mid+1;
            now=ask(i,L);
        }
    }
}
int main()
{
    int t;
    cin>>t;
    int cas=1;
    while(t--)
    {
        int n;
        scanf("%d",&n);
        for(int i=1;i<=n;i++)scanf("%d",&a[i]);
        build(n);
        solve(n);
        int m;
        scanf("%d",&m);
        printf("Case #%d:\n",cas++);
        for(int i=1;i<=m;i++)
        {
            int x,y;
            scanf("%d%d",&x,&y);
            int tmp=ask(x,y);
            printf("%d %lld\n",tmp,mmp[tmp]);
        }
    }
}