首先离散化颜色

设pre[x]表示与x颜色相同的点上一次出现的位置,对于每种颜色开一个set维护

修改时需要修改x、x修改前的后继、x修改后的后继

询问[l,r]等价于询问[l,r]内pre[x]<l的点的权值和

线段树套Treap维护

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<set>
#define N 100010
using namespace std;
typedef long long ll;
int n,m,i,x,y,pre,nxt,a[N],b[N<<1],C,que[N][3],loc[N<<1];char ch;ll ans;
set<int>T[N<<1];set<int>::iterator j,k;
inline void read(int&a){char c;while(!(((c=getchar())>='0')&&(c<='9')));a=c-'0';while(((c=getchar())>='0')&&(c<='9'))(a*=10)+=c-'0';}
inline int lower(int x){
int l=1,r=C,t,mid;
while(l<=r)if(b[mid=(l+r)>>1]<=x)l=(t=mid)+1;else r=mid-1;
return t;
}
struct node{
int val,p;ll v,sum;node*l,*r;
node(){val=v=sum=p=0;l=r=NULL;}
inline void up(){sum=v+l->sum+r->sum;}
}*blank=new(node),*root[N<<2];
inline void Rotatel(node*&x){node*y=x->r;x->r=y->l;x->up();y->l=x;y->up();x=y;}
inline void Rotater(node*&x){node*y=x->l;x->l=y->r;x->up();y->r=x;y->up();x=y;}
void Ins(node*&x,int p,int v){
if(x==blank){
x=new(node);x->val=p;x->l=x->r=blank;x->v=x->sum=v;x->p=std::rand();
return;
}
x->sum+=v;
if(p==x->val){x->v+=v;return;}
if(p<x->val){
Ins(x->l,p,v);
if(x->l->p>x->p)Rotater(x);
}else{
Ins(x->r,p,v);
if(x->r->p>x->p)Rotatel(x);
}
}
ll Ask(node*x,int p){
if(x==blank)return 0;
if(p==x->val)return x->l->sum;
if(p>x->val)return x->v+x->l->sum+Ask(x->r,p);
return Ask(x->l,p);
}
void build(int x,int a,int b){
root[x]=blank;
if(a==b)return;
int mid=(a+b)>>1;
build(x<<1,a,mid),build(x<<1|1,mid+1,b);
}
inline void add(int c,int d,int e){
int x=1,a=1,b=n,mid;
while(a<b){
Ins(root[x],d,e);mid=(a+b)>>1;x<<=1;
if(c<=mid)b=mid;else a=mid+1,x|=1;
}
Ins(root[x],d,e);
}
void ask(int x,int a,int b,int c,int d,int e){
if(c<=a&&b<=d){ans+=Ask(root[x],e);return;}
int mid=(a+b)>>1;
if(c<=mid)ask(x<<1,a,mid,c,d,e);
if(d>mid)ask(x<<1|1,mid+1,b,c,d,e);
}
int main(){
blank->l=blank->r=blank;
for(read(n),i=1;i<=n;i++)read(a[i]),b[++C]=a[i];
for(read(m),i=1;i<=m;i++){
while(!(((ch=getchar())=='U')||(ch=='Q')));
read(que[i][1]),read(que[i][2]);
if(ch=='U')b[++C]=que[i][2];else que[i][0]=1;
}
for(sort(b+1,b+C+1),i=1;i<=C;i++)T[i].insert(0),T[i].insert(n+1);
for(build(i=1,1,n);i<=n;i++){
T[a[i]=lower(a[i])].insert(i);
add(i,loc[a[i]],b[a[i]]);
loc[a[i]]=i;
}
for(i=1;i<=m;i++)if(que[i][0])ans=0,ask(1,1,n,que[i][1],que[i][2],que[i][1]),printf("%lld\n",ans);
else{
x=que[i][1],y=lower(que[i][2]);
k=j=T[a[x]].find(x);k--;pre=*k;k=j;k++;nxt=*k;
add(x,pre,-b[a[x]]);
if(nxt<=n)add(nxt,x,-b[a[x]]),add(nxt,pre,b[a[x]]);
T[a[x]].erase(x);
T[a[x]=y].insert(x);
k=j=T[y].find(x);k--;pre=*k;k=j;k++;nxt=*k;
add(x,pre,b[y]);
if(nxt<=n)add(nxt,pre,-b[y]),add(nxt,x,b[y]);
}
return 0;
}

  

BZOJ2883 : gss2加强版的更多相关文章

  1. bzoj AC倒序

    Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...

  2. 【BZOJ-3514】Codechef MARCH14 GERALD07加强版 LinkCutTree + 主席树

    3514: Codechef MARCH14 GERALD07加强版 Time Limit: 60 Sec  Memory Limit: 256 MBSubmit: 1288  Solved: 490 ...

  3. 斐波拉契数列加强版——时间复杂度O(1),空间复杂度O(1)

    对于斐波拉契经典问题,我们都非常熟悉,通过递推公式F(n) = F(n - ) + F(n - ),我们可以在线性时间内求出第n项F(n),现在考虑斐波拉契的加强版,我们要求的项数n的范围为int范围 ...

  4. vijos P1915 解方程 加强版

    背景 B酱为NOIP 2014出了一道有趣的题目, 可是在NOIP现场, B酱发现数据规模给错了, 他很伤心, 哭得很可怜..... 为了安慰可怜的B酱, vijos刻意挂出来了真实的题目! 描述 已 ...

  5. [NOIP2014] 解方程&加强版 (bzoj3751 & vijos1915)

    大概有$O(m)$,$O(n\sqrt{nm})$,$O(n\sqrt{m})$的3个算法,其中后2个可以过加强版.代码是算法3,注意bzoj的数据卡掉了小于20000的质数. #include< ...

  6. HTML5音乐播放器(最新升级改造加强版)

    最近么,单位里面么老不顺心的,公司一直催要程序员要PHP,然后本宅好不容易推荐了一个,我日嘞,最后待遇变成1.3,吾师最后也同意1.3W,然后还说要考虑... 尼玛,4年多5年不到一点的工作经验,前端 ...

  7. 【转】PowerShell入门(二):PowerShell是Cmd命令行的加强版吗?

    转至:http://www.cnblogs.com/ceachy/archive/2013/01/31/PowerShell_vs_Cmd.html PowerShell是命令行的加强版吗?Power ...

  8. 【BZOJ-4407】于神之怒加强版 莫比乌斯反演 + 线性筛

    4407: 于神之怒加强版 Time Limit: 80 Sec  Memory Limit: 512 MBSubmit: 241  Solved: 119[Submit][Status][Discu ...

  9. 【BZOJ-3545&3551】Peaks&加强版 Kruskal重构树 + 主席树 + DFS序 + 倍增

    3545: [ONTAK2010]Peaks Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 1202  Solved: 321[Submit][Sta ...

随机推荐

  1. unity3d中资源文件从MAX或者MAYA中导出的注意事项

    原地址:http://blog.sina.com.cn/s/blog_6ad33d3501011ekx.html 之前在项目中,没有怎么接触过美术的软件(之前的美术团队很犀利,被他们宠坏了).在自己公 ...

  2. JQuery入门

    JQuery入门 1 jQuery的概述 1.1 jQuery简介 jQuery是一个 JavaScript函数库,它是一个“写的更少,但做的更多”的轻量级 JavaScript 库.jQuery 极 ...

  3. 安装ruby

    (这些文章都是从我的个人主页上粘贴过来的,大家也可以访问我的主页 www.iwangzheng.com) 过程中有点小曲折,我们leader是技术大牛,现在我生命中多了个超高智商处女座man了,还有一 ...

  4. c++中的srand()和rand() 转载 自:http://blog.sina.com.cn/s/blog_624c2c4001012f67.html

    今天看了同事写的小程序,发现了其中出现了srand()和rand()这两个我以前没有用过的函数,当然从名字可以看出肯定能随机数有关,于是网查资料知这两个函数配合一起使用来产生随机数的,哈哈,又长知识了 ...

  5. WPF 路由事件总结

    1.什么是路由事件 已下为MSDN中的定义 功能定义:路由事件是一种可以针对元素树中的多个侦听器(而不是仅针对引发该事件的对象)调用处理程序的事件. 实现定义:路由事件是一个 CLR 事件,可以由 R ...

  6. Nmap备忘单:从探索到漏洞利用(Part3)

    众所周知NMAP是经常用来进行端口发现.端口识别.除此之外我们还可以通过NMAP的NSE脚本做很多事情,比如邮件指纹识别,检索WHOIS记录,使用UDP服务等. 发现地理位置 Gorjan Petro ...

  7. Floyd算法 及其运用

    #include<stdio.h> ][]; ][]; void floyd(int n) { ;k<=n;k++) { ;i<=n;i++) { ;j<=n;j++) ...

  8. 将XML解析成DOM文档

    在支持html5的浏览其中,可以使用标准解析器DOMParser对象进行解析html或者xml等字符串 var data = '<div></div>'; var tmp = ...

  9. django-cms 代码研究(二)bugs?

    djangocms集成到现有项目中后,发现了几个问题: 1. 现有项目的url匹配失效,下面requests请求被交给djangocms处理了 url(r'^admin/', include(admi ...

  10. [flag飞起]

    重度Flag: Rush 未来程序・改 (note: 由于代码存放地址原因(物理)无法本周完成) 轻度Flag: 未来程序・改写完后刷QTREE与主席树 严重Flag: 计算几何...