[Jzoj5606]Yja

题意

你需要找到$n$个点,每个点离原点的距离分别为$R_1,R_2,\ldots,R_n$问$n$个点形成的凸包的最大面积是多少$?$(不要求所有$n$个点都在凸包上)

题解

看到$n\le8$,想这题怕不是爆搜$?$

假设搜出了$m$个点,分别是$a_1\ldots a_m$,为了方便令$a_{m+1}=a_1$

根据多边形面积的定义可以知道我们要最大化的就是$\frac12\sum_{i=1}^mr_{a_i}r_{a_{i+1}}\sin x_i$,其中$\sum_{i=1}^mx_i=2\pi$

这是要求一个有条件最值,考虑用拉格朗日乘数法来试一试

设$f(x_1,\ldots,x_n)=\sum_{i=1}^mr_{a_i}r_{a_{i+1}}\sin x_i,g(x_1,\ldots,x_n)=\sum_{i=1}^mx_i=2\pi$

在这题里面

$$\frac{\partial f+(-\lambda)g}{\partial x_i}=r_{a_i}r_{a_{i+1}}\cos x_i+(-\lambda)=0\\ \Rightarrow \sum_{i=1}^mx_i=\sum_{i=1}^m\arccos\frac{\lambda}{r_{a_i}r_{a_{i+1}}}=2\pi\\ $$
这是 $\arccos$的函数图像

然后考虑到$\arccos$是单调递减的所以可以二分求出当前的$\lambda$

注意$\lambda\in[-min\{r_{a_i}r_{a_{i+1}}\},min\{r_{a_i}r_{a_{i+1}}\}]$

#include<bits/stdc++.h>
#define fp(i,a,b) for(register int i=a,I=b+1;i<I;++i)
#define fd(i,a,b) for(register int i=a,I=b-1;i>I;--i)
#define go(u) for(register int i=fi[u],v=e[i].to;i;v=e[i=e[i].nx].to)
#define file(s) freopen(s".in","r",stdin),freopen(s".out","w",stdout)
template<class T>inline bool cmax(T&a,const T&b){return a<b?a=b,1:0;}
template<class T>inline bool cmin(T&a,const T&b){return a>b?a=b,1:0;}
using namespace std;
const int N=10;
const double eps=1e-8,pi2=2*acos(-1);
typedef int arr[N];
int n,vis[N];double ans,p[N],q[N],r[N];
inline void chk(int m){
    double L,R,T=1e6,mid,tp;p[m+1]=p[1];
    fp(i,1,m)q[i]=p[i]*p[i+1],cmin(T,q[i]);L=-(R=T);
    while(R-L>eps){
        mid=(L+R)*0.5;tp=0;
        fp(i,1,m)tp+=acos(mid/q[i]);
        if(tp>=pi2)L=mid;
        else R=mid;
    }tp=0;if(T-L<eps)return;
    fp(i,1,m)tp+=q[i]*sin(acos(L/q[i]));
    cmax(ans,tp);
}
void dfs(int d){
    if(d>3)chk(d-1);if(d==n+1)return;
    fp(i,1,n-2)if(!vis[i])p[d]=r[i],vis[i]=1,dfs(d+1),vis[i]=0;
}
int main(){
    // #ifndef ONLINE_JUDGE
        file("yja");
    // #endif
    scanf("%d",&n);
    fp(i,1,n)scanf("%lf",r+i);sort(r+1,r+n+1);
    p[1]=r[n],p[2]=r[n-1];dfs(3);
    printf("%.6lf",ans*0.5);
return 0;
}