Codeforces 2C - Commentator problem

最新推荐文章于 2019-12-04 11:07:17 发布

原创最新推荐文章于 2019-12-04 11:07:17 发布 · 360 阅读

0 ·

CC 4.0 BY-SA版权

CodeForces 专栏收录该内容

35 篇文章

订阅专栏

本文介绍了一种使用数学方法求解由三个已知圆组成的三角形内部的共轭点（即内切圆心）的算法。通过初始化内切圆心为三角形质心，并利用梯度下降法进行迭代优化，最终找到使得三个圆到该点的距离比等于其半径比的点。此算法在计算机图形学和几何优化中具有应用价值。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

可以用数学方法做

初试爬山

#include<stdio.h>
#include<cmath>

const int fx[4]={1,-1,0,0},fy[4]={0,0,1,-1};
inline double sqr(double x){return x*x;}
const double eps=1e-6;
double x[3],y[3],r[3],X,Y;

inline double judge(double X,double Y)
{
	static double Tan[3];
	for (int i=0;i<3;i++) Tan[i]=sqrt(sqr(X-x[i])+sqr(Y-y[i]))/r[i];
	double rt=0;
	for (int i=0;i<3;i++) rt+=sqr(Tan[i]-Tan[(i+1)%3]);
	return rt;
}

int main()
{
	for (int i=0;i<3;i++) scanf("%lf%lf%lf",x+i,y+i,r+i);
	X=(x[0]+x[1]+x[2])/3;
	Y=(y[0]+y[1]+y[2])/3;
	for (double st=1;st>eps;)
	{
		bool flag=1;
		for (int i=0;i<4;i++) if (judge(X+fx[i]*st,Y+fy[i]*st)<judge(X,Y))
			X+=fx[i]*st,Y+=fy[i]*st,flag=0;
		if (flag) st/=2;
	}
	if (judge(X,Y)<eps) printf("%.5lf %.5lf\n",X,Y);
}