[计算几何]POJ 1266 三角形的外接圆 圆的参数方程

http://acm.pku.edu.cn/JudgeOnline/problem?id=1266

题目大意是给定一个圆弧上的3个点（先给2个端点，再给定圆弧中间的一个点），要求出覆盖此圆弧的最小矩形（矩形的4个角坐标必须为整数）

先由给定的3个点可以确定这个圆弧所在圆的圆心的半径，用三角形外接圆圆心公式（好长的公式啊）：

x=mx1+(0.5*((mx2-mx1)*(mx2-mx1)+(my2-my1)*(my2-my1))*(my3-my1)

-0.5*((mx3-mx1)*(mx3-mx1)+(my3-my1)*(my3-my1))*(my2-my1))

/((mx2-mx1)*(my3-my1)-(mx3-mx1)*(my2-my1));

y=my1+(0.5*((mx3-mx1)*(mx3-mx1)+(my3-my1)*(my3-my1))*(mx2-mx1)

-0.5*((mx2-mx1)*(mx2-mx1)+(my2-my1)*(my2-my1))*(mx3-mx1))

/((mx2-mx1)*(my3-my1)-(mx3-mx1)*(my2-my1));

下面求出3个点相对于圆心的偏转角度，有了这个角度之后，我们就可以判断出这段弧是一段优弧还是一段劣弧，进而给定一个角度后，我们就可以判断这个角度的点是否在这段弧上。

然后根据圆的性质很容易求得包围这段圆弧的minx,maxx,miny,maxy

PS:又是可恶的浮点误差WA了我N次，利用floor,ceil函数的时候一定要加上误差控制!!!

floor(x+EPS),ceil(x-EPS)

#include <iostream>

#include <algorithm>

#include <cmath>

using namespace std;

#define EPS 1.0e-6

#define pi 3.141592653589793232846

double mx1,my1,mx2,my2,mx3,my3;

bool flag;

double a1,a2,a3;

inline bool zero(double r) {

return fabs(r)<EPS;

}

void getcircle(double mx1,double my1,double mx2,double my2,

double mx3,double my3,double &x,double &y) {

x=mx1+(0.5*((mx2-mx1)*(mx2-mx1)+(my2-my1)*(my2-my1))*(my3-my1)

-0.5*((mx3-mx1)*(mx3-mx1)+(my3-my1)*(my3-my1))*(my2-my1))

/((mx2-mx1)*(my3-my1)-(mx3-mx1)*(my2-my1));

y=my1+(0.5*((mx3-mx1)*(mx3-mx1)+(my3-my1)*(my3-my1))*(mx2-mx1)

-0.5*((mx2-mx1)*(mx2-mx1)+(my2-my1)*(my2-my1))*(mx3-mx1))

/((mx2-mx1)*(my3-my1)-(mx3-mx1)*(my2-my1));

}

double angle(double x,double y,double mx1,double my1) {

double dx=mx1-x,dy=my1-y;

double theta;

if(zero(dx)) {

if(dy>0) return pi*0.5;

return pi*1.5;

}else{

theta=atan(dy/dx);

if(dx<0) theta+=pi;

else if(dy<0) theta+=2*pi;

}

return theta;

}

bool onarc(double angle) {

return ((angle-a1)*(angle-a2)<=0)==flag;

}

int main() {

cin>>mx1>>my1>>mx2>>my2>>mx3>>my3;

double x,y;

getcircle(mx1,my1,mx2,my2,mx3,my3,x,y);

double r=sqrt((x-mx1)*(x-mx1)+(y-my1)*(y-my1));

a1=angle(x,y,mx1,my1);

a2=angle(x,y,mx2,my2);

a3=angle(x,y,mx3,my3);

flag=(a3-a1)*(a3-a2)<=0;

int minx,maxx,miny,maxy;

if(onarc(0)) {

maxx=ceil(x+r-EPS);

}else{

maxx=ceil(max(mx1,mx2)-EPS);

}

if(onarc(pi)) {

minx=floor(x-r+EPS);

}else{

minx=floor(min(mx1,mx2)+EPS);

}

if(onarc(pi*0.5)) {

maxy=ceil(y+r-EPS);

}else{

maxy=ceil(max(my1,my2)-EPS);

}

if(onarc(pi*1.5)) {

miny=floor(y-r+EPS);

}else{

miny=floor(min(my1,my2)+EPS);

}

cout<<(((long long)(maxx-minx))*(maxy-miny))<<endl;

return 0;

}