本题主要靠结论

12pt

爆搜

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

36pt

\(f_{i,j}表示前i个数由状态j转移过来,a_i表示前缀和\)

\(So,f_{i,j}=f_{j,k}+(a_i-a_j)^2\ \ \ \ \ \ \ (a_j-a_k\leq a_i-a_j)\)

时间复杂度\(O(n^3)\)

64pt

我们发现,在\(i\)变大过程中,每个\(j\)对应的\(k\)也只会变大,所以用一个\(g_j\)表示当前\(j\)状态的\(k\)动到哪儿了,再配合后缀最小值均摊转移\(O(1)\)

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

100pt

结论:当\(j\)一定时取\(k\)最大的,并且要合法。

证明

所以:\(f_i=min \{f_j+(a_i-a_j)^2\}\ \ \ \ \ \ \ \ \ \ (a_j-a_{g_j} \leq a_i-a_j)\)

我们把条件变换一下\(\Rightarrow\) \(2*a_j-a_{g_j}\leq a_i\),发现可以用单调队列优化

\(\mathfrak{Talk\ is\ cheap,show\ you\ the\ code.}\)

#include<cstdio>
#include<algorithm>
using namespace std;
# define Type template<typename T>
# define read read1<int>()
Type inline T read1()
{
T t=0;
bool ty=0;
char k;
do k=getchar(),(k=='-')&&(ty=1);while('0'>k||k>'9');
do t=(t<<3)+(t<<1)+(k^'0'),k=getchar();while('0'<=k&&k<='9');
return ty?-t:t;
}
# define fre(k) freopen(k".in","r",stdin);freopen(k".out","w",stdout)
# define ll long long
# define I128 __int128
int s,p[100001],l[100001],r[100001],g[40000001],q[40000001];
# define ret(n) (a[n]*2-a[g[n]])
ll a[40000002];
bool ty;
void into()
{
s=read;ty=read;
if(!ty)for(int i=0;i++^s;)a[i]=read;
else
{
ll x=read,y=read,z=read;
int now=0,b[2];
b[0]=read,b[1]=read;
int m=read;
for(int i=0;i++^m;p[i]=read,l[i]=read,r[i]=read);
for(int i=0;i++^s;)
{
while(p[now]<i)++now;
if(i<=2)a[i]=b[i-1]%(r[now]-l[now]+1)+l[now];
else
{
b[0]^=b[1]^=(b[0]=(y*b[0]+x*b[1]+z)%(1<<30))^=b[1];
a[i]=b[1]%(r[now]-l[now]+1)+l[now];
}
}
}
}
void work()
{
int l=0,r=0;
for(int i=0;i++^s;a[i]+=a[i-1]);
for(int i=0;i++^s;)
{
while(l<r&&ret(q[l+1])<=a[i])++l;
g[i]=q[l];
while(l<r&&ret(q[r])>=ret(i))--r;
q[++r]=i;
}
I128 ans=0;
while(s)ans+=(I128)(a[s]-a[g[s]])*(a[s]-a[g[s]]),s=g[s];
r=0;
do
{
q[++r]=ans%10;
ans=ans/10;
}while(ans);
do
{
printf("%d",q[r]);
}while(--r);
}
int main()
{
//fre("partition");
into();
work();
return 0;
}

最新文章

  1. phonegap android 开发环境搭建
  2. 说说设计模式~适配器模式(Adapter)
  3. js 小数[非]四舍五入
  4. css sprite
  5. struts2视频学习笔记 24-27(国际化)
  6. swift-基础部分
  7. 三种情形容易引起Azure虚拟机重新启动
  8. 【转】jQuery列表拖动排列-jquery list dragsort插件参数和使用方法
  9. OC——NSArray和NSMutableArray
  10. cpu卡,sam卡原理
  11. 反射---hasattr、getattr、setattr、delattr
  12. 优先队列 poj3253 Fence Repair
  13. lambda高级进阶--组合函数
  14. Jquery的ID选择器
  15. 序列化时提示There was an error reflecting type &#39;System.Collections.Generic.List`1
  16. 最短路(Bellman)-Hdu1874畅通工程序
  17. 转:关于将Java编译过的.class文件打成jar可执行文件/JAR详解
  18. ZooKeeper分布式过程协同技术详解2——了解ZooKeeper
  19. android 开发 实现一个ListView套嵌GirdView的滚动布局
  20. go语言指针判等

热门文章

  1. office2016下载安装
  2. disconf的简单使用与远程配置更改为使用本地配置
  3. linux echo命令颜色显示
  4. 源码安装rlwrap 0.43(为了方便使用linux下的sqlplus)
  5. 二十、Python与Mysql交互
  6. 在spring boot中三分钟上手apache顶级分布式链路追踪系统skywalking
  7. JavaScript/JQuery自执行函数
  8. LaTeX公式手册(全网最全)
  9. APIO2019游记
  10. hdu5285-wyh2000 and pupil-(染色法二分图判定)