P4587 [FJOI2016]神秘数

题目大意

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

n<=1e5,m<=1e5，a[i]的总和不超过1e9；

思路

首先推得这样一个结论：

如果我已经能表示出区间[1,maxx]，那么我们考虑下一个加进来的数x

• 如果x>=maxx+2，那么他对答案是没有贡献的，最小不能被表示出来的数还是maxx+1；

• 如果x<=maxx+1，那么对于原来区间的所有数，都可以加上x，此时区间应该为[1,maxx+x]

设当前值域[1,x]

则ans=x+1

若小于等于ans的数的和res≥ans*，则一定有未选的且小于等于ans的数，

令ans=res+1即可。

反之说明答案就是ans，直接break

因为值域很大，可能有1e9，且是区间查询，权值线段树+区间查询=主席树

主席树就是可持久化的前缀权值线段树

#include<bits/stdc++.h>
#define ll long long
using namespace std;

int cnt,root[1000001],n,m;
const int inf=1e9;
struct node{
    int l,r,sum;
}seg[5000001];
//一颗线段树维护[1,sum]这个值域中所有出现过得数的和
//第i棵树表示桶中加入了前i个数后的区间和。
#define mid ((lb+rb)>>1)
void update(int &rt,int lb,int rb,int x){
    seg[++cnt]=seg[rt];rt=cnt;seg[rt].sum+=x;
    if (lb==rb) return;
    (mid>=x?update(seg[rt].l,lb,mid,x):update(seg[rt].r,mid+1,rb,x));
}
//区间查询，查询[l,r]值域，i,j表示限定的区间
int query(int i,int j,int lb,int rb,int l,int r){
    if (!(seg[j].sum-seg[i].sum)||r<lb|rb<l) return 0;
    if (lb>=l&&rb<=r) return (seg[j].sum-seg[i].sum);
    return query(seg[i].l,seg[j].l,lb,mid,l,r)+query(seg[i].r,seg[j].r,mid+1,rb,l,r);
}
int main(){
    scanf("%d",&n);
    for (int i=1;i<=n;i++){
        root[i]=root[i-1];
        int x;
        scanf("%d",&x);
        update(root[i],1,inf,x);
    }
    scanf("%d",&m);
    while (m--){
        int l,r;
        scanf("%d%d",&l,&r);
        int mx=0,pos=0;
        for (int Sum;;){
            Sum=query(root[l-1],root[r],1,inf,mx+1,pos+1);
            if (!Sum) break;
            mx=pos+1;pos+=Sum;
        }
        printf("%d\n",pos+1);
    }
    return 0;
}