傳送門

解題思路

　　如果區間內沒有\(1\)，那么答案就為\(1\)，從這一點繼續歸納。如果區間內有\(x\)個\(1\)，設區間內\([2,x+1]\)的和為\(sum\)，如果\(sum=0\)，那么答案為\(x+1\)，否則\([1,x+sum]\)中的所有數字一定可以被表示，然后這個操作每次使答案至少擴大\(1\)倍，再用一個主席樹維護，時間復雜度\(O(nlognlogA)\)

代碼

#include<iostream> #include<cstdio> #include<cstring> #include<cmath> #include<algorithm>using namespace std; const int N=100005; const int M=N*33; const int inf=1000000000; typedef long long LL;template<class T> void rd(T &x){x=0;char ch=getchar();while(!isdigit(ch)) ch=getchar();while(isdigit(ch)) x=(x<<1)+(x<<3)+ch-'0',ch=getchar(); } int n,m,a[N],rt[N],ls[M],rs[M],sum[M],cnt; LL Sum[M];void build(int &x,int l,int r,int k){x=++cnt;if(l==r) {sum[x]=1;Sum[x]=l;return ;}int mid=(l+r)>>1;if(k<=mid) build(ls[x],l,mid,k);else build(rs[x],mid+1,r,k);Sum[x]=Sum[ls[x]]+Sum[rs[x]]; }void update(int pre,int &x,int l,int r,int k){x=++cnt;ls[x]=ls[pre];rs[x]=rs[pre];if(l==r) {sum[x]=sum[pre]+1;Sum[x]=Sum[pre]+l;return;}int mid=(l+r)>>1;if(k<=mid) update(ls[pre],ls[x],l,mid,k);else update(rs[pre],rs[x],mid+1,r,k);Sum[x]=Sum[ls[x]]+Sum[rs[x]]; }int query_tot(int u,int v,int l,int r,int k){if(l==r) return sum[v]-sum[u];int mid=(l+r)>>1;if(k<=mid) return query_tot(ls[u],ls[v],l,mid,k);else return query_tot(rs[u],rs[v],mid+1,r,k); }LL query_sum(int u,int v,int l,int r,int L,int R){if(L<=l && r<=R) return Sum[v]-Sum[u];int mid=(l+r)>>1;LL ret=0;if(L<=mid) ret+=query_sum(ls[u],ls[v],l,mid,L,R);if(mid<R) ret+=query_sum(rs[u],rs[v],mid+1,r,L,R);return ret; }int main(){rd(n);rd(a[1]);build(rt[1],1,inf,a[1]);for(int i=2;i<=n;i++)rd(a[i]),update(rt[i-1],rt[i],1,inf,a[i]);rd(m);int l,r,now,k,tot,lst;while(m--){rd(l),rd(r);k=0;lst=0;while(1){now=query_sum(rt[l-1],rt[r],1,inf,lst,k+1);if(!now) break;lst=k+2;k=now+k;}printf("%d\n",k+1);}return 0; }

轉載于:https://www.cnblogs.com/sdfzsyq/p/10254219.html

總結

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

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