Description

一个可重复数字集合$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\;\leq\;r)$,求由$a[l],a[l+1],…,a[r]$所构成的可重复数字集合的神秘数。

Input

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

Output

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

Sample Input

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

Sample Output

2
4
8
8
8

HINT

$n,m\;\leq\;10^5,\sum\;a[i]\;\leq\;10^9$

Solution

当$x$为答案时,满足条件$\sum_{a[i]\;\leq\;x}a[i]<x$

枚举$x$,用主席树求和,判断是否合法,若不合法,$x=\sum_{a[i]\;\leq\;x}a[i]$

#include<cmath>
#include<ctime>
#include<queue>
#include<stack>
#include<cstdio>
#include<vector>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define N 100005
#define M 2000005
using namespace std;
typedef long long ll;
struct linetree{
int l,r,s,lc,rc;ll t;
}lt[M];
int a[N],p[N],rt[N],n,m,l,r,siz,cnt,tmp;
inline void bld(int u,int l,int r,int k){
lt[u].l=l;lt[u].r=r;
if(k>=p[l]&&k<=p[r]){
lt[u].s=1;lt[u].t=(ll)(k);
}
if(lt[u].l<lt[u].r){
int mid=lt[u].l+lt[u].r>>1;
lt[u].lc=++cnt;bld(cnt,l,mid,k);
lt[u].rc=++cnt;bld(cnt,mid+1,r,k);
}
}
inline void build(int u,int lst,int l,int r,int k){
lt[u].l=l;lt[u].r=r;
lt[u].s=lt[lst].s+1;lt[u].t=lt[lst].t+(ll)(k);
if(lt[u].l<lt[u].r){
int mid=lt[u].l+lt[u].r>>1;
if(k<=p[mid]){
lt[u].lc=++cnt;
build(cnt,lt[lst].lc,l,mid,k);
if(lt[lst].rc) lt[u].rc=lt[lst].rc;
else{
lt[u].rc=++cnt;
build(cnt,lt[lst].rc,mid+1,r,k);
}
}
else{
lt[u].rc=++cnt;
build(cnt,lt[lst].rc,mid+1,r,k);
if(lt[lst].lc) lt[u].lc=lt[lst].lc;
else{
lt[u].lc=++cnt;
build(cnt,lt[lst].lc,l,mid,k);
}
}
}
}
inline ll ask(int u,int x){
if(!x||!u) return 0;
if(p[lt[u].r]<=x) return lt[u].t;
int mid=lt[u].l+lt[u].r>>1;
if(x<=p[mid]) return ask(lt[u].lc,x);
return lt[lt[u].lc].t+ask(lt[u].rc,x);
}
inline int chk(int u,int x){
if(!x||!u) return 0;
if(p[lt[u].r]<=x) return lt[u].s;
int mid=lt[u].l+lt[u].r>>1;
if(x<=p[mid]) return chk(lt[u].lc,x);
return lt[lt[u].lc].s+chk(lt[u].rc,x);
}
inline void Aireen(){
scanf("%d",&n);
for(int i=1;i<=n;++i){
scanf("%d",&a[i]);p[i]=a[i];
}
sort(p+1,p+1+n);
siz=1;
for(int i=2;i<=n;++i)
if(p[i]!=p[i-1]) p[++siz]=p[i];
rt[1]=++cnt;bld(cnt,1,siz,a[1]);
for(int i=2;i<=n;++i){
rt[i]=++cnt;
build(cnt,rt[i-1],1,siz,a[i]);
}
scanf("%d",&m);
while(m--){
scanf("%d%d",&l,&r);
for(ll i=0,s;;i=s+1ll){
s=ask(rt[r],i)-ask(rt[l-1],i);
if(s>=ask(rt[r],p[siz])-ask(rt[l-1],p[siz])){
printf("%lld\n",s+1ll);break;
}
tmp=upper_bound(p+1,p+1+siz,(int)(s+1ll))-p-1;
if(ask(rt[r],p[tmp])-ask(rt[l-1],p[tmp])<s+1ll){
printf("%lld\n",s+1ll); break;
}
}
}
}
int main(){
freopen("mystic.in","r",stdin);
freopen("mystic.out","w",stdout);
Aireen();
fclose(stdin);
fclose(stdout);
return 0;
}

[bzoj4408][Fjoi2016]神秘数的更多相关文章

  1. (bzoj4408)[FJOI2016]神秘数(可持久化线段树)

    (bzoj4408)[FJOI2016]神秘数(可持久化线段树) bzoj luogu 对于一个区间的数,排序之后从左到右每一个数扫 如果扫到某个数a时已经证明了前面的数能表示[1,x],那么分情况: ...

  2. 题解【bzoj4587 & bzoj4408 [FJOI2016]神秘数】

    Description \(n\) 个数的序列,每次询问一个区间,求最小的一个数使得不能用这个区间中的数之和表示. \(n \leq 10^5, \sum a_i \leq 10^9\) 这两个题一个 ...

  3. 【BZOJ4408】[FJOI2016]神秘数(主席树)

    [BZOJ4408][FJOI2016]神秘数(主席树) 题面 BZOJ 洛谷 题解 考虑只有一次询问. 我们把所有数排个序,假设当前可以表示出的最大数是\(x\). 起始\(x=0\). 依次考虑接 ...

  4. 【LG4587】[FJOI2016]神秘数

    [LG4587][FJOI2016]神秘数 题面 洛谷 题解 首先我们想一想暴力怎么做 对于一段区间\([l,r]\) 我们先将它之间的数升序排序 从左往右扫, 设当前我们可以表示出的数为\([1,x ...

  5. BZOJ4299 & CC FRBSUM:ForbiddenSum & BZOJ4408 & 洛谷4587 & LOJ2174:[FJOI2016]神秘数——题解

    https://www.lydsy.com/JudgeOnline/problem.php?id=4299 https://www.lydsy.com/JudgeOnline/problem.php? ...

  6. [FJOI2016]神秘数(脑洞+可持久化)

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

  7. Luogu P4587 [FJOI2016]神秘数

    一道好冷门的好题啊,算是对于一个小结论和数据结构的一点考验吧 首先看完题目我们发现要从这个神秘数的性质入手,我们观察or手玩可得: 如果有\(x\)个\(1\),那么\([1,x]\)都是可以表示出来 ...

  8. 【BZOJ-4408】神秘数 可持久化线段树

    4408: [Fjoi 2016]神秘数 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 475  Solved: 287[Submit][Status ...

  9. BZOJ 4408 FJOI2016 神秘数 可持久化线段树

    Description 一个可重复数字集合S的神秘数定义为最小的不能被S的子集的和表示的正整数.例如S={1,1,1,4,13},1 = 12 = 1+13 = 1+1+14 = 45 = 4+16 ...

随机推荐

  1. tomcat安装配置.md

    tomcat 安装 安装jdk jdk我使用的是oracle的,如果下载请在oracle的官网上下载.或者你也可以使用openjdk,官网在http://openjdk.java.net/. # ta ...

  2. k-d tree 学习笔记

    以下是一些奇怪的链接有兴趣的可以看看: https://blog.sengxian.com/algorithms/k-dimensional-tree http://zgjkt.blog.uoj.ac ...

  3. dubbo2.5.3 与spring 3.1.x 冲突

    在集成了dubbo2.5.3 的项目中初始化出现 MalformedParameterizedTypeException 检查发现这是因为dubbo2.5.3依赖的springframeworks是2 ...

  4. usb驱动开发12之设备生命线

    函数usb_control_msg完成一些初始化后调用了usb_internal_control_msg之后就free urb.剩下的活,全部留给usb_internal_control_msg去做了 ...

  5. AR 不同 继承映射的问题总结

    在使用AR(Nhibernate) 做ORM时,使用类的继承体系时,它有不同的映射方式,解决的问题不同,带来的问题差异也很大. 1.所有数据 存储在一张表,不同的类使用 DiscriminatorCo ...

  6. Tomcat6配置webdav协议

    Tomcat6默认是支持webdav协议的,只是webapps目录中不带例子而已.要添加支持很简单,在webapps目录中建立webdav目录,在webdav目录中建立WEB-INF目录,加入web. ...

  7. JS多异步之间的协作方案

    场景:使用工具函数downloadAllAsync接收一个URL数组并下载所有文件,结果返回一个存储了文件内容的数组,每个URL对应一个字符串. 好处:downloadAllAsync并不只有清理嵌套 ...

  8. eclipse: workspace出错导致无法启用的解决

    通常我们会在eclipse中创建多个workspace,比如一个用于学习,一个用于工作... ,因为种种原因,时不时会发现eclipse切换workspace后启动失败,提示让你去看workspace ...

  9. JAVA格物致知基础篇:用JAX-RS和Jersey打造RESTful Service

    随着服务器的处理能力越来越强,业务需求量的不断累积,越来越多的公司开始从单一服务器,单一业务承载变成了多服务器,多业务承载的快速扩展的过程中.传统的方法很难满足和应付这种业务量的增长和部署方式的改变. ...

  10. WebApp:如何让安卓的webview缓存webapp的html、js和图片等资源

    一.开发环境     客户端:安卓+webview(vuejs)     服务器端:tomcat 8.0 二.问题     使用安卓原生+web(基于webpack+vuejs)的方式开发了一个安卓应 ...