#include <iostream>
#include <map>
#include <string>
#include <iomanip>

using namespace std;
long double p[500005];
long double dp[500005];
pair<long double,long double>aint[1200000];
int eval(pair<long double,long double>a,pair<long double,long double>b,int unde)
{
    if(a.first*unde+a.second<b.first*unde+b.second)
        return 1;
    return 0;
}
void add(int node,int st,int dr,pair<long double,long double>func)
{
    int mij=(st+dr)/2;
    if(eval(aint[node],func,mij)==1)
        swap(aint[node],func);
    if(st==dr)
    {
        return;
    }
    if(eval(aint[node],func,st)==1)
        add(node*2,st,mij,func);
    else
        add(node*2+1,mij+1,dr,func);
}
long double best;
void qr(int node,int st,int dr,int poz)
{
    best=max(best,aint[node].first*poz+aint[node].second);
    if(st!=dr)
    {
        int mij=(st+dr)/2;
        if(poz<=mij)
            qr(node*2,st,mij,poz);
        else
            qr(node*2+1,mij+1,dr,poz);
    }
}
int main()
{
    int n,i,j,k,l,m,a,b,c,d,z;
    cin>>n;
    for(i=1;i<=n;i++)
    {
        cin>>p[i];
    }
    dp[n]=n*p[n];
    add(1,1,n,{-p[n],dp[n]});
    long double sol=dp[n];
    for(i=n-1;i>=1;i--)
    {
        best=0;
        qr(1,1,n,i);
        dp[i]=p[i]*i+p[i]*best;
        sol=max(sol,dp[i]);
        add(1,1,n,{-p[i],dp[i]});
    }
    cout<<fixed<<setprecision(12)<<sol;
    return 0;
}