淺談主席樹：https://www.cnblogs.com/AKMer/p/9956734.html

題目傳送門：https://www.lydsy.com/JudgeOnline/problem.php?id=4408

我們先把問題簡化，給你一個可重集，求最小的不能被湊出來的正整數。

性質一：假設原集合可以把\([1,num]\)全部湊出來，新加入一個數字\(x\)，就可以湊出\([1,num+x]\)

證明：顯然。

性質二：如果你已經可以湊出\([1,num_1]\)，當前集合里權值在\([1,num_1]\)里的所有數的權值和為\(num_2\)，那么你就可以湊出\([1,num_2]\)。

證明：因為性質一。

性質三：如果你由\([1,num_1]\)得到了\([1,num_2]\)且\(num_2>num_1\)，然后由\([1,num_2]\)得到了\([1,num_3]\)且\(num_3> num_2\)，那么\(num_3\geqslant2*num_1\)。

證明：因為由\(num_2\)轉到\(num_3\)，增加的數字必然在\([num_1,num_2]\)內，所以\(num_3\geqslant num_1+num_1\)。

因為這兩條性質，我們可以將最小的不能湊出的數不斷變大，而且是成倍增長的。先欽點不能湊出的數字為\(1\)，然后根據性質二不斷變大，時間是\(60*logn\)的。

時間復雜度：\(O(60*mlogn)\)

空間復雜度：\(O(nlogn)\)

代碼如下：

#include <cstdio> using namespace std;const int maxn=1e5+5;int n,m,sum; int a[maxn],rt[maxn];int read() {int x=0,f=1;char ch=getchar();for(;ch<'0'||ch>'9';ch=getchar())if(ch=='-')f=-1;for(;ch>='0'&&ch<='9';ch=getchar())x=x*10+ch-'0';return x*f; }struct tree_node {int sum,cnt,ls,rs; };struct chairman_tree {int tot;tree_node tree[maxn*35];void ins(int lst,int &now,int l,int r,int pos) {now=++tot;tree[now]=tree[lst];tree[now].cnt++;tree[now].sum+=pos;if(l==r)return;int mid=(l+r)>>1;if(pos<=mid)ins(tree[lst].ls,tree[now].ls,l,mid,pos);else ins(tree[lst].rs,tree[now].rs,mid+1,r,pos);}int query(int x,int y,int l,int r,int pos) {if(r<=pos)return tree[y].sum-tree[x].sum;int mid=(l+r)>>1,res=tree[tree[y].ls].sum-tree[tree[x].ls].sum;if(pos<=mid)res=query(tree[x].ls,tree[y].ls,l,mid,pos);else res+=query(tree[x].rs,tree[y].rs,mid+1,r,pos);return res;} }T;int main() {n=read();for(int i=1;i<=n;i++) {a[i]=read();sum+=a[i];}for(int i=1;i<=n;i++)T.ins(rt[i-1],rt[i],1,sum,a[i]);m=read();for(int i=1;i<=m;i++) {int l=read(),r=read(),cnt=1;while(cnt<sum) {int tmp=T.query(rt[l-1],rt[r],1,sum,cnt);if(tmp<cnt)break; else cnt=tmp+1;//cnt就是題解里講的num}printf("%d\n",cnt);}return 0; }

轉載于:https://www.cnblogs.com/AKMer/p/9971824.html

總結

以上是生活随笔為你收集整理的BZOJ4408：[FJOI2016]神秘数的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。