码迷,mamicode.com
首页 > 其他好文 > 详细

POJ 3335

时间:2014-08-06 14:32:01      阅读:362      评论:0      收藏:0      [点我收藏+]

标签:blog   os   io   for   ar   div   line   算法   

/*半平面交求核心的增量法:
假设前N-1个半平面交,对于第N个半平面,只需用它来交前N-1个平面交出的多边形。
算法开始时,调整点的方向为顺时针方向,对于是否为顺时针,只需求出其面积,若为正,必为逆时针的。
对于每相邻两点求出一条直线,用该直线去交其半平面,并求出交点及判断原多边形点的方位。
*/ 

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
const int MAXN=110;
const double eps=1e-8;

struct point {
	double x,y;
};
point pts[MAXN],p[MAXN],q[MAXN];
int n,cCnt,curCnt;

int DB(double d){
	if(d>eps) return 1;
	if(d<-eps) return -1;
	return 0;
}

double getArea(point *tmp, int n){
	double ans=0;
	for(int i=1;i<=n;i++)
	ans+=(tmp[i].x*tmp[i+1].y-tmp[i].y*tmp[i+1].x);
	return ans/2;
}

void Adjust(point *ps,int n){
	for(int i = 1; i < (n+1)/2; i ++)
 	swap(ps[i], ps[n-i]);
}

void initial(){
	double area=getArea(pts,n);
	if(DB(area)==1) Adjust(pts,n);
	for(int i=1;i<=n;i++)
	p[i]=pts[i];
	p[n+1]=p[1];
	p[0]=p[n];
	cCnt=n;
}

void getline(point x,point y,double &a,double &b,double &c){
    a = y.y - x.y;
    b = x.x - y.x;
    c = y.x * x.y - x.x * y.y;
}

point intersect(point x,point y,double a,double b,double c){
    double u = fabs(a * x.x + b * x.y + c);
    double v = fabs(a * y.x + b * y.y + c);
    point pt;
    pt.x=(x.x * v + y.x * u) / (u + v);
    pt.y=(x.y * v + y.y * u) / (u + v);
    return  pt;
}

void cut(double a,double b,double c){
	curCnt=0;
	for(int i=1;i<=cCnt;i++){
		if(DB(a*p[i].x+b*p[i].y+c)>=0) q[++curCnt] = p[i];
		else {
			if(DB(a*p[i-1].x + b*p[i-1].y + c )>0){
				q[++curCnt] = intersect(p[i],p[i-1],a,b,c);
			}
			 if(DB(a*p[i+1].x + b*p[i+1].y + c )> 0){
                q[++curCnt] = intersect(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 slove(){
	initial();
	for(int i=1;i<=n;i++){
		double a,b,c;
		getline(pts[i],pts[i+1],a,b,c);
		cut(a,b,c);
	}
}

int main(){
	int t;
	scanf("%d",&t);
	while(t--){
		scanf("%d",&n);
		for(int i=1;i<=n;i++)
		scanf("%lf%lf",&pts[i].x,&pts[i].y);
		pts[n+1]=pts[1];
		slove();
		if(cCnt>=1) printf("YES\n");
		else printf("NO\n");
	}
	return 0;
}

  

POJ 3335,布布扣,bubuko.com

POJ 3335

标签:blog   os   io   for   ar   div   line   算法   

原文地址:http://www.cnblogs.com/jie-dcai/p/3894432.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!