Luogu P5350 序列

ODT

这道题目太毒瘤啦,经过了无数遍的TLE、WA,和RE(TAT),我终于了解了珂朵莉树的强(R)大(E),我会详细的介绍关于TLE,WA和RE的原因。

首先我们看到区间赋值操作和保证数据随机,我们的第一直觉肯定是幸福的珂朵莉树啦,虽然在刻意构造的数据下她的时间复杂度是错误的,但是在随机数据下她的表现十分优秀,甚至可以碾压其他数据结构。

1.首先是split操作,这是ODT的基本操作,可以将完整的区间拆分开来。但应注意的是,

if(it!=s.end()&&it->l==pos)return it;

所以我们需要在n的后面加上一个数防止出锅,

s.insert((node){n+1,n+1,(LL)0});

2.区间加、求和、推平,这些都是ODT模板上的操作,注意一下及时取模就珂以了。

3.让我们来看看新的操作:区间复制:根据zh_dou大佬的讲解,我们知道过期的迭代器是不能使用的,否则直接原地起飞。所以我们的做法就是将我们需要复制的区间先存起来,然后再删除了对应的区间后再将区间复制过去。

void fuzhi(int l1,int r1,int l2,int r2)
{
	IT it1r=split(r1+1),it1l=split(l1);
	int len=0;
	for(IT it=it1l;it!=it1r;++it)
	{
		a[++len].l=it->l;
		a[len].r=it->r;
		a[len].val=it->val;
	}
	IT it2r=split(r2+1),it2l=split(l2);
	s.erase(it2l,it2r);
	for(int i=1;i<=len;++i)
	{
		s.insert(node(a[i].l - l1 + l2,a[i].r - l1 + l2,a[i].val));
	}
}

4.新操作:区间翻转:还是根据迭代器的一套理论,我们应该先存起来在进行操作:

void fanzhuan(int l,int r)
{
	if(l>r)swap(l,r);
	IT it2=split(r+1),it1=split(l);
	int len=0;
	for(IT it=it1;it!=it2;++it)
	{
		a[++len].l=it->l;
		a[len].r=it->r;
		a[len].val=it->val;
	}
	s.erase(it1,it2);
	for(int i=1;i<=len;++i)
	{
		s.insert(node(r-a[i].r+l, r-a[i].l+l, a[i].val));
	}
}

5.毒瘤交换操作:这可能是最毒瘤的一种操作了,RE到我直接自闭。我们可以想像一下我们交换两个数的情景,就会发现这也同样适用于区间的交换,为了 好调bug 简洁明了,我用了两个数组来完成:

void my_swap(int l1,int r1,int l2,int r2)
{
	if(l1>l2){swap(l1,l2);swap(r1,r2);}
	int len1=0,len2=0;
	IT it1r=split(r1+1),it1l=split(l1);
	for(IT it=it1l;it!=it1r;++it)
	{
		a[++len1].l=it->l;
		a[len1].r=it->r;
		a[len1].val=it->val;
	}
	IT it2r=split(r2+1),it2l=split(l2);
	for(IT it=it2l;it!=it2r;++it)
	{
		b[++len2].l=it->l;
		b[len2].r=it->r;
		b[len2].val=it->val;
	}
	it1r=split(r1+1),it1l=split(l1);
	s.erase(it1l,it1r);
	it2r=split(r2+1),it2l=split(l2);
	s.erase(it2l,it2r);
	for(int i=1;i<=len2;++i)s.insert(node(b[i].l - l2 + l1,b[i].r - l2 + l1,b[i].val));
	for(int i=1;i<=len1;++i)s.insert(node(a[i].l - l1 + l2,a[i].r - l1 + l2,a[i].val));
}

然后每个数都要输出的话只用遍历一遍set就珂以啦~(当然你也珂以用来中间调试用)。

void pr()
{
	for(IT it=s.begin();it!=s.end()&&it->r<=n;++it)
	{
		for(int i=it->l;i<=it->r;++i)printf("%lld ",it->val);
	}
}

最后说几点需要注意的事情

1.过期的迭代器千万不要使用,否则自己怎么RE的都不知道,具体的操作就应该是在每次使用迭代器前就应该split一下,防止出锅。

2.遍历的时候最好再开一个迭代器千万不要这样写:

s.erase(it1,it2);
for(;it1!=it2;++it1)
{
	s.insert(node(XXX,XXX,XXX));
}

不要问我为什么我知道(TAT).

3.要开long long,并且还需要及时取模。

4.函数传参不要传错。

5.最好在n+1的位置上加上一个数(珂以是0、233或者666,随你的心意啦)。

6.需要开O2,否则就死掉了QAQ

最后献上我的完整代码

#include<iostream>
#include<cstdio>
#include<set>
#define LL long long
#define IT set<node>::iterator
using namespace std;
int n,q,opt,l1,r1,l2,r2,l,r;
LL val,x;
const int mod=1e9+7,N=500010;
inline int read()
{
	register int x=0,y=1;register char ch=getchar();
	while(ch>'9'||ch<'0'){if(ch=='-')y=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+(ch^48);ch=getchar();}
	return x*y;
}
struct node
{
	int l,r;
	mutable long long val;
	node(int L=0,int R=-1,int V=0):l(L),r(R),val(V){}
	friend bool operator <(const node &a,const node &b){return a.l<b.l;} 
//	int len(){return r-l+1;}
};
node a[N],b[N];
set<node>s;
IT split(int pos)
{
	IT it=s.lower_bound(node(pos));
	if(it!=s.end()&&it->l==pos)return it;
	it--;
	int ll=it->l,lr=it->r;
	LL lv=it->val;
	s.erase(it);
	s.insert(node(ll,pos-1,lv));
	return s.insert(node(pos,lr,lv)).first;
}
LL ask(int l,int r)
{
	IT it2=split(r+1),it1=split(l);
	LL ans=0;
	for(IT it=it1;it!=it2;++it)
		(ans+=(LL)(it->r - it->l + 1)*it->val)%=mod;
	return ans;
}
void add(int l,int r,LL val)
{
	IT it2=split(r+1),it1=split(l);
	for(IT it=it1;it!=it2;++it)(it->val+=val)%=mod;
}
void tuiping(int l,int r,LL val)
{
	IT it2=split(r+1),it1=split(l);
	s.erase(it1,it2);
	s.insert(node(l,r,val));
}
void fuzhi(int l1,int r1,int l2,int r2)
{
	IT it1r=split(r1+1),it1l=split(l1);
	int len=0;
	for(IT it=it1l;it!=it1r;++it)
	{
		a[++len].l=it->l;
		a[len].r=it->r;
		a[len].val=it->val;
	}
	IT it2r=split(r2+1),it2l=split(l2);
	s.erase(it2l,it2r);
	for(int i=1;i<=len;++i)
	{
		s.insert(node(a[i].l - l1 + l2,a[i].r - l1 + l2,a[i].val));
	}
}
void pr()
{
	for(IT it=s.begin();it!=s.end()&&it->r<=n;++it)
	{
		for(int i=it->l;i<=it->r;++i)printf("%lld ",it->val);
	}
}
void my_swap(int l1,int r1,int l2,int r2)
{
	if(l1>l2){swap(l1,l2);swap(r1,r2);}
	int len1=0,len2=0;
	IT it1r=split(r1+1),it1l=split(l1);
	for(IT it=it1l;it!=it1r;++it)
	{
		a[++len1].l=it->l;
		a[len1].r=it->r;
		a[len1].val=it->val;
	}
	IT it2r=split(r2+1),it2l=split(l2);
	for(IT it=it2l;it!=it2r;++it)
	{
		b[++len2].l=it->l;
		b[len2].r=it->r;
		b[len2].val=it->val;
	}
	it1r=split(r1+1),it1l=split(l1);
	s.erase(it1l,it1r);
	it2r=split(r2+1),it2l=split(l2);
	s.erase(it2l,it2r);
	for(int i=1;i<=len2;++i)s.insert(node(b[i].l - l2 + l1,b[i].r - l2 + l1,b[i].val));
	for(int i=1;i<=len1;++i)s.insert(node(a[i].l - l1 + l2,a[i].r - l1 + l2,a[i].val));
}	
void fanzhuan(int l,int r)
{
	if(l>r)swap(l,r);
	IT it2=split(r+1),it1=split(l);
	int len=0;
	for(IT it=it1;it!=it2;++it)
	{
		a[++len].l=it->l;
		a[len].r=it->r;
		a[len].val=it->val;
	}
	s.erase(it1,it2);
	for(int i=1;i<=len;++i)
	{
		s.insert(node(r-a[i].r+l, r-a[i].l+l, a[i].val));
	}
}
int main()
{
	n=read(),q=read();
	for(int i=1;i<=n;++i)
	{
		scanf("%lld",&x);
		s.insert(node(i,i,x));
	}
	s.insert((node){n+1,n+1,(LL)0});
	while(q--)
	{
		opt=read();
		if(opt==1)
		{
			l=read();r=read();
			printf("%lld
",ask(l,r));
		}
		if(opt==2)
		{
			l=read();r=read();scanf("%lld",&val);
			tuiping(l,r,val);
		}
		if(opt==3)
		{
			l=read();r=read();scanf("%lld",&val);
			add(l,r,val);
		}
		if(opt==4)
		{
			l1=read();r1=read();l2=read();r2=read();
			fuzhi(l1,r1,l2,r2);
		}
		if(opt==5)
		{
			l1=read();r1=read();l2=read();r2=read();
			my_swap(l1,r1,l2,r2);
		}
		if(opt==6)
		{
			l=read();r=read();
			fanzhuan(l,r);
		}
	}
	pr();
	return 0;
}
原文地址:https://www.cnblogs.com/wljss/p/11553656.html