CodeForces - 897E Willem, Chtholly and Seniorious(珂朵莉树)

题目链接:点击查看

题目大意:给出一个长度为 n n n 的数列,现在需要执行 m m m 次操作,每次操作分为下列四种情况:

  1. 1   l   r   x 1 \ l \ r \ x 1 l r x: [ l , r ] [l,r] [l,r] 内的位置都加上 x x x
  2. 2   l   r   x 2 \ l \ r \ x 2 l r x: [ l , r ] [l,r] [l,r] 内的数都变为 x x x
  3. 3   l   r   x 3 \ l \ r \ x 3 l r x:求 [ l , r ] [l,r] [l,r] 内第 x x x 小的数
  4. 4   l   r   x   y 4 \ l \ r \ x \ y 4 l r x y:求区间 [ l , r ] [l,r] [l,r] 内的幂次和,对 y y y 取模

数据随机

题目分析:大佬博客

因为都是随机数据,所以直接上珂朵莉树推平就好了

有个奇怪的点是网上讲解珂朵莉树的博客都会加入一个哨兵节点,但我个人感觉是不许要加的,因为所有的运算都是针对于超尾去实现的,而 s t d : : s e t std::set std::set 显然是自带超尾的,所以无需加入哨兵节点防止越界

需要注意的一点细节就是 a i a_i ai​ 可能比较大,在进行快速幂的时候需要先取模,不然会直接爆掉 l o n g   l o n g long \ long long long

代码:

// #pragma GCC optimize(2)
// #pragma GCC optimize("Ofast","inline","-ffast-math")
// #pragma GCC target("avx,sse2,sse3,sse4,mmx")
#include<iostream>
#include<cstdio>
#include<string>
#include<ctime>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<stack>
#include<climits>
#include<queue>
#include<map>
#include<set>
#include<sstream>
#include<cassert>
#include<bitset>
#include<unordered_map>
using namespace std;
typedef long long LL;
typedef unsigned long long ull;
template<typename T>
inline void read(T &x)
{
	T f=1;x=0;
	char ch=getchar();
	while(0==isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
	while(0!=isdigit(ch)) x=(x<<1)+(x<<3)+ch-'0',ch=getchar();
	x*=f;
}
template<typename T>
inline void write(T x)
{
	if(x<0){x=~(x-1);putchar('-');}
    if(x>9)write(x/10);
    putchar(x%10+'0');
}
const int inf=0x3f3f3f3f;
const int N=1e6+100;
struct Node
{
	int l,r;
	mutable LL val;
	Node(int l,int r,LL val):l(l),r(r),val(val){}
	Node(int l):l(l){}
	bool operator<(const Node& t)const{return l<t.l;}
};
set<Node>st;
set<Node>::iterator split(int pos)
{
	auto it=st.lower_bound(Node(pos));
	if(it!=st.end()&&it->l==pos) return it;
	it--;
	int l=it->l,r=it->r;
	LL val=it->val;
	st.erase(it);
	st.insert(Node(l,pos-1,val));
	return st.insert(Node(pos,r,val)).first;
}
void assign(int l,int r,LL val)
{
	auto itr=split(r+1),itl=split(l);
	st.erase(itl,itr);
	st.insert(Node(l,r,val));
}
void add(int l,int r,LL val)
{
	auto itr=split(r+1),itl=split(l);
	for(auto it=itl;it!=itr;it++) it->val+=val;
}
LL kth(int l,int r,int k)
{
	auto itr=split(r+1),itl=split(l);
	vector<pair<LL,int>>node;
	for(auto it=itl;it!=itr;it++) node.emplace_back(it->val,it->r-it->l+1);
	sort(node.begin(),node.end());
	for(auto it:node)
	{
		k-=it.second;
		if(k<=0) return it.first;
	}
	return -1;
}
LL q_pow(LL a,LL b,LL mod)
{
	LL ans=1;
	a%=mod;
	while(b)
	{
		if(b&1) ans=ans*a%mod;
		a=a*a%mod;
		b>>=1;
	}
	return ans;
}
LL query(int l,int r,int x,int mod)
{
	auto itr=split(r+1),itl=split(l);
	LL ans=0;
	for(auto it=itl;it!=itr;it++) ans=(ans+q_pow(it->val,x,mod)*(it->r-it->l+1))%mod;
	return ans;
}
LL seed;
int rnd()
{
	int ret=seed;
	seed=(seed*7+13)%1000000007;
	return ret;
}
int main()
{
#ifndef ONLINE_JUDGE
//  freopen("data.in.txt","r",stdin);
//  freopen("data.out.txt","w",stdout);
#endif
//  ios::sync_with_stdio(false);
	int n,m,vmax;
	read(n),read(m),read(seed),read(vmax);
	for(int i=1;i<=n;i++)
	{
		int num=rnd()%vmax+1;
		st.insert(Node(i,i,num));
	}
	for(int i=1;i<=m;i++)
	{
		int op=rnd()%4+1,l=rnd()%n+1,r=rnd()%n+1,x,y;
		if(l>r) swap(l,r);
		if(op==3) x=rnd()%(r-l+1)+1;
		else x=rnd()%vmax+1;
		if(op==4) y=rnd()%vmax+1;
		if(op==1) add(l,r,x);
		else if(op==2) assign(l,r,x);
		else if(op==3) write(kth(l,r,x)),putchar('\n');
		else if(op==4) write(query(l,r,x,y)),putchar('\n');
	}
    return 0;
}
上一篇:modbus协议之-01-初次见面


下一篇:cnpm安装教程