uva10256如何判断俩个凸包是否相交

本题并没有多大的难度,就是点集分割,然后求凸包,判断俩个凸包是否相交,

重点在于如何判断来个凸包是否相交,

分俩步进行,

1.判断俩个凸包上的任意线段是否规范相交

2.判断各自凸包上的点是否在另一个凸包内部以及边界上

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cctype>
#include<cmath>
#include<string>
#include<set>
#include<map>
#include<vector>
#include<algorithm>
using namespace std;
const double eps=1e-10;
const int N=1e5;
int dcmp(double x)
{
    if(fabs(x)<eps)
        return 0;
    return x>0?1:-1;
}
struct point
{
    double x,y;
    point() {}
    point(double x,double y):x(x),y(y){}
    point operator + (const point &a) const
    {
        return point(x+a.x,y+a.y);
    }
    point operator - (const point &a) const
    {
        return point(x-a.x,y-a.y);
    }
    void in()
    {
         cin>>x>>y;
    }
    bool operator < (const point &a) const
    {
        return x+eps<a.x||(!dcmp(x-a.x)&&y<a.y);
    }
};
double dot(point a,point b)
{
    return a.x*b.x+a.y*b.y;
}
double cross(point a,point b)
{
    return a.x*b.y-a.y*b.x;
}
bool SegmentProperIntersectiom(point a1,point a2,point b1,point b2)
{
    double c1=cross(a2-a1,b1-a1),c2=cross(a2-a1,b2-a1),
           c3=cross(b2-b1,a1-b1),c4=cross(b2-b1,a2-b1);
    return dcmp(c1)*dcmp(c2)<=0&&dcmp(c3)*dcmp(c4)<=0;
}
bool OnSegment(point p,point a1,point a2)
{
    return dcmp(cross(a1-p,a2-p))==0&&dcmp(dot(a1-p,a2-p))<0;
}
int isPointInpolygon(point p,point *ch,int n)
{
    int wn=0;
    for(int i=0;i<n;i++)
    {
        if(OnSegment(p,ch[i],ch[(i+1)%n])) return 1;
        int k=dcmp(cross(ch[(i+1)%n]-ch[i],p-ch[i]));
        int d1=dcmp(ch[i].y-p.y);
        int d2=dcmp(ch[(i+1)%n].y-p.y);
        if(k>0&&d1<=0&&d2>0) wn++;
        if(k<0&&d2<=0&&d1>0) wn--;
    }
    if(wn!=0) return 1;
    return 0;
}
int ConvexHull(point *p,int n,point *ch)
{
    sort(p,p+n);
    int m=0;
    for(int i=0;i<n;i++)
    {
        while(m>1&&cross(ch[m-1]-ch[m-2],p[i]-ch[m-2])<=0) m--;
        ch[m++]=p[i];
    }
    int k=m;
    for(int i=n-2;i>=0;i--)
    {
        while(m>k&&cross(ch[m-1]-ch[m-2],p[i]-ch[m-2])<=0) m--;
        ch[m++]=p[i];
    }
    if(n>1) m--;
    return m;
}
int main()
{
    int m,c;
    while(cin>>m>>c)
    {
        if(m==0&&c==0) break;
        point p1[510],p2[510],ch1[510],ch2[510];
        for(int i=0;i<m;i++)
            p1[i].in();
        for(int i=0;i<c;i++)
            p2[i].in();
        int n1=ConvexHull(p1,m,ch1);
        int n2=ConvexHull(p2,c,ch2);
        int flag=1;
        for(int i=0;i<n1;i++)
            if(isPointInpolygon(ch1[i],ch2,n2)){
                flag=0;
                break;
            }
          //  cout<<flag<<endl;
        for(int i=0;i<n2;i++)
            if(isPointInpolygon(ch2[i],ch1,n1)){
            flag=0;
            break;
            }
         //   cout<<flag<<endl;
        ch1[n1]=ch1[0];
        ch2[n2]=ch2[0];
        for(int i=0;i<n1;i++)
        {
            for(int j=0;j<n2;j++)
            {
                if(SegmentProperIntersectiom(ch1[i],ch1[i+1],ch2[j],ch2[j+1]))
                {
                    flag=0;
                  //  cout<<flag<<endl;
                    break;
                }
            }
            if(flag==0)
                break;
        }
       // cout<<flag<<endl;
        if(flag)
            cout<<"Yes"<<endl;
        else
            cout<<"No"<<endl;
    }
    return 0;
}


  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值