模拟退火果然是一个非常高端的东西，思路神马的全然搞不懂啊~

题目大意：

给出三个圆，求一点到这三个圆的两切线的夹角相等。

解题思路：

对于这个题来说还是有多种思路的 。只是都搞不明确~~   /害羞脸

用模拟退火来解也是一件赌人品的事。由于退火的过程设计的不合理，WA妥妥的。

事实上我也是学了一点点。还不是太明确啊~~

以下是代码：

#include <set>

#include <map>

#include <queue>

#include <math.h>

#include <vector>

#include <string>

#include <stdio.h>

#include <string.h>

#include <stdlib.h>

#include <iostream>

#include <cctype>

#include <algorithm>

#define eps 1e-6

#define pi acos(-1.0)

#define inf 107374182

#define inf64 1152921504606846976

#define lc l,m,tr<<1

#define rc m + 1,r,tr<<1|1

#define zero(a) fabs(a)<eps

#define iabs(x)  ((x) > 0 ? (x) : -(x))

#define clear1(A, X, SIZE) memset(A, X, sizeof(A[0]) * (min(SIZE,sizeof(A))))

#define clearall(A, X) memset(A, X, sizeof(A))

#define memcopy1(A , X, SIZE) memcpy(A , X ,sizeof(X[0])*(SIZE))

#define memcopyall(A, X) memcpy(A , X ,sizeof(X))

#define max( x, y )  ( ((x) > (y)) ?

(x) : (y) )

#define min( x, y )  ( ((x) < (y)) ? (x) : (y) )

using namespace std;

int dir[4][2]={{1,0},{-1,0},{0,1},{0,-1}};

struct circle

{

    double x,y,r;

}cir[3];

double dis(double x,double y,double xx,double yy)

{

    return sqrt((x-xx)*(x-xx)+(y-yy)*(y-yy));

}

double f(double x,double y)

{

    double tmp[3];

    for(int i=0;i<3;i++)

        tmp[i]=dis(x,y,cir[i].x,cir[i].y)/cir[i].r; //视角一半的sin值

    double ans=0;

    for(int i=0;i<3;i++)

        ans+=(tmp[i]-tmp[(i+1)%3])*(tmp[i]-tmp[(i+1)%3]);

    return ans;

}

int main()

{

    double x=0,y=0;

    for(int i=0;i<3;i++)

    {

        scanf("%lf%lf%lf",&cir[i].x,&cir[i].y,&cir[i].r);

        x+=cir[i].x/3;

        y+=cir[i].y/3;

    }

    double step=2;

    while(step>eps)

    {

        double tmp=f(x,y);

        int tag=-1;

        for(int i=0;i<4;i++)

        {

            double cnt=f(x+dir[i][0]*step,y+dir[i][1]*step);

            if(cnt<tmp)

            {

                tmp=cnt;

                tag=i;

            }

        }

        if(tag==-1)

            step/=2;

        else

        {

            x=x+dir[tag][0]*step;

            y=y+dir[tag][1]*step;

        }

    }

    if(f(x,y)<eps)

        printf("%.5lf %.5lf\n",x,y);

    return 0;

}