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

P9403 [POI 2020/2021 R3] Les Bitrables

题意简明,不再阐述。

首先可以对当前两行(假设为第 \(i\)\(i+1\) 行)的情况分类。

  1. \(s_i\leq s_{i+1}\)

此时可以分为三种情况。

一种是从 0 处调 \(x\) 件物品(\(0\leq x\leq s_{i+1}\)),这 \(x\) 件物品显然对应第 \(i+1\)\(p_1,p_2,...,p_x\),然后第 \(i\) 行的后 \(s_i+x-s_{i+1}\) 个物品去到 \(d\)

另一种是从 \(d\) 处调 \(x\) 件物品(\(0\leq x \leq s_{i+1}\)),这 \(x\) 件物品显然对应第 \(i+1\)\(p_{s_{i+1}-x+1},p_{s_{i+1}-x+2},...,p_{s_{i+1}}\),然后第 \(i\) 行的前 \(s_i+x-s_{i+1}\) 个物品去到 0。

还有一种是从 0 处调 \(x\) 件物品,从 \(d\) 处调 \(y\) 件物品(\(0\leq x,0\leq y,x+y+s_i=s_{i+1}\)),分别对应 \(i+1\)\(p_1,p_2,...,p_x\)\(p_{s_{i+1}-y+1},p_{s_{i+1}-y+2},p_{s_{i+1}}\),第 \(i\) 行的物品则对应第 \(i+1\) 行的 \(p_{x+1},p_{x+2},...,p_{s_{i+1}-y}\)

发现最小代价无论在哪种情况中取到其代价随 \(x\) 变化所形成的图像都是单谷的(第三种情况看作 \(y=s_{i+1}-s_{i}-x\)),于是可以三分,对于每一种情况都三分 \(x\) 找到最小代价最后三种情况取最小即可。

  1. \(s_i>s_{i+1}\)

此时也是分三种情况,前两种情况与 \(s_i\leq s_{i+1}\) 时的前两种情况相同,第三种情况变为调前 \(x\) 个元素到 0,调后 \(y\) 个元素到 \(d\)\(0<x,0<y,s_i-x-y=s_{i+1}\))。

最小代价无论在哪种情况中取到其代价随 \(x\) 变化所形成的图像也是单谷的,同上三分求最小代价即可。

时间复杂度 \(O(n\log_3n)\)

#include<bits/stdc++.h>
using namespace std;
#define mp make_pair
#define int long long
#define db double
#define endl '\n'
#define lowbit(x) x&-x
#define intz(x,a) memset(x,a,sizeof(x))
const int N=5e5+5; 
int s[N];vector<int>p[N];
signed main(){int n,d;cin>>n>>d;for(int i=1;i<=n;i++){cin>>s[i];p[i].resize(s[i]+5);for(int j=1;j<=s[i];j++)cin>>p[i][j];}for(int i=1;i<n;i++){if(s[i]<=s[i+1]){int l=0,r=s[i+1],ans=(1ll<<63)-1;while(l<=r){int mid0=l+(r-l)/3,mid1=r-(r-l)/3,sum0=0,sum1=0;for(int j=1;j<=s[i+1];j++)if(j<=mid0)sum0+=p[i+1][j];else if(j<=mid0+s[i])sum0+=abs(p[i+1][j]-p[i][j-mid0]);else sum0+=d-p[i+1][j];for(int j=s[i+1]-mid0+1;j<=s[i];j++)sum0+=min(p[i][j],d-p[i][j]);for(int j=1;j<=s[i+1];j++)if(j<=mid1)sum1+=p[i+1][j];else if(j<=mid1+s[i])sum1+=abs(p[i+1][j]-p[i][j-mid1]);else sum1+=d-p[i+1][j];for(int j=s[i+1]-mid1+1;j<=s[i];j++)sum1+=min(p[i][j],d-p[i][j]);if(sum0>=sum1)l=mid0+1,ans=min(ans,sum1);else r=mid1-1,ans=min(ans,sum0);}l=0,r=s[i+1];while(l<=r){int mid0=l+(r-l)/3,mid1=r-(r-l)/3,sum0=0,sum1=0;for(int j=s[i+1];j;j--)if(j>=s[i+1]-mid0+1)sum0+=d-p[i+1][j];else if(j>=s[i+1]-mid0-s[i]+1)sum0+=abs(p[i+1][j]-p[i][j-s[i+1]+mid0+s[i]]);else sum0+=p[i+1][j];for(int j=1;j<=-s[i+1]+mid0+s[i];j++)sum0+=min(p[i][j],d-p[i][j]);for(int j=s[i+1];j;j--)if(j>=s[i+1]-mid1+1)sum1+=d-p[i+1][j];else if(j>=s[i+1]-mid1-s[i]+1)sum1+=abs(p[i+1][j]-p[i][j-s[i+1]+mid1+s[i]]);else sum1+=p[i+1][j];for(int j=1;j<=-s[i+1]+mid1+s[i];j++)sum1+=min(p[i][j],d-p[i][j]); if(sum0>=sum1)l=mid0+1,ans=min(ans,sum1);else r=mid1-1,ans=min(ans,sum0);}cout<<ans<<endl;}else{int l=0,r=s[i],ans=(1ll<<63)-1;while(l<=r){int mid0=l+(r-l)/3,mid1=r-(r-l)/3,sum0=0,sum1=0;for(int j=1;j<=s[i];j++)if(j<=mid0)sum0+=p[i][j];else if(j<=mid0+s[i+1])sum0+=abs(p[i][j]-p[i+1][j-mid0]);else sum0+=d-p[i][j];for(int j=s[i]-mid0+1;j<=s[i+1];j++)sum0+=min(p[i+1][j],d-p[i+1][j]);for(int j=1;j<=s[i];j++)if(j<=mid1)sum1+=p[i][j];else if(j<=mid1+s[i+1])sum1+=abs(p[i][j]-p[i+1][j-mid1]);else sum1+=d-p[i][j];for(int j=s[i]-mid1+1;j<=s[i+1];j++)sum1+=min(p[i+1][j],d-p[i+1][j]);if(sum0>=sum1)l=mid0+1,ans=min(ans,sum1);else r=mid1-1,ans=min(ans,sum0);}l=0,r=s[i];while(l<=r){int mid0=l+(r-l)/3,mid1=r-(r-l)/3,sum0=0,sum1=0;for(int j=s[i];j;j--)if(j>=s[i]-mid0+1)sum0+=d-p[i][j];else if(j>=s[i]-mid0-s[i+1]+1)sum0+=abs(p[i][j]-p[i+1][j-s[i]+mid0+s[i+1]]);else sum0+=d-p[i][j];for(int j=1;j<=-s[i]+mid0+s[i+1];j++)sum0+=min(p[i+1][j],d-p[i+1][j]);for(int j=s[i];j;j--)if(j>=s[i]-mid1+1)sum1+=d-p[i][j];else if(j>=s[i]-mid1-s[i+1]+1)sum1+=abs(p[i][j]-p[i+1][j-s[i]+mid1+s[i+1]]);else sum1+=d-p[i][j];for(int j=1;j<=-s[i]+mid1+s[i+1];j++)sum1+=min(p[i+1][j],d-p[i+1][j]);if(sum0>=sum1)l=mid0+1,ans=min(ans,sum1);else r=mid1-1,ans=min(ans,sum0);}cout<<ans<<endl;}}return 0;
}
http://www.zskr.cn/news/22066.html

相关文章:

  • P5609 [Ynoi2013] 对数据结构的爱
  • 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属性
  • 《膜拜!适合新手入门的卷积神经网络原理详解教程》读书报告