简单的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. Resin 优化配置

    1.在resin.conf中的session-config内部添加 <session-max>4096</session-max> <session-timeout> ...

  2. Working with macro signatures

    https://docs.kentico.com/k11/macro-expressions/troubleshooting-macros/working-with-macro-signatures ...

  3. 【转】不要使用SBJSON(json-framework)

    原文网址:http://blog.devtang.com/2012/05/05/do-not-use-sbjson/ 不知道为什么,在iOS开发中,有很多人使用 SBJSON (又被称作json-fr ...

  4. 学习java的方式

  5. ManualResetEvent和AutoResetEvent的区别,分享来的

    在讨论这个问题之前,我们先了解这样一种观点,线程之间的通信是通过发信号来进行沟通的.(这不是废话) 先来讨论ManualResetEvent,讨论过程中我会穿插一些AutoResetEvent的内容, ...

  6. logrotate配置文件

    一.logrotate配置文件 /etc/logrotate.conf

  7. LUA 创建文件和文件夹

    创建文件: os.execute('mkdir e:\\aa') 创建文件夹: os.execute("cd.>e:\\wang.ini")

  8. CI中的验证码

    CI中的验证码相对来说使用非常方便直接加载类调用函数以及一些配置,代码如上,比较简单,具体函数可在CI手册的辅助函数参考CAPTCHA辅助函数中查询,CI中的验证码是直接生成验证码图片在你自己创建的C ...

  9. win10下mysql安装过程中遇到的各种坑

    前几天重装系统,又要下回来mysql,但没想到还是遇到了许多麻烦,翻了十多篇博文才搞定,写个总结出来方便以后不要重复踩坑,也给大家参考参考. 1.下载与安装 这个没什么好说的,下载地址网上一大堆,安装 ...

  10. Command "python setup.py egg_info" failed with error code 1 in /tmp/pip-build*解决办法

    easy_install -U setuptools or pip install ipython 亲测有效