博客
关于我
hdu4747 Mex(线段树)
阅读量:230 次
发布时间:2019-03-01

本文共 3306 字,大约阅读时间需要 11 分钟。

题意:

给定长度为n的序列a,

求所有子区间的mex的和.

数据范围:n<=2e5,0<=a(i)<=1e9

解法:

考虑枚举左端点,计算所有右端点的mex,那么对于每个左端点的所有右端点的mex,全部加起来就是答案.考虑递推,令mex[i]表示区间[1,i]的mex,左端点为1的mex[]可以O(1)计算出来.接下来考虑如何用[1,i]的mex[]计算出[2,i]的mex[].[1,i]->[2,i]时,要删掉a[1],考虑删掉a[1]会对哪些mex[i]造成影响,设val=a[1],下一个val出现的位置是pos,即a[pos]=val,因为pos以及pos后面的所有位置,都有val了,删掉无影响,那么[1,pos-1]的mex[]中,大于val的都要变成val,由于mex[]一定是非递减的,那么找到[1,pos-1]中第一个大于val的位置t,那么操作变为[t,pos-1]的区间覆盖,覆盖为val.用线段树做这道题,维护区间max和区间sum,其中区间max用于树上二分找第一个大于val的位置t.因为还需要知道下一个val出现的位置,所以还需要预处理一下nt[].

code:

#include 
using namespace std;#define int long longconst int maxm=2e5+5;int nt[maxm];int a[maxm];int b[maxm];int n;struct Tree{ int ma[maxm<<2]; int sum[maxm<<2]; int laz[maxm<<2]; inline void pp(int node){ sum[node]=sum[node*2]+sum[node*2+1]; ma[node]=max(ma[node*2],ma[node*2+1]); } inline void pd(int node,int l,int r){ if(laz[node]!=-1){ int mid=(l+r)/2; sum[node*2]=(mid-l+1)*laz[node]; sum[node*2+1]=(r-mid)*laz[node]; ma[node*2]=laz[node]; ma[node*2+1]=laz[node]; laz[node*2]=laz[node]; laz[node*2+1]=laz[node]; laz[node]=-1; } } void update(int st,int ed,int val,int l,int r,int node){ //区间覆盖 if(st<=l&&ed>=r){ sum[node]=(r-l+1)*val; ma[node]=val; laz[node]=val; return ; } pd(node,l,r); int mid=(l+r)/2; if(st<=mid)update(st,ed,val,l,mid,node*2); if(ed>mid)update(st,ed,val,mid+1,r,node*2+1); pp(node); } int ask(int st,int ed,int l,int r,int node){ //区间求和 if(st<=l&&ed>=r)return sum[node]; pd(node,l,r); int mid=(l+r)/2; int ans=0; if(st<=mid)ans+=ask(st,ed,l,mid,node*2); if(ed>mid)ans+=ask(st,ed,mid+1,r,node*2+1); return ans; } int ffind(int val,int l,int r,int node){ //树上二分找大于val的第一个位置 if(ma[node]<=val)return n+1; if(l==r)return l; pd(node,l,r); int mid=(l+r)/2; int ans=n+1; ans=ffind(val,l,mid,node*2); if(ans==n+1)ans=ffind(val,mid+1,r,node*2+1); return ans; } void build(int l,int r,int node){ //建树 laz[node]=-1; if(l==r){ sum[node]=b[l]; ma[node]=b[l]; return ; } int mid=(l+r)/2; build(l,mid,node*2); build(mid+1,r,node*2+1); pp(node); }}T;signed main(){ ios::sync_with_stdio(0);cin.tie(0); while(cin>>n&&n){ for(int i=1;i<=n;i++){ cin>>a[i]; } int now=0; map
mark; for(int i=1;i<=n;i++){ mark[a[i]]++; while(mark[now])now++; b[i]=now; } mark.clear(); for(int i=n;i>=1;i--){ if(!mark[a[i]]){ nt[i]=n+1; }else{ nt[i]=mark[a[i]]; } mark[a[i]]=i; } T.build(1,n,1); int ans=0; for(int i=1;i<=n;i++){ ans+=T.ask(i,n,1,n,1); if(i!=n){ int t=T.ffind(a[i],1,n,1); int pos=nt[i]; if(t<=pos-1){ T.update(t,pos-1,a[i],1,n,1); } } } cout<
<

转载地址:http://ptkv.baihongyu.com/

你可能感兴趣的文章
Mysql学习总结(30)——MySQL 索引详解大全
查看>>
Mysql学习总结(31)——MySql使用建议,尽量避免这些问题
查看>>
Mysql学习总结(32)——MySQL分页技术详解
查看>>
Mysql学习总结(33)——阿里云centos配置MySQL主从复制
查看>>
Mysql学习总结(35)——Mysql两千万数据优化及迁移
查看>>
Mysql学习总结(36)——Mysql查询优化
查看>>
Mysql学习总结(37)——Mysql Limit 分页查询优化
查看>>
Mysql学习总结(38)——21条MySql性能优化经验
查看>>
Mysql学习总结(39)——49条MySql语句优化技巧
查看>>
Mysql学习总结(3)——MySql语句大全:创建、授权、查询、修改等
查看>>
Mysql学习总结(40)——MySql之Select用法汇总
查看>>
Mysql学习总结(41)——MySql数据库基本语句再体会
查看>>
Mysql学习总结(42)——MySql常用脚本大全
查看>>
Mysql学习总结(43)——MySQL主从复制详细配置
查看>>
Mysql学习总结(44)——Linux下如何实现mysql数据库每天自动备份定时备份
查看>>
Mysql学习总结(45)——Mysql视图和事务
查看>>
Mysql学习总结(46)——8种常被忽视的SQL错误用法
查看>>
Mysql学习总结(48)——MySql的日志与备份还原
查看>>
Mysql学习总结(49)——从开发规范、选型、拆分到减压
查看>>
Mysql学习总结(4)——MySql基础知识、存储引擎与常用数据类型
查看>>