这题有很多算法,我已知:
分块+二分:O(n*sqrt(n*log(n))
函数式权值分块:O(n*sqrt(n))
带修莫队+权值分块:O(n5/3)
树套树
我打的莫队+分块,跑的死慢
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#include<cstring>
using namespace std;
inline char nc()
{
static char buf[100000],*p1=buf,*p2=buf;
if (p1==p2) { p2=(p1=buf)+fread(buf,1,100000,stdin); if (p1==p2) return EOF; }
return *p1++;
}
inline void read(int &x)
{
char c=nc(),b=1;
for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1;
for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b;
}
namespace Block{
int block,n,tot;
int a[150005],pos[150005];
int sum[1005];
inline void init(int _n){
n=_n;
block=pow(n,2.0/3.0)+1;
for (int i=1;i<=n;i++) pos[i]=(i-1)/block+1;
tot=pos[n];
}
inline void add(int x,int r){
a[x]+=r; sum[pos[x]]+=r;
}
inline int kth(int k){
int i,j;
for (i=1;i<=tot && k>sum[i];i++)
k-=sum[i];
for (j=(i-1)*block+1;j<=min(i*block,n) && k>a[j];j++)
k-=a[j];
return j;
}
inline int rank(int k){
int cnt=0;
for (int i=1;i<=pos[k]-1;i++)
cnt+=sum[i];
for (int i=(pos[k]-1)*block+1;i<k;i++)
cnt+=a[i];
return cnt+1;
}
inline int nxt(int k){
return kth(rank(k)+a[k]);
}
inline int pre(int k){
return kth(rank(k)-1);
}
}
struct Change{
int pos,x,y;
}C[50005];
int clk;
struct Query{
int f;
int l,r,x;
int pos,idx,time;
bool operator < (const Query &B) const{
return pos==B.pos?r<B.r:pos<B.pos;
}
}Q[50005];
int tot;
int n,m,B;
int a[50005],_a[50005];
int sx[150005],icnt;
int ans[100005];
inline int Bin(int x){
return lower_bound(sx+1,sx+icnt+1,x)-sx;
}
inline void deal(int id)
{
if (Q[id].f==1)
ans[Q[id].idx]=Block::rank(Q[id].x);
else if (Q[id].f==2)
ans[Q[id].idx]=sx[Block::kth(Q[id].x)];
else if (Q[id].f==4)
ans[Q[id].idx]=sx[Block::pre(Q[id].x)];
else if (Q[id].f==5)
ans[Q[id].idx]=sx[Block::nxt(Q[id].x)];
}
inline void Solve()
{
for (int i=1;i<=n;i++)
_a[i]=Bin(a[i]);
int l=1,r=0;
for (int i=1;i<=tot;i++)
{
for (int j=Q[i-1].time+1;j<=Q[i].time;j++)
{
_a[C[j].pos]=C[j].y;
if (C[j].pos>=l && C[j].pos<=r)
Block::add(C[j].x,-1),Block::add(C[j].y,1);
}
for (int j=Q[i-1].time;j>Q[i].time;j--)
{
_a[C[j].pos]=C[j].x;
if (C[j].pos>=l && C[j].pos<=r)
Block::add(C[j].y,-1),Block::add(C[j].x,1);
}
for (;l<Q[i].l;l++)
Block::add(_a[l],-1);
for (;l>Q[i].l;l--)
Block::add(_a[l-1],1);
for (;r<Q[i].r;r++)
Block::add(_a[r+1],1);
for (;r>Q[i].r;r--)
Block::add(_a[r],-1);
deal(i);
}
}
int main()
{
int order;
freopen("t.in","r",stdin);
freopen("t.out","w",stdout);
read(n); read(m);
B=pow(n,2.0/3.0)+1;
for (int i=1;i<=n;i++)
read(a[i]),sx[++icnt]=a[i],_a[i]=a[i];
for (int i=1;i<=m;i++)
{
read(order);
if (order==1)
{
Q[++tot].f=order; Q[tot].idx=tot; Q[tot].time=clk;
read(Q[tot].l); read(Q[tot].r); read(Q[tot].x);
sx[++icnt]=Q[tot].x;
}
else if (order==2)
{
Q[++tot].f=order; Q[tot].idx=tot; Q[tot].time=clk;
read(Q[tot].l); read(Q[tot].r); read(Q[tot].x);
}
else if (order==3)
{
++clk;
read(C[clk].pos); read(C[clk].y); C[clk].x=_a[C[clk].pos]; _a[C[clk].pos]=C[clk].y;
sx[++icnt]=C[clk].y;
}
else if (order==4){
Q[++tot].f=order; Q[tot].idx=tot; Q[tot].time=clk;
read(Q[tot].l); read(Q[tot].r); read(Q[tot].x);
sx[++icnt]=Q[tot].x;
}
else if (order==5){
Q[++tot].f=order; Q[tot].idx=tot; Q[tot].time=clk;
read(Q[tot].l); read(Q[tot].r); read(Q[tot].x);
sx[++icnt]=Q[tot].x;
}
}
sort(sx+1,sx+icnt+1);
icnt=unique(sx+1,sx+icnt+1)-sx-1;
Block::init(icnt);
for (int i=1;i<=tot;i++)
{
if (Q[i].f!=2)
Q[i].x=Bin(Q[i].x);
Q[i].pos=(Q[i].l-1)/B+1;
}
sort(Q+1,Q+tot+1);
for (int i=1;i<=clk;i++)
C[i].x=Bin(C[i].x),C[i].y=Bin(C[i].y);
Solve();
for (int i=1;i<=tot;i++)
printf("%d\n",ans[i]);
return 0;
}