Problem:
MathJax.Hub.Config({ tex2jax: { inlineMath: [["$","$"], ["\\(","\\)"]], processEscapes: true }});

Problem E:
HullMarathon

うさぎはフルマラソンという競技が好きである. この競技はチームで行う. チームメンバーは競技開始前に原点に集まる. 競技開始と同時に走り出し, 1 分後に立ち止まる．このとき，チームメンバーの位置の凸包の面積が最も大きなチームが勝ちとなる.

あなたは $N$ 匹のうさぎからなるチームの監督である. $i$ 匹目のうさぎは 1 分で $r_i$ 移動することができる. このチームが最適な戦略をとった場合の, 1 分後の凸包の面積の最大値を求めよ.

Constraints

$N$ will be between 3 and 8, inclusive.

$r_i$ will be an integer between 1 and 1,000, inclusive.

Input

入力は以下の形式で与えられる:

$N$

$r_1$

...

$r_N$

Output

凸包の面積の最大値を表す実数を 1 行に出力せよ. 小数点以下何桁出力してもよいが, 絶対誤差または相対誤差が $10^{-6}$ 以下のとき Accepted になる.

Sample Input 1

4
5
8
58
85

Sample Output 1

2970.000000000

Sample Input 2

6
1
1
1
1
1
1

Sample Output 2

2.598076211

Solution in C++14 (judge verdict: Accepted)

#include <bits/stdc++.h>
using namespace std;

typedef long long ll;
typedef unsigned long long ull;
typedef pair<ll, ll> P;

#define fi first
#define se second
#define repl(i,a,b) for(ll i=(ll)(a);i<(ll)(b);i++)
#define rep(i,n) repl(i,0,n)
#define all(x) (x).begin(),(x).end()
#define dbg(x) cout<<#x"="<<x<<endl
#define mmax(x,y) (x>y?x:y)
#define mmin(x,y) (x<y?x:y)
#define maxch(x,y) x=mmax(x,y)
#define minch(x,y) x=mmin(x,y)
#define uni(x) x.erase(unique(all(x)),x.end())
#define exist(x,y) (find(all(x),y)!=x.end())
#define bcnt __builtin_popcount

#define INF 1e16
#define mod 1000000007
#define PI acos(-1.0)
#define eps 1e-6

int N;
vector<double> r;

int main(){
  cin>>N;
  r.resize(N);
  rep(i,N){
    cin>>r[i];
  }
  sort(all(r));
  double res=0;
  do{
    repl(k,3,N+1){
      double lb=0,ub=PI;
      rep(_,100){
        double X=(lb+ub)/2.0;
        double cos1=cos(X);
        double theta=0;
        bool ok=true;
        rep(i,k){
          double p=((r[0]*r[1])/(r[i]*r[(i+1)%k]))*cos1;
          if(p<-1.0)p=-1.0;
          if(p>+1.0)p=+1.0;
          theta+=acos(((r[0]*r[1])/(r[i]*r[(i+1)%k]))*cos1);
        }
        if(theta<2*PI)lb=X;
        else ub=X;
      }
      double sum=0;
      double cos1=cos(lb);
      double thetasum=0;
      rep(i,k){
        double p=((r[0]*r[1])/(r[i]*r[(i+1)%k]))*cos1;
        if(p<-1.0)p=-1.0;
        if(p>+1.0)p=+1.0;
        sum+=(r[i]*r[(i+1)%k])*sin(acos(p));
        thetasum+=acos(p);
      }
      if(abs(thetasum-2.0*PI)<eps)maxch(res,sum/2.0);
    }
  }while(next_permutation(all(r)));
  printf("%.10f\n",res);
  return 0;
}