「LuoguP1430」 序列取数(区间dp
题目描述
给定一个长为n的整数序列(n<=1000),由A和B轮流取数(A先取)。每个人可从序列的左端或右端取若干个数(至少一个),但不能两端都取。所有数都被取走后,两人分别统计所取数的和作为各自的得分。假设A和B都足够聪明,都使自己得分尽量高,求A的最终得分。
输入输出格式
输入格式:
第一行,一个正整数T,表示有T组数据。(T<=100)
接着T行,每行第一个数为n,接着n个整数表示给定的序列.
输出格式:
输出T行,每行一个整数,表示A的得分
输入输出样例
说明
时限:3s
题解
首先让我们试试暴力思路:
设A的得分为VA,B的得分为VB
那么在(VA-VB)取得最大值时,有VA最大。
证明:VA-VB=VA-(Sum-VA)=2*VA-Sum
设F[L][R]为当前还剩[L][R]时,(先手得分-后手得分)的最大值。
若当前正在处理F[L][R],那么存在三种选择方案:
1.全取。
F[L][R]=∑{ v[i] | L < = i < = R }
2.从左边取一些。(保留L'到R)
F[L][R]=∑{ v[i] | L <= i <= L'-1 }-F[L'][R]
3.从右边取一些。(保留L到R')
F[L][R]=∑{ v[i] | R'+1 <= i <= R }-F[L][R']
关于减去F[L'][R]:
当区间转换到[L',R]时的F值,为转换后的先手-后手,也就是当前意义下的后手-先手,
所以-(后手-先手)=先手-后手。
这样,我们从小到大枚举区间,每次在这三种决策中取一种最优决策,
最后的max(VA-VB)=F[1][n],
max(VA)=(F[1][n]+Sum)/2。
于是我们可以很较为轻松的写出O(n^3)的暴力啦!
/*
qwerta
P1430 序列取数
Unaccepted
40
代码 C++,0.94KB
提交时间 2018-09-19 18:57:13
耗时/内存
19351ms, 4628KB
*/
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
#define R register
inline int read()
{
char ch=getchar();
int x=;bool s=;
while(!isdigit(ch)){if(ch=='-')s=;ch=getchar();}
while(isdigit(ch)){x=x*+ch-'';ch=getchar();}
return s?x:-x;
}
int s[];
int f[][];
//int m[1007][1007];
int main()
{
//freopen("a.in","r",stdin);
int t=read();
while(t--)
{
int n=read();
for(R int i=;i<=n;++i)
s[i]=s[i-]+read();
for(R int len=;len<=n;++len)
for(R int l=;l+len-<=n;++l)
{ int r=l+len-;
f[l][r]=s[r]-s[l-];
for(R int _l=l+;_l<=r;++_l)
f[l][r]=max(f[l][r],s[_l-]-s[l-]-f[_l][r]);
for(R int _r=r-;_r>=l;--_r)
f[l][r]=max(f[l][r],s[r]-s[_r]-f[l][_r]); }
printf("%d\n",(f[][n]+s[n])>>);
}
return ;
}
但是想要通过1e3的数据,O(n^3*t)的时间复杂度肯定是布星的。
所以我们还需要一点儿优化。
考虑状态数是n^2的,雷打不动,所以我们要将贼手伸向状态转移。
以从右取为例:
F[L][R]=max(S[R]-S[R']-F[L][R']) //(L <= R' <= R-1)
=max(S[R]-(S[R']+F[L][R']))
如果我们能知道R'在L到R-1上时,S[R']+F[L][R']的最小值,那么就能O(1)转移了。
设
Min[L][R]=min{S[R']+F[L][R'] | L <= R' <= R }
那么有
Min[L][R-1]=min{S[R']+F[L][R'] | L <= R' <= R-1 }
也就是说
Min[L][R]=min(Min[L][R-1],S[R]+F[L][R])
这样,在循环枚举区间的同时顺便维护一下Min[L][R],就可以实现O(1)转移了。
其实就是把暴力中间那两句话换了种说法而已。
/*
qwerta
P1430 序列取数
Accepted
100
代码 C++,1.63KB
提交时间 2018-09-19 20:05:27
耗时/内存
4910ms, 12444KB
*/
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
#define R register
inline int read()
{
char ch=getchar();
int x=;bool s=;
while(!isdigit(ch)){if(ch=='-')s=;ch=getchar();}
while(isdigit(ch)){x=x*+ch-'';ch=getchar();}
return s?x:-x;
}//快读(这题略卡常)
int s[];
int f[][];
int ml[][];
int mr[][];
//设ml为固定L端时的min值,mr为固定R端时的max值
void write(int x)
{
if(x>)write(x/);
putchar(x%+'');
return;
}//快写
int main()
{
//freopen("a.in","r",stdin);
int t=read();
while(t--)
{
int n=read();
for(R int i=;i<=n;++i)
s[i]=s[i-]+read();//前缀和
for(R int l=;l<=n;++l)
{
f[l][l]=s[l]-s[l-];
ml[l][l]=s[l]+f[l][l];
mr[l][l]=s[l-]-f[l][l];
}//预处理
for(R int len=;len<=n;++len)
for(R int l=,r=len;r<=n;++l,++r)//枚举区间
{
f[l][r]=s[r]-s[l-];
//mr
//for(R int _l=l+1;_l<=r;++_l)
//f[l][r]=max(f[l][r],s[_l-1]-s[l-1]-f[_l][r]);
f[l][r]=max(f[l][r],mr[l+][r]-s[l-]);
//ml
//for(R int _r=r-1;_r>=l;--_r)
//f[l][r]=max(f[l][r],s[r]-s[_r]-f[l][_r]);
f[l][r]=max(f[l][r],s[r]-ml[l][r-]);
//
ml[l][r]=min(ml[l][r-],s[r]+f[l][r]);
mr[l][r]=max(mr[l+][r],s[l-]-f[l][r]);
}
int x=(f[][n]+s[n])>>;
if(x<){putchar('-');write(-x);}
else write(x);
putchar('\n');
//输出答案
}
return ;
}
「LuoguP1430」 序列取数(区间dp的更多相关文章
- 「LuoguP1220」 关路灯(区间dp
题目描述 某一村庄在一条路线上安装了n盏路灯,每盏灯的功率有大有小(即同一段时间内消耗的电量有多有少).老张就住在这条路中间某一路灯旁,他有一项工作就是每天早上天亮时一盏一盏地关掉这些路灯. 为了给村 ...
- 「网络流24题」「LuoguP2774」方格取数问题(最大流 最小割
Description 在一个有 m*n 个方格的棋盘中,每个方格中有一个正整数.现要从方格中取数,使任意 2 个数所在方格没有公共边,且取出的数的总和最大.试设计一个满足要求的取数算法.对于给定的方 ...
- 【Luogu】P1430序列取数(DP)
题目链接 博弈DP太喵了qwq 设f[i][j]表示剩下区间[i,j]要取,先手最大值 明显我们要从这区间里面拿个最大的 就等价于这段区间的前缀和,我们要给对手留下个最小的 就是f[i][j]=sum ...
- LibreOJ #6007. 「网络流 24 题」方格取数 最小割 最大点权独立集 最大流
#6007. 「网络流 24 题」方格取数 内存限制:256 MiB时间限制:1000 ms标准输入输出 题目类型:传统评测方式:文本比较 上传者: 匿名 提交提交记录统计讨论测试数据 题目描述 ...
- Libre 6007 「网络流 24 题」方格取数 / Luogu 2774 方格取数问题 (网络流,最大流)
Libre 6007 「网络流 24 题」方格取数 / Luogu 2774 方格取数问题 (网络流,最大流) Description 在一个有 m*n 个方格的棋盘中,每个方格中有一个正整数.现要从 ...
- 洛谷 P1430 序列取数 解题报告
P1430 序列取数 题目描述 给定一个长为\(n\)的整数序列\((n<=1000)\),由\(A\)和\(B\)轮流取数(\(A\)先取).每个人可从序列的左端或右端取若干个数(至少一个), ...
- loj #2051. 「HNOI2016」序列
#2051. 「HNOI2016」序列 题目描述 给定长度为 n nn 的序列:a1,a2,⋯,an a_1, a_2, \cdots , a_na1,a2,⋯,an,记为 a[1: ...
- 「HNOI2016」序列 解题报告
「HNOI2016」序列 有一些高妙的做法,懒得看 考虑莫队,考虑莫队咋移动区间 然后你在区间内部找一个最小值的位置,假设现在从右边加 最小值左边区间显然可以\(O(1)\),最小值右边的区间是断掉的 ...
- 「SDOI2016」储能表(数位dp)
「SDOI2016」储能表(数位dp) 神仙数位 \(dp\) 系列 可能我做题做得少 \(QAQ\) \(f[i][0/1][0/1][0/1]\) 表示第 \(i\) 位 \(n\) 是否到达上界 ...
随机推荐
- Direct2D教程(三)简单几何图形
从本章开始,我们介绍D2D几何图形. D2D图形分类 Direct2D支持多种类型的几何图形,包括Simple Geometry(简单几何图形) 矩形 圆角矩形 椭圆 Path Geometry(路径 ...
- python 读取共享内存
测试环境 centos7 python3.6.5 首先使用c创建内存,这里的方法是:作为参数读一个二进制数据文件进去,把文件的内容作为共享内存的内容 定义块 #include <stdio.h& ...
- sql分级汇总
--測试数据 create table tb([DB-ID] varchar(10),ENTITY varchar(10),DATE varchar(10),[CUST-NO] int,AMOUNT ...
- MyEclipse的实体关系设计
原文地址:http://www.myeclipsecn.com/learningcenter/database-development/myeclipse-entity-relation-design ...
- java通过http方式下载文件
package com.qiyi; import java.io.BufferedReader; import java.io.IOException; import java.io.InputStr ...
- Spring源码学习(一)资源加载
这里先从最简单的一个Spring例子开始. 下面是Spring的context的配置 <?xml version="1.0" encoding="UTF-8&quo ...
- vuex 介绍
vuex是为vue.js开发的状态管理模式,负责vue的状态管理,状态管理是干啥的呢,举个栗子,比如一个酒店,哪间屋子入住了客人,哪间屋子客人退房了,客人退房后,房间有没有清扫过,这些都需要去记录,以 ...
- svn服务器迁移(windows下)
废话不多说,直接上步骤: 服务端: 1.创建一个备份文件夹 如:D:\svn_bak 2.进入cmd,cd命令到你的svn服务器安装目录的bin文件下,本人的安装目录在 D:\Program File ...
- Struts MVC工作原理(转载)
1.Struts MVC中Model 1 和Model 2简介 我们在开发Web应用时经常提到的一个概念是Model 1/Model 2,那么到底它是什么意思呢?其实它是对采用JSP技术构成Web应用 ...
- MySQL字段名与保留字冲突的问题及解决方法
问题:MySQL字段名与保留字冲突在实际操作是常常出现的.一把会出现下面错误. com.mysql.jdbc.exceptions.jdbc4.MySQLSyntaxErrorException. 解 ...