UVAlive 6131 dp+斜率优化

这道题和06年论文《从一类单调性问题看算法的优化》第一道例题很相似。

题意:给出n个矿的重量和位置,这些矿石只能从上往下运送,现在要在这些地方建造m个heap,要使得,sigma距离*重量最小。

思路:O(n ^ 3)的DP解法是很容易想出来的。

dp[i][j] 表示第i个矿石点是j个heap的最小花费。

dp[i][j] = min(dp[i][j] , dp[k][j - 1] + sigma(sum[i] - sum[k])) 。

其中i , j , k 分别要一重循环,所以复杂度达到10 ^ 9。

这显然是TLE的,所以需要优化。

我们可以来看状态转移方程,dp[i][j] = dp[k][j - 1] +( sum1[i] - sum1[k] ) * a[i] - (sum2[i] - sum2[k]) .其中sum1是1到i的总重量,sum2表示1到i的总重量*距离。

这样,我们就可以进行斜率优化了。

所以这一维就降成O(1)了。那总的复杂度就是O(n ^ 2)。

 

#include <set>
#include <map>
#include <stack>
#include <cmath>
#include <queue>
#include <cstdio>
#include <string>
#include <vector>
#include <iomanip>
#include <cstring>
#include <iostream>
#include <algorithm>
#define Max 2505
#define FI first
#define SE second
#define ll long long
#define PI acos(-1.0)
#define inf 0x7ffffffffffffll
#define LL(x) ( x << 1 )
#define bug puts("here")
#define PII pair<int,int>
#define RR(x) ( x << 1 | 1 )
#define mp(a,b) make_pair(a,b)
#define mem(a,b) memset(a,b,sizeof(a))
#define REP(i,s,t) for( int i = ( s ) ; i <= ( t ) ; ++ i )using namespace std;#define N 1111ll sum1[N] ;
ll sum2[N] ;
ll dp[N][N] ;
int qe[N] ;
ll a[N] , b[N] ;
int n , m ;
//分子
ll getU(int j , int k , int z){return dp[k][j - 1] + sum2[k] - (dp[z][j - 1] + sum2[z]) ;
}
//分母
ll getD(int k , int z){return sum1[k] - sum1[z] ;
}ll getDP(int i ,int j ,int k){return dp[k][j - 1] + (sum1[i] - sum1[k]) * a[i] - (sum2[i] - sum2[k]) ;
}int main() {while(cin >> n >> m){for (int i = 1 ; i <= n ; i ++ )cin >> a[i] >> b[i] ;sum1[0] = sum2[0] = 0 ;for (int i = 1 ; i <= n ; i ++ ){sum1[i] = sum1[i - 1] + b[i] ;sum2[i] = sum2[i - 1] + a[i] * b[i] ;
//            cout << sum1[i] << " " << sum2[i] << endl;}for (int i = 0 ; i <= n ; i ++ )for (int j = 0 ; j <= m ; j ++) dp[i][j] = inf ;dp[0][0] = 0 ;for (int j = 1 ; j <= m ; j ++ ){int head = 0 , tail = 0 ;qe[tail ++ ] = 0 ;for (int i = 1 ; i <= n ; i ++ ){while(head + 1 < tail && getU(j , qe[head + 1] , qe[head]) <= a[i] * getD(qe[head + 1] , qe[head]))head ++ ;dp[i][j] = getDP(i , j , qe[head]) ;while(head + 1 < tail && getU(j , i , qe[tail - 1]) * getD(qe[tail - 1] ,qe[tail - 2]) <=getU(j , qe[tail - 1] , qe[tail - 2]) * getD(i , qe[tail - 1]))tail -- ;qe[tail ++ ] = i ;}}cout << dp[n][m] << endl ;}return 0 ;
}



 

 

转载于:https://www.cnblogs.com/suncoolcat/p/3301700.html

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/360440.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

c语言错误重定义,C语言的重定义错误求解

追问那如果要用到全局变量&#xff0c;要把它写在C文件中&#xff1f;那全局变量的主要作用其实就是跨文件作用域来传递一些东西&#xff0c;写在C文件中&#xff0c;不太好在别的文件中引用吧。用的时候写extern进去&#xff1f;感觉全局变量在C里用起来比较麻烦&#xff0c…

Netty:另一种Web(套接字)服务器

如今&#xff0c; Netty已用于Internet上的各种应用程序中&#xff0c;以处理数千&#xff08;即使不是数百万&#xff09;的聊天对话&#xff0c;包括Minecraft &#xff0c; Twitter和许多其他应用程序在内的多人游戏。 但是&#xff0c;它并没有在开发企业应用程序的企业程序…

arm opcode hook

/**************************************/* 作者:半斤八兩/* 博客:http://cnblogs.com/bjblcracked/* 日期:2013-09-05 11:11/**************************************只是感兴趣&#xff0c;没有其他目的。失误之处敬请诸位大侠赐教!因工作需要,最近一直在研究arm.看到论坛也…

c语言程序整数四则运算,c语言中三个整数随机的四则运算

满意答案hors10722014.01.06采纳率&#xff1a;58% 等级&#xff1a;12已帮助&#xff1a;18274人#include #include #include #define N 10 //随机出10道题目int main(){int num1, num2, num3, count0, result,resultTrue,flag;//result:用户输入结果 resultTrue:正确结果 …

带有调试器的Apache Camel Eclipse工具

大约2个月前&#xff0c; Lars Heineman在 JBoss工具堆栈中的博客中介绍了改进的Apache Camel Eclipse工具。 在即将发布的版本中&#xff0c;他们将Camel调试器与本机Eclipse调试器集成在一起&#xff0c;因此当您使用断点时&#xff0c;您将获得Eclipse调试体验&#xff0c;…

马秀丽C语言程序设计答案pdf,C语言程序设计清华大学马秀丽刘志妩科后习题9答案.doc...

C语言程序设计清华大学马秀丽刘志妩科后习题9答案习 题 9二 编程题main(){int n,*a,i,j,t0;scanf("%d",&n);for(i0;ifor(j0;jscanf("%d",ai*nj);for(i0;ifor(j0;jif (*(ai*nj)!0) t1;if (t0) printf("is xia san jiao fang zhen");else prin…

从Gradle自动将工件提升到Maven Central

快速教程&#xff0c;无需在带有Gradle Nexus Staging插件的Nexus GUI中单击即可将Gradle项目中的工件提升/发布到Maven Central。 介绍 Maven Central &#xff08;又名“中央存储库”&#xff09;&#xff08;可能&#xff09;是世界上最大的一组开源工件&#xff0c;这些工…

c语言指针自定义函数,c语言函数指针定义,指针函数和函数指针的区别

往往&#xff0c;我们一提到指针函数和函数指针的时候&#xff0c;就有很多人弄不懂。下面就由小编详细为大家介绍C语言中函数指针&#xff0c;指针函数和函数指针之间的区别。c语言指针函数定义&#xff1a;函数指针是指向函数的指针变量。 因此“函数指针”本身首先应是指针变…

HarmonyOS的功能及场景应用

一、基本介绍 鸿蒙HarmonyOS主要应用的设备包括智慧屏、平板、手表、智能音箱、IoT设备等。具体来说&#xff0c;鸿蒙系统是一款面向全场景(移动办公、运动健康、社交通信、媒体娱乐等)的分布式操作系统&#xff0c;能够支持手机、平板、智能穿戴、智慧屏、车机等多种终端设备…

java DOM4J 读取XML

最近学习Java&#xff0c;在处理XML文档的时候&#xff0c;查阅相关资料&#xff0c;发现了DOM4J这个jre库&#xff0c;相对C#的XML处理来说&#xff0c;功能还算是跟得上 下面展示一篇我自己写的一个XML读取测试 import java.util.Iterator; import java.io.BufferedReader; i…

c语言如何实现不定参数,C语言中不定参数的实现

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼但是请注意,由於 K&R C 中并不检查参数型态&#xff0c;so 在此是用 ANSIC 来说明&#xff0c;毕竟 ANSI C 是目前所有 C Compiler 都支援的标准&#xff1b;常在网路上看到有人推荐 K&R The C Programming Language 这本书…

SX-BOX试用笔记

1.今天开始试用SX-BOX&#xff0c;我将用它来做自己平时的开发和文档的工作&#xff0c;因为我的职责是一个程序员。 1.屏幕如果支持更大的也清晰的话&#xff0c;我会感觉很舒服。 2.USB设备插入还无法识别。。。 3.现在开始安装开发工具&#xff0c;目前我打算在上面安装Adob…

c语言成绩转换绩点,如何将平时成绩转化为GPA成绩?

新东方留学生资料下载姓名手机号短信验证码图片验证码年级请选择5年级6年级7年级8年级9年级高一高二高三高三复习大一大二大三大四大五研一研二研三博一博二在职课程请选择雅思托福GREGMATSAT其他行政区请选择请选择请选择{"name":{"label":"姓名&quo…

Hibernate查询缓存如何工作

介绍 既然我已经介绍了实体和集合缓存&#xff0c;现在该研究查询缓存的工作原理了。 查询缓存与实体严格相关&#xff0c;它在搜索条件和满足该特定查询过滤器的实体之间绘制关联。 像其他Hibernate功能一样&#xff0c;查询缓存也不像人们想象的那么琐碎。 实体模型 对于我…

Javascript学习笔记1 数论

1.Javascript不用担心内存的回收与对象的销毁&#xff01; 2.Javascript有&#xff1a;infinity、NaN全局变量表示 被0整除的无穷 和 非数字。undefined和null表示 未定义 和 空&#xff0c;undefined和null可以互换&#xff0c;判别二者需用 全等 号&#xff08;不光判断值&…

C语言bcd码减法过程,bcd码的减法运算规则举例.ppt

bcd码的减法运算规则举例第1章&#xff1a;微型计算机基础 本章基本要求&#xff1a; ⑴ 单片微型计算机的含义 ⑵ 各系列单片机的特点 ⑶ 51系列单片机的概念及指标(重点) ⑷ 单片微机工业产品概念 第1章&#xff1a;微型计算机基础 电子计算机是一种能对信息进行加工处理的机…

JavaFX技巧20:有很多需要展示的地方吗? 使用画布!

JavaFX应用程序似乎有两种&#xff1a;第一种使用带有节点和CSS样式的场景图&#xff0c;第二种使用单个画布。 但是&#xff0c;将这两种方法混合使用是完全合法的。 尤其是当您的应用程序必须显示大量详细信息时&#xff0c;您很容易最终创建成千上万个节点。 即使JavaFX的整…

VirtualBox命令更改虚拟硬盘空间

主要是使用VBoxManage命令来操作第一步&#xff1a;打开CMD&#xff0c;进入到virtualbox存放虚拟机的目录中(win7 系统可以直接在文件夹空白处按住shift键右键鼠标选择[在此处打开命令窗口])&#xff0c;输入[vboxmanage list hdds]可以查看所有的虚拟机的信息&#xff0c;这里…

华为手机老是android自动升级,华为手机系统怎么升级 华为手机升级系统的两种方法...

华为手机使用的均为安卓系统&#xff0c;升级方法有多种&#xff0c;各位可以根据条件不同自行选择升级。华为手机升级系统的两种方法方法一&#xff1a;手机联网更新首先把手机联网。点击设置。向下拖动点击最后一项“关于手机”点击“系统软件更新”。点击在线升级即可。友情…

Wildfly,Apache CXF和@SchemaValidation

在过去的几天中&#xff0c;我一直在进行从JBoss 4到Wildfly 8的应用程序迁移。 该应用程序使用了不同的技术&#xff0c;但是我们这里将重点放在XML Web Services JAX-WS上 。 是的&#xff0c;我知道它们已不再流行&#xff0c;但是这些是很久以前开发的&#xff0c;因此需要…