Finding the Radius for an Inserted Circle--2017 ACM-ICPC 亚洲区（南宁赛区）网络赛

1. Finding the Radius for an Inserted Circle

  比赛最后没时间写了，思路大概有，也不知道对不对，然而结束发现输入输出有问题，这个题过了很多人，都是用笛卡尔定理过的，表示第一次听说这个东西，我们的思路是二分，一直以为T组数据，每次输入R与K，然后输出K次答案。。。。

 二分思路：将两个圆心固定在x轴两边，利用等边三角形性质求出第三个圆心在y轴上的坐标，因为所有的点都是很有特点，很容易知道第一个内置圆的圆心即三个圆心所组成三角形的内心，开始思路是这样，因为如果知道了每个圆心，二分其半径即可，后来发现既然每个内置圆的处理方式都是一样，那么直接从第一个内置圆开始，因为求出了 y轴上的圆的圆心，那么第一个内置圆与这个圆的切点即y轴圆心减去R，然后二分半径，切点到内置圆的圆心距离就是半径，所以当前内置圆的圆心可以知道，然后判断到另外两个圆心的距离关系是否相切。具体上代码吧：

const double eps=1e-8;
const double PI=acos(-1.0);
const int INF=1e9+10;
const int MOD=1e9+7;
const int N=1e6+10;
#define quantity sqrt(3.0)
double R;
int k;
struct point
{
    double x,y;
} a,b,c,pre,now;
double dis(point a,point b)//求两点距离
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
bool judge(double mid)
{
    now.y=pre.y-mid;
    double dist=dis(now,a);
    if(dist-R>=mid) return true;
    return false;
}
void solve()
{
    double last=R;
    now.x=pre.x=0,pre.y=c.y-R;
    int flag=0;
    double ans;
    for(int i=1; i<=k; i++)
    {
        double l=0,r=last;
        ans=0;
        if(!flag)
        {
            while(l+eps<=r)
            {
                double mid=(l+r)/2.0;
                if(judge(mid)) l=mid+eps,ans=mid;
                else r=mid-eps;
            }
        }
        last=ans;
        pre.y-=2.0*ans;
        if((long long)last==0) flag=1;
    }
    printf("%d %lld\n",k,((long long)ans));
}
int main()
{
    int t;
    while(~scanf("%d",&t)&&t!=-1)
    {
        scanf("%lf",&R);
        for(int i=0; i<t; i++)
        {
            scanf("%d",&k);
            a.x=-R,a.y=0,b.x=R,b.y=0;
            c.x=0,c.y=quantity*R;
            solve();
        }
    }
    return 0;
}


//1
//152973.6
//4