首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >PKU POJ 1720 SQUARES 解题报告

PKU POJ 1720 SQUARES 解题报告

作者头像
owent
发布2018-08-01 17:26:36
2660
发布2018-08-01 17:26:36
举报
文章被收录于专栏:owentowentowent

题目链接:http://acm.pku.edu.cn/JudgeOnline/problem?id=1720

这题纯计算几何就搞定了,开始我写了个很长很长的代码,但是Wa掉,也不知道是代码那里有疏漏还是精度问题

后来看了一篇解题报告,改为全用long处理,然后根据他的一些方法,代码量也精简了不少

我的原先的方法是判断左边和底边的覆盖,分别计算每个square的覆盖区域和长度

解题报告里也是计算覆盖,不过用斜率代替了我分别计算的左边和底边,同时这里我原先用的是double,但是以斜

率计算可以优化为记录斜率的点,这样就不需要用double,解决了精度问题

然后是判断square的重叠问题,我原先是分成两部分的,左边和底边的在后的被在前的覆盖,报告里这里有更好的

算法,记录边界点,以x轴和y轴坐标和来判断哪个被哪个覆盖

最后只要对每一个square,从上到下找出被覆盖的斜率,然后最后被覆盖的位置相对于圆点的斜率大于square右下

角点的相对于圆点的斜率,那么这个square必然没被覆盖完全,也就是看得到

我的代码:

/**
 * URL:http://acm.pku.edu.cn/JudgeOnline/problem?id=1720
 * Author: OWenT
 * 计算几何
 */
#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
using namespace std;


typedef struct
{
    long x,y;
    long L;
}point;

point zero = {0, 0};
point p[1005];

int cross( point p0, point p1 , point p2)
{
    int t = ( p1.x -p0.x ) * ( p2.y -p0.y ) -( p2.x -p0.x ) * ( p1.y -p0.y);
    if ( t > 0 )
        return 1;
    if ( t < 0 )
        return -1 ;
    return 0 ;
}
bool cmp(point a , point b )
{
    return cross( zero , a, b ) < 0;//按斜率降序排列
}
int main()
{
    long n , i , j , count;
    point tmp1 , tmp2;
    scanf("%ld", &n) ;
    for ( i = 0 ; i < n; i ++)
    {
        scanf("%ld %ld %ld" , &p[i].x, &p[i].y , &p[i].L);
        p[i].y += p[i].L;
    }
    sort(p, p + n, cmp);
    count = 0 ;
    for ( i = 0 ; i < n; i ++)
    {
        tmp1 = p[i] ;
        for ( j = 0 ; j < n; j ++)
        {
            if ( p[i].x + p[i].y <= p[j].x + p[j].y )
                continue;

            tmp2.x = p[j].x + p[j].L ;
            tmp2.y = p[j].y -p[j].L;
            int t1 = cross(zero, p[j] , tmp1 ) ;
            int t2 = cross(zero, tmp2, tmp1 ) ;
            if ( t1 != t2 )//不等的情况必然有覆盖
                tmp1 = tmp2;
            else if(t1 > 0)
                break;//由于按斜率降序排序,第j的点斜率都小于tmp1,那么后面的必然小于
        }
        tmp2.x = p[i].x + p[i].L ;
        tmp2.y = p[i].y - p[i].L ;
        if ( cross(zero, tmp2, tmp1 )> 0 )
            count ++;
    }
    printf("%ld\n", count) ;

    return 0;
}
本文参与 腾讯云自媒体分享计划,分享自作者个人站点/博客。
原始发表:2010-04-03,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 作者个人站点/博客 前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体分享计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档