给定一个三角形求垂心
跟前一篇求圆心的一样。
因为发现了我之前的做法都比较二，于是乎想把整个模板都改改。
我之前的模板的做法仍就很麻烦
我先一点和该点的对边求点斜式。然后两条直线求交点。
想法简单，但是代码长度比较长，
POJ1673
point orthocenter(point p0,point p1,point p2){
double a1,b1,a2,b2,c1,c2;
a1 = p2.x-p1.x; b1=p2.y-p1.y;c1 = 0;
a2 = p2.x-p0.x; b2=p2.y-p0.y;c2 = (p1.x-p0.x)*a2+(p1.y-p0.y)*b2;
double d = a1 * b2 - a2 * b1;
ret.x = p0.x + (c1 * b2 - c2 * b1) / d;
ret.y = p0.y + (a1 * c2 - a2 * c1) / d;
return ret;
}
读入三个点p0,p1,p2返回垂心
根据垂直向量点成为0
(O.x-p0.x)*(p2.x-p1.x)+(O.y-p0.y)*(p2.y-p1.y) = 0; (1)
(O.x-p1.x)*(p2.x-p0.x)+(O.y-p1.y)*(p2.y-p0.y) = 0; (2)
令O.x = x+p0.x;O.y=y+p0.y带入(1)(2)得到
x*(p2.x-p1.x)+y*(p2.y-p1.y) = 0;
x*(p2.x-p0.x)+y*(p2.y-p0.y) = (p2.x-p0.x)*(p1.x-p0.x)+(p2.y-p0.y)*(p1.y-p0.y);
解得(x,y)
O = (x,y)+p0;
#include<iostream>
#include<cstdio>
using namespace std;
struct point{
double x,y;
} p0,p1,p2,ret;
int Test;
point orthocenter(point p0,point p1,point p2){
double a1,b1,a2,b2,c1,c2;
a1 = p2.x-p1.x; b1=p2.y-p1.y;c1 = 0;
a2 = p2.x-p0.x; b2=p2.y-p0.y;c2 = (p1.x-p0.x)*a2+(p1.y-p0.y)*b2;
double d = a1 * b2 - a2 * b1;
ret.x = p0.x + (c1 * b2 - c2 * b1) / d;
ret.y = p0.y + (a1 * c2 - a2 * c1) / d;
return ret;
}
int main(){
freopen("in.txt","r",stdin);
freopen("out.txt","w",stdout);
cin >> Test;
for (int ii = 1;ii<=Test;++ii){
cin>>p0.x>>p0.y>>p1.x>>p1.y>>p2.x>>p2.y;
ret = orthocenter(p0,p1,p2);
printf("%.4lf %.4lf\n",ret.x,ret.y);
}
return 0;
}