优美序列(sequence)

问题描述

Lxy养了N头奶牛，他把N头奶牛用1..N编号，第i头奶牛编号为i。为了让奶牛多产奶，每天早上他都会让奶牛们排成一排做早操。奶牛们是随机排列的。在奶牛排列中，如果一段区间[L,R]中的数从小到大排列后是连续的，他认为这段区间是优美的。比如奶牛排列为：(3, 1, 7, 5, 6, 4, 2)，区间[3,6]是优美的，它包含4,5,6,7连续的四个数，而区间[1,3] 是不优美的。Lxy的问题是：对于给定的一个区间[L,R]（1<=L<=R<=N）, 他想知道，包含区间[L,R]的最短优美区间，比如区间[1,3]的最短优美区间是[1,7]。

输入

第一行为一个整数N，表示奶牛的个数。
第二行为1到N的一个排列，表示奶牛的队伍。
第三行为一个整数M，表示有M个询问。
后面有M行，每行有两个整数L,R表示询问区间。

输出

输出为M行，每行两个整数，表示包含询问区间的最短优美区间。

输入输出样例

样例1

input
7
3 1 7 5 6 4 2
3
3 6
7 7
1 3
output
3 6
7 7
1 7

样例2

input
10
2 1 4 3 5 6 7 10 8 9
5
2 3
3 7
4 7
4 8
7 8
output
1 4
3 7
3 7
3 10
7 10

数据范围

15%的数据满足: 1 <=N,M <= 15；
50%的数据满足：1 <=N,M <= 1000。
100%的数据满足：1 <=N,M <= 100000。

解析

为了给出题人留面子，不让他被骂出个题只有一种解法而且只有满分解法，我先给个80分做法。
预处理一个数组pos[i]，表示数i在原数组的位置(第几个数)。
在询问区间内找一个最小值和最大值，作为新的l和r，新的l和r作为一个询问区间，在pos数组里找一个最小值和最大值，他们再作为新的l和r，形成询问区间，在原数组寻找最小值和最大值，并再次重复刚才的操作……
上述过程可能很无脑很中二，然而这就是模拟题意找出答案的过程。
最大值最小值用st表预处理，一共四个st表，比线段树好写……
那么……
正解是分治，像归并排序那样。
不断求[l,mid] [mid+1,r]的优美序列，直到分成[mid,mid+1]
可以证明是正确的。
事实上分治做的熟练的人可能一眼秒这道题。

代码

80pts

#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
template<class T>
void read(T &res)
{
	res = 0;
	T f = 1;
	char cc = getchar();
	while(cc < '0' || cc > '9')
	{
		if(cc == '-') f = -1;
		cc = getchar();
	}
	while(cc >= '0' && cc <= '9')
	{
		res = res * 10 + cc - '0';
		cc = getchar();
	}
	res *= f;
}
const int maxn = 1e5 + 5;
int n,m;
int seqmin[maxn][50];
int posmin[maxn][50];
int seqmax[maxn][50];
int posmax[maxn][50];
int ansl,ansr;
void init()
{
	for(int j=1;(1<<j) <= n;j++)
    {
        for(int i=0;i+(1<<j)-1 <= n;i++)
        {
            seqmax[i][j] = max(seqmax[i][j-1],seqmax[i+(1<<(j-1))][j-1]); // 1
            seqmin[i][j] = min(seqmin[i][j-1],seqmin[i+(1<<(j-1))][j-1]); // 2
            posmax[i][j] = max(posmax[i][j-1],posmax[i+(1<<(j-1))][j-1]); // 3
            posmin[i][j] = min(posmin[i][j-1],posmin[i+(1<<(j-1))][j-1]); // 4
        }
    }
}
int query(int l,int r,int ju)
{
    int t = int(log((double)(r - l + 1)) / log(2.0));
    if(ju == 1)
    	return max(seqmax[l][t],seqmax[r-(1<<t)+1][t]);
    if(ju == 2)
    	return min(seqmin[l][t],seqmin[r-(1<<t)+1][t]);
	if(ju == 3)
    	return max(posmax[l][t],posmax[r-(1<<t)+1][t]);
    if(ju == 4)
    	return min(posmin[l][t],posmin[r-(1<<t)+1][t]);
}
void change(int inl,int inr)
{
	int templ,tempr,endl,endr;
	templ = query(inl,inr,2);
	tempr = query(inl,inr,1);
	endl = query(templ,tempr,4);
	endr = query(templ,tempr,3);
	if(endl == inl && endr == inr)
	{
		ansl = endl;
		ansr = endr;
		return;
	}
	else
		change(endl,endr);
	
}
int main()
{
	freopen("sequence.in","r",stdin);
	freopen("sequence.out","w",stdout);
	read(n);
	for(int i=1;i<=n;i++)
		read(seqmin[i][0]),seqmax[i][0] = seqmin[i][0];
	for(int i=1;i<=n;i++)
		posmax[seqmin[i][0]][0] = i,posmin[seqmin[i][0]][0] = i;
	init();
	read(m);
	while(m--)
	{
		int l,r;
		read(l);read(r);
		change(l,r);
		printf("%d %d
",ansl,ansr);
	}
	return 0;
}

100pts

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N=100005;
#define pb push_back
#define lc (o<<1)
#define mid (l+r>>1)
#define rc ((o<<1)|1)
struct ask{
    int L,num;
    bool operator <(const ask &u)const{
        return L<u.L;
    }
};
vector<ask> g[N];
multiset<ask> s;
multiset<ask> :: iterator it;
int n,a[N],p[N],mx[N*4],pos[N*4],Q;
int le,ri,w,al[N],ar[N],M,P,tag[N*4];
inline void mt(int o){
    mx[o]=max(mx[lc],mx[rc]);
    pos[o]=(mx[rc]==mx[o]?pos[rc]:pos[lc]);
}
inline void work(int o,int der){ tag[o]+=der,mx[o]+=der;}
inline void pd(int o){ if(tag[o]) work(lc,tag[o]),work(rc,tag[o]),tag[o]=0;}
void build(int o,int l,int r){
    mx[o]=r,pos[o]=r;
    if(l==r) return;
    build(lc,l,mid),build(rc,mid+1,r);
}
void update(int o,int l,int r){
    if(l>=le&&r<=ri){ work(o,w); return;}
    pd(o);
    if(le<=mid) update(lc,l,mid);
    if(ri>mid) update(rc,mid+1,r);
    mt(o);
}
void query(int o,int l,int r){
    if(l>=le&&r<=ri){ if(mx[o]>M) M=mx[o],P=pos[o]; return;}
    pd(o);
    if(ri>mid) query(rc,mid+1,r);
    if(le<=mid) query(lc,l,mid);
}
inline bool can(int cc){
    le=1,ri=it->L,M=P=0;
    query(1,1,n);
    if(M==cc){ al[it->num]=P,ar[it->num]=cc; return 1;}
    return 0;
}
inline void solve(){
    build(1,1,n),w=1;
    ask inf=(ask){233333,0};
    for(int i=1;i<=n;i++){
        for(ask x:g[i]) s.insert(x);
        ri=p[a[i]-1],le=1;
        if(ri&&ri<i) update(1,1,n);
        ri=p[a[i]+1],le=1;
        if(ri&&ri<i) update(1,1,n);
        for(;s.size();s.erase(it)){
            it=--s.lower_bound(inf);
            if(!can(i)) break;
        }
    }
}
int main(){
  freopen("sequence.in","r",stdin);
  freopen("sequence.out","w",stdout);
     
    scanf("%d",&n);
    for(int i=1;i<=n;i++) {
    	scanf("%d",&a[i]);
    	p[a[i]]=i;
    };
    scanf("%d",&Q);
    for(int i=1,l,r;i<=Q;i++){
        scanf("%d %d",&l,&r);
        g[r].pb((ask){l,i});
    }
    solve();
    for(int i=1;i<=Q;i++)  printf("%d %d
",al[i],ar[i]);
    return 0;
}