zoukankan      html  css  js  c++  java
  • joj 2536

     计算几何求解最小外接矩形的问题,这道题是暑假训练时做的一道题,当时不会做,选拔赛在即,突然想看看自己进步了多少,便做了一下,想法很容易就有了(说明进步了。。),但总是wa,逼我买了前辈的代码(178point啊),但还是找不到哪错了,最后终于发现最后四个数的初始化有问题,天,我已经初始化为整形最大值了(数据为double类型),还是不够大,天知道他给的测试数据有多变态。。。。话说马上选拔赛了,总感觉心里很不安,哎,尽人事,听天命吧。

    #include<math.h>
    #include<iostream>
    #include<stdio.h>
    using namespace std;
    const double PI = 3.14159265;
    struct POINT
    {
         double x;
         double y;
    };
    POINT bian[100];
    POINT rotate(POINT o,double alpha,POINT p)
    {
       POINT tp;
       p.x -=o.x;
       p.y -=o.y;
       tp.x=p.x*cos(alpha) - p.y*sin(alpha)+o.x;
       tp.y=p.y*cos(alpha) + p.x*sin(alpha)+o.y;
       return tp;
    }
    POINT  circle_center(POINT  pt[3]) 

       double  x1,x2,x3,y1,y2,y3; 
       double  x  =  0; 
       double  y  =  0; 
       x1  =  pt[0].x; 
       x2  =  pt[1].x; 
       x3  =  pt[2].x; 
       y1  =  pt[0].y; 
       y2  =  pt[1].y; 
       y3  =  pt[2].y; 
       x=((y2-y1)*(y3*y3-y1*y1+x3*x3-x1*x1)-(y3-y1)*(y2*y2-y1*y1+x2*x2-x1*x1))/(2*(x3-x1)*(y2-y1)-2*((x2-x1)*(y3-y1)));
       y=((x2-x1)*(x3*x3-x1*x1+y3*y3-y1*y1)-(x3-x1)*(x2*x2-x1*x1+y2*y2-y1*y1))/(2*(y3-y1)*(x2-x1)-2*((y2-y1)*(x3-x1)));
       POINT center;
       center.x  =  x  ; 
       center.y  =  y  ; 
       return center;
    }
    int main()
    {
        int n,i,iii=0;
        while(scanf("%d",&n),n)
        {  
            iii++;
            POINT p[3];
            for(i=0;i<3;i++)
                scanf("%lf%lf",&p[i].x,&p[i].y);
            POINT center=circle_center(p);
            bian[0]=p[0];
            double m=(2*PI)/n;
            for(i=1;i<n;i++)
                bian[i]=rotate(center,m,bian[i-1]);
            double xmax=bian[0].x,ymax=bian[0].y,xmin=bian[0].x,ymin=bian[0].y;
            for(i=0;i<n;i++)
            {
                if(xmax<bian[i].x)
                   xmax=bian[i].x;
                if(ymax<bian[i].y)
                   ymax=bian[i].y;
                if(xmin>bian[i].x)
                   xmin=bian[i].x;
                if(ymin>bian[i].y)
                   ymin=bian[i].y;
            }
            double area=(ymax-ymin)*(xmax-xmin);
            printf("Polygon %d: %.2lf\n",iii,area);
        }
        return 0;
    }

  • 相关阅读:
    phalcon之视图缓存
    Java NIO框架Netty教程(一) – Hello Netty
    setsockopt的作用
    支持向量机通俗导论(理解SVM的三层境地)
    quartz中的corn表达式(转)
    Applet 数字签名技术全然攻略
    SJTU 3001. 二哥的幸运
    OGRE之跳出漫长的编译等待
    VB.NET 数组的定义 动态使用 多维数组
    【Python】用Python的“结巴”模块进行分词
  • 原文地址:https://www.cnblogs.com/dchipnau/p/4985960.html
Copyright © 2011-2022 走看看