【BZOJ 1563】 （四边形优化、决策单调性）

1563: [NOI2009]诗人小G

Description

1 #include<cstdio>
2 #include<cstdlib>
3 #include<cstring>
4 #include<iostream>
5 #include<algorithm>
6 using namespace std;
7 #define LL long long
8 #define LD long double
9 #define Maxn 100010
10 const long double INF=1e18;
11
12 char s[50];
13 LD a[Maxn],f[Maxn],len[Maxn],sm[Maxn],L;
14 int lt[Maxn],nt[Maxn],st[Maxn],P;
15
16 LD qpow(LD x,int b)
17 {
18     if(x<0) x=-x;
19     LD ans=1.0;
20     while(b)
21     {
22         if(b&1) ans*=x;
23         x*=x;
24         b>>=1;
25     }
26     return ans;
27 }
28
29 LD cal(int i,int j)
30 {
31     return f[j]+qpow(sm[i]-sm[j]+(LD)i-(LD)j-1.0-L,P);
32 }
33
34 bool check(int mid,int x,int y)
35 {
36     return cal(mid,x)>=cal(mid,y);
37 }
38
39 int n;
40 int ffind(int l,int r,int x,int y)
41 {
42     int ans=n+1;
43     while(l<=r)
44     {
45         int mid=(l+r)>>1;
46         if(check(mid,x,y)) ans=mid,r=mid-1;
47         else l=mid+1;
48     }
49     return ans;
50 }
51
52 int main()
53 {
54     int T;
55     scanf("%d",&T);
56     while(T--)
57     {
58         cin>>n>>L>>P;sm[0]=0;
59         for(int i=1;i<=n;i++)
60         {
61             scanf("%s",s);len[i]=(LD)strlen(s);
62             sm[i]=sm[i-1]+len[i];
63         }
64         for(int i=0;i<=n;i++) f[i]=INF+1,st[i]=n+1,lt[i]=i-1,nt[i]=i+1;
65         f[0]=0;st[0]=1;
66         int now=0;
67         for(int i=1;i<=n;i++)
68         {
69             while(st[nt[now]]<=i) now=nt[now];
70             f[i]=cal(i,now);
71             while(st[lt[i]]>i)
72             {
73                 if(check(st[lt[i]],lt[i],i))
74                 {
75                     lt[i]=lt[lt[i]];
76                     nt[lt[i]]=i;
77                 }
78                 else break;
79             }
80             st[i]=ffind(st[lt[i]],n,lt[i],i);
81             if(st[i]>n) nt[lt[i]]=nt[i],lt[nt[i]]=lt[i];
82         }
83         if(f
>INF) printf("Too hard to arrange\n");
84         else cout<<(LL)f
<<endl;
85         printf("--------------------\n");
86     }
87     return 0;
88 }