2020秦皇岛C、 Cameraman

最新推荐文章于 2022-05-11 20:55:35 发布

nof_uck

最新推荐文章于 2022-05-11 20:55:35 发布

阅读量332

点赞数 1

分类专栏： ACM 文章标签： 1024程序员节

本文链接：https://blog.csdn.net/weixin_43571920/article/details/109257245

版权

ACM 专栏收录该内容

52 篇文章 1 订阅

订阅专栏

本文详细介绍了如何解决一个二维空间中找到最大夹角的问题。通过枚举相邻点并计算它们与目标点之间的距离，再根据角度范围进行判断，最终确定最大夹角。代码实现中涉及了极坐标转换、三角函数应用和排序等技术，对于理解算法和数学在编程中的应用具有参考价值。

摘要由CSDN通过智能技术生成

思路

很明显，alex在的位置肯定是bob的位置，因为这样形成的夹角才最大。
然后枚举相邻的点，求一下之间的距离即可。

#pragma GCC optimize(2)
#include<bits/stdc++.h>
using namespace std;
const int man = 2e5+10;
#define IOS ios::sync_with_stdio(0)
#define ull unsigned ll
#define uint unsigned
#define pai pair<int,int>
#define pal pair<ll,ll>
#define IT iterator
#define pb push_back
#define fi first
#define se second
#define For(i,j,k) for (int i=(int)(j);i<=(int)(k);++i)
#define Rep(i,j,k) for (int i=(int)(j);i>=(int)(k);--i)
#define endl '\n'
#define ll long long
const ll mod = 1e9+7;
int x,y,w,h;
double a[6],d[man];
double get_d(pai a){
	return atan2(a.second-y,a.first-x);
}
double pi = acos(-1);
double get(double d){
	double ans = 0;
	if(d>=a[1]&&d<a[2]){
		if(d<=-pi/2.0)ans = x - tan(fabs(d+pi/2.0))*y;
		else ans = x + tan(fabs(pi/2.0+d))*y;
	}else if(d>=a[2]&&d<a[3]){
		ans = w + y;
		if(d<=0)ans -= tan(fabs(d))*(w-x);
		else ans += tan(d)*(w-x);
	}else if(d>=a[3]&&d<a[4]){
		ans = w + h + w - x;
		if(d<=pi/2.0)ans -= tan(pi/2.0-d)*(h-y);
		else ans += tan(d - pi/2.0)*(h-y);
	}else if(d>=a[4]||d<a[1]){
		ans = 2.0*w + h;
		if(d>=a[4])ans += h - y - tan(pi-d)*x;
		else ans += h - y + tan(pi+d)*x;
	}
	return ans;
}

signed main() {
	#ifndef ONLINE_JUDGE
		//freopen("in.txt", "r", stdin);
		//freopen("out.txt","w",stdout);
	#endif
	int t,T = 0;
	scanf("%d",&t);
	while(t--){
		int n;
		scanf("%d%d",&w,&h);
		scanf("%d%d",&x,&y);
		scanf("%d",&n);
		a[1] = get_d(make_pair(0,0));
		a[2] = get_d(make_pair(w,0));
		a[3] = get_d(make_pair(w,h));
		a[4] = get_d(make_pair(0,h));
		// for(int i = 1;i <= 4;++i){
		// 	cout<<a[i]<<endl;
		// }
		for(int i = 1;i <= n;++i){
			pai tp;
			scanf("%d%d",&tp.first,&tp.second);
			d[i] = get(get_d(tp));
			//cout<<"i:"<<i<<" "<< d[i]<<endl;
		}
		sort(d+1,d+1+n);
		double ans = 2.0*(w+h);
		ans -= (d[n] - d[1]);
		for(int i = 2;i <= n;++i){
			ans = max(ans,d[i] - d[i-1]);
		}
		printf("Case #%d: %.8lf\n",++T,ans);
	}
	return 0;
}