【BZOJ1563】【NOI2009】诗人小G(dp+决策单调性)
2017-03-15 21:37
411 查看
Description
Input
Output
对于每组数据,若最小的不协调度不超过1018,则第一行一个数表示不协调度若最小的不协调度超过1018,则输出”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
总共10个测试点,数据范围满足:
测试点 T N L P
1 ≤10 ≤18 ≤100 ≤5
2 ≤10 ≤2000 ≤60000 ≤10
3 ≤10 ≤2000 ≤60000 ≤10
4 ≤5 ≤100000 ≤200 ≤10
5 ≤5 ≤100000 ≤200 ≤10
6 ≤5 ≤100000 ≤3000000 2
7 ≤5 ≤100000 ≤3000000 2
8 ≤5 ≤100000 ≤3000000 ≤10
9 ≤5 ≤100000 ≤3000000 ≤10
10 ≤5 ≤100000 ≤3000000 ≤10
所有测试点中均满足句子长度不超过30。
题解:%%%PoPoQQQ大神!!!%%%byvoid大神!!!
首先,可以列出dp方程,这样可得30分(因为有个讨厌的指数)。
列出了dp方程但是只能拿暴力分数的时候,就要想怎么优化了(除非dp根本是错的)。决策单调性我不会证(我是连导数都不会的蒟蒻),但是打表还是可以看出来的。
由于列出的方程是这样的:
我们可以看出这是一个1D1D动态规划,那么决策区间就是连续的段落,于是我们维护一个上凸壳,每次更新的时候用二分就好了。
注意数据比较大,用long double 算完转long long。
代码如下:
Input
Output
对于每组数据,若最小的不协调度不超过1018,则第一行一个数表示不协调度若最小的不协调度超过1018,则输出”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
总共10个测试点,数据范围满足:
测试点 T N L P
1 ≤10 ≤18 ≤100 ≤5
2 ≤10 ≤2000 ≤60000 ≤10
3 ≤10 ≤2000 ≤60000 ≤10
4 ≤5 ≤100000 ≤200 ≤10
5 ≤5 ≤100000 ≤200 ≤10
6 ≤5 ≤100000 ≤3000000 2
7 ≤5 ≤100000 ≤3000000 2
8 ≤5 ≤100000 ≤3000000 ≤10
9 ≤5 ≤100000 ≤3000000 ≤10
10 ≤5 ≤100000 ≤3000000 ≤10
所有测试点中均满足句子长度不超过30。
题解:%%%PoPoQQQ大神!!!%%%byvoid大神!!!
首先,可以列出dp方程,这样可得30分(因为有个讨厌的指数)。
列出了dp方程但是只能拿暴力分数的时候,就要想怎么优化了(除非dp根本是错的)。决策单调性我不会证(我是连导数都不会的蒟蒻),但是打表还是可以看出来的。
由于列出的方程是这样的:
我们可以看出这是一个1D1D动态规划,那么决策区间就是连续的段落,于是我们维护一个上凸壳,每次更新的时候用二分就好了。
注意数据比较大,用long double 算完转long long。
代码如下:
#include<iostream> #include<stdio.h> #include<algorithm> #include<string.h> #include<math.h> #define ll long double #define inf 9000000000000000000 #define MAX 1000000000000000000LL using namespace std; int t,n,l,p,top; ll sum[100005],f[100005],from[100005]; char ch[100005][35]; struct nod { int l,r,p; nod(){} nod(int a,int b,int c):l(a),r(b),p(c){} }q[100005]; ll read() { ll x=0; char c=getchar(); while(c<'0' || c>'9') c=getchar(); while(c<='9' && c>='0'){x=x*10+c-'0';c=getchar();} return x; } ll pow(ll x) { if(x<0) x=-x; ll ans=1; for(int i=1;i<=p;i++) ans*=x; return ans; } ll cal(int j,int i) { return f[j]+pow(sum[i]-sum[j]+(i-j-1)-l); } int find(nod a,int b) { int l=a.l,r=a.r; while(l<=r) { int mid=l+r>>1; if(cal(a.p,mid)<cal(b,mid)) l=mid+1; else r=mid-1; } return l; } void dp() { int hd=1,tl=0; q[++tl]=nod(0,n,0); for(int i=1;i<=n;i++) { if(hd<=tl && i>q[hd].r) hd++; f[i]=cal(q[hd].p,i);from[i]=q[hd].p; if(hd>tl || cal(i,n)<=cal(q[tl].p,n)) { while(hd<=tl && cal(i,q[tl].l)<=cal(q[tl].p,q[tl].l)) tl--; if(hd>tl) q[++tl]=nod(i,n,i); else { int t=find(q[tl],i); q[tl].r=t-1; q[++tl]=nod(t,n,i); } } } } int main() { t=read(); while(t--) { n=read(),l=read(),p=read(); for(int i=1;i<=n;i++) scanf("%s",ch[i]); for(int i=1;i<=n;i++) sum[i]=sum[i-1]+strlen(ch[i]); dp(); if(f >MAX) printf("Too hard to arrange\n"); else printf("%lld\n",(long long)(f )); printf("--------------------\n"); } return 0; }
相关文章推荐
- [BZOJ1563][NOI2009]诗人小G(dp+决策单调性)
- [BZOJ1563][NOI2009]诗人小G(决策单调性优化DP)
- BZOJ 1563: [NOI2009]诗人小G 决策单调性DP
- bzoj 1563 [NOI2009]诗人小G 四边形不等式 决策单调dp
- BZOJ_1563_[NOI2009]诗人小G_决策单调性
- bzoj1563: [NOI2009]诗人小G【决策单调性优化dp】
- [bzoj1563][NOI2009]诗人小G(决策单调性优化)
- [NOI2009]诗人小G(决策单调性优化dp)
- 决策单调性Ⅰ:四边形不等式(bzoj 1563: [NOI2009]诗人小G)
- ★【动态规划】【决策单调性优化】【NOI2009】诗人小G
- 【BZOJ 1563】 (四边形优化、决策单调性)
- [BZOJ 1563][NOI 2009]诗人小G(四边形优化DP)
- bzoj2369 && 2687 -- 决策单调性优化DP
- bzoj 2216: [Poi2011]Lightning Conductor(DP决策单调性)
- 【BZOJ1563】【NOI2009】诗人小G
- Bzoj:[Poi2011]Lightning Conductor:决策单调性优化DP详解
- bzoj4709 -- 决策单调性优化DP
- bzoj1563 [NOI2009]诗人小G
- BZOJ.1010.[HNOI2008]玩具装箱toy(DP 斜率优化/单调队列 决策单调性)
- [决策单调性 分治||单调栈 DP] BZOJ 2739 最远点