题目大意:给定一个长度为n的序列,至多将序列分成m+1段,每段序列都有权值,权值为序列内两个数两两相乘之和。m<=n<=1000.

分析:令w[i,j]表示区间[i,j]中两两乘积之和,f[i][j]表示前j个数分成i段的最小值。

f[i][j]=f[i-1][k]+w[k+1,j]
w[k+1,j]可以转换为w[1,j]-w[1,k]-sum[k]*(sum[j]-sum[k])
其中sum[j]表示前j个数的前缀和。
f[i][j]=f[i-1][k]+w[j]-w[k]-sum[k]*(sum[j]-sum[k])
令y=f[i-1][k]+w[j]-w[k]+sum[k]^2,x=sum[k],b=sum[j],g=f[i][j],则有:

y-bx=g

此为直线方程,b为定值,要求g最小,即为直线的截距最小。平面上有若干点(x,y),这些点是由各个决策点产生的。而将直线从下往上平移,它接触到的第一个点即为最佳决策点。因为斜率b是上升的,所以,下一阶段的直线方程斜率更高,于是最佳决策点一定形成了下凸包序列。

#include<iostream>
#include<cstdio>
#include<cstring>
#define MAXN 1005
#define LL long long int
LL f[MAXN][MAXN],w[MAXN],sum[MAXN];
#define FZ(i,p) (f[i-1][p]-w[p]+sum[p]*sum[p])
int n,m,num[MAXN];
int que[MAXN],head,tail;
#define MAXZ (1LL<<22)
bool turnup(int i,int p1,int p2,int p3) //p1>p2>p3
{
LL y1=FZ(i,p1);
LL x1=sum[p1];
LL y2=FZ(i,p2);
LL x2=sum[p2];
LL y3=FZ(i,p3);
LL x3=sum[p3];
if((x2-x3)*(y1-y2)>(x1-x2)*(y2-y3))return ;
else return ;
}
int main()
{
while(scanf("%d%d",&n,&m)&&(n||m))
{
memset(sum,,sizeof sum);
memset(w,,sizeof w);
m++;
for(int i=;i<=n;i++)
{scanf("%d",&num[i]);
sum[i]=sum[i-]+num[i];
w[i]=w[i-]+sum[i-]*num[i];
}
for(int i=;i<=n;i++)
f[][i]=w[i];
for(int i=;i<=m;i++)
{
//f[i-1][i-1]=0;
head=tail=;
que[tail++]=i-;
for(int j=i;j<=n;j++)
{
while(head<tail-&&FZ(i,que[head+])-FZ(i,que[head])<sum[j]*(sum[que[head+]]-sum[que[head]]))head++;
int k=que[head];
f[i][j]=f[i-][k]+w[j]-w[k]-sum[k]*(sum[j]-sum[k]);
while(head<tail-&&turnup(i,j,que[tail-],que[tail-])==)
tail--;
que[tail++]=j;
}
}
/*
for(int i=1;i<=m;i++)
{for(int j=1;j<=n;j++)
printf("%I64d ",f[i][j]);
printf("\n");
}
*/
printf("%I64d\n",f[m][n]);
}
}

本题也可以用平行四边形优化。

f[i][j]=f[i-1][k]+w[k+1,j]

w[k+1,j]很明显满足区间包含性质和平行四边形性质,所以f[i][j]也满足平行四边形性质,所以设s[i][j]表示f[i][j]的最佳决策点。s[i-1][j]<=s[i][j]<=s[i][j+1]。

从常识角度来思考,也是比较好想的。

如果数组元素不变,分段数增加,肯定分界点会往两边扩展,若分段数减少,则分界点会往中间靠拢,所以s[i][j]<=s[i+1][j]

如果分段数不变,数的个数增加(右边补充进来),分界点应该往右微调;如果数的个数减少(从右边剔除),分界点应该往左微调,所以有s[i][j-1]<=s[i][j]

所以s[i-1][j]<=s[i][j]<=s[i][j+1]

 #include<iostream>
#include<cstring>
#include<cstdio>
#define MAXN 1005
#define MAXM 1005
#define LL long long int
#define min(a,b) ((a)<(b)?(a):(b))
LL f[MAXN][MAXM],sum[MAXN],w[MAXN],s[MAXN][MAXM];
int num[MAXN],t,n,m;
using namespace std;
void pre()
{
for(int i=;i<=n;i++)
sum[i]=sum[i-]+num[i];
for(int i=;i<=n;i++)
w[i]=w[i-]+sum[i-]*num[i];
}
int main()
{
while(scanf("%d%d",&n,&m)&&(n||m))
{
m++;
for(int i=;i<=n;i++)
scanf("%d",&num[i]);
memset(f,0x3f,sizeof f);
pre();
for(int i=;i<=n;i++)
{for(int j=i;j<=n;j++)
s[i][j]=i-;
}
for(int i=;i<=n;i++)
s[i][m+]=i-;
for(int i=;i<=n;i++)
f[i][i]=;
for(int j=;j<=n;j++)
f[j][]=w[j]; for(int i=;i<=n;i++)
for(int j=min(m,i);j>=;j--)
{
for(int k=s[i][j+];k>=s[i-][j];k--)
{
LL temp=f[k][j-]+w[i]-w[k]-sum[k]*(sum[i]-sum[k]);
if(f[i][j]>temp)
{f[i][j]=temp;
s[i][j]=k;
}
}
}
printf("%I64d\n",f[n][m]);
}
}

HDU2829的更多相关文章

  1. HDU2829 Lawrence —— 斜率优化DP

    题目链接:https://vjudge.net/problem/HDU-2829 Lawrence Time Limit: 2000/1000 MS (Java/Others)    Memory L ...

  2. HDU2829 Lawrence(斜率优化dp)

    学了模板题之后上网搜下斜率优化dp的题目,然后就看到这道题,知道是斜率dp之后有思路就可以自己做不出来,要是不事先知道的话那就说不定了. 题意:给你n个数,一开始n个数相邻的数之间是被东西连着的,对于 ...

  3. HDU-2829 Lawrence (DP+四边形不等式优化)

    题目大意:有n个敌方军火库呈直线排列,每个军火库有一个值vi,并且任意相邻的两个库之间都有通道相连.对于任意一条连起来的军火库链,它对我方的威胁可以用函数w(i,j)表示为:w(i,j)=vi*sum ...

  4. hdu2829 Lawrence

    题目链接:戳我 朴素DP:\(dp[i][j]=dp[i-1][k]+cost[k+1][j]\) 其中dp[i][j]表示炸第i次的时候,处理到前j个的最小值是多少.cost[i][j]表示的是i, ...

  5. hdu2829 四边形优化dp

    Lawrence Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total S ...

  6. [HDU2829] Lawrence [四边形不等式优化dp]

    题面: 传送门 思路: 依然是一道很明显的区间dp 我们设$dp\left[i\right]\left[j\right]$表示前$j$个节点分成了$i$块的最小花费,$w\left[i\right]\ ...

  7. 斜率优化板题 HDU2829 Lawrence

    题目大意:给定一个长度为nnn的序列,至多将序列分成m+1m+1m+1段,每段序列都有权值,权值为序列内两个数两两相乘之和.求序列权值和最小为多少? 数据规模:m<=n<=1000.m&l ...

  8. 单调队列 && 斜率优化dp 专题

    首先得讲一下单调队列,顾名思义,单调队列就是队列中的每个元素具有单调性,如果是单调递增队列,那么每个元素都是单调递增的,反正,亦然. 那么如何对单调队列进行操作呢? 是这样的:对于单调队列而言,队首和 ...

  9. DP总结 ——QPH

    常见优化 单调队列 形式 dp[i]=min{f(k)} dp[i]=max{f(k)} 要求 f(k)是关于k的函数 k的范围和i有关 转移方法 维护一个单调递增(减)的队列,可以在两头弹出元素,一 ...

随机推荐

  1. oc小总结

    oc的一些总结 下面几个问题是oc中需要掌握的内容 1.如何掌握一个方法的方法名 2.一个对象调用一个autorelease,什么时候释放 3.字典和数组,集合都有什么特点 4.如何定义一个类 5.类 ...

  2. Matlab位运算笔记

    本文为转载其他地方的文章; MATLAB函数 1.matlab函数bitset 设置数的某一位二进制位为1. <Simulink与信号处理> 使用方法 C = bitset(A,bit) ...

  3. ES6-函数扩展

  4. jsp页面 date转化成string

    1.先引入fmt标签<%@ taglib prefix="fmt" uri="http://java.sun.com/jsp/jstl/fmt"%> ...

  5. 当我把电脑500G硬盘格式化了的时候,在想些什么

    所以我就想重新尝试一下怎么装多系统,bootloader,MBR到底是怎么回事儿. 简要记录一下: 首先,在可以作为启动盘的U盘里面,备份一份可以用的WinXP镜像+校园网客户端安装文件+驱动精灵安装 ...

  6. 网站不能访问(httperrLog【Timer_MinBytesPerSecond】【Timer_ConnectionIdle】)(转载)

    在\LogFiles\HTTPERR的日志(C:\Windows\System32\LogFiles\HTTPERR)中发现了大量Timer_MinBytesPerSecond,Timer_Conne ...

  7. php5.2.6+apache2.2.15配置

    首先下载软件,忘记php下载地址了,apache是官网. 文件名 httpd-2.2.15-win32-x86-openssl-0.9.8m-r2.msi php-5.2.6-win32-instal ...

  8. [转]ssh常用用法小结

    ssh常用用法小结 1.连接到远程主机: 命令格式 : ssh name@remoteserver 或者 ssh remoteserver -l name 说明:以上两种方式都可以远程登录到远程主机, ...

  9. Linux系统编程@多线程编程(一)

    多线程编程 涉及操作系统原理概念 时间片 进程状态 上下文: 对进程来说,就是进程的执行环境,具体就是各个变量和数据,包括所有的寄存器变量.打开的文件.内存信息等. 进程的写时复制:由于一般 fork ...

  10. 论文笔记之:DeepCAMP: Deep Convolutional Action & Attribute Mid-Level Patterns

    DeepCAMP: Deep Convolutional Action & Attribute Mid-Level Patterns CVPR 2016 本文提出一种 分割图像 patch 的 ...