解析

本题的决策单调性可以说是显然
但是本题是同维度（其实只有一维）自左向右转移，分治的写法是不能奏效的
所以我们使用决策点调性的另一种实现方法：二分栈
具体来说，维护[1,n]各自的最优转移点
一开始的转移点全是0，把三元组(0,1,n)(0,1,n)(0,1,n)push入栈，分别表示转移点和左右端点
然后从1扫到n
扫到i的时候，利用i当前的最优决策点进行转移
然后尝试更新后面的最优决策点
不断取出栈顶，如果栈顶区间在 l 的位置还是不如当前的转移，显然这个区间没有用了，弹出栈顶
然后对于无法使之弹出的区间，二分出其与自己在何时会发生最优转移的变化，如果这个变化的位置pl<=n，把(i,pl,n)(i,pl,n)(i,pl,n)入栈，并把之前的栈顶区间右端点改为pl-1
然后就解决啦！

注意事项

至少是我易错的

二分的时候当calc相等的时候不能草率！可能是都返回了inf，所以要按照这个时候那个绝对值的正负分类讨论

while(st<ed){int mid=(st+ed)>>1;ll a=calc(j,mid),b=calc(i,mid);if(a<b||(a==b&&sum[mid]-sum[i]-l-1>0)) ed=mid;else st=mid+1;
}

代码

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define il inline
#define debug(a) fprintf(stderr,a)
const int N=1e5+100;
const int M=3e6+100;
const int mod=998244353;
inline ll read(){ll x=0,f=1;char c=getchar();while(!isdigit(c)){if(c=='-') f=-1;c=getchar();}while(isdigit(c)){x=x*10+c-'0';c=getchar();}return x*f;
}
int n,p,l;
ll f[N],sum[N];
ll ksm(ll x,ll k){if(x<0) x=-x;ll res=1;while(k){if(k&1) res*=x;x=x*x;k>>=1;}return res;
}
inline ll calc(int i,int j){//printf("  calc:%d %d jd=%lld>%lf\n",i,j,abs(sum[j]-sum[i]-l-1),(1.0*log(1e18-f[i])/log(p)));if(abs(sum[j]-sum[i]-l-1)&&(f[i]>1e18||p>1.0*log(2e18-f[i])/log(abs(sum[j]-sum[i]-l-1)))) return 2e18+1;else return f[i]+ksm(sum[j]-sum[i]-l-1,p);
}
#define pr pair<int,int>
#define mkp make_pair
pr q[N];//id st
int st,ed;
char s[N][33];
int pre[N];
int find(int i,int j){int st=j+1,ed=n+1;while(st<ed){int mid=(st+ed)>>1;ll a=calc(j,mid),b=calc(i,mid);if(a<b||(a==b&&sum[mid]-sum[i]-l-1>0)) ed=mid;else st=mid+1;}//printf("find:i=%d j=%d st=%d\n",i,j,st);return st;
}
int a[N],b[N],tot;
int main(){#ifndef ONLINE_JUDGEfreopen("a.in","r",stdin);freopen("a.out","w",stdout);#endifint T=read();while(T--){n=read();l=read();p=read();for(int i=1;i<=n;i++){scanf(" %s",s[i]+1);sum[i]=sum[i-1]+strlen(s[i]+1)+1;//printf("i=%d sum=%d\n",i,sum[i]);}st=ed=1;q[1]=mkp(0,1);int flag=0;for(int i=1;i<=n;i++){while(st<ed&&q[st+1].second<=i) st++;int now=q[st].first;pre[i]=now;f[i]=calc(now,i);//printf("i=%d now=%d f=%lld\n",i,now,f[i]);while(st<=ed){int pl=find(q[ed].first,i);if(pl<=q[ed].second) ed--;else{if(pl<=n) q[++ed]=mkp(i,pl);break;}}if(st>ed) q[++ed]=mkp(i,i+1);}//debug("ok");if(f[n]>1e18) printf("Too hard to arrange\n");else{printf("%lld\n",f[n]);tot=0;int pl=n;while(pl){b[++tot]=pl;pl=pre[pl];a[tot]=pl+1;}for(int i=tot;i>=1;i--){for(int j=a[i];j<=b[i];j++){printf("%s",s[j]+1);if(j!=b[i]) putchar(' ');}printf("\n");}}printf("--------------------\n");}return 0;
}
/*
3 1
3 1 33 2
1 1 2
3 1 3
*/

洛谷P1912：诗人小G（二分栈、决策单调性）相关推荐

P1912 [NOI2009]诗人小G
P1912 [NOI2009]诗人小G 思路: 平行四边形不等式优化dp 因为f(j, i) = abs(sum[i]-sum[j]+i-j-1-l)^p 满足平行四边形不等式 j < i f( ...
决策单调性Ⅰ：四边形不等式（bzoj 1563: [NOI2009]诗人小G）
题目描述: 给出n个数字和常数L,你可以任意合并相邻的两个数字a[x]和a[x+1],并得出一个新的数a[x]+a[x+1]+1,一通合并后得到一个有若干个数的序列,这个序列的不协调值为∑|(a[i] ...
【NOI 2009】诗人小G
[NOI 2009]诗人小G Problem Description 小 \(G\) 是一个出色的诗人,经常作诗自娱自乐.但是,他一直被一件事情所困扰,那就是诗的排版问题. 一首诗包含了若干个句子,对 ...
洛谷P1462 通往奥格瑞玛的道路二分答案+最短路SPFA
洛谷P1462 通往奥格瑞玛的道路二分答案+最短路SPFA 二分交费最多的一次的钱数然后只将符合要求的边加入图中如果到终点的最短路大于等于血量或者直接起点不能到达终点那么说明不符合要求需要 ...
【洛谷4005】小Y和地铁（搜索）
[洛谷4005]小Y和地铁(搜索) 题面洛谷有点长. 题解首先对于需要被链接的两个点,样例中间基本上把所有的情况都给出来了. 但是还缺了一种从下面绕道左边在从整个上面跨过去在从右边绕到下面来的情 ...
NOI2009 诗人小G
题目小G是一个出色的诗人,经常作诗自娱自乐.但是,他一直被一件事情所困扰,那就是诗的排版问题. 一首诗包含了若干个句子,对于一些连续的短句,可以将它们用空格隔开并放在一行中,注意一行中可以放的句子数 ...
洛谷 1373 dp 小a和uim之大逃离良心题解
洛谷 1373 dp 这题还不算太难,,当初看的时候不是很理解题意,以为他们会选择两条不同的路径,导致整体思路混乱传送门其实理解题意和思路之后还是敲了不短的时间,一部分身体原因再加上中午休息不太好 ...
一个诗人的一生——诗人小G的人生
看到群里有人说周末作文是以一个诗人的一生为题写作文-于是就突发脑洞写了这么个东西. 哇发出来好羞耻啊(x 诗与人生某一天,在某一个地方,小G诞生了. 小G原本是有自己的名字的,但是在某一天之后,他就 ...
【洛谷1361】小M的作物（最小割）
传送门洛谷 Solution 这是一个比较实用的套路,很多题目都有用,而且这个套路难以口胡出来. 考虑把每一个附加贡献重新建一个点,然后向必需的点连边,流量为val. 然后直接种植的从源点向这个点连 ...

洛谷P1912：诗人小G（二分栈、决策单调性）

解析

注意事项

代码

洛谷P1912：诗人小G（二分栈、决策单调性）相关推荐

最新文章

热门文章