题意:已知起点、圆、矩形,要求计算从起点开始,经过圆(和圆上任一点接触即可),到达矩形的路径的最短距离。(可以穿过园)。

分析:没什么好的方法,凭感觉圆上的每个点对应最短距离,应该是一个凸函数,用三分来解。不过应该是分成两部分,用两次三分来解。具体原因不明,通过实验只能得出三分必须是针对一个凸函数,很明显,从0~2*pi不是一个凸函数,但同理,0~pi,pi~2*pi也不一定是个凸函数。个人认为,网络上流传的[0,pi][pi,2*pi]这种分法,并不存在合理性。继续思考= =

    另外,直接暴力枚举也是能过的。

错误:分析圆上的点与矩形的位置关系写挫了,虽然现在还是研究不明白到底出了什么问题,应该是边界吧。

 double rec(double x,double y)
{
check();
if(x-xi<eps&&y-yj>eps)
return len(x,y,xi,yj);
else if(x-xi<eps&&y-yi<eps)
return len(x,y,xi,yi);
else if(x-xj>eps&&y-yj>eps)
return len(x,y,xj,yj);
else if(x-xj>eps&&y-yi<eps)
return len(x,y,xj,yj); if(x-xi>eps&&x-xj<eps)
if(y-yj>eps)
return fabs(y-yj);
else if(y-yi<eps)
return fabs(yi-y);
else if(y-yi>eps&&y-yj<eps)
if(x-xi<eps)
return fabs(xi-x);
else if(x-xj>eps)
return fabs(x-xj);
}

错误的写法

 #include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<algorithm>
using namespace std; const double pi=acos(-1.0);
const double eps=1e-; double sx,sy;
double rx,ry,r;
double xi,yi,xj,yj; double degree()
{
return atan((ry-sy)/(rx-sx));
} double len(double x,double y,double a,double b)
{
return sqrt((x-a)*(x-a)+(y-b)*(y-b));
} void change(double& x,double &y)
{
double p;
p=x;
x=y;
y=p;
} void check()
{
if(xi>xj)
change(xi,xj);
if(yi>yj)
change(yi,yj);
} double rec(double x,double y)
{
double lenx=,leny=;
check();
if(x<xi)
lenx=xi-x;
else if(x>xj)
lenx=x-xj;
if(y<yi)
leny=yi-y;
else if(y>yj)
leny=y-yj;
return sqrt(lenx*lenx+leny*leny);
} double f(double m)
{
double dx=rx+r*cos(m);
double dy=ry+r*sin(m); double len1=len(dx,dy,sx,sy);
double len2=rec(dx,dy); return len1+len2;
} int main()
{
int n;
while(~scanf("%lf%lf",&sx,&sy))
{
if(fabs(sx)<eps&&fabs(sy)<eps)
return ;
scanf("%lf%lf%lf",&rx,&ry,&r);
scanf("%lf%lf%lf%lf",&xi,&yi,&xj,&yj); double l=degree(),r=l+pi;
while(fabs(r-l)>eps)
{
double m1=l+(r-l)/;
double m2=r-(r-l)/;
if(f(m1)<f(m2))
r=m2;
else
l=m1;
}
double l1=l;
l=degree()+pi;r=l+pi;
while(fabs(r-l)>eps)
{
double m1=l+(r-l)/;
double m2=r-(r-l)/;
if(f(m1)<f(m2))
r=m2;
else
l=m1;
}
double l2=l;
printf("%.2f\n",min(f(l1),f(l2)));
}
return ;
}

暴力枚举,奇葩的精度。不过相比较而言,这个算法根据有可信度。

 #include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<algorithm>
using namespace std; const double INF=1e10;
const double pi=acos(-1.0);
const double eps=1e-; double sx,sy;
double rx,ry,r;
double xi,yi,xj,yj; double len(double x,double y,double a,double b)
{
return sqrt((x-a)*(x-a)+(y-b)*(y-b));
} void change(double& x,double &y)
{
double p;
p=x;
x=y;
y=p;
} void check()
{
if(xi>xj)
change(xi,xj);
if(yi>yj)
change(yi,yj);
} double rec(double x,double y)
{
double lenx=,leny=;
check();
if(x<xi)
lenx=xi-x;
else if(x>xj)
lenx=x-xj;
if(y<yi)
leny=yi-y;
else if(y>yj)
leny=y-yj;
return sqrt(lenx*lenx+leny*leny);
} double f(double m)
{
double dx=rx+r*cos(m);
double dy=ry+r*sin(m); double len1=len(dx,dy,sx,sy);
double len2=rec(dx,dy); return len1+len2;
} int main()
{
int n;
while(~scanf("%lf%lf",&sx,&sy))
{
if(fabs(sx)<eps&&fabs(sy)<eps)
return ;
scanf("%lf%lf%lf",&rx,&ry,&r);
scanf("%lf%lf%lf%lf",&xi,&yi,&xj,&yj); double m=INF;
for(double i=;i<=;i+=0.01)
{
m=min(m,f(i/*pi));
}
printf("%.2f\n",m);
}
return ;
}

hdu 4454 Stealing a Cake(计算几何:最短距离、枚举/三分)的更多相关文章

  1. hdu 4454 Stealing a Cake(三分之二)

    pid=4454" target="_blank" style="">题目链接:hdu 4454 Stealing a Cake 题目大意:给定 ...

  2. hdu 4454 Stealing a Cake (三分)

    Stealing a Cake Time Limit: 5000/2000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) ...

  3. HDU 4454 Stealing a Cake(枚举角度)

    题目链接 去年杭州现场赛的神题..枚举角度..精度也不用注意.. #include <iostream> #include <cstdio> #include <cstr ...

  4. HDU 4454 Stealing a Cake --枚举

    题意: 给一个点,一个圆,一个矩形, 求一条折线,从点出发,到圆,再到矩形的最短距离. 解法: 因为答案要求输出两位小数即可,精确度要求不是很高,于是可以试着爆一发,暴力角度得到圆上的点,然后求个距离 ...

  5. hdu 4454 Stealing a Cake 三分法

    很容易想到三分法求解,不过要分别在0-pi,pi-2pi进行三分. 另外也可以直接暴力枚举…… 代码如下: #include<iostream> #include<stdio.h&g ...

  6. hdu 4454 Stealing a Cake

    简单的计算几何: 可以把0-2*pi分成几千份,然后找出最小的: 也可以用三分: #include<cstdio> #include<cmath> #include<al ...

  7. hdu 4454 Stealing a Cake(三分法)

    给定一个起始点,一个矩形,一个圆,三者互不相交.求从起始点->圆->矩形的最短距离. 自己画一画就知道距离和会是凹函数,不过不是一个凹函数.按与水平向量夹角为圆心角求圆上某点坐标,[0, ...

  8. HDU 4454 - Stealing a Cake(三分)

    我比较快速的想到了三分,但是我是从0到2*pi区间进行三分,并且漏了一种点到边距离的情况,一直WA了好几次 后来画了下图才发现,0到2*pi区间内是有两个极值的,每个半圆存在一个极值 以下是代码 #i ...

  9. HDU - 4454: Stealing a Cake (圆上三分)

    pro:给定一个蛋糕,一个矩阵房子,一只蚂蚁.最开始三者两两相离,问蚂蚁触摸到蛋糕后再触摸矩阵的最短距离.结果保留两位小数,坐标的绝对值<1e4: sol:由于坐标不大,而且精度要求不高,不难想 ...

随机推荐

  1. Codeforces 897 B.Chtholly's request-思维题(处理前一半)

      B. Chtholly's request   time limit per test 2 seconds memory limit per test 256 megabytes input st ...

  2. 小程序收集formid跳转后收集不到

    为了突破微信小程序模板消息限制,需收集到足够的推送码,即每次提交表单时获取到的formId.一个formId代表着开发者有向当前用户推送模板消息的一次权限.当表单组件中的属性report-submit ...

  3. gedit插件配置

    Ubuntu用户 sudo apt-get install gedit-plugins Fedora用户 yum install gedit-plugins 使用gEdit搭配terminal来写程序 ...

  4. Linux系统入门命令100条 转

    https://www.howtoforge.com/linux-commands/ 2017-04-27 RiboseYim 睿哥杂货铺 Author : Himanshu Arora 原文:htt ...

  5. 打包工具 使用帮助 inno setup

    http://wenku.baidu.com/link?url=0VRJ8n9am1KgVAAqwz-AU1htXamo7Vh0d4QIdGG6_LcTrZBdb7lRim8Jx6M8KaLJDQm1 ...

  6. CM3大礼包

  7. 批量删除Redis中的key

    bin/redis-cli -h 192.168.46.151 -p 6379 keys "rulelist*" | xargs bin/redis-cli  -h 192.168 ...

  8. Elasticsearch教程(六) elasticsearch Client创建

    Elasticsearch  创建Client有几种方式. 首先在 Elasticsearch  的配置文件 elasticsearch.yml中.定义cluster.name.如下: cluster ...

  9. [ElasticSearch]Java API 之 词条查询(Term Level Query)

    1. 词条查询(Term Query)  词条查询是ElasticSearch的一个简单查询.它仅匹配在给定字段中含有该词条的文档,而且是确切的.未经分析的词条.term 查询 会查找我们设定的准确值 ...

  10. angular controller的一些用法

    最近公司的项目是es6+angular.其中的代码格式还在逐步摸索中.感谢今天同事每天帮我解惑. 今天简单梳理一下controller的一些用法 之前看书所熟知的都是 这是最普通的一种 //html ...