问题描述
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;
}