首页 > 代码库 > SPOJ TTM
SPOJ TTM
这道题一眼看去就是一个可持久化线段树,但是是区间修改,由于wyx说此题复杂度是O(nlogn)的,我就没写树套树,然后就自己yy了一个离线做法。
我们考虑直接模拟这个过程,对于一个B操作,我们直接将之前的操作的影响清除,这样每个操作最多会被计算2次,但是有一个问题就是H操作不太好弄,我们可以先离线下来,然后维护一下H操作的这个时间t被经过了多少次,这样我们就知道这个H操作是我们第几次模拟后的查询,这样就可以通过一个普通的线段树来维护了。
这道题好像正解是区间修改的主席树,我之前只会写套上树状数组的,不过感觉这个没什么用...就先不学了吧..
1 #include<cstdio> 2 #include<cstring> 3 #include<algorithm> 4 #include<vector> 5 using namespace std; 6 #define maxn 100005 7 typedef long long LL; 8 struct tree 9 { 10 int l,r; 11 LL sum,lazy; 12 }t[maxn*4]; 13 struct ask 14 { 15 int opt,l,r,d,t; 16 }q[maxn]; 17 int n,m,T,cnt[maxn],st[maxn],top,cur[maxn],a[maxn]; 18 LL ans[maxn]; 19 vector<pair<int,int> > g[maxn]; 20 21 inline int read(void) 22 { 23 int x=0,f=1; 24 char ch=getchar(); 25 while (ch>‘9‘||ch<‘0‘) 26 { 27 if (ch==‘-‘) f=-1; 28 ch=getchar(); 29 } 30 while (ch>=‘0‘&&ch<=‘9‘) 31 { 32 x=x*10+ch-‘0‘; 33 ch=getchar(); 34 } 35 return x*f; 36 } 37 38 void build(int x,int l,int r) 39 { 40 t[x].l=l;t[x].r=r; 41 if (l==r) 42 { 43 t[x].sum=a[l]; 44 return; 45 } 46 int mid=(l+r)>>1; 47 build(2*x,l,mid); 48 build(2*x+1,mid+1,r); 49 t[x].sum=t[2*x].sum+t[2*x+1].sum; 50 } 51 52 void update(int x) 53 { 54 t[2*x].sum+=t[x].lazy*(t[2*x].r-t[2*x].l+1); 55 t[2*x+1].sum+=t[x].lazy*(t[2*x+1].r-t[2*x+1].l+1); 56 t[2*x].lazy+=t[x].lazy; 57 t[2*x+1].lazy+=t[x].lazy; 58 t[x].lazy=0; 59 } 60 61 void change(int x,int l,int r,int z) 62 { 63 if (t[x].l==l&&t[x].r==r) 64 { 65 t[x].sum+=(LL)z*(r-l+1); 66 t[x].lazy+=z; 67 return; 68 } 69 if (t[x].lazy) update(x); 70 int mid=(t[x].l+t[x].r)>>1; 71 if (l>mid) change(2*x+1,l,r,z); 72 else if (r<=mid) change(2*x,l,r,z); 73 else 74 { 75 change(2*x,l,mid,z); 76 change(2*x+1,mid+1,r,z); 77 } 78 t[x].sum=t[2*x].sum+t[2*x+1].sum; 79 } 80 81 LL query(int x,int l,int r) 82 { 83 if (t[x].l==l&&t[x].r==r) return t[x].sum; 84 if (t[x].lazy) update(x); 85 int mid=(t[x].l+t[x].r)>>1; 86 if (l>mid) return query(2*x+1,l,r); 87 else if (r<=mid) return query(2*x,l,r); 88 else return query(2*x,l,mid)+query(2*x+1,mid+1,r); 89 } 90 91 int main() 92 { 93 n=read();m=read(); 94 for (int i=1;i<=n;i++) a[i]=read(); 95 for (int i=1;i<=m;i++) 96 { 97 char s[3]; 98 scanf("%s",s); 99 if (s[0]==‘C‘) 100 { 101 q[i].opt=1;cnt[++T]++; 102 q[i].l=read();q[i].r=read();q[i].d=read(); 103 q[i].t=T; 104 } 105 if (s[0]==‘Q‘) 106 { 107 q[i].opt=2;q[i].l=read(); 108 q[i].r=read();g[T].push_back(make_pair(cnt[T],i)); 109 q[i].t=T; 110 } 111 if (s[0]==‘H‘) 112 { 113 q[i].opt=3;q[i].l=read();q[i].r=read(); 114 int t=read(); 115 g[t].push_back(make_pair(cnt[t],i)); 116 q[i].t=t; 117 } 118 if (s[0]==‘B‘) 119 { 120 q[i].opt=4;q[i].d=read();q[i].t=T; 121 T=q[i].d; 122 } 123 } 124 build(1,1,n);T=0; 125 for (int i=1;i<=m;i++) 126 if ((q[i].opt==3||q[i].opt==2)&&q[i].t==0) 127 ans[i]=query(1,q[i].l,q[i].r); 128 memset(cnt,0,sizeof cnt); 129 for (int i=1;i<=m;i++) 130 { 131 if (q[i].opt==2||q[i].opt==3) continue; 132 if (q[i].opt==1) 133 { 134 cnt[++T]++; 135 change(1,q[i].l,q[i].r,q[i].d); 136 st[++top]=i; 137 int siz=g[T].size(); 138 while (cur[T]<siz&&cnt[T]==g[T][cur[T]].first) 139 { 140 int id=g[T][cur[T]].second; 141 ans[id]=query(1,q[id].l,q[id].r); 142 cur[T]++; 143 } 144 } 145 if (q[i].opt==4) 146 { 147 while (top>0&&q[st[top]].t>q[i].d) 148 change(1,q[st[top]].l,q[st[top]].r,-q[st[top]].d),top--; 149 T=q[i].d; 150 } 151 } 152 for (int i=1;i<=m;i++) 153 if (q[i].opt==2||q[i].opt==3) printf("%lld\n",ans[i]); 154 return 0; 155 }
SPOJ TTM
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。