UVA 10609 数学几何递推

最新推荐文章于 2024-07-15 22:10:26 发布

fwm_94

最新推荐文章于 2024-07-15 22:10:26 发布

阅读量593

点赞数

分类专栏： UVA 文章标签： uva 算法 algorithm

本文链接：https://blog.csdn.net/u012287002/article/details/38526715

版权

UVA 专栏收录该内容

1 篇文章 0 订阅

订阅专栏

用题目给的方法把一个线段分形，最后输出每个点的坐标。

大体思路很简单，就是找到两个4等分点，然后递归处理每一段。直接用解析几何的方法推导会比较繁杂，这里使用两个公式。用定比分点公式求出2个4等分点，然后就是已知两点，求等边三角形的另一点。平移坐标系，用旋转矩阵把一个点转60度（控制旋转方向，这样另一个点就有唯一解），再把坐标系平移回去就能求出新点的坐标。

#include <cstdio>
#include <cmath>
#include <algorithm>
#include <vector>
#define eps 1e-8
using namespace std;
const double PI=acos(-1.0);
int sgn(double x)
{
	if(fabs(x)<eps)
	   	return 0;
	else if(x>0)
		return 1;
	else
		return -1;
}
struct node
{
	double x,y;
	bool operator <(const node &rhs)const
	{
		if(sgn(x-rhs.x)!=0)
			return sgn(x-rhs.x)<0;
		else
			return sgn(y-rhs.y)<0;
	}
};
double T;
vector<node>ans;
void solve(double x1,double y1,double x2,double y2,bool flag)
{
	node tmp;
	if(flag)
	{
		tmp.x=x1;tmp.y=y1;ans.push_back(tmp);
	}
	tmp.x=x2;tmp.y=y2;ans.push_back(tmp);
	double x3,y3,x4,y4,k,x,y,nx,ny,len;
	k=1.0/3;
	x3=(x1+k*x2)/(k+1);
	y3=(y1+k*y2)/(k+1);
	k=3;
	x4=(x1+k*x2)/(k+1);
	y4=(y1+k*y2)/(k+1);
	len=sqrt((x3-x4)*(x3-x4)+(y3-y4)*(y3-y4));
	if(sgn(len-T)<0)
		return;
	nx=x4-x3;ny=y4-y3;
	x=cos(PI*5/3)*nx-sin(PI*5/3)*ny;
	y=sin(PI*5/3)*nx+cos(PI*5/3)*ny;
	x+=x3;y+=y3;
	solve(x3,y3,x,y,1);
	solve(x,y,x4,y4,0);
	return;
}
int main()
{
	double x1,x2,y1,y2;int cot=1;
	while(scanf("%lf%lf%lf%lf%lf",&x1,&y1,&x2,&y2,&T)!=EOF)
	{
		if(sgn(T-1)<0)
			return 0;
		ans.clear();
		solve(x2,y2,x1,y1,1);
		sort(ans.begin(),ans.end());
		printf("Case %d:\n",cot++);
		int n=ans.size();
		printf("%d\n",n);
		for(int i=0;i<n;i++)
			printf("%.5f %.5f\n",ans[i].x,ans[i].y);
	}
}