codeforces 497 d Gears

最新推荐文章于 2025-01-20 21:44:42 发布

whai362

最新推荐文章于 2025-01-20 21:44:42 发布

阅读量914

点赞数

CC 4.0 BY-SA版权

本文链接：https://blog.csdn.net/whai362/article/details/42099259

计算几何专栏收录该内容

5 篇文章

订阅专栏

本文介绍了一种判断两个绕不同中心旋转的多边形是否相交的方法。通过将问题简化为点绕中心旋转与线段绕另一中心旋转相交的问题，并利用两层循环进行判断。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目链接

点击打开链接

题意：给出点P,Q,和两个多边形，这两个多边形分别绕P,Q以相同的角速度旋转，问这两个多边形能不能相交。

限制：

0<=|x|,|y|<=1e4

多边形边边数小于1000

题解：

先考虑一个线段绕P旋转，一个点绕Q旋转的情况，在什么情况下点和线段相交。

题目就变成1000条线段和1000个点能不能相交，两层循环直接判。

/*cf283div1d
  题意：给出点P,Q,和两个多边形，这两个多边形分别绕P,Q以相同的角速度旋转，问这两个多边形能不能相交。
  限制：
  0<=|x|,|y|<=1e4
  多边形边边数小于1000
  题解：
  先考虑一个线段绕P旋转，一个点绕Q旋转的情况，在什么情况下点和线段相交。
  题目就变成1000条线段和1000个点能不能相交，两层循环直接判。
 */
#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
#define LL __int64
struct Point{
	double x,y;
	Point(){}
	Point(int _x,int _y){ x=_x; y=_y; }
	void scan(){
		scanf("%lf%lf",&x,&y);
	}
	Point operator + (Point a){
		return Point(x+a.x,y+a.y);
	}
	Point operator - (Point a){
		return Point(x-a.x,y-a.y);
	}
	double operator * (Point a){
		return x*a.x+y*a.y;
	}
	double operator ^ (Point a){
		return x*a.y-y*a.x;
	}
	Point rotate(){
		return Point(-y,x);
	}
	double len(){
		return *this * *this;
	}
};
int n,m;
Point P,Q;
const int N=1010;
Point A[N],B[N];
bool ans=false;
double R;
int sgn(double x){
	if(x==0) return 0;
	if(x<0) return -1;
	return 1;
}
double getMinD(Point O,Point C,Point D){
	double minD = min( (O-C).len(),(O-D).len() );
	Point OC=C-O;
	Point OD=D-O;
	Point OH=(D-C).rotate();
	if( sgn(OH^OC) * sgn(OH^OD) !=-1 )
		return minD;
	double s=fabs(OC^OD);
	double h2=s*s/(C-D).len();
	return h2;
}
bool gao(Point O,Point C,Point D){
	double maxD=max( (O-C).len(),(O-D).len() );
	double minD=getMinD(O,C,D);
	if(minD <= R && R <= maxD)
		return true;
	return false;
}
int main(){
	P.scan();
	scanf("%d",&n);
	for(int i=0;i<n;++i)
		A[i].scan();
	Q.scan();
	scanf("%d",&m);
	for(int i=0;i<m;++i)
		B[i].scan();
	Point PQ=Q-P;
	R=PQ.len();
	for(int i=0;i<n;++i){
		Point O=A[i]+PQ;
		for(int j=0;j<m;++j){
			ans |= gao(O,B[j],B[(j+1)%m]);
		}
	}
	for(int i=0;i<m;++i){
		Point O=B[i]-PQ;
		for(int j=0;j<n;++j){
			ans |= gao(O,A[j],A[(j+1)%n]);
		}
	}
	if(ans)
		puts("YES");
	else
		puts("NO");
	return 0;
}