题目背景
🐸在🐸模仿大赛中获得第二名,🐸非常开心,于是问了你一个很短的问题。
题面
给定长为 n n n 的正整数序列 { a i } \{a_i\} {ai}, q q q 次询问 [ l , r ] [l,r] [l,r] 的最短子段 [ l ′ , r ′ ] [l',r'] [l′,r′] ,满足所有在 [ l , r ] [l,r] [l,r] 中出现的数也在 [ l ′ , r ′ ] [l',r'] [l′,r′] 中出现。你只需要输出 [ l ′ , r ′ ] [l',r'] [l′,r′] 的长度 r ′ − l ′ + 1 r'-l'+1 r′−l′+1 。
1 ≤ n , q , a i ≤ 2 × 1 0 6 1\leq n,q,a_i\leq 2\times 10^6 1≤n,q,ai≤2×106 ,不强制在线,5s,1024MB。
题解
要找出现所有数的区间,我们可以求出现数的种类最多的情况下,最小的区间。
存在一种扫描线的 n 2 n^2 n2 做法,
我们把询问按照左端点放到序列上,从右往左扫,每次扫到位置 i i i 时,枚举区间 [ i , j ] [i,j] [i,j] ,算出 [ i , j ] [i,j] [i,j] 中不同的数的种类数 f ( i , j ) f(i,j) f(i,j) ,再用它更新 [ j , n ] [j,n] [j,n] 内的右端点。
考虑怎么优化。
首先, f ( i , j ) f(i,j) f(i,j) 可以用线段树维护,每次往左加入了一个数 x x x 后,把右边上一个等于 x x x 的位置删掉(实际上这种位置作为 j j j 已经产生不了贡献了),查询区间内的个数就是 f ( i , j ) f(i,j) f(i,j) 。
然后是每次的更新,我们可以对全局加一个更新的懒标记。区间的懒标记 ( x , y ) (x,y) (x,y) 表示这个区间前面有 x x x 种数,有效左端点在 y y y 。
我们惊喜地发现懒标记是可以合并的,只需要以 x x x 为第一关键字, y y y 为第二关键字保留最大值就行。当我们决定删掉一个位置,在这之前先把传到底的懒标记拿出来更新就好了。
尽量别用更劣的数据结构,常数非常重要。
维护右端点的答案可以随便用个树状数组。
时间复杂度 O ( n log n ) O(n\log n) O(nlogn) 。需要卡常。
CODE
#include<map>
#include<set>
#include<cmath>
#include<ctime>
#include<queue>
#include<stack>
#include<random>
#include<bitset>
#include<vector>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<unordered_map>
#pragma GCC optimize(2)
using namespace std;
#define MAXN 2000005
#define LL long long
#define ULL unsigned long long
#define ENDL putchar('\n')
#define DB double
#define lowbit(x) (-(x) & (x))
#define FI first
#define SE second
#define PR pair<int,int>
#define UIN unsigned int
int xchar() {
static const int maxn = 1000000;
static char b[maxn];
static int pos = 0,len = 0;
if(pos == len) pos = 0,len = fread(b,1,maxn,stdin);
if(pos == len) return -1;
return b[pos ++];
}
#define getchar() xchar()
inline LL read() {
LL f = 1,x = 0;int s = getchar();
while(s < '0' || s > '9') {if(s<0)return -1;if(s=='-')f=-f;s = getchar();}
while(s >= '0' && s <= '9') {x = (x<<1) + (x<<3) + (s^48);s = getchar();}
return f*x;
}
void putpos(LL x) {if(!x)return ;putpos(x/10);putchar((x%10)^48);}
inline void putnum(LL x) {
if(!x) {putchar('0');return ;}
if(x<0) putchar('-'),x = -x;
return putpos(x);
}
inline void AIput(LL x,int c) {putnum(x);putchar(c);}
int n,m,s,o,k;
int a[MAXN],nx[MAXN];
int qr[MAXN],as[MAXN];
vector<int> bu[MAXN];
PR c[MAXN];
inline void MX(PR &x,PR y) {if(y > x) x = y;}
inline void addc(int x,PR y) {while(x<=n)MX(c[x],y),x+=lowbit(x);}
inline PR summ(int x) {PR s;while(x>0)MX(s,c[x]),x&=(x-1);return s;}
struct it{
int ct,r; PR lz;
it(){ct=r=0;lz={0,0};}
it(int x) {ct=1;r=x;lz={0,0};}
}tre[MAXN<<2];
int M,bt;
int Max(int a,int b) {return a>b ? a:b;}
void maketree(int n) {M=1;bt=0;while(M<n+2) M<<=1,bt++;}
inline void upd(int s) {
if(s >= M) return ;
tre[s].ct = tre[s<<1].ct + tre[s<<1|1].ct;
tre[s].r = Max(tre[s<<1].r,tre[s<<1|1].r);
}
inline void pushdown(int s) {
if(s >= M) return ;
MX(tre[s<<1].lz,tre[s].lz);
MX(tre[s<<1|1].lz,{tre[s].lz.FI+tre[s<<1].ct,tre[s].lz.SE});
tre[s].lz = {0,0};
return ;
}
inline void pushr(int x) {
int s = M+x; for(int i = bt;i > 0;i --) pushdown(s>>i);
if(tre[s].ct && tre[s].lz.SE) addc(x,{tre[s].lz.FI+1,tre[s].lz.SE - x - 1});
}
inline void addtree(int x,it y) {
pushr(x);
int s = M+x;tre[s] = y;s >>= 1;
while(s) upd(s),s >>= 1;
}
inline PR findtree(int l,int r) {
int cn = 0,rr = 0,s = M+l-1,t = M+r+1;
for(int i = bt;i > 0;i --) pushdown(s>>i),pushdown(t>>i);
for(;(s>>1) ^ (t>>1);s >>= 1,t >>= 1) {
if(!(s&1)) cn += tre[s^1].ct,rr = Max(rr,tre[s^1].r);
if(t & 1) cn += tre[t^1].ct,rr = Max(rr,tre[t^1].r);
} return {cn,rr};
}
int main() {
freopen("frog.in","r",stdin);
freopen("frog.out","w",stdout);
n = read();
for(int i = 1;i <= n;i ++) {
a[i] = read();
}
maketree(n);
m = read();
for(int i = 1;i <= m;i ++) {
s = read(); qr[i] = read();
as[i] = qr[i] - s + 1;
bu[s].push_back(i);
}
for(int i = n;i > 0;i --) {
int x = a[i];
if(nx[x]) addtree(nx[x],it());
nx[x] = i; addtree(i,it(i));
MX(tre[1].lz,{0,i});
for(auto &j:bu[i]) {
int r = qr[j];
PR pr = findtree(i,r);
if(pr.SE) pushr(pr.SE);
PR s = summ(r);
if(s.FI == pr.FI) as[j] = min(as[j],-s.SE);
}
}
for(int i = 1;i <= m;i ++) {
AIput(as[i],'\n');
}
return 0;
}