poj 1673

      之前看了题目分类，才事先知道是要求三角形的垂心，不然真的很难看出是垂心。这里有个比较好的方法来证出该点就是垂心：

 

http://hi.baidu.com/pojoflcc/blog/item/9452a4dade5c996cd1164e22.html

 

我用的方法是有向面积。

 

以下是代码：

 

1. #include<cstdio>
   #include<iostream>
   using namespace std;
   const double eps=1e-8;
   
2. struct point
   {
    double x,y;
   };
4. double det(double x1,double y1,double x2,double y2)
   {
    return(x1*y2-x2*y1);
   }
6. int main()
   {
7.  int T;
    scanf("%d",&T);
    while(T--)
    {
     point a,b,c;
     scanf("%lf%lf%lf%lf%lf%lf",&a.x,&a.y,&b.x,&b.y,&c.x,&c.y);
     double a1,a2,b1,b2,c1,c2;
     double x,y;
     b1=b.y-a.y;
     a1=b.x-a.x;
     c1=a.y*c.y-b.y*c.y+a.x*c.x-b.x*c.x;
     b2=c.y-a.y;
     a2=c.x-a.x;
     c2=a.y*b.y-b.y*c.y+a.x*b.x-b.x*c.x;
     x=det(c2,b2,c1,b1)/det(a1,b1,a2,b2)+eps;
     y=det(c1,a1,c2,a2)/det(a1,b1,a2,b2)+eps;
     printf("%.4lf %.4lf/n",x,y);
    }
    return 0;
   }