hdu4720-(Naive and Silly Muggles)

    补题中，这道题数据有点弱。

    题目大意：给一个三角形的三个顶点的坐标，再给出一个点的坐标(x,y)，问最后这个点(x,y)是否落在由三角形构成的最小的圆里面（包括圆的边）。题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4720。

    思路：这道题队友写的，后来我补题的时候才发现他的思路不太对，但却ac了，他直接求三角形的重心算出圆的半径，这样也行。这道题应该求三角形的外心，可以由三个坐标直接得外心的坐标，公式可以自己推。

    推理过程如下：设外心坐标为(a,b)，三角形三个坐标分别为(x1,y1)，(x2,y2)，(x3,y3)，则

    (a-x1)*(a-x1)+(b-y1)*(b-y1)=(a-x2)*(a-x2)+(b-y2)*(b-y2)

    (a-x1)*(a-x1)+(b-y1)*(b-y1)=(a-x3)*(a-x3)+(b-y3)*(b-y3)

    两条式子联立求出a、b。

    外心求出的是外切圆，对于钝角三角形需要另外算，因为钝角三角形的外切圆不是最小圆。对于钝角三角形而言，最小圆的圆心是对边的中点。还好输入样例错误了，不然没发现。

以下是ac代码：

#include<iostream>
#include<sstream>
#include<malloc.h>
#include<map>
#include<list>
#include<bitset>
#include<cctype>
#include<iomanip>
#include<utility>
#include<stdlib.h>
#include<functional>
#include<stdio.h>
#include<stack>
#include<queue>
#include<deque>
#include<algorithm>
#include<string>
#include<string.h>
#include<math.h>
#include<set>
#include<vector>
using namespace std;

int main(){
    int t;
    scanf("%d",&t);
    for(int i=1;i<=t;i++){
        double x1,x2,x3,y1,y2,y3;
        double x,y;
        double a,b,r;
        scanf("%lf%lf%lf%lf%lf%lf%lf%lf",&x1,&y1,&x2,&y2,&x3,&y3,&x,&y);
        if((x2-x1)*(x3-x1)+(y2-y1)*(y3-y1)<0){
            a=(x2+x3)/2;
            b=(y2+y3)/2;
            r=(a-x2)*(a-x2)+(b-y2)*(b-y2);
        }
        else if((x1-x2)*(x3-x2)+(y1-y2)*(y3-y2)<0){
            a=(x1+x3)/2;
            b=(y1+y3)/2;
            r=(a-x1)*(a-x1)+(b-y1)*(b-y1);
        }
        else if((x1-x3)*(x2-x3)+(y1-y3)*(y2-y3)<0){
            a=(x1+x2)/2;
            b=(y1+y2)/2;
            r=(a-x1)*(a-x1)+(b-y1)*(b-y1);
        }
        else{
            a=((x1*x1+y1*y1-x2*x2-y2*y2)*(y1-y3)-(x1*x1+y1*y1-x3*x3-y3*y3)*(y1-y2))/(2*((y1-y3)*(x1-x2)-(y1-y2)*(x1-x3)));
            b=((x1*x1+y1*y1-x2*x2-y2*y2)*(x1-x3)-(x1*x1+y1*y1-x3*x3-y3*y3)*(x1-x2))/(2*((x1-x3)*(y1-y2)-(x1-x2)*(y1-y3)));
            r=(a-x1)*(a-x1)+(b-y1)*(b-y1);
        }
        printf("Case #%d: ",i);
        if((a-x)*(a-x)+(b-y)*(b-y)<=r)
            printf("Danger\n");
        else
            printf("Safe\n");
    }
    return 0;
}