洛谷P7907 [Ynoi2005] rmscne
数据结构好题
首先将询问离线,扫描线扫答案区间的左端点。
设和 \(l\) 颜色相同的下一个位置为 \(x\) 。
那么对于左端点 \(\leq l\),\(l \leq\) 右端点 $<x $ 的询问, \(l\) 再往右移动的话一定是不合法的,找到这些询问并且删除即可。
设所有出现位置全部在 \((l,x)\) 外的颜色中在 \(x\) 后最小的位置为 \(y\)。
对于右端点 \(>=y\) 的询问,他的答案区间如果以 \(l\) 为端点一定不优,因为如果以 \(l\) 为端点他一定会包含 \(y\),也就包含 \(x\),不如往右移,不用考虑。
所以我们只需要维护右端点在 \(x,y\) 之间的区间就可以了,这可以二维数点来维护。
但是这样会出现问题,就是会修改到之间已经删除掉的不合法的区间。
只要在删除的时候把右端点改成答案区间的右端点就好了。
代码
#include<bits/stdc++.h>
#define N 2001001
#define MAX 2001
#define eps 1e-10
using namespace std;
typedef int ll;
typedef double db;
const db PI=acos(-1);
const ll mod=64123,inv2=(mod+1)/2,inf=1e9;
inline char nc(){
static char buf[100000],*p1=buf,*p2=buf;
return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline ll read(){
char ch=nc();int sum=0;
while(!(ch>='0'&&ch<='9'))ch=nc();
while(ch>='0'&&ch<='9')sum=sum*10+ch-48,ch=nc();
return sum;
}
ll n,a[N],m,nxt[N],las[N],ans[N],pd[N],ls[N],rs[N],headv[N],nexv[N],totv,headg[N],nexg[N],totg,totl,nexl[N],v[N],g[N];
inline void vdd(ll x,ll y)
{
v[++totv]=y;
nexv[totv]=headv[x];
headv[x]=totv;
return;
}
inline void gdd(ll x,ll y)
{
g[++totg]=y;
nexg[totg]=headg[x];
headg[x]=totg;
return;
}
struct point
{
ll l,r;
}q[N];
bool vis[N];
struct cc
{
ll x,l,r,v;
};
struct node
{
ll val,tag;
}seg[N<<2];
cc lis[N];
inline void build(ll pos,ll l,ll r)
{
seg[pos].val=seg[pos].tag=inf;
if(l!=r)
{
ll mid=l+r>>1;
build(pos<<1,l,mid);
build(pos<<1|1,mid+1,r);
}
return;
}
inline void add(ll pos,ll num)
{
seg[pos].val=seg[pos].val<num?seg[pos].val:num;
return;
}
ll s,t,num;
inline void upgrade(ll pos,ll l,ll r)
{
if(l>t||r<s)
return;
else if(l>=s&&r<=t)
return add(pos,num);
if(num>=seg[pos].val)
return;
ll mid=l+r>>1;
upgrade(pos<<1,l,mid);
upgrade(pos<<1|1,mid+1,r);
return;
}
ll res=inf;
inline void query(ll pos,ll l,ll r,ll p)
{
if(l>p||r<p)
return;
res=min(res,seg[pos].val);
if(l==r)
return;
ll mid=l+r>>1;
if(p<=mid)
query(pos<<1,l,mid,p);
else
query(pos<<1|1,mid+1,r,p);
return;
}
struct ct
{
ll siz;
}st[N<<2];
inline ct operator +(ct x,ct y)
{
return ct{x.siz+y.siz};
}
inline void upd(ll pos,ll l,ll r,ll p,ll num)
{
ll now=pd[p];
while(now)
{
st[now].siz+=num;
now>>=1;
}
return;
}
struct ncp
{
ll num;
}sg[N<<2];
inline ncp operator +(ncp x,ncp y)
{
return ncp{x.num+y.num};
}
inline void update(ll pos,ll l,ll r,ll p)
{
ll now=pd[p];
while(now)
{
sg[now].num++;
now>>=1;
}
return;
}
inline ll findpre(ll pos,ll l,ll r,ll p)
{
if(l>p||!st[pos].siz)
return 0;
else if(l==r)
return l;
ll mid=l+r>>1;
ll tmp=findpre(pos<<1|1,mid+1,r,p);
if(tmp)
return tmp;
return findpre(pos<<1,l,mid,p);
}
inline ll findnxt(ll pos,ll l,ll r,ll p)
{
if(r<p||!st[pos].siz)
return n+1;
else if(l==r)
return l;
ll mid=l+r>>1;
ll tmp=findnxt(pos<<1,l,mid,p);
if(tmp!=n+1)
return tmp;
return findnxt(pos<<1|1,mid+1,r,p);
}
ll cur;
inline void digui(ll pos,ll l,ll r)
{
if(!sg[pos].num)
return;
else if(l==r)
{
ll tmp=findpre(1,1,n,l);
for(int k=headv[l];k;k=nexv[k])
{
int j=v[k];
ans[j]=ans[j]<tmp-num+1?ans[j]:tmp-num+1;
q[j].r=tmp;
}
headv[l]=0;
sg[pos].num=0;
return;
}
ll mid=l+r>>1;
digui(pos<<1|1,mid+1,r);
digui(pos<<1,l,mid);
sg[pos]=sg[pos<<1]+sg[pos<<1|1];
return;
}
inline void work(ll pos,ll l,ll r)
{
if(!sg[pos].num)
return;
if(l>t||r<s)
return;
else if(l>=s&&r<=t)
return digui(pos,l,r);
ll mid=l+r>>1;
work(pos<<1|1,mid+1,r);
work(pos<<1,l,mid);
sg[pos]=sg[pos<<1]+sg[pos<<1|1];
return;
}
inline void prework(ll pos,ll l,ll r)
{
if(l==r)
pd[l]=pos;
else
{
ll mid=l+r>>1;
prework(pos<<1,l,mid);
prework(pos<<1|1,mid+1,r);
}
return;
}
static char buf[100000],*pp=buf;
template<class T>inline void pc(T c)
{
if(pp-buf==100000)fwrite(buf,1,100000,stdout),pp=buf;
*pp++=c;
}
inline void fsh(){fwrite(buf,1,pp-buf,stdout);pp=buf;}
char cs[N],tot;
inline void write(ll x)
{
tot=0;
cs[++tot]='\n';
while(x)
cs[++tot]=(x%10)^48,x/=10;
while(tot)
pc(cs[tot--]);
}
signed main()
{
n=read();
for(int i=1;i<=n;i++)
a[i]=read();
for(int i=0;i<N;i++)
las[i]=n+1;
for(int i=n;i;i--)
{
nxt[i]=las[a[i]];
las[a[i]]=i;
}
m=read();
prework(1,1,n);
for(int i=1;i<=n;i++)
{
if(vis[a[i]])
continue;
vis[a[i]]=true;
upd(1,1,n,i,1);
}
for(int i=1;i<=m;i++)
{
ans[i]=inf;
q[i].l=read();
q[i].r=read();
gdd(q[i].l,i);
}
for(int i=1;i<=n;i++)
{
for(int k=headg[i];k;k=nexg[k])
{
int j=g[k];
vdd(q[j].r,j);
update(1,1,n,q[j].r);
}
s=i,t=nxt[i]-1,num=i;
work(1,1,n);
ll tmp=findnxt(1,1,n,nxt[i]);
lis[++totl]=cc{i,nxt[i],tmp-1,findpre(1,1,n,nxt[i])-i+1};
if(nxt[i]!=n+1)
upd(1,1,n,nxt[i],1);
}
totv=0;
for(int i=1;i<=n;i++)
headv[i]=0;
for(int i=1;i<=totl;i++)
vdd(lis[i].x,i);
build(1,1,n);
ll cnt=0;
for(int i=n;i;i--)
{
for(int k=headv[i];k;k=nexv[k])
{
int j=v[k];
s=lis[j].l,t=lis[j].r,num=lis[j].v;
upgrade(1,1,n),cnt++;
}
for(int k=headg[i];k;k=nexg[k])
{
int j=g[k];
res=inf;
query(1,1,n,q[j].r);
ans[j]=min(ans[j],res);
}
}
for(int i=1;i<=m;i++)
write(ans[i]);
fsh();
quick_exit(0);
}