hihocoder第16周（RMQ-ST算法）

题目链接：

huangjing

思路：新学会的一种算法，RMQ（Rangle Minimum Query）从名字来看，觉得就是查询最小值的，哈哈，

大白上有仔细的讲解。dp[i][j]=min(dp[i][j-1],dp[i+(1<<(j-1))][j-1])，具体的图请看大白，真是神奇的方法。。

预处理的复杂度为O(N*logN),查询就是O(1)的复杂度。

方法二：线段树解法  复杂度为O(Q*logN).

题目：

题目1 : RMQ-ST算法

时间限制: 10000ms

单点时限: 1000ms

内存限制: 256MB

描述

小Hi和小Ho在美国旅行了相当长的一段时间之后，终于准备要回国啦！而在回国之前，他们准备去超市采购一些当地特产——比如汉堡（大雾）之类的回国。

但等到了超市之后，小Hi和小Ho发现者超市拥有的商品种类实在太多了——他们实在看不过来了！于是小Hi决定向小Ho委派一个任务：假设整个货架上从左到右拜访了N种商品，并且依次标号为1到N，每次小Hi都给出一段区间[L, R]，小Ho要做的是选出标号在这个区间内的所有商品重量最轻的一种，并且告诉小Hi这个商品的重量，于是他们就可以毫不费劲的买上一大堆东西了——多么可悲的选择困难症患者。

（虽然说每次给出的区间仍然要小Hi来进行决定——但是小Hi最终机智的选择了使用随机数生成这些区间！但是为什么小Hi不直接使用随机数生成购物清单呢？——问那么多做什么！）

提示一：二分法是宇宙至强之法！（真的么？）

提示二：线段树不也是二分法么？

输入

每个测试点（输入文件）有且仅有一组测试数据。

每组测试数据的第1行为一个整数N，意义如前文所述。

每组测试数据的第2行为N个整数，分别描述每种商品的重量，其中第i个整数表示标号为i的商品的重量weight_i。

每组测试数据的第3行为一个整数Q，表示小Hi总共询问的次数。

每组测试数据的第N+4~N+Q+3行，每行分别描述一个询问，其中第N+i+3行为两个整数Li, Ri，表示小Hi询问的一个区间[Li, Ri]。

对于100%的数据，满足N<=10^6，Q<=10^6, 1<=Li<=Ri<=N，0<weight_i<=10^4。

输出

对于每组测试数据，对于每个小Hi的询问，按照在输入中出现的顺序，各输出一行，表示查询的结果：标号在区间[Li, Ri]中的所有商品中重量最轻的商品的重量。

样例输入

10
7334
1556
8286
1640
2699
4807
8068
981
4120
2179
5
3 4
2 8
2 4
6 8
7 10

样例输出

1640
981
1556
981
981

代码1：

#include<cstdio>
#include<cmath>
#include<algorithm>
#include<iostream>
using namespace std;

const int maxn=1000000+10;


int a[maxn],dp[maxn][20],n,q,ans[maxn];

void RMQ_init()
{
    for(int i=1;i<=n;i++)  dp[i][0]=a[i];
    for(int j=1;(1<<j)<=n;j++)
        for(int i=1;i+(1<<(j-1))-1<=n;i++)
           dp[i][j]=min(dp[i][j-1],dp[i+(1<<(j-1))][j-1]);
}

int RMQ_Query(int l,int r)
{
    int k=0;
    while(1<<(k+1)<=r-l+1)  k++;
    return min(dp[l][k],dp[r-(1<<k)+1][k]);
}

int main()
{
    int l,r,k;
    while(~scanf("%d",&n))
    {
        for(int i=1;i<=n;i++)
            scanf("%d",&a[i]);
        RMQ_init();
        scanf("%d",&q);
        for(int i=1;i<=q;i++)
        {
            scanf("%d%d",&l,&r);
            ans[i]=RMQ_Query(l,r);
        }
        for(int i=1;i<=q;i++)
          printf("%d\n",ans[i]);
    }
    return 0;
}

代码2：

#include<cstdio>
#include<iostream>
#include<algorithm>
using namespace std;

const int maxn=1000000+10;

int tree[maxn<<2],n,q,ans[maxn];

void push_up(int dex)
{
    tree[dex]=min(tree[dex<<1],tree[dex<<1|1]);
}

void buildtree(int dex,int l,int r)
{
    if(l==r)
    {
        scanf("%d",&tree[dex]);
        return;
    }
    int mid=(l+r)>>1;
    buildtree(dex<<1,l,mid);
    buildtree(dex<<1|1,mid+1,r);
    push_up(dex);
}

int Query(int dex,int L,int R,int l,int r)
{
    if(l<=L&&r>=R)
        return tree[dex];
    int mid=(L+R)>>1;
    if(l>mid)  return Query(dex<<1|1,mid+1,R,l,r);
    else if(r<=mid)  return Query(dex<<1,L,mid,l,r);
    else return min(Query(dex<<1,L,mid,l,r),Query(dex<<1|1,mid+1,R,l,r));
}

int main()
{
    int l,r;
    while(~scanf("%d",&n))
    {
         buildtree(1,1,n);
         scanf("%d",&q);
         for(int i=1;i<=q;i++)
         {
            scanf("%d%d",&l,&r);
            ans[i]=Query(1,1,n,l,r);
         }
        for(int i=1;i<=q;i++)
            printf("%d\n",ans[i]);
    }
    return 0;
}