题目描述

一个可重复数字集合S的神秘数定义为最小的不能被S的子集的和表示的正整数。例如S={1,1,1,4,13},
1 = 1
2 = 1+1
3 = 1+1+1
4 = 4
5 = 4+1
6 = 4+1+1
7 = 4+1+1+1
8无法表示为集合S的子集的和,故集合S的神秘数为8。
现给定n个正整数a[1]..a[n],m个询问,每次询问给定一个区间[l,r](l<=r),求由a[l],a[l+1],…,a[r]所构成的可重复数字集合的神秘数。

输入

第一行一个整数n,表示数字个数。
第二行n个整数,从1编号。
第三行一个整数m,表示询问个数。
以下m行,每行一对整数l,r,表示一个询问。

输出

对于每个询问,输出一行对应的答案。

样例输入

5
1 2 4 9 10
5
1 1
1 2
1 3
1 4
1 5

样例输出

2
4
8
8
8


题解

主席树的一道神题

我们先想暴力怎么做:把一段区间的数取出来,排个序,从小到大选择。如果$a1$~$a_{i-1}$能够表示$1~x$,此时加入$a_i$,如果$a_i\le x+1$,那么就可以表示$x+a_i$,否则x就是答案。

试着优化一下这个过程:设$a_{i-1}=k$,$a_i=y$,1~i-1的神秘数为ans=x+1,那么显然$ans=\sum\limits_{t=1}^{i-1}a_t$。此时如果存在k+1~ans的数就可以更新ans。更具体地,如果k+1~ans内的数的和为s,那么ans+=s;而ans为1~k的数的和+1,故ans的新值应该赋为1~ans的数的和。

说了这么多废话有什么用?我们可以发现每次ans的增量都大于等于前一次的ans,所以这个过程的时间复杂度应该为$O(\log a)$。

而事实上我们并不能把区间拿出来排序,所以需要使用数据结构,上一个主席树就好了。

时间复杂度为$O(n\log^2n)$

#include <cstdio>
#include <algorithm>
#define N 100010
using namespace std;
int v[N] , a[N] , root[N] , ls[N << 5] , rs[N << 5] , sum[N << 5] , tot;
void insert(int p , int l , int r , int x , int &y)
{y = ++tot , sum[y] = sum[x] + a[p];if(l == r) return;int mid = (l + r) >> 1;if(p <= mid) rs[y] = rs[x] , insert(p , l , mid , ls[x] , ls[y]);else ls[y] = ls[x] , insert(p , mid + 1 , r , rs[x] , rs[y]);
}
int query(int p , int l , int r , int x , int y)
{if(r <= p) return sum[y] - sum[x];int mid = (l + r) >> 1;if(p <= mid) return query(p , l , mid , ls[x] , ls[y]);else return query(p , mid + 1 , r , rs[x] , rs[y]) + sum[ls[y]] - sum[ls[x]];
}
int main()
{int n , m , i , x , y , ans , tmp;scanf("%d" , &n);for(i = 1 ; i <= n ; i ++ ) scanf("%d" , &a[i]) , v[i] = a[i];sort(a + 1 , a + n + 1);for(i = 1 ; i <= n ; i ++ ) v[i] = lower_bound(a + 1 , a + n + 1 , v[i]) - a;for(i = 1 ; i <= n ; i ++ ) insert(v[i] , 1 , n , root[i - 1] , root[i]);a[n + 1] = 1 << 30;scanf("%d" , &m);while(m -- ){scanf("%d%d" , &x , &y) , ans = 1;while((tmp = query(upper_bound(a + 1 , a + n + 2 , ans) - a - 1 , 1 , n , root[x - 1] , root[y])) >= ans)ans = tmp + 1;printf("%d\n" , ans);}return 0;
}

转载于:https://www.cnblogs.com/GXZlegend/p/7115254.html

【bzoj4408】[Fjoi 2016]神秘数 主席树相关推荐

  1. BZOJ44084299[Fjoi 2016]神秘数——主席树

    题目描述 一个可重复数字集合S的神秘数定义为最小的不能被S的子集的和表示的正整数.例如S={1,1,1,4,13}, 1 = 1 2 = 1+1 3 = 1+1+1 4 = 4 5 = 4+1 6 = ...

  2. BZOJ4408: [Fjoi 2016]神秘数

    考虑处理单组询问 一开始ans=1,每次统计<=ans的数的和sum,若sum>=ans,就用sum+1更新ans 最坏情况下,ans的更新是1,2,3,5,8-.是斐波拉契数列 所以an ...

  3. Bzoj 4408: [Fjoi 2016]神秘数 可持久化线段树,神题

    4408: [Fjoi 2016]神秘数 Time Limit: 10 Sec  Memory Limit: 128 MB Submit: 177  Solved: 128 [Submit][Stat ...

  4. 4408: [Fjoi 2016]神秘数

    4408: [Fjoi 2016]神秘数 Time Limit: 10 Sec  Memory Limit: 128 MB Submit: 452  Solved: 273 [Submit][Stat ...

  5. BZOJ 4408: [Fjoi 2016]神秘数(可持久化线段树)

    题目描述 http://www.lydsy.com/JudgeOnline/problem.php?id=4408 题目大意:求最小不能被一段区间中某些数的和表示的数.(还是看题面吧) 思路 可持久化 ...

  6. ●BZOJ 4408 [Fjoi 2016]神秘数

    题链: http://www.lydsy.com/JudgeOnline/problem.php?id=4408 题解: 主席树 首先,对于一些数来说, 如果可以我们可以使得其中的某些数能够拼出 1- ...

  7. bzoj 4408: [FJOI2016]神秘数 主席树

    先看一下对于给定的一段如何暴力求最大值,首先将这一段排序,然后如果存在某一个数,这个数比它前面的数的前缀和至少大2,那么答案就是它前面那个数的前缀和+1. 那么假设现在处理了前面较小的一些数之后的答案 ...

  8. LUOGU P4587 [FJOI2016]神秘数(主席树)

    传送门 解题思路 如果区间内没有\(1\),那么答案就为\(1\),从这一点继续归纳.如果区间内有\(x\)个\(1\),设区间内\([2,x+1]\)的和为\(sum\),如果\(sum=0\),那 ...

  9. [FJOI 2016]bzoj 4408 神秘数 - 线段树

    题目大意:给你一列数,多次询问用一个区间的数字形成一个可重集合,最小的不能被表示为其一个子集的数字是多少. 题解:考虑给你一个可重集合你怎么算:从小到大排序,假设用前x个数字不能表示的最小都数字是an ...

最新文章

  1. 直接送显示器!屏幕大了,打起代码才带劲!
  2. springCloud之配置中心学习
  3. 屏蔽鼠标右键、Ctrl+N、Shift+F10、F11、F5刷新、退格键
  4. matlab rbf函数_基于径向基函数(RBF)的无网格伪谱法与程序实现(2)——微分矩阵...
  5. Python 析构方法__del__
  6. Razor.js,基于JavaScript的Razor实现
  7. 关于8位AD_DA转换芯片的采样率问题
  8. 如何用jar命令对java工程进行打包
  9. python拉格朗日插值法_【统计学】拉格朗日插值法的一种python实现方式
  10. 记录CodeForces第一次比赛经历
  11. LRU算法的实现(STL+模拟)
  12. Python3 - 深入学习Docker容器管理
  13. Protel 99 SE 应用技术问答
  14. 高能!一大波奇葩挖矿方式来袭~
  15. html5 自动设置全屏,HTML5 body如何设置自适应全屏 HTML5 body设置自适应全屏代码
  16. STM32低功耗模式
  17. 7-45 连年 7-69 求婚
  18. Swiper 轮播插件
  19. Allegro自动备份PCB设计文件的方法
  20. 用英雄联盟的方式讲解JavaScript设计模式!

热门文章

  1. Python 知识点笔记一常用数据类型
  2. CTFshow 反序列化 web257
  3. OpenCV中图像窗口的鼠标事件
  4. [YTU]_2442( C++习题 矩阵求和--重载运算符)
  5. 3.3 matlab用switch语句实现选择结构
  6. matlab程序和程序文件
  7. 程序员编程艺术第一章(第二节)
  8. php数组为什么其他语言,PHP语言特性和各版本的差异
  9. Dask.distributed创建client时出现报错
  10. 运行SSIS包的几种方式