poj 2540 Hotter Colder 半平面交

最新推荐文章于 2020-11-23 18:39:51 发布
最新推荐文章于 2020-11-23 18:39:51 发布
阅读量701 收藏
点赞数
分类专栏： ACM_计算几何 半平面交 文章标签： poj 计算几何
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/yskyskyer123/article/details/52294620
版权

题目链接

题意：

在一个(0,0) ,(10,0),(10,10),(0,10)四个点围成的正方形房屋内藏着一个人B，现在A蒙着眼从(0,0)出发，之后经过一些点停下，B会告诉他这个位置和上一次比是距自己近了，还是远了。

每一步输出B可能所在区域的面积。

解法：

首先一点，如果是”Same”,那么之后都输出0.00 。

每一步半平面交之后求凸多边形面积。

对于两个连续的停留点A、B，首先确定中点midP，然后 v⃗ =AB−→− ，之后根据 (x,y)∗(y,−x)=0 ，将 v⃗ =(v.y,−v.x) ,这就求出了中垂线方向上的向量。

我们假定是Colder，那么A点必然在 v⃗  左边，所以直接判断A点 v⃗  是否在左边，如果不是调换 v⃗  的方向，那么现在A点必然在 v⃗  左边，如果这时发现是Hotter,那么再次调换 v⃗  的方向。
这样midP和 v⃗  就构成了一条射线(一个半平面)。

代码：

#include<cstdio>
#include<string>
#include<cstring>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<vector>
using namespace std;

#define all(x) (x).begin(), (x).end()
#define for0(a, n) for (int (a) = 0; (a) < (n); (a)++)
#define for1(a, n) for (int (a) = 1; (a) <= (n); (a)++)
#define mes(a,x,s)  memset(a,x,(s)*sizeof a[0])
#define mem(a,x)  memset(a,x,sizeof a)
#define ysk(x)  (1<<(x))
typedef long long ll;
typedef pair<int, int> pii;
const int INF =0x3f3f3f3f;
const int maxn= 100   ;
const double PI=acos(-1.0);
const double eps=1e-10;
int dcmp(double x)
{
    if(fabs(x)<eps)  return 0;
    else return x<0?-1:1;
}
struct Point
{
    double x,y;
    Point(double x=0,double y=0):x(x),y(y) {};
    bool operator ==(const Point B)const {return dcmp(x-B.x)==0&&dcmp(y-B.y)==0;}

    bool operator<(const Point& b)const
    {
        return dcmp(x-b.x)<0|| dcmp(x-b.x)==0 &&dcmp(y-b.y)<0;
    }
};
typedef Point Vector;
Vector operator -(Vector A,Vector B) {return Vector(A.x-B.x,A.y-B.y); }


double Cross(Vector A,Vector B)//叉乘
{
    return A.x*B.y-A.y*B.x;
}

double Dot(Vector A,Vector B)//点乘
{
    return A.x*B.x+A.y*B.y;
}

Vector operator +(Vector A,Vector B) {return Vector(A.x+B.x,A.y+B.y); }
Vector operator *(Vector A,double p) {return Vector(A.x*p,A.y*p); }
Vector operator /(Vector A,double p) {return Vector(A.x/p,A.y/p); }
Vector operator -(Vector A)  {return  Vector(-A.x,-A.y);}

struct Line
{
    Point p,p2;
    Vector v;
    double ang;
    Line(){}
    Line(Point a,Vector v):p(a),v(v){ang=atan2(v.y,v.x); }//点线式
    void twoPointIntial(Point p,Point p2)//两点式
    {
        this->p=p;
        this->p2=p2;
        v=  p2-p;
        ang=atan2(v.y,v.x);
    }
    Point point(double t)//参数方程求参数t对应点
    {
        return p+v*t;
    }
    bool operator<(const Line & L)const
    {
        return ang<L.ang;
    }
};


Point GetIntersection(Line a,Line b)
{
    Vector u=a.p-b.p;
    double t=Cross(b.v,u)/Cross(a.v,b.v);
    return a.p+a.v*t;
}

bool OnRight(Line L,Point p)
{
    return dcmp(Cross(L.v,p-L.p))<0;
}


int HalfplaneIntersection(Line* L,int n,Point * poly)
{
    sort(L,L+n);
    int first,last;
    Point *p=new Point[n];
    Line *q=new Line[n];
    q[first=last=0]=L[0];
    for(int i=1;i<n;i++)
    {
        while(first<last&&OnRight(L[i],p[last-1] ) )  last--;
        while(first<last&&OnRight(L[i],p[first]))  first++;
        q[++last]=L[i];
        if(fabs(Cross(q[last].v,q[last-1].v) )<eps )
        {
            last--;
            if(!OnRight(q[last],L[i].p) )  q[last]=L[i];
        }
        if(first<last) p[last-1]=GetIntersection(q[last-1],q[last] );
    }
    while(first<last&&OnRight(q[first],p[last-1])) last--;
    if(last-first<=1)  return 0;//小于3个，不存在
    p[last]=GetIntersection(q[last],q[first]);

    int m=0;
    for(int i=first;i<=last;i++) poly[m++]=p[i];
    return m;

}

double PolygonArea(Point *p,int n)
{
    double area=0;
    for(int i=1;i<n-1;i++)
    {
        area+=Cross(p[i]-p[0],p[i+1]-p[0]);
    }
    return area/2;
}

Point getMidPoint(Point a,Point b)
{
    return Point( (a.x+b.x)/2,(a.y+b.y)/2 );
}
int n,np;
Line L[maxn+10];
Point p[maxn+10],poly[maxn+10];
string s;



int main()
{
   std::ios::sync_with_stdio(false);
   double x,y;
   bool zero=0;n=0;
    L[n++].twoPointIntial(Point(0,0),Point(10,0));
    L[n++].twoPointIntial(Point(10,0),Point(10,10));
    L[n++].twoPointIntial(Point(10,10),Point(0,10));
    L[n++].twoPointIntial(Point(0,10),Point(0,0));
    np=1;
   while(cin>>p[np].x>>p[np].y)
   {
        cin>>s;
        if(s=="Same") zero=1;

        Vector v=p[np]-p[np-1];
        v=Vector(v.y,-v.x);
        Point midP=getMidPoint(p[np],p[np-1]);

        if(dcmp(Cross(v,p[np-1]-midP ) )<0 )  v=-v;
        if(s=="Hotter")  v=-v;

        L[n++]=Line( midP,v);

        if(zero)
        {
            puts("0.00");
        }
        else
        {
            int m=HalfplaneIntersection(L,n,poly);
            double s=PolygonArea(poly,m);
            if(dcmp(s)==0 )  puts("0.00");
            else printf("%.2f\n",s);
        }
        np++;
   }
   return 0;
}




#include<cstdio>
#include<string>
#include<cstring>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<vector>
using namespace std;

#define all(x) (x).begin(), (x).end()
#define for0(a, n) for (int (a) = 0; (a) < (n); (a)++)
#define for1(a, n) for (int (a) = 1; (a) <= (n); (a)++)
#define mes(a,x,s)  memset(a,x,(s)*sizeof a[0])
#define mem(a,x)  memset(a,x,sizeof a)
#define ysk(x)  (1<<(x))
typedef long long ll;
typedef pair<int, int> pii;
const int INF =0x3f3f3f3f;
const int maxn= 100   ;
const double PI=acos(-1.0);
const double eps=1e-10;
int dcmp(double x)
{
    if(fabs(x)<eps)  return 0;
    else return x<0?-1:1;
}
struct Point
{
    double x,y;
    Point(double x=0,double y=0):x(x),y(y) {};
    bool operator ==(const Point B)const {return dcmp(x-B.x)==0&&dcmp(y-B.y)==0;}

    bool operator<(const Point& b)const
    {
        return dcmp(x-b.x)<0|| dcmp(x-b.x)==0 &&dcmp(y-b.y)<0;
    }
};
typedef Point Vector;
Vector operator -(Vector A,Vector B) {return Vector(A.x-B.x,A.y-B.y); }


double Cross(Vector A,Vector B)//叉乘
{
    return A.x*B.y-A.y*B.x;
}

double Dot(Vector A,Vector B)//点乘
{
    return A.x*B.x+A.y*B.y;
}

Vector operator +(Vector A,Vector B) {return Vector(A.x+B.x,A.y+B.y); }
Vector operator *(Vector A,double p) {return Vector(A.x*p,A.y*p); }
Vector operator /(Vector A,double p) {return Vector(A.x/p,A.y/p); }
Vector operator -(Vector A)  {return  Vector(-A.x,-A.y);}

struct Line
{
    Point p,p2;
    Vector v;
    double ang;
    Line(){}
    Line(Point a,Vector v):p(a),v(v){ang=atan2(v.y,v.x); }//点线式
    void twoPointIntial(Point p,Point p2)//两点式
    {
        this->p=p;
        this->p2=p2;
        v=  p2-p;
        ang=atan2(v.y,v.x);
    }
    Point point(double t)//参数方程求参数t对应点
    {
        return p+v*t;
    }
    bool operator<(const Line & L)const
    {
        return ang<L.ang;
    }
};


Point GetIntersection(Line a,Line b)
{
    Vector u=a.p-b.p;
    double t=Cross(b.v,u)/Cross(a.v,b.v);
    return a.p+a.v*t;
}

bool OnRight(Line L,Point p)
{
    return dcmp(Cross(L.v,p-L.p))<0;
}


int HalfplaneIntersection(Line* L,int n,Point * poly)
{
    sort(L,L+n);
    int first,last;
    Point *p=new Point[n];
    Line *q=new Line[n];
    q[first=last=0]=L[0];
    for(int i=1;i<n;i++)
    {
        while(first<last&&OnRight(L[i],p[last-1] ) )  last--;
        while(first<last&&OnRight(L[i],p[first]))  first++;
        q[++last]=L[i];
        if(fabs(Cross(q[last].v,q[last-1].v) )<eps )
        {
            last--;
            if(!OnRight(q[last],L[i].p) )  q[last]=L[i];
        }
        if(first<last) p[last-1]=GetIntersection(q[last-1],q[last] );
    }
    while(first<last&&OnRight(q[first],p[last-1])) last--;
    if(last-first<=1)  return 0;//小于3个，不存在
    p[last]=GetIntersection(q[last],q[first]);

    int m=0;
    for(int i=first;i<=last;i++) poly[m++]=p[i];
    return m;

}

double PolygonArea(Point *p,int n)
{
    double area=0;
    for(int i=1;i<n-1;i++)
    {
        area+=Cross(p[i]-p[0],p[i+1]-p[0]);
    }
    return area/2;
}

Point getMidPoint(Point a,Point b)
{
    return Point( (a.x+b.x)/2,(a.y+b.y)/2 );
}
int n,np;
Line L[maxn+10];
Point p[maxn+10],poly[maxn+10];
string s;



int main()
{
   std::ios::sync_with_stdio(false);
   double x,y;
   bool zero=0;n=0;
    L[n++].twoPointIntial(Point(0,0),Point(10,0));
    L[n++].twoPointIntial(Point(10,0),Point(10,10));
    L[n++].twoPointIntial(Point(10,10),Point(0,10));
    L[n++].twoPointIntial(Point(0,10),Point(0,0));
    np=1;
   while(cin>>p[np].x>>p[np].y)
   {
        cin>>s;
        if(s=="Same") zero=1;

        Vector v=p[np]-p[np-1];
        v=Vector(v.y,-v.x);
        Point midP=getMidPoint(p[np],p[np-1]);

        if(dcmp(Cross(v,p[np-1]-midP ) )<0 )  v=-v;
        if(s=="Hotter")  v=-v;

        L[n++]=Line( midP,v);

        if(zero)
        {
            puts("0.00");
        }
        else
        {
            int m=HalfplaneIntersection(L,n,poly);
            double s=PolygonArea(poly,m);
            if(!m||dcmp(s)==0 )  puts("0.00");
            else printf("%.2f\n",s);
        }
        np++;
   }
   return 0;
}
PoemK
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
poj2451半平面交
10-17
计算几何半平面交，poj2451 一题是模版题，这里的代码可以作为模版使用，速度比较快
acm之半平面交，很少有的好资源
09-28
例如，在“Hotter and Colder”游戏中，根据A对B位置的反馈，我们可以用半平面来定义可能的P点所在区域，并通过每次反馈后更新这个区域，计算面积。而在“Nice Milk”问题中，SRbGa想要最大化面包蘸牛奶的面积，这就...
POJ 2540 Hotter Colder
weixin_34041003的博客
07-08 61
http://poj.org/problem?id=2540 题意:给你每次行走的路径，而且告诉你每次离一个点光源是远了还是近了，要求每次光源可能存在的位置的面积。 思路:如果出现"same"，说明光源在中垂线上，面积为0，否则我们考虑远了或者近了，实际上就是在路径中两点连成直线的中垂线就是半平面直线，近了在这条直线的远侧，远了在这条直线的近侧。 1 #include<cstd...
POJ2540：半平面交
FlyingJack
04-11 130
POJ2540 题解 如果比之前更远，那么在两点的中垂线靠近之前的那边，反之另一边。 我们可以找出a和b的终点，在把确定向量ab，根据情况旋转90或270。 如果出现same那就不再计算 了。 代码 #include <iostream> #include <cstdio> #include <algorithm> #include <cstr...
poj2540 半平面交求可行性区域解线性规划
kth_xixi
11-23 130
题意： 现有一个[0,0]到[10,10]的方格子，已知这里面有一点火源，但你不知道其确切位置你初始在(0,0)每次移动到下一个点并告知你距离火源和上一次比较是更近还是更远或者Same。对于每次移动你需要给出目前火源可能在的区域面积。 思路: 每次找到lst与cur点的中垂线，显然hotter的一边计入答案与原先剩下的面积半平面交 求中垂线： 中点为mid，法向量Fa，中垂线向量(mid,end) end =mid+Fa； Vector QAQ; QAQ.x = cur.x -...
poj 2540 Hotter Colder(线性规划可行区域)
志当存高远
09-07 1245
Hotter Colder Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 1836   Accepted: 764 Description The children's game Hotter Colder is played as follows. Player
POJ2540-Hotter Colder（半平面交）
mowayao的专栏
08-20 1238
Hotter Colder Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 2523   Accepted: 1045 Description The children's game Hotter Colder is played as follows. Playe
半平面相关题解1
08-03
本文将讨论半平面交及其在解决各种几何问题中的应用。 半平面交是指在二维平面上，由一组有向直线所定义的半平面（即直线一侧的区域）的交集。这个问题的关键在于确定这些半平面在特定区域内（例如题目中的正方形[0...
POJ.rar_poj java_poj1048
最新发布
09-24
【标题】"POJ.rar_poj java_poj1048" 涉及的知识点主要围绕编程竞赛中的“约瑟夫环”问题，这里是一个加强版，使用Java语言进行解决。 【描述】"POJ1048，加强版的约瑟夫问题 难度中等" 提示我们，这个问题是编程...
poj.rar_poj
09-23
标题中的"poj.rar_poj"暗示了这是一个与POJ（Programming Online Judge）相关的压缩文件，POJ是一个在线编程竞赛平台，主要供程序员们练习和提交算法解决方案。在这个压缩包中，很可能包含了用户在POJ上参与编程挑战...
【半平面交】poj2540 Hotter Colder
每天心塞一点点
04-16 645
写了三天的半平面交总觉得自己姿势有问题。。 #include #include #include #include #define eps 1e-8 using namespace std; struct Point { double x, y; Point(double _x = 0. , double _y = 0. ) : x(_x), y(
poj2540——半平面交（线性规划问题）
mumei的博客
09-19 299
题目链接：http://poj.org/problem?id=2540 The children's game Hotter Colder is played as follows. Player A leaves the room while player B hides an object somewhere in the room. Player A re-enters at positi...
uva 319 - Pendulum（几何）
不慌不忙、不急不躁
08-19 749
题目链接：uva 319 - Pendulum 注意高度不能高过水平线，一种周期是绕某点一直转圈，一种周期是返回起点。 #include #include #include #include #include using namespace std; const double pi = 4 * atan(1); const double eps = 1e-8; i
poj 2540
[CHN]liurui
07-31 359
题目概述一个房间，左下角坐标(0,0)，右上角坐标(10,10)，A把某物藏于某处，B从(0,0)出发，每次走到一个位置，A会说距离藏物点更近了(Hotter)，更远了(Colder)或距离没变(Same)，给定每次B到的位置(x,y)和A喊的话，求可能藏物的区域的面积 A每次只会说Hotter,Colder或Same时限1000ms/3000ms输入每行两个浮点数x,y，代表B到的位置，以及一个
poj 2540 Hotter Colder （半平面交+线性规划）
clover_hxy的博客
01-04 379
Hotter Colder Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 3255   Accepted: 1304 Description The children's game Hotter Colder is played as follows. Playe
POJ 2540 半平面交
u013112097的专栏
07-23 353
半平面交，套模板。
poj 2540 && uva 10084 Hotter Colder（半平面交）
丁丁丁丁丁
08-11 429
题目大意：“更冷更热”是一个游戏：两个人在一个左下角与右上角左边为（0,0)，（10,10)的正方形房间里，甲闭上眼睛，乙在房间里藏一个东西。然后甲猜这个东西在哪里，第一次必须猜（0,0)。第二次开始每猜一个位置乙都要回答“Hotter”， “Colder”， “Same”。  “Hotter”表示新猜的点比上次猜的点要近，“Colder"表示新猜的点比上次猜的点要远，“Same"表示两次猜得一
poj 2540 求半平面交
Glory 的 博客
09-27 197
题意 ： 给你n 个点和一个标志，告诉你离这个点远了还是近了，让你求出那个点的可能位置。题解 ： 不难发现其实每一个点对应这一条垂直平分线然后看看这个区域有多大就可以了。#include <iostream> #include <algorithm> #include <cstring> #include <cstdio> #include <cmath> #include <string> usi

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值