| WZJ的数据结构（负三十三） |
| 难度级别：C； 运行时间限制：7000ms； 运行空间限制：262144KB； 代码长度限制：2000000B |
试题描述 |
请你设计一个数据结构，完成以下功能： 给定一个大小为N的整数组A，要求你回答执行N次操作。操作分两种： 操作1：每次操作给你l,r,v三个参数，求Al至Ar中值<=v的个数。 操作2：每次操作给你l,r,v三个参数，将Al至Ar所有数的值设为v。 |
输入 |
第一行为一个正整数N。 第二行为N个整数Ai。 接下来N行每行4个正整数t,l,r,v。若t=2表示操作1，t=1表示操作2。 |
输出 |
对每个操作1输出结果。 |
输入示例 |
6 1 2 2 1 2 3 2 1 4 2 2 1 4 1 1 2 4 1 2 1 4 1 1 3 5 2 2 2 5 2 |
输出示例 |
4 2 4 4 |
其他说明 |
1<=N<=200000 1<=Ai,v<=10^9 |
XYZ大爷论文题“一个简单的序列问题”。
一上午的青春全部奉献给这道题了（3棵线段树调得我要报警）。
先将区间修改转化成带删除的区间插入，用一棵线段树来做。
再用线段树分治将问题转化成静态问题，用一棵线段树和2个nlogn的邻接表。
最后排序后转化成区间加与区间查询的问题，用一棵线段树来做。
说起来真是简单呢！
#include<cstdio>
#include<cctype>
#include<queue>
#include<cmath>
#include<cstring>
#include<algorithm>
#define rep(i,s,t) for(int i=s;i<=t;i++)
#define dwn(i,s,t) for(int i=s;i>=t;i--)
#define ren for(int i=first[x];i!=-1;i=next[i])
using namespace std;
inline int read() {
int x=,f=;char c=getchar();
for(;!isdigit(c);c=getchar()) if(c=='-') f=-;
for(;isdigit(c);c=getchar()) x=x*+c-'';
return x*f;
}
const int maxn=;
struct Query {
int tp,l,r,x,begin,end,id;
bool operator < (const Query& ths) const {return begin<ths.begin;}
}Q[maxn],Q2[maxn*];
struct Solver {
int id,l,r,x;
bool operator < (const Solver& ths) const {return x<ths.x||(x==ths.x&&id<ths.id);}
}A[maxn*];
int n,ans[maxn],m;
int setv[maxn*],is[maxn*];
void pushdown(int o) {
int lc=o<<,rc=lc|;
if(setv[o]) {
setv[lc]=setv[rc]=setv[o];
is[lc]=is[rc]=;
setv[o]=;
}
}
int lastv,L,R;
void find(int o,int l,int r,int v) {
if(!is[o]) return;
if(setv[o]) {
if(setv[o]==lastv) R=r;
else {
if(lastv!=-) Q2[++m]=(Query){,L,R,Q[lastv].x,lastv,v-,Q[lastv].id};
lastv=setv[o];L=l;R=r;
}
setv[o]=;
}
else if(l<r) {
int mid=l+r>>,lc=o<<,rc=lc|;
find(lc,l,mid,v);find(rc,mid+,r,v);
}
is[o]=;
}
void update(int o,int l,int r,int ql,int qr,int v) {
if(ql<=l&&r<=qr) find(o,l,r,v),setv[o]=v;
else {
pushdown(o);
int mid=l+r>>,lc=o<<,rc=lc|;
if(ql<=mid) update(lc,l,mid,ql,qr,v);
if(qr>mid) update(rc,mid+,r,ql,qr,v);
}
is[o]=;
}
int first1[maxn*],next1[maxn*],id1[maxn*],ToT1;
int first2[maxn*],next2[maxn*],id2[maxn*],ToT2,N;
void query(int o,int l,int r,int pos,int v) {
N=max(N,o);
id1[++ToT1]=v;next1[ToT1]=first1[o];first1[o]=ToT1;
if(l==r) return;
int mid=l+r>>,lc=o<<,rc=lc|;
if(pos<=mid) query(lc,l,mid,pos,v);
else query(rc,mid+,r,pos,v);
}
void modify(int o,int l,int r,int ql,int qr,int v) {
N=max(N,o);
if(ql<=l&&r<=qr) {
id2[++ToT2]=v;next2[ToT2]=first2[o];first2[o]=ToT2;
}
else {
int mid=l+r>>,lc=o<<,rc=lc|;
if(ql<=mid) modify(lc,l,mid,ql,qr,v);
if(qr>mid) modify(rc,mid+,r,ql,qr,v);
}
}
int addv[maxn*],sumv[maxn*];
void pushdown2(int o,int l,int r) {
int mid=l+r>>,lc=o<<,rc=lc|;
if(addv[o]) {
addv[lc]+=addv[o];addv[rc]+=addv[o];
sumv[lc]+=addv[o]*(mid-l+);sumv[rc]+=addv[o]*(r-mid);
addv[o]=;
}
}
void update2(int o,int l,int r,int ql,int qr,int v) {
int mid=l+r>>,lc=o<<,rc=lc|;
if(ql<=l&&r<=qr) addv[o]+=v;
else {
pushdown2(o,l,r);
if(ql<=mid) update2(lc,l,mid,ql,qr,v);
if(qr>mid) update2(rc,mid+,r,ql,qr,v);
}
if(l<r) sumv[o]=sumv[lc]+sumv[rc];
else sumv[o]=;
if(addv[o]) sumv[o]+=addv[o]*(r-l+);
}
int query2(int o,int l,int r,int ql,int qr) {
if(ql<=l&&r<=qr) return sumv[o];
pushdown2(o,l,r);
int mid=l+r>>,lc=o<<,rc=lc|,ans=;
if(ql<=mid) ans+=query2(lc,l,mid,ql,qr);
if(qr>mid) ans+=query2(rc,mid+,r,ql,qr);
return ans;
}
int tmp[maxn];
int main() {
n=read();
rep(i,,n) Q[Q[i].id=i].tp=,Q[i].l=i,Q[i].r=i,Q[i].x=read();
rep(i,n+,*n) Q[Q[i].id=i].tp=read(),Q[i].l=read(),Q[i].r=read(),Q[i].x=read();
n<<=;
Q[n+].tp=;Q[n+].l=;Q[n+].r=n>>;
rep(i,,n+) {
if(Q[i].tp==) {
lastv=-;
update(,,n>>,Q[i].l,Q[i].r,i);
if(lastv!=-) Q2[++m]=(Query){,L,R,Q[lastv].x,lastv,i-,Q[lastv].id};
}
else Q2[++m]=Q[i],Q2[m].begin=Q2[m].end=i;
}
rep(i,,m) {
if(Q2[i].tp==) modify(,,n,Q2[i].begin,Q2[i].end,i);
else query(,,n,Q2[i].begin,i);
}
rep(i,,N) if(first1[i]&&first2[i]) {
int cnt=;
for(int j=first1[i];j;j=next1[j]) A[++cnt]=(Solver){Q2[id1[j]].id,Q2[id1[j]].l,Q2[id1[j]].r,Q2[id1[j]].x};
for(int j=first2[i];j;j=next2[j]) A[++cnt]=(Solver){,Q2[id2[j]].l,Q2[id2[j]].r,Q2[id2[j]].x};
sort(A+,A+cnt+);
rep(j,,cnt) {
if(!A[j].id) update2(,,n>>,A[j].l,A[j].r,);
else ans[A[j].id]+=query2(,,n>>,A[j].l,A[j].r);
}
rep(j,,cnt) if(!A[j].id) update2(,,n>>,A[j].l,A[j].r,-);
}
rep(i,,n) if(Q[i].tp==) printf("%d\n",ans[i]);
return ;
}
8.6重写了一遍，30min写完，WA了两处233.
错误：
void pushdown(int o) {
int lc=o<<,rc=lc|;
if(setv[o]) {
setv[lc]=setv[rc]=;//setv[o]
is[lc]=is[rc]=;
setv[o]=;
}
}
int query(int o,int l,int r,int ql,int qr) {
if(ql<=l&&r<=qr) return sumv[o];
//pushdown(o,l,r);
int mid=l+r>>,lc=o<<,rc=lc|,res=;
if(ql<=mid) res+=query(lc,l,mid,ql,qr);
if(qr>mid) res+=query(rc,mid+,r,ql,qr);
return res;
}8.10啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊标程被分块虐了！！！！！！！！！！！！！！！！！！！！！！！！！！！！
#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<queue>
#include<cstring>
#define PAU putchar(' ')
#define ENT putchar('\n')
using namespace std;
const int maxn=+,maxb=,inf=-1u>>;
int z[maxn],st[maxb],en[maxb],size,n,A[maxn],B[maxn],set[maxb],siz[maxb];
inline int read(){
int x=,sig=;char ch=getchar();
for(;!isdigit(ch);ch=getchar())if(ch=='-')sig=;
for(;isdigit(ch);ch=getchar())x=*x+ch-'';
return sig?x:-x;
}
inline void write(int x){
if(x==){putchar('');return;}if(x<)putchar('-'),x=-x;
int len=,buf[];while(x)buf[len++]=x%,x/=;
for(int i=len-;i>=;i--)putchar(buf[i]+'');return;
}
void build(int b){
for(int i=st[b];i<=en[b];i++)B[i]=A[i];sort(B+st[b],B+en[b]+);return;
}
void down(int b){
if(set[b]!=inf){for(int i=st[b];i<=en[b];i++)A[i]=B[i]=set[b];set[b]=inf;}return;
}
void sett(int L,int R,int cv){
down(z[L]);for(int i=L;i<=R;i++)A[i]=cv;build(z[L]);return;
}
void setb(int L,int R,int cv){
for(int b=L;b<=R;b++)set[b]=cv;return;
}
void seto(int ql,int qr,int cv){
if(z[ql]==z[qr])sett(ql,qr,cv);
else sett(ql,en[z[ql]],cv),setb(z[ql]+,z[qr]-,cv),sett(st[z[qr]],qr,cv);return;
}
int queryb(int L,int R,int cv){
int ans=;
for(int b=L;b<=R;b++){
if(set[b]!=inf){
if(set[b]<=cv)ans+=siz[b];
}else ans+=upper_bound(B+st[b],B+en[b]+,cv)-B-st[b];
}return ans;
}
int queryt(int L,int R,int cv){
int ans=;
if(set[z[L]]!=inf){
if(set[z[L]]<=cv)ans+=(R-L+);
}else for(int i=L;i<=R;i++)if(A[i]<=cv)ans++;
return ans;
}
int queryo(int ql,int qr,int cv){
if(z[ql]==z[qr])return queryt(ql,qr,cv);
else return queryt(ql,en[z[ql]],cv)+queryb(z[ql]+,z[qr]-,cv)+queryt(st[z[qr]],qr,cv);
}
void init(){
n=read();size=sqrt((double)n*0.9);int tp,ql,qr,cv;
for(int i=;i<=n;i++){
A[i]=B[i]=read();
z[i]=(i-)/size+;
if(!st[z[i]])st[z[i]]=i;
en[z[i]]=i;
}
for(int b=;b<=z[n];b++)build(b),set[b]=inf,siz[b]=en[b]-st[b]+;
for(int i=;i<=n;i++){
tp=read();ql=read();qr=read();cv=read();
if(tp==)write(queryo(ql,qr,cv)),ENT;
else seto(ql,qr,cv);
}
return;
}
void work(){
return;
}
void print(){
return;
}
int main(){init();work();print();return ;}