bzoj1563 [NOI2009]诗人小G(dp+决策单调性+二分)

首先有很显然的 O ( n 2 ) dp
f [ i ] = m i n { f [ j ] + | s [ i ] s [ j ] + i j 1 L | p }
30分。
然后考虑p=2的,化一下式子可以发现是个斜率优化。
20分。
考虑4,5两个L非常小的点,可以贪心,每行如果放了不止一个句子,长度还大于了2L,那么一定不优,因此可以限制每行的长度,复杂度 O ( N L )
20分。

通过算法二的斜率优化,我们猜想这个dp也是有决策单调性的。
(要不怎么做呢对吧qaq)
于是我们可以单调队列+二分维护最优决策区间即可。
复杂度 O ( n l o g n )

当然我们还是有大神证明了四边形不等式的:portal

#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define inf 1e18
#define N 100010
#define ld long double
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-') f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,L,p,s[N],path[N];
ld f[N];
char a[N][35];
struct Icefox{
    int l,r,id;
    Icefox(int _l=0,int _r=0,int _id=0){l=_l;r=_r;id=_id;}
}q[N];
inline ld ksm(ld x,int k){
    if(x<0) x=-x;ld res=1;for(;k;k>>=1,x=x*x) if(k&1) res=res*x;return res;
}
inline ld calc(int i,int j){
    return f[j]+ksm(s[i]-s[j]+i-j-1-L,p);
}
inline int ask(int l,int r,int k1,int k2){
    while(l<=r){
        int mid=l+r>>1;
        if(calc(mid,k1)<=calc(mid,k2)) l=mid+1;
        else r=mid-1;
    }return l-1;
}
inline void print(int x){
    if(path[x]) print(path[x]);
    for(int i=path[x]+1;i<=x;++i){
        for(int j=1;a[i][j];++j) putchar(a[i][j]);
        putchar(i==x?'\n':' ');
    }
}
int main(){
//  freopen("a.in","r",stdin);
    int tst=read();
    while(tst--){
        n=read();L=read();p=read();
        for(int i=1;i<=n;++i) scanf("%s",a[i]+1),s[i]=s[i-1]+strlen(a[i]+1);
        int qh=1,qt=0;q[++qt]=Icefox(1,n,0);
        for(int i=1;i<=n;++i){
            if(qh<=qt&&q[qh].r<i) ++qh;
            if(qh<=qt) ++q[qh].l;
            f[i]=calc(i,q[qh].id);path[i]=q[qh].id;
            if(qh>qt||calc(n,i)<calc(n,q[qt].id)){
                while(qh<=qt&&calc(q[qt].l,i)<=calc(q[qt].l,q[qt].id)) --qt;
                if(qh<=qt){
                    int t=ask(q[qt].l,q[qt].r,q[qt].id,i);
                    q[qt].r=t;q[++qt]=Icefox(t+1,n,i);
                }else q[++qt]=Icefox(i+1,n,i);
            }
        }
        if(f[n]>inf) puts("Too hard to arrange");
        else printf("%lld\n",(ll)f[n]);//print(n);
        puts("--------------------");
    }return 0;
}

猜你喜欢

转载自blog.csdn.net/icefox_zhx/article/details/80904671