Feng Shui POJ

题目链接：https://cn.vjudge.net/problem/POJ-3384

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include <iostream>#define eps 1e-10
using namespace std;
struct point{double x,y;
}points[110],p[110],q[110];
int n,curcnt,cCnt,m;
void guizheng(){for(int i=1;i<=n;++i){q[i]=points[n-i+1];}for(int i=1;i<=n;++i){points[i]=q[i];}
}
double juli(point p1, point p2){return (p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y);
}
void getline(point p1,point p2,double &a,double &b,double &c){a=p2.y-p1.y;b=p1.x-p2.x;c=p2.x*p1.y-p2.y*p1.x;
}
void init(){for(int i=1;i<=n;++i){p[i]=points[i];}p[n+1]=p[1];p[0]=p[n];cCnt=n;
}
point getinter(point p1,point p2,double a,double b,double c){double u=fabs(a*p1.x+b*p1.y+c);double v=fabs(a*p2.x+b*p2.y+c);point temp;temp.x=(v*p1.x+u*p2.x)/(u+v);temp.y=(v*p1.y+u*p2.y)/(u+v);return temp;
}
void cut(double a,double b,double c){curcnt=0;for(int i=1;i<=cCnt;++i){if(a*p[i].x+b*p[i].y+c>=0)q[++curcnt]=p[i];else{if(a*p[i-1].x+b*p[i-1].y+c>0){q[++curcnt]=getinter(p[i],p[i-1],a,b,c);}if(a*p[i+1].x+b*p[i+1].y+c>0)q[++curcnt]=getinter(p[i],p[i+1],a,b,c);}}for(int i=1;i<=curcnt;++i){p[i]=q[i];}p[curcnt+1]=q[1];p[0]=p[curcnt];cCnt=curcnt;
}
void solve(){points[n+1]=points[1];init();for(int i=1;i<=n;++i){point ta,tb,tc;tc.x=points[i+1].y-points[i].y;tc.y=points[i].x-points[i+1].x;double k=m/sqrt(tc.x*tc.x+tc.y*tc.y);tc.x*=k;tc.y*=k;ta.x=points[i].x+tc.x;ta.y=points[i].y+tc.y;tb.x=points[i+1].x+tc.x;tb.y=points[i+1].y+tc.y;double a,b,c;getline(ta,tb,a,b,c);cut(a,b,c);}double max=-1;int r1,r2;for(int i=1;i<=cCnt;++i){for(int j=i+1;j<=cCnt;++j){double temp= juli(p[i], p[j]);if(temp>max){max=temp;r1=i;r2=j;}}}printf("%.6lf %.6lf %.6lf %.6lf\n",p[r1].x,p[r1].y,p[r2].x,p[r2].y);
}
int main()
{int i,j,k;while(~scanf("%d%d",&n,&m)){for(i=1;i<=n;++i){scanf("%lf%lf",&points[i].x,&points[i].y);}solve();}return 0;
}

Feng Shui POJ - 3384相关推荐

Feng Shui POJ - 3384 [半平面交]
Feng Shui POJ - 3384 题意:n个顶点的凸包,放入2个半径为r的圆,可以重叠,要求面积最大,输出2个圆的圆心坐标(保留4位小数) 思路:找出圆心的可行域(内推r,求半平面交),再求核 ...
poj 3384 Feng Shui (Half Plane Intersection)
3384 -- Feng Shui 构造半平面交,然后求凸包上最远点对. 这题的题意是给出一个凸多边形区域,要求在其中放置两个半径为r的圆(不能超出凸多边形区域),要求求出两个圆心,使得多边形中没有被 ...
POJ - 3384 Feng Shui（半平面交）
链接 Feng Shui 题意将两个半径为 rrr 的圆放入一个多边形中,两个圆占据最大面积时圆心坐标是多少: 思路在多边形中放圆,最大圆的圆心一定在多边形内核中: 将多边形的每条边都内推 rrr ...
CodeForces - Feng Shui(半平面交)
题目链接:http://codeforces.com/gym/101650/attachments Time Limit: 2000MS Memory Limit: 65536K Descriptio ...
poj 3384 Feng Shui 半平面交
http://poj.org/problem?id=3384 给定一个多边形,在多边形内放有两个相同的圆,使两个圆尽可能多的覆盖多边形.输出最终两个圆心的位置. 最优的放置方法必定是圆内切于两条边,那 ...
POJ 3384 Feng Shui(半平面交)
题意 : 给你一个凸多边形,让你在其中找两个圆,使得圆的覆盖面积最大. 这个题目和 poj 3525 有点类似,那个题目是一个圆,想到两者的联系,可以发现两个圆覆盖面积最大就是重叠面积最小,怎样使得重 ...
POJ 3384 Feng Shui
坑爹的题,,就是先向内推进r,然后半平面交得出圆心范围. 然后旋转卡壳求最远点对这道题首先样例给的就不对,不知道什么情况然后一直WA,后来看了discuss,试了下这组数据: 4 1 0 0 0 ...
[POJ3384]Feng Shui（半平面交+凸包）
题目: 我是超链接题意: 在一个凸多边形中放两个半径固定的圆,要求输出使覆盖面积最大的(可重叠)两个圆圆心题解: Emmm...套路缩小凸多边形,形成的凸包是圆心可以在的位置,然后暴力找最远点! ...
Using the five elements of fashion colors to create the best feng shui in auspicious Home
Romania, just a legend Licheng Peng Galaxy SOHO opening-day sales 46 Pan Shiyi Students start the On ...

Feng Shui POJ - 3384

Feng Shui POJ - 3384相关推荐

最新文章

热门文章