hdu 4643(计算几何)
题意:容易理解
分析:切换的地点为两个基站所在直线的中垂线与两座城市所在直线的交点。
代码实现:
#include <cstdio>
#include <cmath>
#include <algorithm>
#define maxn 60
#define eps 1e-7
using namespace std;
int dcmp(double x) //控制精度
{
if(fabs(x)<eps) return ;
else return x<?-:;
}
double toRad(double deg) //角度转弧度
{
return deg/180.0*acos(-1.0);
}
struct Point
{
double x,y;
Point(){}
Point(double x,double y):x(x),y(y) {}
void input()
{
scanf("%lf %lf",&x,&y);
}
};
typedef Point Vector;
Vector operator+( Vector A, Vector B ) //向量加
{
return Vector( 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 );
}
bool operator<(const Point& A, const Point& B ) //两点比较
{
return dcmp( A.x - B.x ) < || ( dcmp( A.x - B.x ) == && dcmp( A.y - B.y ) < );
}
bool operator==( const Point& a, const Point& b ) //两点相等
{
return dcmp( a.x - b.x ) == && dcmp( a.y - b.y ) == ;
}
struct Line
{
Point s,e;
Vector v;
Line() {}
Line(Point s,Point v,int type)://法向量式
s(s),v(v){}
Line(Point s,Point e):s(s),e(e)//两点式
{v=e-s;} };
double Dot(Vector A,Vector B)//向量点乘
{
return A.x*B.x+A.y*B.y;
}
double Length(Vector A)//向量模
{
return sqrt(Dot(A,A));
}
double Angle(Vector A,Vector B)//向量夹角
{
return acos(Dot(A,B)/Length(A)/Length(B));
}
double Cross(Vector A,Vector B)//向量叉积
{
return A.x*B.y-A.y*B.x;
}
double Area2(Point A,Point B,Point C )//向量有向面积
{
return Cross(B-A,C-A);
}
double Dist(Point A,Point B)
{
return Length(A-B);
}
Vector Rotate(Vector A, double rad)//向量逆时针旋转
{
return Vector(A.x*cos(rad)-A.y*sin(rad),A.x*sin(rad)+A.y*cos(rad));
}
Vector Normal(Vector A)//向量单位法向量
{
double L=Length(A);
return Vector(-A.y/L,A.x/L);
}
Point GetLineIntersection(Line l1,Line l2)//两直线交点
{
Point P=l1.s;
Vector v=l1.v;
Point Q=l2.s;
Vector w=l2.v;
Vector u=P-Q;
double t=Cross(w,u)/Cross(v,w);
return P+v*t;
}
double DistanceToLine(Point P,Line L)//点到直线的距离
{
Point A,B;
A=L.s,B=L.e;
Vector v1=B-A,v2=P-A;
return fabs(Cross(v1,v2))/Length(v1);
}
double DistanceToSegment(Point P, Line L)//点到线段的距离
{
Point A,B;
A=L.s,B=L.e;
if(A==B) return Length(P-A);
Vector v1=B-A,v2=P-A,v3=P-B;
if (dcmp(Dot(v1,v2))<) return Length(v2);
else if (dcmp(Dot(v1,v3))>) return Length(v3);
else return fabs(Cross(v1,v2)) / Length(v1);
}
Point GetLineProjection(Point P,Line L)// 点在直线上的投影
{
Point A,B;
A=L.s,B=L.e;
Vector v=B-A;
return A+v*(Dot(v,P-A)/Dot(v,v));
}
bool OnSegment(Point p,Line l)//点在线段上包括端点
{
Point a1=l.s;
Point a2=l.e;
return dcmp(Cross(a1-p,a2-p))==&&dcmp(Dist(p,a1)+Dist(p,a2)-Dist(a1,a2))==;
}
bool Paralled(Line l1,Line l2)//直线平行
{
return dcmp(Cross(l1.e-l1.s,l2.e-l2.s))==;
}
bool SegmentProperIntersection(Line l1,Line l2)//线段相交
{
if(Paralled(l1,l2))
{
return false;
}
Point t=GetLineIntersection(l1,l2);
if(OnSegment(t,l1))
{
return true;
}
return false;
}
int ConvexHull(Point *p,int n,Point *ch) //求凸包
{
sort(p,p+n);
int m=;
for ( int i = ; i < n; ++i )
{
while ( m > && Cross( ch[m - ] - ch[m - ], p[i] - ch[m - ] ) <= ) --m;
ch[m++] = p[i];
}
int k = m;
for ( i = n - ; i >= ; --i )
{
while ( m > k && Cross( ch[m - ] - ch[m - ], p[i] - ch[m - ] ) <= ) --m;
ch[m++] = p[i];
}
if ( n > ) --m;
return m;
}
double PolygonArea(Point *p,int n) //多边形有向面积
{
double area=;
for (int i=;i<n-;++i)
area+=Cross(p[i]-p[],p[i+]-p[]);
return area/2.0;
}
///*******************************************************/// int N,M,Q;
Point city[maxn],base[maxn]; int solve(int a,int b)
{
Line train=Line(city[a],city[b]);
int res=;
int i,j,k;
for(i=; i<M; i++)
{
for(j=i+; j<M; j++)
{
Line vert=Line((base[i]+base[j])/,Normal(base[i]-base[j]),);
if(Paralled(train,vert))
{
// printf("caonima\n");
continue;
}
Point t=GetLineIntersection(train,vert);
if(!OnSegment(t,train)) continue;
double d1=Dist(t,base[i]); int flag=;
for(k=;k<M;k++)
{
if(k==i||k==j) continue;
double d2=Dist(t,base[k]);
if(dcmp(d2-d1)<)
{
flag=;
break;
}
}
if(flag)
{
res++;
}
}
}
return res;
}
int main()
{
//freopen("input.txt","r",stdin);
int i,j;
while(scanf("%d %d",&N,&M)==)
{
for(i=; i<N; i++)
{
city[i].input();
}
for(i=; i<M; i++)
{
base[i].input();
}
scanf("%d",&Q);
while(Q--)
{
int a,b;
scanf("%d %d",&a,&b);
a--,b--;
printf("%d\n",solve(a,b));
}
}
return ;
}
hdu 4643(计算几何)的更多相关文章
- hdu 4643 GSM 计算几何 - 点线关系
/* hdu 4643 GSM 计算几何 - 点线关系 N个城市,任意两个城市之间都有沿他们之间直线的铁路 M个基站 问从城市A到城市B需要切换几次基站 当从基站a切换到基站b时,切换的地点就是ab的 ...
- HDU 4643 GSM 简单计算几何
今天比赛的时候略坑, admin告诉我询问Q的个数不超过n^2, 赛后敲了个 O(Q*m^3)的复杂度,但这个复杂度常数比较低,可能在除以个小常数, 300ms过了,真心无语,数据应该水了吧,比赛的时 ...
- HDU 4643 GSM (2013多校5 1001题 计算几何)
GSM Time Limit: 4000/2000 MS (Java/Others) Memory Limit: 65535/32768 K (Java/Others)Total Submiss ...
- hdu 2108:Shape of HDU(计算几何,判断多边形是否是凸多边形,水题)
Shape of HDU Time Limit: 3000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Tota ...
- HDU 2202 计算几何
最大三角形 Time Limit: 5000/2000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submi ...
- *HDU 2108 计算几何
Shape of HDU Time Limit: 3000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Tota ...
- HDU 5784 (计算几何)
Problem How Many Triangles (HDU 5784) 题目大意 给定平面上的n个点(n<2000),询问可以组成多少个锐角三角形. 解题分析 直接统计锐角三角形较困难,考虑 ...
- hdu 4720 计算几何简单题
昨天用vim练了一道大水题,今天特地找了道稍难一点的题.不过也不是很难,简单的计算几何而已.练习用vim编码,用gdb调试,结果居然1A了,没调试...囧... 做法很简单,无非就是两种情况:①三个巫 ...
- HDU 6205[计算几何,JAVA]
题目链接[http://acm.hdu.edu.cn/showproblem.php?pid=6206] 题意: 给出不共线的三个点,和一个点(x,y),然后判断(x,y)在不在这三个点组成的圆外. ...
随机推荐
- JavaPersistenceWithHibernate第二版笔记Getting started with ORM-001用JPA和Hibernate实现HellowWorld(JTA、Bitronix)
一.结构 二.model层 1. package org.jpwh.model.helloworld; import javax.persistence.Entity; import javax.pe ...
- SPRING IN ACTION 第4版笔记-第九章Securing web applications-011-把敏感信息请求转为https(requiresChannel())
1.把包含敏感信息的请求转为https请求,则较为安全,但如何只把有需要安全的请求转为https,而不是不加分辩就把所有请求都转为https呢?可以用requiresChannel() @Overri ...
- RxJava+Retrofit+MVP构建的App——聚合资讯
RtfRxMVP 聚合资讯APP,提供热点资讯,天气预报以及笑话精选服务,使用 Retrofit + RxJava + MVP 构建代码. Hello U 这是我的一个练习项目,第一次尝试运用 MVP ...
- iOS中检测硬件和传感器
首先要知道,你需要查看所需的硬件或传感器是否存在,而不是假设设备有哪些功能.举个例子,你不能假设只有iPhone才有麦克风,而应该使用API来查看麦克风是否存在.下面这段代码的第一个优势在于,它能自动 ...
- Linux下SSH各配置项解释
关于ssh 设置的相关总结(ssh最大连接数.ssh连接时长.安全性配置等) 以redhat6.3为例 ssh配置文件在: /etc/ssh/sshd_config 可以打开查看相应配置,默认情况下只 ...
- 【周全考虑】CodeForces 245B——Internet Address
来源:点击打开链接 看上去很简单的一道题,可是错的次数却不少. 题目要求是将一串字母转化成网址——形如格式http(ftp)://xxx.ru/xxxx的样子,看上去很简单,可是还是很容易出错. 刚开 ...
- Mac下无法拷贝文件到移动硬盘
Mac下无法拷贝文件到移动硬盘? 是移动硬盘的文件格式的问题. Mac系统无法识别 NTFS 格式的文件. 将移动硬盘格式化为 exFAT 格式的. 别担心,exFAT 格式的硬盘在Windows下也 ...
- ccnu-线段树联系-单点更新2-B
B - 单点更新2 Time Limit:3000MS Memory Limit:32768KB 64bit IO Format:%I64d & %I64u Submit Status Des ...
- Evaluate Reverse Polish Notation(堆栈)
Evaluate the value of an arithmetic expression in Reverse Polish Notation. Valid operators are +, -, ...
- windows平台时间函数性能比较QueryPerformanceCounter,GetTickCount,ftime,time,GetLocalTime,GetSystemTimeAsFileTime
http://gmd20.blog.163.com/blog/static/168439232012113111759514/ 执行 10000000 次, 耗时 2258,369 微秒 Qu ...