hdu 4747 mex 线段树+思维
http://acm.hdu.edu.cn/showproblem.php?pid=4747
题意:
我们定义mex(l,r)表示一个序列a[l]....a[r]中没有出现过得最小的非负整数, 然后我们给出一个长度为n的序列,求他所有的连续的子序列的mex(l,r)的和。
思路:
首先因为n的最大值就是2*10^5 所有我们字需要考虑200000之内的数就好了,然后O(2*n)可以求出(1,1),(1,2), (1,3),(1,4) ... (1,n)来 mex是不减的。
然后我们考虑将第一个数拿走我们就能够得到(2,2),(2,3) ......(2,n) , 如何求他们?下边给出图解:
下边是粘贴别人的,感觉有个例子很好理解。
例: 1, 6,0,2,3,1,4,3
初始mex 0, 0,2,3,4,4,5,5 mex[1,r]
删除1后 0, 0,1,1,1,4,5,5 mex[2,r]
……
当删除第一个1后,红色的mex不变!,紫色的mex值变为1,橙色的mex值不变,删除点的mex置0
因此,用线段树维护一个单调不递增队列,每次求和。查找位置时用二分。线段树延时标记即可。
ps:我这里二笔了一下,题意一下大家,lazy一定要出事化为-1,因为这里面有置0的操作。我因此wa好多次。。
#include <iostream>
#include <cstdio>
#include <cmath>
#include <vector>
#include <cstring>
#include <algorithm>
#include <string>
#include <set>
#include <functional>
#include <numeric>
#include <sstream>
#include <stack>
#include <map>
#include <queue> #define CL(arr, val) memset(arr, val, sizeof(arr)) #define lc l,m,rt<<1
#define rc m + 1,r,rt<<1|1
#define pi acos(-1.0)
#define ll __int64
#define L(x) (x) << 1
#define R(x) (x) << 1 | 1
#define MID(l, r) (l + r) >> 1
#define Min(x, y) (x) < (y) ? (x) : (y)
#define Max(x, y) (x) < (y) ? (y) : (x)
#define E(x) (1 << (x))
#define iabs(x) (x) < 0 ? -(x) : (x)
#define OUT(x) printf("%I64d\n", x)
#define keyTree (chd[chd[root][1]][0])
#define Read() freopen("din.txt", "r", stdin)
#define Write() freopen("dout.txt", "w", stdout); #define M 100007
#define N 200017 using namespace std; int dx[4]={-1,1,0,0};
int dy[4]={0,0,-1,1}; const int inf = 0x7f7f7f7f;
const int mod = 1000000007; const double eps = 1e-8; int mex[N],a[N],next[N];
int vis[N];
ll val[4*N], lz[4*N];
int n; void pushup(int rt)
{
val[rt] = val[rt<<1] + val[rt<<1|1];
}
void pushdown(int rt,int m)
{
if (lz[rt] != -1)
{
lz[rt<<1] = lz[rt<<1|1] = lz[rt];
val[rt<<1] = lz[rt] * (m - (m>>1));
val[rt<<1|1] = lz[rt] * (m>>1);
lz[rt] = -1;
}
}
void build(int l, int r, int rt)
{
lz[rt] = -1;
val[rt] = 0;
if (l == r)
{
val[rt] = mex[l];
return ;
}
int m = (l + r) >> 1;
build(lc); build(rc);
pushup(rt);
}
void update(int L, int R, ll sc, int l, int r, int rt)
{
if (l >= L && r <= R)
{
lz[rt] = sc;
val[rt] = sc*(r - l + 1);
return ;
}
pushdown(rt,r - l + 1);
int m = (l + r) >> 1;
if (L <= m) update(L,R,sc,lc);
if (R > m) update(L,R,sc,rc);
pushup(rt);
}
ll query(int pos, int l,int r, int rt)
{
if (l == r) return val[rt];
int m = (l + r) >> 1;
pushdown(rt, r - l + 1);
if (pos <= m) return query(pos,lc);
else return query(pos,rc);
}
int BSR(int l, int r, int v)
{
int ans = -1;
while (l <= r)
{
int mid = (l + r) >> 1;
if (query(mid,1,n,1) > v)
{
ans = mid;
r = mid - 1;
} else l = mid + 1;
}
return ans;
}
int main()
{
while (~scanf("%d",&n))
{
if (!n) break;
CL(vis,0); CL(next,0);
for (int i = 1; i <= n; ++i)
{
scanf("%d",&a[i]);
a[i] = min(a[i],200001);
if (vis[a[i]]) next[vis[a[i]]] = i;
vis[a[i]] = i;
}
for (int i = 1; i <= n; ++i) if (!next[i]) next[i] = n + 1;
CL(vis,0); int last = 0;
for (int i = 1; i <= n; ++i)
{
vis[a[i]] = 1;
while (true)
{
if (!vis[last])
{
mex[i] = last;
break;
}
last++;
}
}
build(1,n,1); ll ans = 0;
for (int i = 1; i <= n; ++i)
{
ans += val[1];
if (i == n) continue;
int l = i + 1;
int r = next[i] - 1;
int pos = BSR(l,r,a[i]);
if (pos != -1) update(pos, r, a[i], 1, n, 1);
update(i, i, 0, 1, n, 1);
}
printf("%I64d\n",ans);
}
return 0;
}
hdu 4747 mex 线段树+思维的更多相关文章
- hdu 4747 Mex( 线段树? 不,区间处理就行(dp?))
Mex Time Limit: 15000/5000 MS (Java/Others) Memory Limit: 65535/65535 K (Java/Others)Total Submis ...
- HDU 4747 Mex ( 线段树好题 + 思路 )
参考:http://www.cnblogs.com/oyking/p/3323306.html 相当不错的思路,膜拜之~ 个人理解改日补充. #include <cstdio> #incl ...
- hdu 4747【线段树-成段更新】.cpp
题意: 给出一个有n个数的数列,并定义mex(l, r)表示数列中第l个元素到第r个元素中第一个没有出现的最小非负整数. 求出这个数列中所有mex的值. 思路: 可以看出对于一个数列,mex(r, r ...
- hdu 3333 离线线段树 + 思维/树状数组 /在线主席树
#include<iostream> #include<cstdio> #include<string> #include<cmath> #includ ...
- hdu 4031 attack 线段树区间更新
Attack Time Limit: 5000/3000 MS (Java/Others) Memory Limit: 65768/65768 K (Java/Others)Total Subm ...
- hdu 4288 离线线段树+间隔求和
Coder Time Limit: 20000/10000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total Su ...
- hdu 3016 dp+线段树
Man Down Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total S ...
- HDU 4747 Mex 递推/线段树
题目链接: acm.hdu.edu.cn/showproblem.php?pid=4747 Mex Time Limit: 15000/5000 MS (Java/Others)Memory Limi ...
- HDU 4747 Mex (2013杭州网络赛1010题,线段树)
Mex Time Limit: 15000/5000 MS (Java/Others) Memory Limit: 65535/65535 K (Java/Others)Total Submis ...
随机推荐
- netty之LengthFieldBasedFrameDecoder解码器
官方api:http://netty.io/4.1/api/io/netty/handler/codec/LengthFieldBasedFrameDecoder.html package com.e ...
- mongodb拆库分表脚本
脚本功能: 1. 将指定的报告文件按照指定的字段.切库切表策略切分 2. 将切分后的文件并发导入到对应的Mongodb中 3. 生成日志文件和done标识文件 使用手册: -h 打印帮助信息,并 ...
- CH5E09 能量相连【区间DP】
5E09 能量项链 0x5E「动态规划」练习 描述 在Mars星球上,每个Mars人都随身佩带着一串能量项链.在项链上有N颗能量珠.能量珠是一颗有头标记与尾标记的珠子,这些标记对应着某个正整数.并且, ...
- webView返回不刷新
[[NSUserDefaults standardUserDefaults] setInteger:2 forKey: @"WebKitCacheModelPreferenceKey&quo ...
- EOS Dapp开发(1)-基于Docker的开发环境搭建
随着EOS主网的上线,相信基于EOS的Dapp开发会越来越多,查阅了很多资料相关的开发资料都不是很多,只能自己摸索,按照网上仅有的几篇教程,先git clonehttps://github.com/E ...
- mysql 数据操作 多表查询 子查询 带EXISTS关键字的子查询
带EXISTS关键字的子查询 EXISTS关字键字表示存在. EXISTS 判断某个sql语句的有没有查到结果 有就返回真 true 否则返回假 False 如果条件成立 返回另外一条sql语句的返 ...
- Linux资源使用配置文件 /etc/security/limits.conf
Linux资源使用配置文件 /etc/security/limits.conf http://www.linuxidc.com/Linux/2012-05/59489.htm Linux就这个范儿P5 ...
- 使用Atom预览markdown
1.打开任意.md文件(markdown源文件)2.windows : ctrl + shift + pmac : command + shift + p这条命令跟Sublime Text是一样的,打 ...
- EasyUI Tooltip 提示框
通过 $.fn.tooltip.defaults 重写默认的 defaults. 当用户移动鼠标指针在某个元素上时,出现提示信息窗口用来显示额外信息.提示内容可以包含任何来自页面的或者通过 ajax ...
- Spark 参数配置的几种方法
1.Spark 属性Spark应用程序的运行是通过外部参数来控制的,参数的设置正确与否,好与坏会直接影响应用程序的性能,也就影响我们整个集群的性能.参数控制有以下方式:(1)直接设置在SparkCon ...