线段树专题-等差子序列

感谢

感谢孙耀峰的线段树PPT,使我获益匪浅.

题目来源

BZOJ−2124BZOJ-2124BZOJ−2124

题意

给出长度为nnn的1−n1-n1−n的排列AAA
问是否存在一组1≤p1≤p2≤...≤pl≤n,l≥31 \le p_1 \le p_2 \le ... \le p_l \le n,l \ge 31≤p1≤p2≤...≤pl≤n,l≥3
使得Ap1,Ap2,...,AplA_{p1},A_{p2},...,A_{pl}Ap1,Ap2,...,Apl构成等差序列
数据范围n≤104n \le 10^4n≤104

题解

只要能形成长度为333的等差子序列,就直接输出YesYesYes即可.

即需要寻找i,j,ki,j,ki,j,k满足Ai=Aj−t,Ak=Aj+t,且t>0.A_i = A_j - t,A_k = A_j + t,且t > 0.Ai=Aj−t,Ak=Aj+t,且t>0.

根据转化后的式子,我们只需要枚举jjj,判断是否存在这样的ttt即可.

考虑到AAA是一个排列,排列中的数两两不同,因此我们可以考虑开一个数据结构维护一下内容:

当枚举到AjA_jAj的时候,我可以很快的查询任何一个AiA_iAi在AjA_jAj的左边还是右边,(记如果i<ji<ji<j则AjA_jAj应该标记为000表示在左边,若i>ji>ji>j则AjA_jAj应该标记为111表示在右边).

这样的话,如果010101串[Aj−t,Aj−1][A_j-t,A_j-1][Aj−t,Aj−1]与[Aj+1,Aj+t][A_j+1,A_j+t][Aj+1,Aj+t]之间有一位不相同(也就是说存在关于AjA_jAj的对称位置的两个数在AjA_jAj两侧,说明可以形成等差子序列).

因此只有当以上两个01串,翻转过来完全相同的情况下,才不存在以AjA_jAj作为中心的等差序列.

如何判断010101串相同呢?答案是哈希.

建立一颗线段树,维护的是010101串的正向和反向哈希值.

从小到大枚举Aj,[1,n]A_j,[1,n]Aj,[1,n],然后更新线段树,并判断是否有以AjA_jAj为对称中心的等差序列.

代码

#include <iostream>
#include <cstring>
const int N = 10007;
typedef unsigned long long ull;
int T,n;
int a[N];ull pow[N];
struct work{ull hash0,hash1;int len;
};#define pr(x) std::cout << #x << ": " << x << std::endl
struct segtree{ull hashs[N * 4][2];void maintain(int rt,int l,int r) {int mid = (l + r) / 2;hashs[rt][0] = hashs[rt*2][0] + hashs[rt*2+1][0] * pow[mid - l + 1];hashs[rt][1] = hashs[rt*2+1][1] + hashs[rt*2][1] * pow[r - mid];}void change(int rt,int l,int r,int pos,int val) {if(l == r) {hashs[rt][0] = hashs[rt][1] = val;return ;}   int mid = (l + r) / 2;if(pos > mid) change(rt*2+1,mid+1,r,pos,val);elsechange(rt*2,l,mid,pos,val);maintain(rt,l,r);}work query(int rt,int l,int r,int ql,int qr) {if(r < ql || qr < l) {return (work){0,0,0};} else if(ql <= l && r <= qr) {return (work){hashs[rt][0],hashs[rt][1],r - l + 1};}int mid = (l + r) / 2;work wa = query(rt*2,l,mid,ql,qr);work wb = query(rt*2+1,mid+1,r,ql,qr);return (work){wa.hash0 + wb.hash0 * pow[wa.len],wb.hash1 + wa.hash1 * pow[wb.len],wa.len + wb.len};}}seg;int main()
{pow[0] = 1;for(int i = 1;i < N;++i) pow[i] = 3 * pow[i-1];std::ios::sync_with_stdio(false);std::cin >> T;while(T--) {int f = 0;memset(&seg,0,sizeof(seg));std::cin >> n;for(int i = 1;i <= n;++i) {std::cin >> a[i];seg.change(1,1,n,a[i],1);if(a[i] == 1 || a[i] == n) continue;int len = std::min(a[i]-1,n - a[i]);work wa = seg.query(1,1,n,a[i]-len,a[i]-1);work wb = seg.query(1,1,n,a[i]+1,a[i]+len);if(wa.hash0 != wb.hash1) {f = 1;}}if(!f) std::cout << "N" << std::endl;elsestd::cout << "Y" << std::endl;}return 0;
}

线段树专题-等差子序列 BZOJ-2124相关推荐

线段树专题-黑白棋盘 BZOJ-1453
线段树专题-黑白棋盘题目来源 BZOJ−1453BZOJ-1453BZOJ−1453 题意 QQQ次操作每次操作给出(x,y)(x,y)(x,y),将(x,y)(x,y)(x,y)个格子颜色取反 ...
2021牛客暑期多校训练营7 B-xay loves monotonicity（线段树+不降子序列）
P4198 楼房重建线段树维护以某点为开头的最长不下降子序列使用下面calc函数能够计算线段树u维护的区间中,以x为开头最长不下降子序列的个数. calc需要维护区间最值. template< ...
2018 UESTC 线段树专题
A - 一棵简单的线段树 A[1...n]初始全为0. 1. 给两个数p 和 x(1≤p≤n),单点更新 A[p] <- x 2. 给两个数L和R (1≤L<R≤n), L到R区间里这几 ...
蒟蒻君的刷题日记Day12（线段树专题T4）：P8082 [COCI2011-2012#4] KEKS 线段树版题解
解题思路看题解区的大佬们用的都是单调栈,本蒟蒻献上一篇线段树题解. 整个数最大,首先位数是确定的,则肯定优先考虑高位大小. 大体思路就是从前向后依次求出每一位的值(好像是废话). 对于第 iii 位 ...
线段树专题 A(单点更新)
#include <iostream> #include<cstdio> #include<cstring> using namespace std; const ...
HDU 4417 Super Mario（线段树||树状数组+离线操作之线段树篇）
Mario is world-famous plumber. His "burly" figure and amazing jumping ability reminded in ...
【POJ 2482】 Stars in Your Window（线段树+离散化+扫描线）
[POJ 2482] Stars in Your Window(线段树+离散化+扫描线) Time Limit: 1000MS Memory Limit: 65536K Total Submiss ...
【BZOJ3821/UOJ46】玄学（二进制分组，线段树）
[BZOJ3821/UOJ46]玄学(二进制分组,线段树) 题面 BZOJ UOJ 题解呜,很好的题目啊QwQ. 离线做法大概可以线段树分治,或者直接点记录左右两次操作时的结果,两个除一下就可以直接 ...
【BZOJ4991】我也不知道题目名字是什么（线段树）
[BZOJ4991]我也不知道题目名字是什么(线段树) 题面 BZOJ 题解对于线段树维护的区间维护以下东西: 区间左(右)端开始(结束)的最长(短)子串的长度左端右端的值,以及当前区间内的答案 ...

线段树专题-等差子序列 BZOJ-2124

线段树专题-等差子序列

感谢

题目来源

题意

题解

代码

线段树专题-等差子序列 BZOJ-2124相关推荐

最新文章

热门文章