SGU 120 Arhipelago(计算几何)

原创 2013年12月02日 12:46:32

题意:给出一个正N边形的两个顶点,求正N边形的所有顶点。

思路:因为基本没写过几何,所以写得十分蛋疼,改了好久才过……这题主要是求中心点,只要求出中心点,那么通过旋转,就可以求出所有的点了。首先求出给出的两个顶点与中心点的夹角d=abs(n1-n2)*(2π/N)。然后把这条线转到中心上,然后就可以搞了……


代码:


#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<map>
#include<queue>
#include<set>
#include<stack>
#include<cmath>
#include<vector>
#define inf 0x3f3f3f3f
#define Inf 0x3FFFFFFFFFFFFFFFLL
#define eps 1e-9
#define pi acos(-1.0)
using namespace std;
typedef long long ll;
const int maxn=200+10;
struct Point
{
    double x,y;
    Point (double x=0,double y=0) :x(x),y(y){}
};
typedef Point Vector;
Point pt[maxn];
int N,N1,N2;
int dcmp(double x)
{
    if(abs(x)<eps) return 0;
    return x>0?1:-1;
}
double ocmp(double x)
{
    if(abs(x)<eps) return 0;
    return x;
}
double Len(Point a,Point b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
//x1=xcos(a)+ysin(a)
//y1=ycos(a)-xsin(a)
Vector rotate(Vector vc,double a)
{
    Vector res;
    res.x=vc.x*cos(a)+vc.y*sin(a);
    res.y=vc.y*cos(a)-vc.x*sin(a);
    return res;
}
Point getcp(Point p1,Point p2,double d)
{
    Point cp;
    int m=abs(N1-N2);
    double a=d*m,b;
    Point mid;
    mid.x=(p1.x+p2.x)/2.0;
    mid.y=(p1.y+p2.y)/2.0;
    double tmp=pi-a;
    if(dcmp(tmp)>=0) b=tmp/2.0;
    else b=(-tmp)/2.0;
    double D=Len(p1,p2)/(2.0*cos(b));
    if((N%2==0)&&(m==N/2)) return mid;
    Vector vc=Vector(p2.x-p1.x,p2.y-p1.y);
    if(m>N/2) b=2.0*pi-b;
    vc=rotate(vc,b);
    cp.x=p1.x+vc.x*D/Len(p1,p2);
    cp.y=p1.y+vc.y*D/Len(p1,p2);
    return cp;
}
int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    Point p1,p2,cp;
    cin>>N>>N1>>N2;
    cin>>p1.x>>p1.y;
    cin>>p2.x>>p2.y;
    if(N1>N2) {swap(N1,N2);swap(p1,p2);}
    double d=2.0*pi/N;
    cp=getcp(p1,p2,d);
    int id=N1;
    pt[N1]=p1;
    Vector vc=Vector(p1.x-cp.x,p1.y-cp.y);
    Vector vt;
    for(int i=1;i<N;++i)
    {
        id++;
        if(id==N+1) id=1;
        vt=rotate(vc,d*i);
        pt[id].x=cp.x+vt.x;
        pt[id].y=cp.y+vt.y;
    }
    pt[N2]=p2;
    for(int i=1;i<=N;++i)
      printf("%.6lf %.6lf\n",ocmp(pt[i].x),ocmp(pt[i].y));
    return 0;
}

相关文章推荐

SGU 120. Archipelago 计算几何 正n边形个顶点坐标

题目链接点这儿 QAQ这题做的真痛苦。。。 题意就是

sgu265:Wizards(计算几何)

题目大意:       ~~~~~~一个空间直角坐标系内,对一个点有三种操作:       ~~~~~~1.1.平移(x,y,z)⇒(x+a,y+b,z+c)(x,y,z)\Rightarrow ...

SGU 129 Inheritance(计算几何)

题意:给出n个点围成的多边形,求m个线段在多边形区域内的长度。 思路:这题又把我坑了一天Orz……果然不适合做几何题么。。。首先根据题中给出的条件,由于从多边形中的一个点到多边形中的另一个点不经过多...
  • qian99
  • qian99
  • 2014年01月11日 22:59
  • 787

sgu110-计算几何

110. Dungeon time limit per test: 0.25 sec.  memory limit per test: 4096 KB The mission of sp...

sgu198:Get Out!(计算几何+spfa判负环)

题意: 平面上有一堆圆,再给出一个圆P,问P能否在其他圆的缝隙中穿过到达无穷远处? 分析: 把每一个圆加上圆P的半径,圆P就成了一个点。 我们把每两个相交的圆之间连两条有向...

sgu 110 Dungeon 三维计算几何

题意: 思路:

SGU 124. Broken line 计算几何 判断点是否在多边形内

题目链接点这儿

SGU110 计算几何 Computational Geometry

题意:三维空间N个球体,给一个光束的方向,求它在这些球体间的反射情况(球面反射)。 Problem: N spheres are in a 3D-space. Give you a directio...

SGU 253 计算几何 判定点是否在凸包内

/******************************************************************************* 心血来潮上SGU敲了道计算几何~...

sgu120

题目 120. Archipelagotime limit per test: 0.25 sec. memory limit per test: 4096 KBArchipelago Ber-Isla...
  • zxozxo4
  • zxozxo4
  • 2014年12月09日 22:08
  • 271
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:SGU 120 Arhipelago(计算几何)
举报原因:
原因补充:

(最多只允许输入30个字)