POJ 3368

题目给了不降序列,让我们去区间询问区间最长连续相等序列。
我们这里对于例子
10 3
-1 -1 1 1 1 1 3 10 10 10
做出这样的二元对。
(-1,2),(1,4),(3,1),(10,4)
这样一来,对于二元对第二个数字即出现次数,可以用rmq来进行维护,表示求二元对中,第二个数字最大的那个。即是我们需要求的最大序列长度。
然而有问题的是,如果查询,1,3,这样的区间,(-1,2)和(1,4)这两个二元对都没有包含满,而且可以看出来用rmq查询这两个区间是错误的,那么我们可以预处理一下,使得问题转换成为,rmq模型。
比如 我们给的区间在四个二元对当中,其中最左边的二元对没有包含满,最右边的二元对没有包含满,那么我们先除掉这两个区间(因为这两个区间可以直接算出来最大长度),剩下两个二元对就可以rmq查询了。
注意结尾,是只输入一个n=0,就结束了,并不是n=q=0才结束。
如果用rmq进行查询的时候,找最大k使得i+k在区间(i,j)里面的时候,用log去找k的时候,记得注意强制类型转换,不过还是直接按照lrj的紫书的模板些就好了,毕竟时间复杂度也就增加了log级别的而已,而且是肯定不容易错误。
rmq类型以后用dp写,不容易出错,用线段树去推很容易就出错了。
注意stl容器不要直接写在while里面,用宏定义会减少很多时间,防止tle。

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define up(i,a,b)  for(int i=a;i
#define dw(i,a,b)  for(int i=a;i>b;i--)
#define upd(i,a,b) for(int i=a;i<=b;i++)
#define dwd(i,a,b) for(int i=a;i>=b;i--)
//#define local
typedef long long ll;
const double esp = 1e-6;
const double pi = acos(-1.0);
const long long INF = 0x3f3f3f3f;
using namespace std;
typedef pair<int, int> pir;
int n, q;
int mod[100015];
int dp[100010][100];
int l[100010];
int r[100010];
int st[100010];
int mp[100010];
void init(int len)
{
	upd(i, 1, len)
	{
		dp[i][0] = st[i];
	}
	for (int j = 1; (1 << j) <= len; j++)
	{
		for (int i = 1; i + (1 << j) - 1 <= len; i++)
		{
			dp[i][j] = max(dp[i][j - 1], dp[i + (1 << (j - 1))][j - 1]);
		}
	}
}
int rmq_querry(int x, int y)
{
	if (x > y)return 0;//为了后面代码写的方便
	int k = 0;
	//while ((1 << (k + 1)) <= y - x + 1)k++;
	k = int(log((double)y - x + 1)/log(2.0));//注意 int要给括号给后面全部
	return max(dp[x][k], dp[y - (1 << k) + 1][k]);
}
int main()
{
	while (scanf("%d",&n)!=EOF)
	{
		if (n == 0)break;
		scanf("%d", &q);
		int len = 0;
		int qq = 1;
		mod[0] = 1e9;
		upd(i, 1, n)
		{
			scanf("%d", &mod[i]);
			if (mod[i] != mod[i - 1])
			{
				l[++len] = i;
				r[len] = i;
				qq= 1;
				st[len] = qq;
			}
			else
			{
				r[len] = i;
				st[len] = ++qq;
			}
			mp[i] = len;
		}//预处理,模板
		init(len);
		/*upd(i, 1, len)
			cout << st[i];*/
		int x, y;
		while (q--)
		{
			scanf("%d %d", &x, &y);
			if (mp[x] == mp[y])
			{
				printf("%d\n", y - x + 1);
			}
			else
			{
				int llen = r[mp[x]] - x + 1;
				int rlen = y - l[mp[y]] + 1;
				printf("%d\n", max(max(llen, rlen), rmq_querry(mp[x] + 1, mp[y] - 1)));//这里直接全部放进去判断就好了
			}
		}
	}
	return 0;
}

你可能感兴趣的:(rmq)