博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
HDU 6619 Horse 斜率优化dp
阅读量:5291 次
发布时间:2019-06-14

本文共 1724 字,大约阅读时间需要 5 分钟。

#include
#define fi first#define se second#define INF 0x3f3f3f3f#define LNF 0x3f3f3f3f3f3f3f3f#define fio ios::sync_with_stdio(false);cin.tie(0);cout.tie(0)#define pqueue priority_queue#define NEW(a,b) memset(a,b,sizeof(a))const double pi=4.0*atan(1.0);const double e=exp(1.0);const int maxn=1e4+8;typedef long long LL;typedef unsigned long long ULL;const LL mod=1e9+7;const ULL base=1e7+7;const int maxp=26+5;using namespace std;LL a[maxn],f[maxn];LL g[maxn],gg[maxn];LL dp[maxn][58];LL q[58][maxn];LL l[58],r[58];LL gdp(int i,int j,int k) {
return dp[j][k]-gg[j]-g[j]*(i-j)+gg[i];}LL gfz(int j,int h,int k) {
return dp[j][k]-gg[j]+g[j]*(j)-(dp[h][k]-gg[h]+g[h]*(h));}LL gfm(int j,int h) {
return g[j]-g[h];}int main(){ int t; scanf("%d",&t); int n,m,k; while(t--){ scanf("%d%d%d",&n,&k,&m); for(int i=1;i<=n;i++){ scanf("%lld",&a[i]); g[i]=g[i-1]+a[i]; gg[i]=g[i]+gg[i-1]; } for(int i=1;i<=k;i++){ l[i]=r[i]=1; q[i][l[i]]=0; } for(int i=1;i<=n;i++){ for(int j=1;j<=i&&j<=k+1;j++){ while(l[j-1]
<=i*gfm(q[j-1][l[j-1]+1],q[j-1][l[j-1]])) l[j-1]++; dp[i][j]=gdp(i,q[j-1][l[j-1]],j-1); while(l[j]
<=gfm(i,q[j][r[j]])*gfz(q[j][r[j]],q[j][r[j]-1],j)) r[j]--; q[j][++r[j]]=i; } } for(int i=1;i<=n;i++){ f[i]=a[i]*(n-i+1); } sort(f+1,f+1+n); LL ans=0; for(int i=n;i>=n-m+1;i--){ ans+=f[i]; } ans-=dp[n][k+1]; printf("%lld\n",ans); }}

 

转载于:https://www.cnblogs.com/Profish/p/11297891.html

你可能感兴趣的文章
最长公共子串问题(LCS)
查看>>
TortoiseSVN is locked in another working copy
查看>>
PHP魔术方法之__call与__callStatic方法
查看>>
ubuntu 安装后的配置
查看>>
Html学习_简易个人网页制作
查看>>
angular中ng-bind指令小案例
查看>>
jqery总结
查看>>
Lodop获取客户端主网卡ip地址是0.0.0.0
查看>>
VSCODE更改文件时,提示:EACCES: permission denied的解决办法(mac电脑系统)
查看>>
web前端之路,js的一些好书(摘自聂微东 )
查看>>
【模板】对拍程序
查看>>
微信小程序开发初体验
查看>>
dos批处理(bat)运行exe
查看>>
关键字
查看>>
Pycharm安装Markdown插件
查看>>
上传图片并预览
查看>>
哈夫曼编码_静态库
查看>>
【转】redo与undo
查看>>
C#更新程序设计
查看>>
常用Request对象获取请求信息
查看>>