BZOJ 1563 诗人小G
Description
Input
Output
对于每组数据,若最小的不协调度不超过\(10^{18}\),则第一行一个数表示不协调度若最小的不协调度超过\(10^{18}\),则输出"\(Too\;hard\;to\;arrange\)"(不包含引号)。每个输出后面加"\(--------------------\)"。
Sample Input
4
4 9 3
brysj,
hhrhl.
yqqlm,
gsycl.
4 9 2
brysj,
hhrhl.
yqqlm,
gsycl.
1 1005 6
poet
1 1004 6
poet
Sample Output
108
\(--------------------\)
32
\(--------------------\)
Too hard to arrange
\(--------------------\)
1000000000000000000
\(--------------------\)
【样例说明】
前两组输入数据中每行的实际长度均为\(6\),后两组输入数据每行的实际长度均为\(4\)。一个排版方案中每行相邻两个句子之间的空格也算在这行的长度中(可参见样例中第二组数据)。每行末尾没有空格。
Hint
1D1D动态规划裸题。
dp方程:$$f_{i}=min(f_{j}+(pre_{i}-pre_{j}+1)^{p})$$
方便起见,我们在\(pre_{i}\)上加个\(1\),于是dp方程变为$$f_{i}=min(f_{j}+(pre_{i}-pre_{j})^{p})$$
这个方程很明显地满足单调性(令对于\(f_{i}\)的转移\(k\)优于\(j\)(\(k > j\)),则对于\(v>i\)的\(f_{v}\)的转移也有\(k\)优于\(j\))。
斜率优化优化很明显当\(p \ne 2\)是行不通的。所以我们转向\(O(nlogn)\)的1D1D的动态规划。
1D1D动态规划是用二分单调栈来实现的。原理便是决策单调性,对于每个已经确定的\(f_{i}\),看其能更新的那一段后缀为那一段。在单调栈中进行二分,与之前的决策进行比较。不懂可以参考一下代码,代码应该好懂:
#include<cstring>
#include<cstdio>
#include<cstdlib>
using namespace std;
#define limit (1e18)
#define maxn 100010
#define maxm 40
int N,L,P,pre[maxn],top;
char s[maxm];
long double f[maxn];
struct node { int l,r,key; }stack[maxn];
inline long double qsm(int a,int b)
{
long double ret = 1;
while (b--) ret *= 1.0*a;
return ret;
}
inline long double calc(int a,int b)
{
return f[b]+qsm(abs(pre[a]-pre[b]-L),P);
}
inline int find(int a)
{
int l = 1,r = top,mid;
while (l <= r)
{
mid = (l + r) >> 1;
if (stack[mid].l<=a&&stack[mid].r>=a) return stack[mid].key;
if (a < stack[mid].l) r = mid - 1;
else l = mid + 1;
}
}
inline void updata(int now)
{
int l = 1,r;
while (top)
{
if (calc(stack[top].l,stack[top].key) >= calc(stack[top].l,now))
--top;
else
{
l = stack[top].l,r = stack[top].r;
while (l <= r)
{
int mid = (l + r) >> 1;
if (calc(mid,stack[top].key) >= calc(mid,now)) r = mid - 1;
else l = mid + 1;
}
stack[top].r = r;
break;
}
}
if (l <= N) stack[++top] = (node){l,N,now};
}
inline void dp()
{
f[0] = 0;
stack[top = 1] = (node) {1,N,0};
for (int i = 1;i <= N;++i)
{
int key = find(i);
f[i] = calc(i,key);
updata(i);
}
}
int main()
{
freopen("1563.in","r",stdin);
freopen("1563.out","w",stdout);
int T; scanf("%d",&T);
while (T--)
{
scanf("%d %d %d\n",&N,&L,&P);
L++;
for (int i = 1;i <= N;++i)
{
scanf("%s",s);
pre[i] = strlen(s)+1+pre[i-1];
}
dp();
if (f[N] > limit) printf("Too hard to arrange\n");
else printf("%.0Lf\n",f[N]);
printf("--------------------\n");
}
fclose(stdin); fclose(stdout);
return 0;
}
BZOJ 1563 诗人小G的更多相关文章
- [BZOJ 1563] [NOI 2009] 诗人小G(决策单调性)
[BZOJ 1563] [NOI 2009] 诗人小G(决策单调性) 题面 一首诗包含了若干个句子,对于一些连续的短句,可以将它们用空格隔开并放在一行中,注意一行中可以放的句子数目是没有限制的.小 G ...
- 1563: [NOI2009]诗人小G
1563: [NOI2009]诗人小G https://lydsy.com/JudgeOnline/problem.php?id=1563 分析: 直接转移f[i]=f[j]+cost(i,j),co ...
- C++之路进阶——codevs2933(诗人小G)
2933 诗人小G 2009年NOI全国竞赛 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 大师 Master 题目描述 Description 小G是一个出色的诗人 ...
- LG1912 [NOI2009]诗人小G
题意 题目描述 小G是一个出色的诗人,经常作诗自娱自乐.但是,他一直被一件事情所困扰,那就是诗的排版问题. 一首诗包含了若干个句子,对于一些连续的短句,可以将它们用空格隔开并放在一行中,注意一行中可以 ...
- bzoj1563: [NOI2009]诗人小G 决策单调性(1D1D)
目录 题目链接 题解 代码 题目链接 bzoj1563: [NOI2009]诗人小G 题解 \(n^2\) 的dp长这样 \(f_i = min(f_j + (sum_i - sum_j - 1 - ...
- 【Luogu1912】【NOI2009】诗人小G(动态规划)
[Luogu1912][NOI2009]诗人小G(动态规划) 题面 洛谷 题解 原来\(NOI\)这么多神仙题... 考虑一个极其明显的\(dp\) 设\(f[i]\)表示前\(i\)个句子产生的最小 ...
- BZOJ1563/洛谷P1912 诗人小G 【四边形不等式优化dp】
题目链接 洛谷P1912[原题,需输出方案] BZOJ1563[无SPJ,只需输出结果] 题解 四边形不等式 什么是四边形不等式? 一个定义域在整数上的函数\(val(i,j)\),满足对\(\for ...
- [NOI2009]诗人小G --- DP + 决策单调性
[NOI2009]诗人小G 题目描述: 小G是一个出色的诗人,经常作诗自娱自乐. 但是,他一直被一件事情所困扰,那就是诗的排版问题. 一首诗包含了若干个句子,对于一些连续的短句,可以将它们用空格隔开并 ...
- NOI 2009A 诗人小G
NOI 2009A 诗人小G 诗人小G [问题描述] 小G是一个出色的诗人,经常作诗自娱自乐.但是,他一直被一件事情所困扰,那就是诗的排版问题. 一首诗包含了若干个句子,对于一些连续的短句,可以将它们 ...
随机推荐
- 关于java里小数点的保留
关于java里小数点的保留 1.先给大家看一个代码. import java.util.*;import java.text.*; public class A { public static voi ...
- header的用法小结(转)
php header()函数的具体作用是向客户端发送一个原始 HTTP 标头[Http Header]到客户端. 标头 (header) 是服务器以 HTTP 协义传 HTML 资料到浏览器前所送出的 ...
- spring mvc DispatcherServlet详解之一--request通过HandlerMaping获取控制器Controller过程
整个spring mvc的架构如下图所示: 现在来讲解DispatcherServletDispatcherServlet的第一步:获取控制器. HandlerMapping HandlerMappi ...
- Customizing the Test Runner
There are several situations where you want to customize Robolectric's test runner to perform some o ...
- int? 参数是这个的时候 是可以传入null的 而int的就不行
such as pager.CurrentPageIndex = (page != null ? (int)page : 1);
- C# 里窗体里(windows form)怎么播放音乐
在.NET的winform里面,没有托管的音乐播放器,API只能播放WAV格式,对于MP3等形式的音频文件,就要依赖于 MediaPlayer里,嘿嘿 使用的方法: 在toolbox上点右键,选择“选 ...
- HTML5 WebAudioAPI-实例(二)
简单播放实例1: var url='../content/audio/海阔天空.mp3'; if (!window.AudioContext) { alert('您的浏览器不支持AudioContex ...
- jdbc - Insert 'Date' value in PreparedStatement
“preparedStatement.setDate()”方法接受的是 'java.sql.Date' 类型的参数,而我们一般格式化日期所使用的是'java.util.Date'中的'SimpleDa ...
- Tensor神经网络进行知识库推理
本文是我关于论文<Reasoning With Neural Tensor Networks for Knowledge Base Completion>的学习笔记. 一.算法简介 网络的 ...
- Html禁止粘贴 复制 剪切
oncopy="return false;" onpaste="return false;" oncut="return false;"