[SHTSC 2014] 信号增幅仪
最小覆盖圆算法。看着题解半蒙半抄的搞过去了…
主要参考以下
http://blog.csdn.net/acdreamers/article/details/9406735
http://blog.csdn.net/lthyxy/article/details/6661250
http://blog.himdd.com/archives/2666
其中有一个求外心的过程是错的…害我调了好久…还把x和y搞反…可是就算是错的我居然过了80%的数据…
仍然几个疑问:
(1)最小覆盖圆算法时间复杂度的证明
(2)椭圆的中心未定的情况下,为什么以原点为中心旋转&缩放即可?(我几何不好QAQ)
- program amplifier;
- uses math;
- type point=record
- x,y:extended;
- end;
- const maxn=;
- pi=3.14159265358979;
- eps=0.00000001;
- var n,i,p:longint;
- a,r:extended;
- o:point;
- v:array[..maxn+] of point;
- procedure rotate(i:longint);
- var x,y:extended;
- begin
- x:=v[i].x*cos(-a)-v[i].y*sin(-a);
- y:=v[i].x*sin(-a)+v[i].y*cos(-a);
- v[i].x:=x/p;
- v[i].y:=y;
- end;
- function circumcenter(p,q,r:point):point;
- var a,b,c,d,e,f:extended;
- begin
- a:=q.x-p.x;
- b:=q.y-p.y;
- c:=-(q.x*q.x+q.y*q.y-p.x*p.x-p.y*p.y)/;
- d:=r.x-p.x;
- e:=r.y-p.y;
- f:=-(r.x*r.x+r.y*r.y-p.x*p.x-p.y*p.y)/;
- circumcenter.y:=(-c*d+f*a)/(b*d-e*a);
- circumcenter.x:=(-c*e+f*b)/(a*e-b*d);
- end;
- function dis(a,b:point):extended;
- begin
- dis:=sqrt(sqr(a.x-b.x)+sqr(a.y-b.y));
- end;
- procedure min_cover_circle;
- var i,j,k:longint;
- begin
- o:=v[];r:=;
- for i:= to n do
- begin
- if dis(v[i],o)>r+eps then
- begin
- o:=v[i];r:=;
- for j:= to i do
- if dis(v[j],o)>r+eps then
- begin
- o.x:=(v[i].x+v[j].x)/;
- o.y:=(v[i].y+v[j].y)/;
- r:=dis(v[j],o);
- for k:= to j do
- if dis(v[k],o)>r+eps then
- begin
- o:=circumcenter(v[i],v[j],v[k]);
- r:=dis(o,v[i]);
- end;
- end;
- end;
- end;
- end;
- begin
- //assign(input,'amplifier20.in');reset(input);
- //assign(output,'amplifier20.out');rewrite(output);
- readln(n);
- for i:= to n do
- readln(v[i].x,v[i].y);
- readln(a);a:=a/*pi;
- readln(p);
- for i:= to n do
- rotate(i);
- if n= then
- begin
- writeln(dis(v[],v[])/::);
- //close(input);close(output);
- halt;
- end;
- min_cover_circle;
- writeln(r::);
- //close(input);close(output);
- end.
amplifier
其实说起来也不是很难,增量法我在考场居然也想到了,就是那时没想到可以以原点为中心,算法也写的有问题- -
来一发考场逗比错误代码,有空分析下错哪里了0 0
- program amplifier;
- uses math;
- const pi=3.141592653589;
- type vtype=record
- x,y:real;
- end;
- var n,p,i,n1,n2,n3:longint;
- xx,yy,xt,yt,a,aa,x0,y0:real;
- v:array[..] of vtype;
- procedure cac(x1,y1,x2,y2,x3,y3:real);
- var a1,a2,b1,b2,c1,c2:real;
- begin
- if x1*y2+x2*y3+x3*y1-x3*y2-x2*y1-x1*y3= then exit;
- a1:=*(x1-x3);b1:=*p*p*(y1-y3);c1:=(x1*x1-x3*x3)+p*p*(y1*y1-y3*y3);
- a2:=*(x1-x2);b2:=*p*p*(y1-y2);c2:=(x1*x1-x2*x2)+p*p*(y1*y1-y2*y2);
- xt:=(b2*c1-b1*c2)/(a1*b2-a2*b1);
- yt:=(a1*c2-a2*c1)/(a1*b2-a2*b1);
- a:=sqrt(sqr(x1-xt)+p*p*sqr(y1-yt));
- end;
- function ini(t:integer;x0,y0,a:real):boolean;
- var tt:real;
- begin
- tt:=sqr(v[t].x-x0)+p*p*sqr(v[t].y-y0);
- if tt<=a*a then exit(true) else exit(false);
- end;
- procedure work2;
- begin
- xx:=(v[n1].x+v[n2].x)/;
- yy:=(v[n1].y+v[n2].y)/;
- a:=sqrt(sqr(v[n1].x-xx)+p*p*sqr(v[n1].y-yy));
- writeln((a/p)::);
- end;
- function line(n1,n2,n3:integer):boolean;
- var x1,x2,x3,y1,y2,y3:real;
- begin
- x1:=v[n1].x;x2:=v[n2].x;x3:=v[n3].x;
- y1:=v[n1].y;y2:=v[n2].y;y3:=v[n3].y;
- if x1*y2+x2*y3+x3*y1-x3*y2-x2*y1-x1*y3= then exit(true) else exit(False);
- end;
- procedure workline;
- begin
- if ((v[n1].x<v[n2].x) and (v[n2].x<v[n3].x)) or ((v[n3].x<v[n2].x) and (v[n2].x<v[n1].x)) then
- begin
- n2:=n3;n3:=n3+;exit;
- end;
- if ((v[n2].x<v[n1].x) and (v[n1].x<v[n3].x)) or ((v[n3].x<v[n1].x) and (v[n1].x<v[n2].x)) then
- begin
- n1:=n2;n2:=n3;n3:=n3+;exit;
- end;
- if ((v[n1].x<v[n3].x) and (v[n3].x<v[n2].x)) or ((v[n2].x<v[n3].x) and (v[n3].x<v[n1].x)) then
- begin
- n3:=n3+; exit;
- end;
- end;
- begin
- assign(input,'amplifier.in');reset(input);
- assign(output,'amplifier.out');rewrite(output);
- readln(n);
- for i:= to n do
- readln(v[i].x,v[i].y);
- readln(a);readln(p);
- for i:= to n do
- begin
- x0:=v[i].x;y0:=v[i].y;
- v[i].x:=cos(-a/*pi)*x0-sin(-a/*pi)*y0;
- v[i].y:=sin(-a/*pi)*x0+cos(-a/*pi)*y0;
- end;
- v[n+].x:=;v[n+].y:=;
- if n= then writeln();
- if n= then
- begin
- n1:=;n2:=;
- work2;
- end;
- if n>= then
- begin
- n1:=;n2:=;n3:=;
- while (line(n1,n2,n3)) and (n3<>n+) do workline;
- if n3=n+ then
- begin
- work2;
- exit;
- end;
- cac(v[n1].x,v[n1].y,v[n2].x,v[n2].y,v[n3].x,v[n3].y);
- xx:=xt;yy:=yt;aa:=a;
- for i:=n3+ to n do
- begin
- if not ini(i,xx,yy,aa) then
- begin
- cac(v[n1].x,v[n1].y,v[n2].x,v[n2].y,v[i].x,v[i].y);
- if ini(n3,xt,yt,a) then
- begin
- n3:=i;xx:=xt;yy:=yt;aa:=a;continue;
- end;
- cac(v[n1].x,v[n1].y,v[i].x,v[i].y,v[n3].x,v[n3].y);
- if ini(n2,xt,yt,a) then
- begin
- n2:=n3;n3:=i;xx:=xt;yy:=yt;aa:=a;continue;
- end;
- cac(v[i].x,v[i].y,v[n2].x,v[n2].y,v[n3].x,v[n3].y);
- if ini(n1,xt,yt,a) then
- begin
- n1:=n2;n2:=n3;n3:=i;xx:=xt;yy:=yt;aa:=a;continue;
- end;
- end;
- end;
- writeln((a/p)::);
- end;
- close(input);close(output);
- end.
amplifier-wrong
[SHTSC 2014] 信号增幅仪的更多相关文章
- BZOJ 3564: [SHOI2014]信号增幅仪 最小圆覆盖
3564: [SHOI2014]信号增幅仪 题目连接: http://www.lydsy.com/JudgeOnline/problem.php?id=3564 Description 无线网络基站在 ...
- [LOJ 2190] 「SHOI2014」信号增幅仪
[LOJ 2190] 「SHOI2014」信号增幅仪 链接 链接 题解 坐标系直到 \(x\) 轴与椭圆长轴平行 点的坐标变换用旋转公式就可以了 因为是椭圆,所以所有点横坐标除以 \(p\) 然后最小 ...
- 【bzoj3564】 [SHOI2014]信号增幅仪
题目描述: 无线网络基站在理想状况下有效信号覆盖范围是个圆形.而无线基站的功耗与圆的半径的平方成正比. 现给出平面上若干网络用户的位置,请你选择一个合适的位置建设无线基站.... 就在你拿起键盘准备开 ...
- BZOJ3564 : [SHOI2014]信号增幅仪
先把所有点绕原点逆时针旋转(360-a)度,再把所有点横坐标除以放大倍数p,最后用随机增量法求最小圆覆盖即可. 时间复杂度期望$O(n)$ #include<cstdio> #includ ...
- BZOJ 3564 信号增幅仪
题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=3564 题意:给出平面上n个点,画出一个椭圆,椭圆的长轴是短轴的p倍,且长轴的方向为x轴逆时 ...
- [BZOJ 3564] [SHOI2014] 信号增幅仪 【最小圆覆盖】
题目链接:BZOJ - 3564 题目分析 求最小椭圆覆盖,题目给定了椭圆的长轴与 x 轴正方向的夹角,给定了椭圆长轴与短轴的比值. 那么先将所有点旋转一个角度,使椭圆长轴与 x 轴平行,再将所有点的 ...
- BZOJ3564 信号增幅仪
http://www.lydsy.com/JudgeOnline/problem.php?id=3564 思路:先旋转坐标系,再缩进x坐标,把椭圆变成圆,然后做最小圆覆盖. 还有,为什么用srand( ...
- BZOJ 3564: [SHOI2014]信号增幅仪(随机增量法)
如果是个圆的话好办,如果是拉成椭圆呢?直接压回去!!! 然后随机增量法就行了 CODE: #include<cstdio> #include<iostream> #includ ...
- 2018.10.15 bzoj3564: [SHOI2014]信号增幅仪(坐标处理+最小圆覆盖)
传送门 省选考最小圆覆盖? 亦可赛艇(你们什么都没看见) 在大佬的引领下成功做了出来. 就是旋转坐标使椭圆的横轴跟xxx轴平行. 然后压缩横坐标使得其变成一个圆. 然后跑最小覆盖圆就可以了. 注意题目 ...
随机推荐
- 在 ASP.NET MVC 中使用 HTML Helpers 的那些事
在 ASP.NET MVC 中使用 HTML Helpers 方法,可以返回得到标准的 HTML 标签,就像 <input>.<button> 或者 <img> 等 ...
- html select的事件 方法 属性
事件 onactivate 当对象设置为活动元素时触发. onafterupdate 当成功更新数据源对象中的关联对象后在数据绑定对象上触发. onbeforeactivate 对象要被设置为当前元素 ...
- 在ionic/cordova中使用极光推送插件(jpush)
Stpe1:创建一个项目(此处使用的是tab类型的项目,创建方式可参照我前一篇如何离线创建Ionic1项目) Stpe2:修改项目信息 打开[config.xml]修改下图内容:
- 回调函数及数组中sort()方法实现排序的原理
1.回调函数:把一个方法A当一个参数值传递到另外一个函数B中,在B执行的过程当中我们随时根据需求让A方法执行: 什么是回调 :它是异步编程基本的方法,需要异步处理的时候一般采用后续传递的方式,将后 ...
- A*算法 -- 八数码问题和传教士过河问题的代码实现
前段时间人工智能的课介绍到A*算法,于是便去了解了一下,然后试着用这个算法去解决经典的八数码问题,一开始写用了挺久时间的,后来试着把算法的框架抽离出来,编写成一个通用的算法模板,这样子如果以后需要用到 ...
- linux下使用多线程编写的聊天室
自从开始学linux网络编程后就想写个聊天室,一开始原本打算用多进程的方式来写,可是发觉进程间的通信有点麻烦,而且开销也大,后来想用多线程能不能实现呢,于是便去看了一下linux里线程的用法,实际上只 ...
- 多个div同时居中的写法
多个div在某个div的中间,他们个数不一定但是需要在那个父级div中显示(和margin:0 auto一样的效果) <!DOCTYPE html><html lang=" ...
- svn-多个项目版本库和自动同步更新post-commit
由于项目测试需求,需要远程服务器上使用svn做版本控制. 需求: 1,项目test1,项目test2,各自独立版本库,各自独立用户权限,便于项目管理 2,同步提交,本地svn提交至版本库后,服务器上的 ...
- 用原生js获取class
<!DOCTYPE html> <html> <head> <meta http-equiv="Content-Type" content ...
- IIS发布网站出现“未能加载文件或程序集“System.Data.SQLite”或它的某一个依赖项。”的解决方法
未能加载文件或程序集“System.Data.SQLite”或它的某一个依赖项.试图加载格式不正确的程序. 说明: 执行当前 Web 请求期间,出现未经处理的异常.请检查堆栈 ...