博客
关于我
hdu4747 Mex(线段树)
阅读量:229 次
发布时间: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学习总结(68)——MYSQL统计每天、每周、每月、每年数据 SQL 总结
查看>>
Mysql学习总结(69)——Mysql EXPLAIN 命令使用总结
查看>>
Mysql学习总结(6)——MySql之ALTER命令用法详细解读
查看>>
Mysql学习总结(70)——MySQL 优化实施方案
查看>>
Mysql学习总结(71)——MySQL 重复记录查询与删除总结
查看>>
Mysql学习总结(71)——数据库介绍(MySQL安装 体系结构、基本管理)再回顾
查看>>
Mysql学习总结(72)——MySQL 开发者开发,设计规范再总结
查看>>
Mysql学习总结(73)——MySQL 查询A表存在B表不存在的数据SQL总结
查看>>
Mysql学习总结(74)——慢SQL!压垮团队的最后一根稻草!
查看>>
Mysql学习总结(75)——并发量大、数据量大的互联网业务数据库设计军规
查看>>
Mysql学习总结(76)——MySQL执行计划(explain)结果含义总结
查看>>
Mysql学习总结(77)——温故Mysql数据库开发核心原则与规范
查看>>
Mysql学习总结(78)——MySQL各版本差异整理
查看>>
Mysql学习总结(79)——MySQL常用函数总结
查看>>
Mysql学习总结(7)——MySql索引原理与使用大全
查看>>
Mysql学习总结(80)——统计数据库的总记录数和库中各个表的数据量
查看>>
Mysql学习总结(81)——为什么MySQL不推荐使用uuid或者雪花id作为主键?
查看>>
Mysql学习总结(82)——MySQL逻辑删除与数据库唯一性约束如何解决?
查看>>
Mysql学习总结(83)——常用的几种分布式锁:ZK分布式锁、Redis分布式锁、数据库分布式锁、基于JDK的分布式锁方案对比总结
查看>>
Mysql学习总结(84)—— Mysql的主从复制延迟问题总结
查看>>