关于查询区间最小没出现的自然数的cdq方法的可行性探讨

关于查询区间最小没出现的自然数的cdq方法的可行性探讨_第1张图片这道题显然有一个可持久化线段树的做法。

首先我们意识到这个极广的值域没有什么用处。

我们首先想到必然存在一个   x \ x  x使答案为   x \ x  x或者   [ 1 , x − 1 ] \ [1,x-1]  [1,x1]中的最小的不存在的数字。

所以首先我们想到找到这个   x \ x  x然后把没有必要的数字全部删掉,这样剩下的数字必然再   [ 1 , n ] \ [1,n]  [1,n]中。

之后我们想到将   a i \ a_{i}  ai   i \ i  i点存在转换为   a i \ a_{i}  ai   [ 1 , i − 1 ] , [ i + 1 , n ] \ [1,i-1],[i+1,n]  [1,i1],[i+1,n]中不存在,设   i \ i  i的两个值为   q = i − 1 , w = i + 1 \ q=i-1,w=i+1  q=i1,w=i+1,这样任意询问   [ l , r ] \ [l,r]  [l,r]就是   q ≥ r \ q \ge r  qr或者   w ≤ l \ w \le l  wl的数字中的最小值。这显然是一个一维偏序问题,那么我们显然可以将每个询问看做两个,直接桶排序即可。复杂度   O ( n ) \ O(n)  O(n)

实际上这种方法仅限于值的两两不同。

如果相同的值是存在的,那么我们依然可以将问题改变。显然如果同时存在两个值相同   a i = a j , i < j \ a_{i}=a_{j},i ai=aj,i<j,那么显然为   a i \ a_{i}  ai   [ 1 , i − 1 ] , [ i + 1 , j − 1 ] , [ j + 1 , n ] \ [1,i-1],[i+1,j-1],[j+1,n]  [1,i1],[i+1,j1],[j+1,n]不存在,我们将其拆为三个操作,   q 1 = 1 , w 1 = i − 1 , q 2 = i + 1 , w 2 = j − 1 , q 3 = j + 1 , w 3 = n \ q_{1}=1,w_{1}=i-1,q_{2}=i+1,w_{2}=j-1,q_{3}=j+1,w_{3}=n  q1=1w1=i1,q2=i+1,w2=j1,q3=j+1,w3=n。对于每个询问   [ l , r ] \ [l,r]  [l,r],即   q ≤ l \ q \le l  ql   w ≥ r \ w \ge r  wr的最小值。那么我们先按照   q \ q  q或者   l \ l  l排序,cdq中按照   r \ r  r或者   w \ w  w,然后求解即可。   O ( n log ⁡ 2 n ) \ O(n \log_{2}{n})  O(nlog2n)

#include
using namespace std;
inline int read()
{
    char ch = getchar(); int x = 0;
    while (!isdigit(ch)) ch = getchar();
    while (isdigit(ch)){ x = x*10+ch-'0'; ch = getchar(); }
    return x;
}
int n,a[300300],b[300300],m,xxx,tot=0,head,tail,ans[200200];
struct nobe
{
	int op,l,r,w,id;
}q[600600],tmp[600600];
inline bool cmp(nobe a,nobe b)
{
	return (a.l^b.l) ? (a.l<b.l) : (a.op<b.op);
}
inline void cdq(int l,int r)
{
	if(l==r) return ;
	int mid=(l+r)>>1;
	cdq(l,mid);
	cdq(mid+1,r);
	int i=l,j=mid+1,ccnt=0,sum=999999999;
	while(j<=r)
	{
		while((q[i].r>=q[j].r)&&(i<=mid))
		{
			tmp[++ccnt]=q[i];
			if(q[i].op&1)
			{
				sum=min(sum,q[i].w);
			}
			++i;
		}
		tmp[++ccnt]=q[j];
		if(q[j].op^1)
		{
			ans[q[j].id]=min(ans[q[j].id],sum);
		}
		++j;
	}
	while(i<=mid)
	{
		tmp[++ccnt]=q[i];
		++i;
	}
	i=1;
	while(i<=ccnt)
	{
		q[l+i-1]=tmp[i];
		++i;
	}
}
int main()
{
	memset(b,0,sizeof(b));
	n=read();
	m=read();
	int i=1;
	while(i<=n)
	{
		a[i]=read();
		if(a[i]<=n) ++b[a[i]];
		++i;
	}
	i=-1;
	while(i<=n)
	{
		if(b[i+1]) ++i;
		else break;
	}
	xxx=i+1;
	i=0;
	while(i<xxx)
	{
		b[i]=0;
		++i;
	}
	i=1;
	while(i<=n)
	{
		if(a[i]<xxx)
		{
			if((i^1)&&(b[a[i]]^i-1))
			{
				q[++tot]=(nobe){1,b[a[i]]+1,i-1,a[i],0};
			}
			b[a[i]]=i;
		}
		++i;
	}
	i=0;
	while(i<xxx)
	{
		if(b[i]^n) q[++tot]=(nobe){1,b[i]+1,n,i,0};
		++i;
	}
	i=1;
	while(i<=m)
	{
		int l,r;
		l=read();
		r=read();
		ans[i]=xxx;
		q[++tot]=(nobe){2,l,r,0,i};
		++i; 
	}
	sort(q+1,q+tot+1,cmp);
	cdq(1,tot);
	i=1;
	while(i<=m)
	{
		printf("%d\n",ans[i]);
		++i;
	}
	return 0;
}

实际上只出现一次的数字依然可以桶排序来解决,将大大优化常数。

你可能感兴趣的:(分治)