简单的RMQ,可我怎么写都WA。不明白,找了一个和我相似的贴过了,要赶着去外婆家。

#include <iostream>
#include <algorithm>
#include <cstring>
#include <string>
#include <cstdio>
#include <cmath>
#include <queue>
#include <vector>
#include <map>
#include <set>
#include <stack>
#define eps 1e-5
#define MAXN 255
#define MAXM 111111
#define INF 1000000000
using namespace std;
int mx[MAXN][MAXN][8], mi[MAXN][MAXN][8];
int n, b, q, a, c;
void rmqinit()
{
int l = int(log(double(n)) / log(2.0)) ;
for(int k = 1; k <= n ; k++)
for(int j = 1; j <= l; j++)
for(int i = 1; i + (1 << (j - 1))- 1 <= n; i++)
{
mx[k][i][j] = max(mx[k][i][j - 1], mx[k][i + (1 << (j - 1 ))][j - 1]) ;
mi[k][i][j] = min(mi[k][i][j - 1], mi[k][i + (1 << (j - 1 ))][j - 1]) ;
}
}
int rmqmax(int lx, int ly, int rx, int ry) // lx, ly为左上角的点 rx ry为右下角的点
{
int l = int(log(double(ry - ly + 1)) / log(2.0));
int ret = -1;
for(int k = lx; k <= rx ; k++)
ret = max(ret, max(mx[k][ly][l], mx[k][ry - (1 << l) + 1][l]));
return ret;
}
int rmqmin(int lx, int ly, int rx, int ry) // lx, ly为左上角的点 rx ry为右下角的点
{
int l = int(log(double(ry - ly + 1)) / log(2.0));
int ret = INF;
for(int k = lx; k <= rx ; k++)
ret = min(ret, min(mi[k][ly][l], mi[k][ry - (1 << l) + 1][l]));
return ret;
}
int main()
{
while(scanf("%d%d%d", &n, &b, &q) != EOF)
{
for(int i = 1; i <= n; i++)
for(int j = 1; j <= n; j++)
{
scanf("%d",&a);
mx[i][j][0] = mi[i][j][0] = a ;
}
rmqinit();
while(q--)
{
scanf("%d%d", &a, &c) ;
int rx = a + b - 1;
if(rx > n) rx = n;
int ry = c + b - 1;
if(ry > n) ry = n;
printf("%d\n", rmqmax(a, c, rx, ry) - rmqmin(a, c, rx, ry)) ;
}
}
return 0;
}

  

MINE:

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <climits>
#include <string.h>
#include <queue>
#include <cmath>
#include <vector>
using namespace std; int num[255][255];
int f1[255][255][30];
int f2[255][255][30];
int n,b,q;
const int inf=1000000000; int rmq_max(int p,int i, int j) {
int k = (int)(log(double(j-i+1)) / log(2.0)), t1;
t1 = max(f1[p][i][k], f1[p][j - (1<<k) + 1][k]);
return t1;
}
int rmq_min(int p,int i, int j) {
int k = (int)(log(double(j-i+1)) / log(2.0)), t2;
t2 = min(f2[p][i][k], f2[p][j - (1<<k) + 1][k]);
return t2;
} int main(){
int l,u;
while(scanf("%d%d%d",&n,&b,&q)!=EOF){
int k = (int) (log((double)n) / log(2.0));
for(int i=0;i<n;i++){
for(int j=0;j<n;j++)
scanf("%d",&num[i][j]);
for(int j = 0; j < n; j++) {
f1[i][j][0] = num[i][j];
f2[i][j][0] = num[i][j];
}
for(int p = 1; p <= k; p++) {
for(int t = 0; t + (1 << p) - 1 < n; t++) {
int m = t + (1 << (p - 1));
f1[i][t][p] = max(f1[i][t][p-1], f1[i][m][p-1]);
f2[i][t][p] = min(f2[i][t][p-1], f2[i][m][p-1]);
}
}
}
int maxn=-1,minn=inf;
for(int i=1;i<=q;i++){
scanf("%d%d",&u,&l);
l--;u--;
for(int p=u;p<(u+b);p++){
maxn=max(maxn,rmq_max(p,l,l+b-1));
minn=min(minn,rmq_min(p,l,l+b-1));
}
printf("%d\n",maxn-minn);
}
}
return 0;
}

  

POJ 2019的更多相关文章

  1. POJ 2019 Cornfields [二维RMQ]

    题目传送门 Cornfields Time Limit: 1000MS   Memory Limit: 30000K Total Submissions: 7963   Accepted: 3822 ...

  2. poj 2019 二维rmq *

    题目大意:给出一个N*N矩形,每个格子上有一个价值.询问一个b*b的矩形在左上角的位置(x,y),(x+b-1,y+b-1)这一部分的最大值-最小值是多少. 模板题 #include <stdi ...

  3. [POJ 2019] Cornfields

    Cornfields Time Limit: 1000MS   Memory Limit: 30000K Total Submissions: 5516   Accepted: 2714 Descri ...

  4. 二维 ST POJ 2019

    题目大意:给你一个n*n的矩阵,每次给你一个点(x,y),以其为左上角,宽度为b的矩阵中最小的数值和最大数值的差是多少?  一共k个询问. 思路:简单的二维st. 定义dp(i,j,k,L)表示以(i ...

  5. POJ 2019 Cornfields(二维RMQ)

    相比以前的RMQ不同的是,这是一个二维的ST算法 #include<iostream> #include<cstring> #include<cstdio> #in ...

  6. Cornfields POJ - 2019(二维RMQ板题)

    就是求子矩阵中最大值与最小值的差... 板子都套不对的人.... #include <iostream> #include <cstdio> #include <sstr ...

  7. POJ 2019 Cornfields (二维RMQ)

    Cornfields Time Limit: 1000MS   Memory Limit: 30000K Total Submissions: 4911   Accepted: 2392 Descri ...

  8. POJ 2019 Cornfields 二维线段树的初始化与最值查询

    模板到不行.. 连更新都没有.. .存个模板. 理解留到小结的时候再写. #include <algorithm> #include <iostream> #include & ...

  9. ACM第一阶段学习内容

    一.知识目录 字符串处理 ................................................................. 3 1.KMP 算法 .......... ...

随机推荐

  1. 最长公共子序列 nlogn

    先来个板子 #include<bits/stdc++.h> using namespace std; , M = 1e6+, mod = 1e9+, inf = 1e9+; typedef ...

  2. bzoj1081: [SCOI2005]超级格雷码(dfs)

    1081: [SCOI2005]超级格雷码 题目:传送门 题解: 又是一道水题... 因为之前做过所以知道规律: 如n=2 B=3: 00 10 20    21 11 01    02 12 22 ...

  3. oracle中关于删除表purge语句和闪回语句的基本使用

    语法: drop table ... purge; 例子:drop table test purge; purge是直接删除表,不保留到回收站,10G开始默认drop表式改名移动到回收站; 闪回(fl ...

  4. k8s 安装并试用Istio service mesh

    本文根据官网的文档整理而成,步骤包括安装istio 0.5.1并创建一个bookinfo的微服务来测试istio的功能. 文中使用的yaml文件可以在kubernetes-handbook的manif ...

  5. Hadoop MapReduce编程 API入门系列之自定义多种输入格式数据类型和排序多种输出格式(十一)

    推荐 MapReduce分析明星微博数据 http://git.oschina.net/ljc520313/codeexample/tree/master/bigdata/hadoop/mapredu ...

  6. 使用Android ADT最新开发工具后,新建项目出现appcompat v7 他是什么?

    做Android开发的朋友最近会发现,更新ADT至22.6.0版本之后,创建新的安装项目,会出现appcompat_v7的内容.并且是创建一个新的内容就会出现.这到底是怎么回事呢?原来appcompa ...

  7. Java电子书下载地址

    http://www.itpub.net/search.php?searchid=1660&orderby=lastpost&ascdesc=desc&searchsubmit ...

  8. Android java处理保留小数点后几位

    方式一: 四舍五入  double   f   =   111231.5585;  BigDecimal   b   =   new   BigDecimal(f);  double   f1   = ...

  9. 文件类型总结 MIME

    来源网上https://www.cnblogs.com/zhongcj/archive/2008/11/03/1325293.html {".3gp", "video/3 ...

  10. trigger事件就是继承某一个类的事件.

    <html><head><script type="text/javascript" src="/jquery/jquery.js" ...