当前位置: 首页 > news >正文

P5609 [Ynoi2013] 对数据结构的爱

线段树维护 \(c_i\) 表示最小初值使根据题意经过节点表示区间 \([l,r]\) 过程中共减去 \(i\cdot p\),区间 \([l,r]\) 操作过程中最多减去 \((r-l+1)\cdot p\)

区间合并就是

\[tr_{u,c_{x+y}}=\min\{\max\{tr_{ls,{c_x}},tr_{rs,{c_y}}-sum_{ls}+x\cdot p\}\} \]

,其中 \(x,y\) 的合法条件是,满足初值为 \(tr_{ls,c_x}\) 的最大值即 \(tr_{ls,c_{x+1}}-1\),在操作完 \(ls\) 表示区间后大于等于 \(tr_{rs,c_y}\)\(tr_{ls,c_{x+1}}-1+sum_{ls}-x\cdot p \ge tr_{rs,c_y}\) 等价 \(tr_{ls,c_{x+1}}+sum_{ls}-x\cdot p > tr_{rs,c_y}\)

\(tr_{ls,c_x}\) 即满足左区间减 \(x\cdot p\) 的最小初值,
\(tr_{rs,c_{x+1}}\) 即满足右区间减 \(y\cdot p\) 的最小初值与操作完左区间的值的差值加上满足左区间条件的最小初值 \(tr_{ls,c_x}\)\(tr_{rs,c_y}-(tr_{ls,c_x}+sum_{ls}-x\cdot p)+tr_{ls,c_x}=tr_{rs,c_y}-sum_{ls}+x\cdot p\)

这样的话时间复杂度接近 \(O(qn^2)\)(?)。

\(tr_{ls,c_{x+1}}-1+sum_{ls}-x\cdot p \ge tr_{rs,c_y}\Leftrightarrow tr_{ls,c_{x+1}}-1 \ge tr_{rs,c_y}-sum_{ls}+x\cdot p\Leftrightarrow \max(tr_{ls,c_x},tr_{rs,c_y}-sum_{ls}+x\cdot p)<tr_{ls,c_x+1}\),则 \(\max(tr_{ls,c_x},tr_{rs,c_y}-sum_{ls}+x\cdot p)<\max(tr_{ls,c_{x+1}},tr_{rs,c_{y-1}}-sum_{ls}+(x+1)\cdot p)\) 则对于更新 \(c_{x+y},x\) 取最小值时则更新值最小,所以 \(c_{x+y}\) 只用 \(x,y\) 合法且 \(x\) 最小的状态转移,所以对于 \(x,y\) 的枚举范围只是所有 \(y\) 使得 \(x-1,y\) 不合法但 \(x,y\) 合法即可。

貌似不用证明 \(tr_{ls,c_{x+1}}-1+sum_{ls}-x\cdot p\) 具有单调性,因为 \(tr_{rs,c_y}\) 显然具有单调性,那么对于 任意 \(x\) 合法的 \(y\) 必然是一个前缀,而 \(y\) 只用枚举对于 \(x-1\) 不合法,\(x\) 合法的状态,所以只需从 \(x-1\) 第一个不合法的状态往后枚举即可,所以可以用双指针解决。

时间复杂度 \(O(n\log n+m\log^2n)\)

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define endl '\n'
#define max(x,y) (x)>(y)?(x):(y)
#define min(x,y) (x)<(y)?(x):(y)
const int N=1e6+5;int a[N];ll p;
inline ll read(){ll x=0,f=1;char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar();return x*f;
}
struct node{int l,r,lf,rt;ll sum;}tr[N<<2];int tot;ll q[N*30];
inline void pushup(int u){int mid=tr[u].l+tr[u].r>>1;tr[u].sum=tr[u<<1].sum+tr[u<<1|1].sum;for(int i=0,j=0;i<=mid-tr[u].l+1;++i)for((j>tr[u].r-mid?--j:0);j<=tr[u].r-mid;++j){if(q[tr[u<<1].lf+i+1]-i*p+tr[u<<1].sum<=q[tr[u<<1|1].lf+j]){--j;break;}q[tr[u].lf+i+j]=max(q[tr[u<<1].lf+i],q[tr[u<<1|1].lf+j]+i*p-tr[u<<1].sum);}
}
void build(int u,int l,int r){tr[u]={l,r,0};tr[u].lf=tot+1;for(int i=0;i<=r-l+3;i++)q[++tot]=1e18;tr[u].rt=tot;q[tr[u].lf]=-1e18;if(l==r){q[tr[u].lf+1]=p-a[l],tr[u].sum=a[l];return;}int mid=l+r>>1;build(u<<1,l,mid),build(u<<1|1,mid+1,r);pushup(u);
}
ll query(int u,int l,int r,ll v){if(tr[u].l<l||tr[u].r>r){int mid=tr[u].l+tr[u].r>>1;if(l<=mid)v=query(u<<1,l,r,v);if(r>mid)v=query(u<<1|1,l,r,v);return v;}return v+tr[u].sum-p*(upper_bound(q+tr[u].lf,q+1+tr[u].rt,v)-q-tr[u].lf-1);
}
int main(){int n,m;ll lst=0;n=read(),m=read(),p=read();for(int i=1;i<=n;++i)a[i]=read();build(1,1,n);//cout<<q[tr[1].lf+n]<<endl;while(m--){int l,r;ll x;l=read()^lst,r=read()^lst,x=read()^lst;lst=query(1,l,r,x);printf("%lld\n",lst);lst=(lst%n+n)%n;}return 0;
}
http://www.zskr.cn/news/22064.html

相关文章:

  • STM32 代码
  • 剑指offer-35、数组中的逆序对
  • 2025 年太阳能厂家最新推荐:全场景系统企业综合实力榜,含热水 / 发电 / 光伏热等领域优质品牌测评
  • 完整教程:AI应用生成平台:数据库、缓存与存储
  • 2025 年电缆桥架生产厂家最新推荐排行榜:含北方 / 河北 / 瓦楞 / 防火 / 模压 / 镀锌桥架品牌及合作案例盘点
  • 2025 年胰岛素泵厂家最新推荐排行榜:国产实力厂家技术、口碑及全场景适配方案全景解析软针植入/平衡式留置针/无异物感胰岛素泵厂家推荐
  • 进程的内存管理
  • 深入理解Java内存模型与volatile关键字:从理论到实践
  • 312、金缕衣
  • AVG Clear:彻底卸载AVG产品的专业工具
  • 深入解析:安卓 WPS Office v18.21.0 国际版
  • 科学背景如何赋能云计算业务战略
  • .netframework中自带的dll
  • 【GitHub每日速递 251016】23k star,Daytona:90ms内极速运行AI代码,安全弹性基础设施来袭!
  • MVCC、幻读、间隙锁与临键锁(三)
  • 生物聚酯塑料回收技术创新与商业应用
  • # ️ MySQL vs PostgreSQL架构深度对比分析报告
  • # 韩国数据中心大火:647套系统因缺失双活集体宕机22小时
  • LlamaIndex API Example
  • 图文并茂展示CSS li 排版大合集,总有一款是你刚好需要的
  • The lamentable decline of reading
  • [FT.COM]The world should prepare for the looming quantum era
  • 函数的类型注释器
  • DshanPI-A1 RK3576 gmrender-resurrect B站投屏
  • 组件级异步加载与预加载策略
  • SAM系列论文浅析
  • 为什么一部电影,一本书一看就喜欢
  • 02020601 Web API01-顶级语句、全局using指令、可空类型、record类型(自动重写ToString、Equals)、init和private属性
  • 《膜拜!适合新手入门的卷积神经网络原理详解教程》读书报告
  • Linux系统下对拍