正題

題目鏈接:https://ac.nowcoder.com/acm/contest/11161/E

---

題目大意

給出$n$個數(shù)字的一個序列，$m$個操作。

給出$l,r,k$，求一個最大的$x$使得${\sum }_{i=l}^{r}max\{a_i?x,0\}\ge k$

單點修改

---

解題思路

帶修的比較麻煩，用帶修莫隊的話需要平衡一下時間復雜度，可以用分塊來做。

這樣修改是$O(1)$的，但是詢問的話樸素的想法是二分然后統(tǒng)計，這個$O(m\sqrt{n}\log n)$顯然是過不了的。

但是如果改為一個個塊從后往前跳確定答案在哪個塊，然后在塊里枚舉就好了。

時間復雜度$O(m{n}^{\frac{2}{3}}+m\sqrt{n}\log n)$

---

code

#include<cstdio> #include<cstring> #include<algorithm> #include<cmath> #define ll long long using namespace std; const ll N=1e5+10; struct node{ll l,r,k,id,t; }q[N]; ll n,m,T,Q,a[N],L[N],R[N],pos[N],p[N],c[N]; ll qnt,mnt,s[N],cnt[N],v[N],ans[N]; bool cmp(node x,node y){if(x.l/T!=y.l/T)return x.l<y.l;if(x.r/T!=y.r/T)return x.r<y.r;return x.t<y.t; } void Add(ll x,ll f){s[pos[x]]+=x*f;cnt[pos[x]]+=f;v[x]+=f;return; } ll Query(ll k){if(!k)return 100000;ll pt,sum=0,ct=0;for(pt=Q;pt>=1;pt--){sum+=s[pt];ct+=cnt[pt];if(sum-ct*R[pt-1]>=k){sum-=s[pt];ct-=cnt[pt];break;}}if(!pt)return -1;for(ll i=R[pt];i>=L[pt];i--){sum+=v[i]*i;ct+=v[i];if(sum-ct*(i-1)>=k)return i-1;}return -1; } signed main() {scanf("%lld%lld",&n,&m);for(ll i=1;i<=n;i++)scanf("%lld",&a[i]);T=pow(n*m,1.0/3.0);Q=316;for(ll i=1;i<=Q;i++)L[i]=R[i-1]+1,R[i]=i*Q;++Q;L[Q]=R[Q-1]+1;R[Q]=1e5;for(ll i=1;i<=Q;i++)for(ll j=L[i];j<=R[i];j++)pos[j]=i;for(ll i=1;i<=m;i++){ll op;scanf("%lld",&op);if(op==0){++qnt;q[qnt].id=qnt;q[qnt].t=mnt;scanf("%lld%lld%lld",&q[qnt].l,&q[qnt].r,&q[qnt].k);}else ++mnt,scanf("%lld%lld",&p[mnt],&c[mnt]);}sort(q+1,q+1+qnt,cmp);ll l=1,r=0,t=0;for(ll i=1;i<=qnt;i++){while(l<q[i].l)Add(a[l],-1),l++;while(l>q[i].l)l--,Add(a[l],1);while(r<q[i].r)r++,Add(a[r],1);while(r>q[i].r)Add(a[r],-1),r--; while(t<q[i].t){t++;if(l<=p[t]&&p[t]<=r)Add(a[p[t]],-1),Add(c[t],1);swap(a[p[t]],c[t]);}while(t>q[i].t){swap(a[p[t]],c[t]);if(l<=p[t]&&p[t]<=r)Add(a[p[t]],1),Add(c[t],-1);t--;}ans[q[i].id]=Query(q[i].k);}for(ll i=1;i<=qnt;i++)printf("%lld\n",ans[i]);return 0; }

總結

以上是生活随笔為你收集整理的牛客挑战赛48E-速度即转发【带修莫队,分块】的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。