计算几何【套模板,推荐模板】 Separating Pebbles UVALive - 7461

13 篇文章 0 订阅
2 篇文章 0 订阅

题目链接


= =生平第一次计算几何就GG?

╮(╯▽╰)╭                因为生成线段的时候为了省时选择的第二层循环起始错误而错误。。。

不过我又一次见识到了板子的重要性。。。

这题如果最后多试试还是可以过的,当然没有如果。

推荐模板:可下载

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
const int INF = 1E200;
const int EP=1E-10;
const int MAXV=300;
const double PI = 3.14159265;
using namespace std;
struct POINT
{
    double x;
    double y;
    POINT(double a=0, double b=0)
    {
        x=a;
        y=b;
    }
};
struct LINESEG
{
    POINT s;
    POINT e;
    LINESEG(POINT a, POINT b)
    {
        s=a;
        e=b;
    }
    LINESEG() {}
};
struct LINE
{
    double a;
    double b;
    double c;
    LINE(double d1=1, double d2=-1, double d3=0)
    {
        a=d1;
        b=d2;
        c=d3;
    }
};
double multiply(POINT sp, POINT ep, POINT op)
{
    return ((sp.x-op.x)*(ep.y-op.y)-(ep.x-op.x)*(sp.y-op.y));
}
bool intersect(LINESEG u, LINESEG v)
{
    return ((max(u.s.x, u.e.x)>=min(v.s.x, v.e.x))&&
            (max(v.s.x, v.e.x)>=min(u.s.x, u.e.x))&&
            (max(u.s.y, u.e.y)>=min(v.s.y, v.e.y))&&
            (max(v.s.y, v.e.y)>=min(u.s.y, u.e.y))&&
            (multiply(v.s, u.e, u.s)* multiply(u.e, v.e, u.s)>=0)&&
            (multiply(u.s, v.e, v.s)* multiply(v.e, u.e, v.s)>=0));
}
int main()
{
    int K, n;
    scanf("%d",&K);
    while(K--)
    {
        POINT p[300],q[300];
        LINESEG aLineseg[62501], bLineseg[62501];
        scanf("%d",&n);
        int op=0,oq=0;
        for(int i=0; i<n; i++)
        {
            int x,y,c;
            scanf("%d%d%d",&x,&y,&c);
            if(c==0)
            {
                p[op].x = x;
                p[op++].y = y;
            }
            else
            {
                q[oq].x = x;
                q[oq++].y  = y;
            }
        }
        int flag=1;
        for(int i=0;i<op;i++)
        {
            for(int j=0;j<oq;j++)
            {
                if(p[i].x==q[j].x && p[i].y == q[j].y)
                    {flag=0;break;}
            }
            if(flag==0) break;
        }
        if(!flag) 
        {printf("1\n");continue;}
        int a1=0, b1=0;
        for(int i=0; i<op; i++)
        {
            for(int j=0; j<op; j++)                    //就是这里
            {
                aLineseg[a1++] = LINESEG(p[i],p[j]);
            }
        }
        for(int i=0; i<oq; i++)
        {
            for(int j=0; j<oq; j++)                    //和这里   应该是因为储存的方式 这是个有向线段
            {
                bLineseg[b1++] = LINESEG(q[i],q[j]);
            }
        }
        for(int i=0; i<a1; i++)
        {
            for(int j=0; j<b1; j++)
            {
                if(intersect(aLineseg[i],bLineseg[j]))
                {
                    flag=0;
                    break;
                }
            }
        }
        if(flag) printf("1\n");
        else printf("0\n");
    }
    return 0;
}



评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值