hdu6097 Mindis 2017多校第六场1002 圆的反演

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;
}