【树套树】【bzoj 3236】: [Ahoi2013]作业

http://www.lydsy.com/JudgeOnline/problem.php?id=3236

BIT套SBT

Treap没救了，无论怎么优化常数都过不了，于是乎随便改成了SBT。。。。。。

开始SBT的sz我用的s维护结果狂RE。。。。。。

最后90s卡过。。。。。比莫队慢。。。。。。。

被虐了2h。。。。。。

早知道写个莫队。。。。。。

莫队明明是O(nsqrt(n)logn)我的是O(nlognlogn)为什么会慢想不通啊。。。。。。。。。。。

一定是我太弱了。。。。。。。。
Orz 5s过的神犇。。。。。。。。

。。。。。。。。。。。。。。。。。。。。

以后有时间试一下http://hi.baidu.com/greencloud/item/0e25878061444bde5e0ec1f7

//#define _TEST _TEST#include <cstdio>#include <cstring>#include <cstdlib>#include <iostream>#include <cmath>#include <algorithm>using namespace std;/************************************************Code By willinglive    Blog:http://willinglive.cf************************************************/#define rep(i,l,r) for(int i=l,___t=(r);i<=___t;i++)#define per(i,r,l) for(int i=r,___t=(l);i>=___t;i--)#define MS(arr,x) memset(arr,x,sizeof(arr))#define LL long long#define INE(i,u,e) for(int i=head[u];~i;i=e[i].next)inline const int getint(){    int r=0,k=1;char c=getchar();    for(;c<'0'||c>'9';c=getchar())if(c=='-')k=-1;    for(;c>='0'&&c<='9';c=getchar())r=r*10+c-'0';    return k*r;}/////////////////////////////////////////////////int n,m;int A[100010];inline int rnd(){return (rand()<<16)|rand();}struct data{int id,l,r,a,b;}q[1000010];bool cmp(data a,data b){return a.r<b.r;}int ans1[1000010],ans2[1000010];int last[100010];/////////////////////////////////////////////////namespace SBT{/////////////////////////////////#define LS T[o].l#define RS T[o].rint root[100010];struct data{int l,r,s,sz,cnt,x,w,sum;}T[1000010];int sz;inline void update(int o){T[o].s=T[LS].s+T[RS].s+T[o].cnt;T[o].sz=T[LS].sz+T[RS].sz+1;T[o].sum=T[LS].sum+T[RS].sum+T[o].w;}inline void l_rot(int &o){int t=RS;RS=T[t].l;T[t].l=o;T[t].s=T[o].s;T[t].sz=T[o].sz;T[t].sum=T[o].sum;update(o);o=t;}inline void r_rot(int &o){int t=LS;LS=T[t].r;T[t].r=o;T[t].s=T[o].s;T[t].sz=T[o].sz;T[t].sum=T[o].sum;update(o);o=t;}inline void maintain(int &o,bool b){if(!b){if(!LS) return;if(T[T[LS].l].sz>T[RS].sz) r_rot(o);else if(T[T[LS].r].sz>T[RS].sz) l_rot(LS),r_rot(o);else return;}else{if(!RS) return;if(T[T[RS].r].sz>T[LS].sz) l_rot(o);else if(T[T[RS].l].sz>T[LS].sz) r_rot(RS),l_rot(o);else return;}    maintain(LS,0); maintain(RS,1);    maintain(o,0);  maintain(o,1);}inline void insert(int &o,int x){if(o==0){o=++sz; T[o].sz=T[o].s=T[o].cnt=1; T[o].x=x; T[o].w=0; T[o].sum=0;return;}T[o].s++; T[o].sz++;if(x<T[o].x){insert(LS,x);maintain(o,0);}else if(x>T[o].x){insert(RS,x);maintain(o,1);}else T[o].cnt++;}inline int query1(int o,int x){int res=0;while(o){if(x<T[o].x) o=LS;else if(x>T[o].x) res+=T[LS].s+T[o].cnt,o=RS;else return res+T[LS].s+T[o].cnt;}return res;}inline int query2(int o,int x){int res=0;while(o){if(x<T[o].x) res+=T[RS].s+T[o].cnt,o=LS;else if(x>T[o].x) o=RS;else return res+T[RS].s+T[o].cnt;}return res;}inline void add(int o,int clr,int x){while(o){T[o].sum+=x;if(clr<T[o].x) o=LS;else if(clr>T[o].x) o=RS;else T[o].w+=x,o=0;}}inline int query3(int o,int x){int res=0;while(o){if(x<T[o].x) o=LS;else if(x>T[o].x) res+=T[LS].sum+T[o].w,o=RS;else return res+T[LS].sum+T[o].w;}return res;}inline int query4(int o,int x){int res=0;while(o){if(x<T[o].x) res+=T[RS].sum+T[o].w,o=LS;else if(x>T[o].x) o=RS;else return res+T[RS].sum+T[o].w;}return res;}}/////////////////////////////////namespace BIT{/////////////////////////////////inline void insert(int o,int x){for(;o<=n;o+=o&-o)    SBT::insert(SBT::root[o],x);}inline int query1(int o,int a,int b){using namespace SBT;int s=0;for(;o>0;o-=o&-o)    s+=SBT::query2(root[o],a) + SBT::query1(root[o],b);    return s;}inline int query2(int o,int a,int b){using namespace SBT;int s=0;for(;o>0;o-=o&-o)    s+=SBT::query3(root[o],b) + SBT::query4(root[o],a) - SBT::query3(root[o],n);    return s;}inline void add(int o,int clr,int x){for(;o<=n;o+=o&-o)    SBT::add(SBT::root[o],clr,x);}}//////////////////////////////////////////////////////////////////////////////////void input(){using namespace BIT;    scanf("%d%d",&n,&m);    int x;    rep(i,1,n)    insert(i,A[i]=getint());    rep(i,1,m)    q[i].id=i, q[i].l=getint(), q[i].r=getint(), q[i].a=getint(), q[i].b=getint();   sort(&q[1],&q[m+1],cmp);}void solve(){using namespace BIT;int cur=1;int id,l,r,a,b;rep(i,1,n){int &lst=last[A[i]];        if(lst==0) add(i,A[i],1);        else add(lst,A[lst],-1), add(i,A[i],1);        last[A[i]]=i;        for(;q[cur].r==i && cur<=m;cur++){id=q[cur].id; l=q[cur].l-1; r=q[cur].r; a=q[cur].a; b=q[cur].b;int s1=query1(r,a,b)-r;int s2=query1(l,a,b)-l;            ans1[id]=s1-s2;                        s1=query2(r,a,b);            s2=query2(l,a,b);            ans2[id]=s1-s2;}}rep(i,1,m) printf("%d %d\n",ans1[i],ans2[i]);}/////////////////////////////////////////////////int main(){    #ifndef _TEST    freopen("std.in","r",stdin); freopen("std.out","w",stdout);    #endif    input(),    solve();    return 0;}