http://acm.hdu.edu.cn/showproblem.php?pid=6097
题意:P,Q两点在一个半径为r,圆心在原点的圆内,且与圆心等距。求 圆上一点D,使PD+QD最小。输出最小的PD+QD。
题解:官方题解:http://bestcoder.hdu.edu.cn/blog/
很不幸不总是中垂线上的点取到最小值,考虑点在圆上的极端情况。
做P点关于圆的反演点P',OPD与ODP'相似,相似比是|OP| : r。
Q点同理。
极小化PD+QD可以转化为极小化P'D+Q'D。
当P'Q'与圆有交点时,答案为两点距离,否则最优值在中垂线上取到。
时间复杂度 O(1)
根据圆的反演,将P,Q与d的距离转化为反演点与d的距离。代码实现有点杂乱,主要就是分两种情况:1、反演点连线与圆没有交点,D取P,Q中垂线,2、取反演点连线与圆的交点。代码:
#include<bits/stdc++.h>
#define debug(a) cout<<a<<endl
#define mem(a,b) memset(a,b,sizeof(a))
#define LL long long
#define lson l,mid,root<<1
#define rson mid+1,r,root<<1|1
#define MIN_INT (-2147483647-1)
#define MAX_INT 2147483647
#define MAX_LL 9223372036854775807i64
#define MIN_LL (-9223372036854775807i64-1)
using namespace std;
const int N = 100000 + 5;
const int mod = 1000000000 + 7;
const double eps = 1e-6;
struct node{
double x,y;
};
int main(){
int t;
node p,q;
node P,Q,DD;
double r,k,d;
scanf("%d",&t);
while(t--){
scanf("%lf",&r);
scanf("%lf%lf%lf%lf",&p.x,&p.y,&q.x,&q.y);
d=sqrt(p.x*p.x+p.y*p.y);
k=r*r/(d*d);
P.x=p.x*k,P.y=p.y*k,Q.x=q.x*k,Q.y=q.y*k;
DD.x=(P.x+Q.x)/2,DD.y=(P.y+Q.y)/2;
double dd=sqrt(DD.x*DD.x+DD.y*DD.y);
if(dd<=r){//取反演点连线与圆的交点
double len=sqrt((P.x-Q.x)*(P.x-Q.x)+(P.y-Q.y)*(P.y-Q.y));
printf("%.7f\n",len*d/r);
}
else{//反演点连线与圆没有交点,D取P,Q中垂线
node dd;
dd.x=(p.x+q.x)/2,dd.y=(p.y+q.y)/2;
d=sqrt((dd.x-p.x)*(dd.x-p.x)+(dd.y-p.y)*(dd.y-p.y));
double len=r-sqrt(dd.x*dd.x+dd.y*dd.y);
printf("%.7f\n",2*sqrt(len*len+d*d));
}
}
return 0;
}