poj 3845 Fractal

最新推荐文章于 2017-12-22 08:45:55 发布

原创最新推荐文章于 2017-12-22 08:45:55 发布 · 501 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

学生时代刷题归档专栏收录该内容

146 篇文章

订阅专栏

本文介绍了一种基于二维空间的几何变换算法，通过旋转和平移操作实现路径上的点坐标计算。利用结构化的点对象进行向量运算，包括点乘、叉乘等，并通过角度获取与比例缩放实现了复杂路径上特定比例点的定位。

妈蛋。。。WA到死啊。。最后才发现原来是自己把一个地方的下表写错了。。去年买了个包的。。。。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<math.h>
#include<string>
#include<algorithm>
#define max_n 1000000
using namespace std;
const double eps=1e-12;
const double pi=acos(-1.0);
struct P
{
	double x,y;
	P(){}
	P(double x,double y):x(x),y(y){
	}
	P operator +(P p){
		return P(x+p.x,y+p.y);
	}
	P operator -(P p){
		return P(x-p.x,y-p.y);
	}
	P operator *(double d){
		return P(x*d,y*d);
	}
	double det(P p){
		return x*p.y-y*p.x;
	}
	double dot(P p){
		return x*p.x+y*p.y;
	}
};
double dis(P a,P b)
{
	double c=(a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
	return sqrt(c);
}
double dblcmp(double k)
{
	if(fabs(k)<eps)
		return 0;
	return k>0?1:-1;
}
P p[max_n];
double g[max_n];
int t,n,d;
double f;
double get_ang(P a,P b,P c)
{
	double A=(b-a).det(c-a);
	if(dblcmp(A)==0)
		A=0;
	double B=(b-a).dot(c-a);
	if(dblcmp(B)==0)
		B=0;
	return atan2(A,B);
}
P get_rotation(P a,P b,double angle,double s)
{
	double cosl=cos(-angle)*s;
	double sinl=sin(-angle)*s;
	b.x -= a.x; b.y -= a.y;  
    P c;  
    c.x = b.x * cosl - b.y * sinl + a.x;  
    c.y = b.x * sinl + b.y * cosl + a.y;  
    return c;  
}
void solve()
{
	int t;
	double len=0,res=0;
	double dx,dy;
	for(int i=0;i<n-1;i++)
		len+=dis(p[i],p[i+1]);
	g[0]=0;
	for(int i=1;i<n;i++)
	{
		res+=dis(p[i],p[i-1]);
		g[i]=res/len;
	}
	for(int i=1;i<d;i++)
	{
		t=0;
		while(t<n && dblcmp(g[t]-f)<0)t++;
		if(dblcmp(g[t]-f)==0)
		{
			printf("(%.7lf,%.7lf)\n",p[t].x,p[t].y);
			return ;
		}
		P tmp;
		dx=p[t-1].x-p[0].x;
		dy=p[t-1].y-p[0].y;
		tmp.x=p[n-1].x+dx;
		tmp.y=p[n-1].y+dy;
		double angle=get_ang(p[t-1],p[t],tmp);
		double s=dis(p[t-1],p[t])/dis(p[0],p[n-1]);
		for(int j=1;j<n;j++)
			p[j]=get_rotation(p[0],p[j],angle,s);
		for(int j=0;j<n;j++)
			p[j].x+=dx,p[j].y+=dy;
		/*for(int j=0;j<n;j++)
			printf("%lf %lf\n",p[j].x,p[j].y);
		system("pause");*/
		f=(f-g[t-1])/(g[t]-g[t-1]);
	}
	t=0;
	while(t<n && dblcmp(g[t]-f)<0)t++;
	if(dblcmp(g[t]-f)==0)
	{
		printf("(%.7lf,%.7lf)\n",p[t].x,p[t].y);
		return ;
	}
	f=(f-g[t-1])/(g[t]-g[t-1]);
	dx=p[t-1].x+(p[t].x-p[t-1].x)*f;
	dy=p[t-1].y+(p[t].y-p[t-1].y)*f;
	printf("(%.7lf,%.7lf)\n",dx,dy);
}
int main()
{
	scanf("%d",&t);
	while(t--)
	{
		scanf("%d",&n);
		for(int i=0;i<n;i++)
			scanf("%lf%lf",&p[i].x,&p[i].y);
		scanf("%d%lf",&d,&f);
		solve();
		system("pause");
	}
	return 0;
}