HDU 5531
题目大意:
给定一个n边形的顶点
以每个顶点为圆心画圆(半径可为0)
每个顶点的圆要和它相邻顶点的圆相切(不相邻的可相交)
求所有圆的最小面积总和并给出所有圆的半径
设半径为r1 r2 ... rn,顶点距离为L1 L2 ... Ln
当顶点数为奇数时 由
r1+r2=L1
r2+r3=L2
......
rn+r1=Ln
可得 r1+r1=L1-L2+L3-......+Ln
只要中间所有半径 r >= 0 那么r1就有唯一解
当顶点数为偶数时 由
r1+r2=L1
r2+r3=L2
......
rn+r1=Ln
可得 0=L1-L2+L3-......-Ln
即 L1+L3+...+Ln-1=L2+L4+...+Ln 时 才有解
由 上式可得
r2=L1-r1
r3=L2-L1+r1
......
r1=Ln-......+r1
也就是每个 r 都可以转换为
第偶数个顶点时 r=x+r1 或
第奇数个顶点时 r=x-r1 的形式
已知 r>=0 得
第偶数个顶点时 x+r1>=0 即 r1>=-x
第奇数个顶点时 x-r1>=0 即 r1<=x
那么由所有顶点的约束 就能得到r1的范围 [L,R]
圆的面积总和 S=PI*(r1^2+r2^2+r3^2+...+rn^2)
而每个r^2都可以转换为 r^2=(r1+x)^2 或 r^2=(r1-x)^2
如: r2^2=(r1-L1)^2, r3^2=(r1+L2-L1)^2
即 递推时 x=Li-x
第偶数个顶点时 ri^2=(r1-x)^2
第奇数个顶点时 ri^2=(r1+x)^2
那么 S=PI*(a*r1^2+b*r1+c)
如: r2^2=(r1-L1)^2, r3^2=(r1+L2-L1)^2
则 S1=PI*(1*r1^2+0*r1+0) (a=1,b=0,c=0,x=0)
S2=PI*(S1+r2^2)=PI*(S1+(r1-x)^2)
=PI*(S1+r1^2-2*x*r1+x*x) (a=a+1,b=b+(-1)*2*x,c=c+x*x,x=L1-x)
S3=PI*(S2+r3^2)=PI*(S2+(r1+x)^2)
=PI*(S1+r1^2+2*x*r1+x*x) (a=a+1,b=b+(1)*2*x,c=c+x*x,x=L2-x)
即
第偶数个顶点时 a+=1,b+=(-1)*2*x,c+=x*x,x=Li-x
第奇数个顶点时 a+=1,b+=1*2*x,c+=x*x,x=Li-x
又 a*r1^2+b*r1+c=0 时存在极值点 1-b/(2*a)
再与r1的范围比较一下 就能得到r1 从而得到S的最小值
或者直接利用[L,R] 用三分法求解
#include <bits/stdc++.h>
#define ll long long
#define INF 0x3f3f3f3f
using namespace std; const double eps = 1e-;
const double PI = acos(-1.0);
double add(double a,double b) {
if(abs(a+b)<eps*(abs(a)+abs(b))) return ;
return a+b;
}
struct P {
double x,y;
P(){};
P(double _x,double _y):x(_x),y(_y){}
P operator - (P p) {
return P(add(x,-p.x),add(y,-p.y));
}
double dot(P p) {
return add(x*p.x,y*p.y);
}
}p[];
int n;
double ans;
double r[], len[];
double lenV(P p) {
return sqrt(p.dot(p));
}
double getArea(double r0) {
double ans=; r[]=r0;
for(int i=;i<n;i++) {
if(r[i]<-eps) return ;
ans+=r[i]*r[i];
r[i+]=len[i]-r[i];
//printf("len[i]%lf r[i]%lf\n",len[i],r[i]);
}
return ans*PI;
} bool solve() {
ll sum=;
for(int i=;i<n;i++) {
len[i]=lenV(p[(i+)%n]-p[i]);
sum+= i% ? -len[i]:len[i];
}
if(n&) {
double r0=sum/2.0;
if(r0<-eps) return ;
ans=getArea(r0);
//printf("ans = %lf\n",ans);
} else {
if(sum) return ;
ll a,b,c;
a=b=c=;
ll t=,L=,R=INF,k=;
for(int i=;i<n;i++) {
//printf("%lld %lld %lld\n",a,b,c);
a++;
b+=2.0*k*t;
c+=t*t;
if(k==) L=max(L,-t);
else R=min(R,t);
k*=-;
t=len[i]-t;
}
if(L>R) return ;
double mid=-b/2.0/a;
mid=max(mid,(double)L);
mid=min(mid,(double)R);
ans=getArea(mid);
}
if(ans<eps) return ;
return ;
} int main()
{
int t; scanf("%d",&t);
while(t--) {
scanf("%d",&n);
for(int i=;i<n;i++)
scanf("%lf%lf",&p[i].x,&p[i].y);
if(solve()) {
printf("%.2f\n",ans);
for(int i=;i<n;i++)
printf("%.2f\n",r[i]);
} else printf("IMPOSSIBLE\n");
} return ;
}
HDU 5531的更多相关文章
- HDU 5531 Rebuild (2015长春现场赛,计算几何+三分法)
Rebuild Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 262144/262144 K (Java/Others)Total S ...
- 2015ACM/ICPC亚洲区长春站 E hdu 5531 Rebuild
Rebuild Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 262144/262144 K (Java/Others)Total S ...
- 2015长春 HDU 5531 Rebuild
题意:n个顶点组成的多边形能否形成正多边形? #include <cstdio> #include <cstring> #include <cmath> #incl ...
- HDU 5531 Rebuild
2015 ACM/ICPC 长春现场赛 E题 三分. 如果节点个数是奇数,那么直接列方程可以求解,因为,如果第一个圆半径变大,必然导致最后一个圆的半径变大, 所以,节点是奇数的时候,要么无解,要么只有 ...
- HDOJ 2111. Saving HDU 贪心 结构体排序
Saving HDU Time Limit: 3000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total ...
- 【HDU 3037】Saving Beans Lucas定理模板
http://acm.hdu.edu.cn/showproblem.php?pid=3037 Lucas定理模板. 现在才写,noip滚粗前兆QAQ #include<cstdio> #i ...
- hdu 4859 海岸线 Bestcoder Round 1
http://acm.hdu.edu.cn/showproblem.php?pid=4859 题目大意: 在一个矩形周围都是海,这个矩形中有陆地,深海和浅海.浅海是可以填成陆地的. 求最多有多少条方格 ...
- HDU 4569 Special equations(取模)
Special equations Time Limit:1000MS Memory Limit:32768KB 64bit IO Format:%I64d & %I64u S ...
- HDU 4006The kth great number(K大数 +小顶堆)
The kth great number Time Limit:1000MS Memory Limit:65768KB 64bit IO Format:%I64d & %I64 ...
随机推荐
- Jmeter json和正在提取及断言
Jmeter json和正在提取及断言 一.提取所需要断言的内容: 响应数据如下:加入需要提取userid为10477的值 { "pub": { , ", " ...
- Tomcat运行错误示例
tomcat运行错误示例 当出现这种错误时,如果是自己配置的tomcat,需要找/conf/server.xml文件. 如果是使用的eclipse tomcat 插件,需要在你的工作空间 找到 如下文 ...
- 安装debian总结以及编译linux内核
1. 安装debian 使用unetbootin(http://unetbootin.sourceforge.net/)来创建启动盘,并且下载debian的基本包. 将磁盘进行压缩操作,并且保留出一个 ...
- Flink DataStream API
Data Sources 源是程序读取输入数据的位置.可以使用 StreamExecutionEnvironment.addSource(sourceFunction) 将源添加到程序.Flink 有 ...
- JAVA集合--Iterator接口
本文首发于cartoon的博客 转载请注明出处:https://cartoonyu.github.io/cartoon-blog 上一篇文章中我在集合元素的遍历中已经有涉及到I ...
- 68.最大k乘积问题 (15分)
C时间限制:3000 毫秒 | C内存限制:3000 Kb题目内容:设I是一个n位十进制整数.如果将I划分为k段,则可得到k个整数.这k个整数的乘积称为I的一个k乘积.试设计一个算法,对于给定的I和 ...
- jmeter-测试webservice接口
测试webservice接口(soap类型接口) 一.webservice协议的本质 一个经过封装的post类型的HTTP请求 Web service一般就是用SOAP协议通过HTTP来调用它,其实他 ...
- IPC$渗透使用
简介 IPC$(Internet Process Connection)是共享"命名管道"的资源,它是为了让进程间通信而开放的命名管道,通过提供可信任的用户名和口令,连接双方可以建 ...
- linux 查看服务器序列号(S/N)
[root@node1~]# dmidecode -t 查看支持的参数 dmidecode: option requires an argument -- 't' Type number or key ...
- js排他功能示例
<!DOCTYPE html><html lang="en"><head> <meta charset="UTF-8" ...