#叉积#zoj 1041 poj 1106 ssl 1232 雷达覆盖问题 Transmitters

原创于 2019-01-04 20:05:48 发布 · 273 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#zoj 1041 #poj 1106 #ssl 1232 #Transmitters #雷达覆盖问题

计算几何专栏收录该内容

11 篇文章

订阅专栏

博客探讨了如何利用叉积来解决雷达覆盖问题，包括zoj 1041, poj 1106, ssl 1232等题目，通过分析和代码展示如何确定雷达旋转时能覆盖到的最大点数。" 6749033,821594,使用boost::multi_index_container和boost::shared_ptr管理数据,"['C++', '容器', '内存管理', '游戏开发', '数据结构']

题目

以雷达心为圆心的半圆形雷达覆盖范围有多个点雷达可旋转，求最多覆盖数（含在边界的）

分析

首先如果超过了范围，首先肯定是不可以的，再而可以枚举一个合法的点，那么这个点与雷达心所在的直线上正好把圆形分成两个半圆，那如何区分呢，就可以用叉积来区分对于该直线，每个点对于合法的点是处于顺时针还是逆时针，这样就可以求出答案

代码

#include <cstdio>
#define rr register
#define max(a,b) (((a)>(b))?(a):(b))
#define diss(i) (xx-x[i])*(xx-x[i])+(yy-y[i])*(yy-y[i])
using namespace std;
int xx,yy,n,x[1001],y[1001]; double rad;
signed main(){
	while (scanf("%d%d%lf%d",&xx,&yy,&rad,&n)==4){
		if (rad<0) return 0;
		rr int ans=0;
		for (rr int i=1;i<=n;++i){
			scanf("%d%d",&x[i],&y[i]);
			if (diss(i)>rad*rad) --n,--i;//不合法
		}
		for (rr int i=1;i<=n;++i){
			rr int lef=0,rig=0;
			for (rr int j=1;j<=n;++j){
				rr int m=(x[i]-xx)*(y[j]-yy)-(x[j]-xx)*(y[i]-yy);
				lef+=m>=0; rig+=m<=0;//分成两个部分
			}
			ans=max(ans,max(lef,rig));//取最大值
		}
		printf("%d\n",ans);
	}
}