有
n−1
个限制,为
xi∗xi+1≥ai
,这里要求
xi
都是正实数。
要求
∑xi
n≤2000
解肯定是分成若干不短于2的段 每段内部取等号且和最小
那么我们用
fi,j
表示
[i,j]
是最后一段取等号的
如果
[i,j]
和
[i+1,k]
取等号时的
xi
和
xi+1
满足
xi∗xi+1>ai
就可以转移到
fi+1,k
朴素做是
O(n3)
的
在
i
处按照
O(n2logn)
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
inline char nc(){
static char buf[100000],*p1=buf,*p2=buf;
return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline void read(int &x){
char c=nc(),b=1;
for (;!(c>='0' && c<='9');c=nc()) if (c=='0') b=-1;
for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b;
}
const int N=2005;
int n,a[N];
double lx[N][N],rx[N][N],val[N][N];
double f[N][N];
typedef pair<double,double> abcd;
abcd tmp[N]; int tot;
double g[N];
int main(){
int T,Case=0;
freopen("t.in","r",stdin);
freopen("t.out","w",stdout);
read(T);
while (T--){
read(n);
for (int i=1;i<n;i++) read(a[i]);
for (int i=1;i<=n;i++){
double cur=1,a=0,b=0;
for (int j=i;j<=n;j++){
if (~(j-i)&1) a+=cur; else b+=cur;
val[i][j]=2.0*sqrt(a*b);
lx[i][j]=sqrt(b/a);
if (~(j-i)&1) rx[i][j]=cur*lx[i][j]; else rx[i][j]=cur/lx[i][j];
cur=::a[j]/cur;
}
}
for (int i=1;i<=n;i++) for (int j=i;j<=n;j++) f[i][j]=1e30;
for (int i=1;i<=n;i++) f[1][i]=val[1][i];
for (int i=1;i<n;i++){
tot=0;
for (int j=1;j<=i;j++)
tmp[++tot]=abcd(rx[j][i],f[j][i]);
sort(tmp+1,tmp+tot+1);
g[tot+1]=1e30;
for (int j=tot;j;j--)
g[j]=min(g[j+1],tmp[j].second);
for (int k=i+2;k<=n;k++){
double tem=a[i]/lx[i+1][k];
int iter=lower_bound(tmp+1,tmp+tot+1,abcd(tem-1e-8,-1e30))-tmp;
f[i+1][k]=val[i+1][k]+g[iter];
}
}
double ans=1e30;
for (int i=1;i<=n;i++)
ans=min(ans,f[i][n]);
printf("Case #%d: %.5lf\n",++Case,ans);
}
return 0;
}