江苏网站建设哪家好,免费咨询医生在线男科,如何更改地图上的店名,给别人做网站收钱违法吗正题
题目链接:https://ac.nowcoder.com/acm/contest/11161/E 题目大意
给出nnn个数字的一个序列#xff0c;mmm个操作。
给出l,r,kl,r,kl,r,k#xff0c;求一个最大的xxx使得∑ilrmax{ai−x,0}≥k\sum_{il}^rmax\{a_i-x,0\}\geq k∑ilrmax{ai−x,0}≥k单点修改 解题思…正题
题目链接:https://ac.nowcoder.com/acm/contest/11161/E 题目大意
给出nnn个数字的一个序列mmm个操作。
给出l,r,kl,r,kl,r,k求一个最大的xxx使得∑ilrmax{ai−x,0}≥k\sum_{il}^rmax\{a_i-x,0\}\geq k∑ilrmax{ai−x,0}≥k单点修改 解题思路
带修的比较麻烦用带修莫队的话需要平衡一下时间复杂度可以用分块来做。
这样修改是O(1)O(1)O(1)的但是询问的话朴素的想法是二分然后统计这个O(mnlogn)O(m\sqrt n\log n)O(mnlogn)显然是过不了的。
但是如果改为一个个块从后往前跳确定答案在哪个块然后在块里枚举就好了。
时间复杂度O(mn23mnlogn)O(mn^{\frac{2}{3}}m\sqrt n\log n)O(mn32mnlogn) code
#includecstdio
#includecstring
#includealgorithm
#includecmath
#define ll long long
using namespace std;
const ll N1e510;
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.ly.l;if(x.r/T!y.r/T)return x.ry.r;return x.ty.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,sum0,ct0;for(ptQ;pt1;pt--){sums[pt];ctcnt[pt];if(sum-ct*R[pt-1]k){sum-s[pt];ct-cnt[pt];break;}}if(!pt)return -1;for(ll iR[pt];iL[pt];i--){sumv[i]*i;ctv[i];if(sum-ct*(i-1)k)return i-1;}return -1;
}
signed main()
{scanf(%lld%lld,n,m);for(ll i1;in;i)scanf(%lld,a[i]);Tpow(n*m,1.0/3.0);Q316;for(ll i1;iQ;i)L[i]R[i-1]1,R[i]i*Q;Q;L[Q]R[Q-1]1;R[Q]1e5;for(ll i1;iQ;i)for(ll jL[i];jR[i];j)pos[j]i;for(ll i1;im;i){ll op;scanf(%lld,op);if(op0){qnt;q[qnt].idqnt;q[qnt].tmnt;scanf(%lld%lld%lld,q[qnt].l,q[qnt].r,q[qnt].k);}else mnt,scanf(%lld%lld,p[mnt],c[mnt]);}sort(q1,q1qnt,cmp);ll l1,r0,t0;for(ll i1;iqnt;i){while(lq[i].l)Add(a[l],-1),l;while(lq[i].l)l--,Add(a[l],1);while(rq[i].r)r,Add(a[r],1);while(rq[i].r)Add(a[r],-1),r--; while(tq[i].t){t;if(lp[t]p[t]r)Add(a[p[t]],-1),Add(c[t],1);swap(a[p[t]],c[t]);}while(tq[i].t){swap(a[p[t]],c[t]);if(lp[t]p[t]r)Add(a[p[t]],1),Add(c[t],-1);t--;}ans[q[i].id]Query(q[i].k);}for(ll i1;iqnt;i)printf(%lld\n,ans[i]);return 0;
}
