1. 首页
2. RMQ
3. 本页

HDU 5875 Function (RMQ+二分)

题目链接：点我~~

题意：给出a1~an,和 q 个询问 [i,j] 表示询问 aimod…aj 的值。

思路：每次找到下一个不超过当前余数的数。用RMQ记录区间最小的数。因为每次余数减少至少一半，复杂度应该是 O(q*logn*logAi).

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int, int> PI;
typedef pair< PI, int> PII;
const double eps=1e-5;
const double pi=acos(-1.0);
const int mod=1e9+7;
const int INF=0x3f3f3f3f;
#define lson l, m, rt<<1
#define rson m+1, r, rt<<1|1
const int MAXN = 100100;

int d[MAXN<<2];
int dp[MAXN][30];
int mm[MAXN];
int a[MAXN];

void intRMQ(int n,int b[])
{
    mm[0]=-1;
    for(int i=1; i<=n; ++i)
    {
        mm[i]=((i&(i-1))==0)?mm[i-1]+1:mm[i-1];
        dp[i][0]=b[i];
    }
    for(int j=1; j<=mm[n]; ++j)  //for(int j=1;(1<<j)<=n;++j)
        for(int i=1; i+(1<<j)-1<=n; ++i)
        {
            dp[i][j]=min(dp[i][j-1],dp[i+(1<<(j-1))][j-1]);
        }

}

int rmq(int x,int y)
{
    int k=mm[y-x+1];
    return min(dp[x][k],dp[y-(1<<k)+1][k]);
}

int binfind(int l,int r,int x)
{
    while(l<r)
    {
        int mid=l+((r-l)>>1);
        if(rmq(l,mid)<=x)
            r=mid;
        else
            l=mid+1;
    }
    if(a[l]<=x)
    {
        return l;
    }
    return -1;
}

int main()
{
    int n,t,q;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&n);
        for(int i=1; i<=n; i++)
            scanf("%d",&a[i]);
        intRMQ(n,a);
        scanf("%d",&q);
        for(int i=1; i<=q; i++)
        {
            int l,r;
            scanf("%d%d",&l,&r);
            int ans=a[l];
            while(l!=r)
            {
                int pos=binfind(l+1,r,ans);
                if(pos==-1)
                    break;
                l=pos;
                ans%=a[pos];
            }
            printf("%d\n",ans);
        }
    }
    return 0;

}