【BZOJ4059】Non-boring sequences

Solution

  记序列为\(a\),计算出与\(a_i\)相等的前一个元素的位置\(pre_i\),以及后一个元素的位置\(nex_i\),显然,对于那些左端点处于\((pre_i,i]\)以及右端点处于\([i,nex_i)\)的区间都可以认为是合法的。
  
  那么我们可以将每个区间\([l,r]\)抽象成一个二维平面的点\((l,r)\),每一个元素可以使得一部分区间合法,可以抽象为一个横坐标范围为\([pre_i+1,i]\)且纵坐标范围为\([i,nex_i-1]\)的矩形。对所有矩形进行求交,如果矩形覆盖了所有的合法点(也就是一个上三角形),那么序列就是合法的,否则不合法。
  
​  但是这样跑的非常慢!在BZOJ上AC的扫描线代码在我们OJ上是完全被卡常的.....有没有更加简便的做法?
  
  考虑\(check(l,r)\)表示所有左右端点处于\([l,r]\)的区间是否合法。则答案就是\(check(1,n)\)
  
​  枚举\(i\in[l,r]\),一旦找到一个\(i\)使得\(pre_i<l\)\(nex_i>r\),就停下。如果找不到,显然序列\([l,r]\)本身就不合法,直接返回\(false\)。此时我们直接可以得知,左端点处于\([l,i]\)且右端点处于\([i,r]\)的区间全部是合法的!还剩下左右断电都处于\([l,i)\)\((i,r]\)的区间未检查,返回\(check(l,i-1)\&\&check(i+1,r)\)即可。
  
​  但是枚举这一步的复杂度我们没有保证,最坏总复杂度会达到\(O(n^2)\)。我们需要一个思想:从两头向中间同时推进枚举,直到遇到第一个所需点为止。
  
  时间复杂度是\(T(n)=\max\{T(i)+T(n-i)+\min(n,n-i)\}=O(n\lg n)\)
  
  为什么呢?感性地讲,如果把递归步骤倒过来看,就是一个启发式合并!关键就在\(min(n,n-i)\)这里,每一层递归贡献的复杂度恰好是关键点与边缘的距离,而我们每次都找最靠近边缘的一个关键点,相当于启发式里面的对较小的部分进行操作的思想一样。
  
  然后就做完了,这种两端向中间枚举的思想很值得学习和思考。
  
  
  

Code

  

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long ll;
const int N=200005;
int n,a[N];
int dizlis[N],dizcnt;
int pre[N],nex[N],mark[N];
void Diz(){
	dizcnt=n;
	memcpy(dizlis,a,(dizcnt+1)*sizeof(int));
	sort(dizlis+1,dizlis+1+dizcnt);
	dizcnt=unique(dizlis+1,dizlis+1+dizcnt)-dizlis-1;
	for(int i=1;i<=n;i++)
		a[i]=lower_bound(dizlis+1,dizlis+1+dizcnt,a[i])-dizlis;
}
bool check(int l,int r){
	if(l>r) return true;
	int p1=l,p2=r;
	while(p1<=p2){
		if(pre[p1]<l&&nex[p1]>r) 
			return check(l,p1-1)&&check(p1+1,r);
		p1++;
		if(pre[p2]<l&&nex[p2]>r) 
			return check(l,p2-1)&&check(p2+1,r);
		p2--;
	}
	return false;
}
int Main(){
	scanf("%d",&n);
	for(int i=1;i<=n;i++) scanf("%d",&a[i]);
	Diz();
	for(int i=1;i<=dizcnt;i++) mark[i]=0;
	for(int i=1;i<=n;i++)
		pre[i]=mark[a[i]],mark[a[i]]=i;
	for(int i=1;i<=dizcnt;i++) mark[i]=n+1;
	for(int i=n;i>=1;i--)
		nex[i]=mark[a[i]],mark[a[i]]=i;
	puts(check(1,n)?"non-boring":"boring");
	return 0;
}
int main(){
	int T;
	scanf("%d",&T);
	while(T--) Main();
	return 0;
}
posted @ 2018-06-21 10:48  RogerDTZ  阅读(106)  评论(0编辑  收藏  举报