简单的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. Hadoop-2.2.0中文文档——Apache Hadoop 下一代 MapReduce (YARN)

    MapReduce在hadoop-0.23中已经经历了一次全然彻底地大修.就是如今我们叫的MapReduce 2.0 (MRv2) or YARN. MRv2的基本思想是把JobTracker分成两个 ...

  2. php简单表格函数

    php简单表格函数 代码 <?php //ctrl+shift+j /** * @param unknown $rows * @param unknown $cols * @param stri ...

  3. poj 2031--Building a Space Station(prim)

    Building a Space Station Time Limit: 1000MS   Memory Limit: 30000K Total Submissions: 6064   Accepte ...

  4. Head First 设计模式 —— 装饰器模式

    1. 装饰器模式与继承 与装饰器模式相比,继承更容易造成类爆炸: 装饰器模式:利用组合取代继承:

  5. JavaScript:Browser 对象

    ylbtech-JavaScript:Browser 对象 1.  Window 对象返回顶部 1. Window 对象 Window 对象 Window 对象表示浏览器中打开的窗口. 如果文档包含框 ...

  6. Node.js:NPM 使用介绍

    ylbtech-Node.js:NPM 使用介绍 1.返回顶部 1. NPM 使用介绍 NPM是随同NodeJS一起安装的包管理工具,能解决NodeJS代码部署上的很多问题,常见的使用场景有以下几种: ...

  7. [JavaEE] 20141228_Java类文章搜集

    http://www.blogjava.net/jiangshachina 博客园java频道 Maven入门--概念与实例(原) Maven入门--较复杂的实例(原) Maven插件使用收集(原) ...

  8. ubantu下NiFi单节点安装部署

    第一步,首先下载安装包:http://nifi.apache.org/download.html,博主下载的是1.4.0版本,直接下载的是编译后的文件. 第二步:将压缩包上传到服务器相应目录下,并且解 ...

  9. php简单测试slim框架的功能

    php简单测试slim框架的功能 监听主路径/ $app->get( '/', function () { $template = <<<EOT<!DOCTYPE htm ...

  10. Blender Python UV 学习

    Blender Python UV 学习 1. bmesh面转换 bm = bmesh.from_edit_mesh(bpy.context.edit_object.data) bm.faces.en ...