bzoj1563(单调决策性DP)

这个题一看数字就很爆炸啊。。不过先不管他。。

先对长度+1，然后取个前缀和设为a

设d[i]为取到i个句子且当前行以i句结尾的最小代价

d[i]=max{d[j]+|a[i]-a[j]-L-1|^p}

显然f(x)=|x|^p是个下凸函数，即x越大导数越大。。

然后考虑2个决策点j<k<i

如果d[j]+|a[i]-a[j]-L-1|^p>d[k]+|a[i]-a[k]-L-1|^p，那么此后j就不会比k更优了。。

所以用单调队列存递增的决策点，用上面的式子去队尾就可以了。。

然而数非常爆炸，在比较的时候很麻烦，于是用long double(double都不行= =)来做大数比较，虽然有精度误差但是问题不大。。

复杂度O(Tnlognlogp)，常数巨大。。

 

 

/**
 *　　　　　　　　┏┓　　 　┏┓
 * 　　　　　　　┏┛┗━━━━━━━┛┗━━━┓
 * 　　　　　　　┃　　　　　　　┃ 　
 * 　　　　　　　┃　　　━　　 　┃
 * 　　　　　　　┃　＞　　　＜　┃
 * 　　　　　　　┃　　　　　　　┃
 * 　　　　　　　┃...　⌒　... 　┃
 * 　　　　　　　┃　　　　　　　┃
 * 　　　　　　　┗━┓　　　┏━┛
 * 　　　　　　　　　┃　　　┃　Code is far away from bug with the animal protecting　　　　　　　　　　
 * 　　　　　　　　　┃　　　┃   神兽保佑,代码无bug
 * 　　　　　　　　　┃　　　┃　　　　　　　　　　　
 * 　　　　　　　　　┃　　　┃  　　　　　　
 * 　　　　　　　　　┃　　　┃
 * 　　　　　　　　　┃　　　┃　　　　　　　　　　　
 * 　　　　　　　　　┃　　　┗━━━┓
 * 　　　　　　　　　┃　　　　　　　┣┓
 * 　　　　　　　　　┃　　　　　　　┏┛
 * 　　　　　　　　　┗┓┓┏━┳┓┏┛
 * 　　　　　　　　　　┃┫┫　┃┫┫
 * 　　　　　　　　　　┗┻┛　┗┻┛
 */
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>
#include<cmath>
#include<map>
#include<stack>
#include<set>
#include<bitset>
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,l,r) for(int i=l;i>=r;i--)
#define link(x) for(edge *j=h[x];j;j=j->next)
#define mem(a) memset(a,0,sizeof(a))
#define ll long long
#define eps 1e-8
#define succ(x) (1LL<<(x))
#define lowbit(x) (x&(-x))
#define sqr(x) ((x)*(x))
#define mid (x+y>>1)
#define NM 100005
#define nm 200005
#define pi 3.1415926535897931
using namespace std;
const ll inf=1e18;
ll read(){
    ll x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
    return f*x;
}








int n,L,p,qh,qt,q[NM],f[NM],v[NM],a[NM];
char _s[50];
long double d[NM];
long double qpow(long double x,int t){return t?qpow(sqr(x),t>>1)*(t&1?x:1):1;}

int main(){
    freopen("data.in","r",stdin);
    int _=read();while(_--){
	mem(d);mem(f);mem(v);
	n=read();L=1+read();p=read();
	inc(i,1,n){scanf("%s",_s);a[i]=a[i-1]+1+strlen(_s);}
	q[qh=qt=1]=0;
	inc(i,1,n){
	    f[i]=max(f[i],f[i-1]);
	    while(qh<=qt&&f[i]!=q[qh])qh++;
	    d[i]=d[q[qh]]+qpow(abs(a[i]-a[q[qh]]-L),p);
	    int s=i+1;
	    while(qh<=qt){
		s=n+1;
#define fun(t) d[t]+qpow(abs(a[mid]-a[t]-L),p)
		for(int x=i+1,y=n;x<=y;)
		    if(fun(q[qt])>fun(i))s=mid,y=mid-1;else x=mid+1;
		if(s<=v[q[qt]])qt--;else break;
	    }
	    if(s<=n)f[s]=i,v[i]=s,q[++qt]=i;
	}
	//inc(i,1,n)printf("%d ",a[i]-a[i-1]);putchar('\n');
	//inc(i,1,n)printf("%d ",f[i]);putchar('\n');
	//inc(i,1,n)printf("%lld ",d[i]);putchar('\n');
	if(d[n]>inf)printf("Too hard to arrange\n");else printf("%lld\n",(ll)d[n]);
	printf("--------------------\n");
    }
    return 0;
}

 

 

 

1563: [NOI2009]诗人小G

Time Limit: 100 Sec  Memory Limit: 64 MB
Submit: 3017  Solved: 1017
[Submit][Status][Discuss]

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。
 

Source

 

[Submit][Status][Discuss]