uva12097

小刀刺大熊

于 2022-04-03 14:52:23 发布

阅读量307

点赞数

CC 4.0 BY-SA版权

分类专栏：练习题目文章标签： c++

本文链接：https://blog.youkuaiyun.com/seanbill/article/details/123937682

练习题目专栏收录该内容

80 篇文章

订阅专栏

这是一个关于计算几何的问题，程序通过二分查找方法求解。输入包含多个测试用例，每个用例包括圆的数量n和至少需要覆盖的圆的个数f，接着输入n个圆的半径平方与π的乘积。程序最终找到最小半径，使得至少f+1个圆被覆盖。

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

#include <iostream>
#include <istream>
#include <sstream>
#include <vector>
#include <stack>
#include <list>
#include <map>
#include <set>
#include <deque>
#include <queue>
#include <cstring>
#include <unordered_map>
#include <unordered_set>
#include <algorithm>
#include <numeric>
#include <chrono>
#include <ctime>
#include <cmath>
#include <cctype>
#include <string>
#include <cstdio>
#include <iomanip>


#include <thread>
#include <mutex>
#include <condition_variable>
#include <functional>
#include <iterator>
using namespace std;
const int maxn = 1e4 + 7;
const double PI = acos(-1.0);
int n, f;
double input[maxn];


bool ok(double mid) {
	int sum = 0;
	for (int i = 0; i < n; i++)
		sum += floor(input[i] / mid);
	return sum >= f + 1;
}
int main()
{
	int t;
	cin >> t;
	while (t--) {
		cin >> n >> f;
		double l = 0, r = 0,t;
		for (int i = 0; i < n; i++) {
			cin >> t;
			input[i] = PI * t * t;
			r = max(r, input[i]);
		}
		while (r - l > 1e-5) {
			double mid = (l + r) / 2;
			if (ok(mid)) {
				l = mid;
			}
			else r = mid;
		}
		cout << fixed << setprecision(4) << l << endl;
	}
	return 0;
}