某考试 T2 yja
2.1 Description
在平面上找 n 个点, 要求这 n 个点离原点的距离分别为 r1, r2, ..., rn.
最大化这 n 个点构成的凸包面积, 凸包上的点的顺序任意.
2.2 Input Format
第一行一个整数 n.
接下来一行 n 个整数依次表示 ri .
2.3 Output Format
输出一个实数表示答案, 要求绝对误差或相对误差 ≤ 10−6.
2.4 Sample
2.4.1 Input
4
5
8
58
85
2.4.2 Output
2970
2.5 Constraints
对于前 20% 的数据, n ≤ 3;
对于前 40% 的数据, n ≤ 4;
对于另 20% 的数据, r1 = r2 = ... = rn;
对于 100% 的数据, 1 ≤ n ≤ 8, 1 ≤ ri ≤ 1000.
考试的时候看了一会题就感觉是个拉格朗日乘数法的题,但是我为什么往叉积求面积上想了23333,于是就出来了n个限制然后还得三分套三分套三分.......
于是在想如何优化求偏导都是零的过程中时间过去了2.5h 。。。(我好菜啊我现在退役算了)。
正解其实就是拉格朗日乘数法,但是求面积是用(1/2) * a * b * sin(a和b的夹角) ,于是这样就只有一个总限制,那就是所有角的和=2*π。
并且利用所有变量的偏导等于0,我们可以得到 λ = r1 * r2 *cos(r1和r2夹角) = .... ,因为最优情况肯定没有超过π的角,而cos在[0,π]上单调,所以我们二分一下λ 就可以得到最后的答案。
当然因为不一定所有点都在凸包上出现,并且每个点出现的顺序不一定,所以我们还要枚举选ri前几大的点(贪心)和它们的顺序(全排列一下就好了)。
#include<bits/stdc++.h>
#define ll long long
#define D double
using namespace std;
const int maxn=10;
const D pi=acos(-1);
D miu,l,r,thita[maxn];
int n,R[maxn]; inline D calc(){
D tot=0;
for(int i=1;i<n;i++) thita[i]=acos(miu/(R[i]*R[i+1])),tot+=thita[i];
thita[n]=acos(miu/(R[n]*R[1])),tot+=thita[n];
return tot;
} inline void solve(){
D ans=0,now,le,ri; sort(R+1,R+n+1);
for(;n>=3;){
le=-R[1]*R[2],ri=-le; while(1){
l=le,r=ri; while(r-l>=1e-8){
miu=(l+r)/2;
if(calc()>=2*pi) l=miu;
else r=miu;
} if(fabs(calc()-2*pi)<=1e-6){
now=0;
for(int i=1;i<n;i++) now+=sin(thita[i])*(D)R[i]*(D)R[i+1];
now+=sin(thita[n])*(D)R[n]*(D)R[1];
ans=max(ans,now);
} if(!next_permutation(R+1,R+n+1)) break;
} sort(R+1,R+n+1),n--;
for(int i=1;i<=n;i++) R[i]=R[i+1];
} printf("%.11lf\n",ans/2);
} int main(){
scanf("%d",&n);
for(int i=1;i<=n;i++) scanf("%d",R+i);
solve();
return 0;
}
某考试 T2 yja的更多相关文章
- 9.13 考试 T2 区间
删区间 题意: 给出一个长度为
- 某考试T2 frog
题目背景 无 题目描述 数轴上有 n 只青蛙,分别编号为 1 到 n.青蛙 i 的初始位置的坐标为 xi. 它们准备进行如下形式的移动:每轮包括 m 次跳跃,第 i 次跳跃由青蛙 ai(1 < ...
- 某考试 T2 Tree
2 树 2.1 题目描述 给一棵n 个节点的树,节点分别编号为0 到n - 1.你可以通过如下的操作来修改这棵树:首先先删去树上的一条边,此时树会分裂为两个连通块,然后在两个连通块之间加上一条新的边使 ...
- 题解 2020.10.24 考试 T2 选数
题目传送门 题目大意 见题面. 思路 本来以为zcx.pxj变强了,后来发现是SPJ出问题了...考试的时候感觉有点人均啊...结果自己还是只想出来一半. 我们假设 \(f(x)=(\lfloor\f ...
- 2019.3.7考试T2 离线数论??
$ \color{#0066ff}{ 题目描述 }$ 一天,olinr 在 luogu.org 刷题,一点提交,等了一分钟之后,又蛙又替. olinr 发动了他的绝招,说:"为啥啊???&q ...
- 某考试 T2 orzcyr
非常nice的一道行列式的题目. 考虑如果没有路径不相交这个限制的话,那么这个题就是一个行列式:设 a[i][j] 为从编号第i小的源点到编号第j小的汇点的路径条数,那么矩阵a[][]的行列式就是的答 ...
- 2019.2.26考试T2 矩阵快速幂加速DP
\(\color{#0066ff}{题解 }\) 可以发现, 数据范围中的n特别小,容易想到状压 可以想到类似于状压DP的思路,按列进行转移 那么应该有3维,\(f[i][j][k]\)代表到第i列, ...
- 2019.2.10考试T2, 多项式求exp+生成函数
\(\color{#0066ff}{ 题目描述 }\) 为了减小文件大小,这里不写一堆题目背景了. 请写一个程序,输入一个数字N,输出N个点的森林的数量.点有标号. 森林是一种无向图,要求图中不能存在 ...
- 某考试 T2 Seg
Seg [问题描述]数轴上有n条线段,第i条线段的左端点是a[i],右端点是b[i].Bob发现1~2n共2n个整数点,每个点都是某条线段的端点.这些线段有如下两类特点:1 x y,表示第x条线段和第 ...
随机推荐
- 【转】windows server 2012 安装 VC14(VC2015) 安装失败解决方案
系统环境如下:cmd命令行-输入 systeminfo 如下图 - The VC14 builds require to have the Visual C++ Redistributable for ...
- (转)用@Resource注解完成属性装配
http://blog.csdn.net/yerenyuan_pku/article/details/52858878 前面我们讲过spring的依赖注入有两种方式: 使用构造器注入. 使用属性set ...
- 闲着蛋疼没事干,写个Mac端的Kcptun Client管理器
原理: 执行一行脚本 输入服务器地址,端口,密码等做了图形化编辑 可以控制Kcptun是否正在运行 App已上传github https://github.com/nicky2k8/KcptunCli ...
- MySQL-01 MySQL数据库安装指南
学习要点 MySQL数据库的安装和设置 下载mysql mysql官网:https://www.mysql.com/downloads/ 主要版本: Oracle MySQL Cloud Servic ...
- cobbler 无人值守-安装
环境准备 准备两台主机,如centos6和centos7 centos7当作server服务器 关闭selinux 关闭防火墙 安装 cobbler包光盘里是没有的,要配置epel源,这里就说怎么配置 ...
- RPM Package Manager
本文大部分内容来自鸟哥的Linux私房菜,并且由作者根据自己的学习情况做了一些更改,鸟哥原文链接 1. 程序的安装方式 源代码安装:利用厂商释出的Tarball 来进行软件的安装,每次安装程序都需要检 ...
- 特殊权限和facl
目 录 第1章 FACL访问控制 1 1.1 FACL简介 1 1.2 getfacl命令查看acl权限 1 1.3 setfacl设置文件的cal权限 1 1.4 批量添加a ...
- 关于web页面优化
简单汇总了一下web的优化方案(主要的前端优化策略) 减少http请求次数 文件合并(js.css.图片):ps:多个图片合并之后,总体积会变小 内联图片,即data:URL scheme,但容易导致 ...
- php S3
转载自:http://www.cnblogs.com/wangxusummer/p/6398772.html <?php /* * To change this license header, ...
- MySQL数据库安装与Navicat Premium 12 安装与破解
一.文件下载: MySQL:官网:https://www.mysql.com/downloads/(现在最新的是5.7版) 下载路径:"Downloads" ==>> ...