POJ3122 Pie

题意

n个半径为r[i]的圆，要求切出F+1个面积相等的块（对形状没有要求），求块的面积最大是多少。

思路

令x表示切出来的块的面积。
发现我们很难通过已知信息求得x，但我们可以很容易的判断x是否是一个合法的答案。同时，若x为最优解，则一定切不出F+1个比x面积更大的块；一定能切出F+1个比x面积更小的块。满足二分性。
则解法为：二分块的大小，对每个圆能切出的块的数量进行求和记为sum，比较sum和F+1对二分区间进行转移。

代码

#include <stdio.h>
#include <string.h>
#define max(a,b) (a>b?a:b)
#include <iostream>
#include <iomanip>
using namespace std;

const double pi=31415926.5358979323;
int n,m;
long long s[10001],max_s;
int check(long long x)
{
	if (x==0) return 1;
	int sum=0;
	for (int i=1;i<=n;++i) sum+=s[i]/x;
	if (sum>=m) return 1;
	return 0;
}
int main()
{
	int T;
	for (scanf("%d",&T);T--;)
	{
		scanf("%d%d",&n,&m);
		m++;
		for (int i=1;i<=n;++i) 
			scanf("%lld",&s[i]),
			s[i]=s[i]*s[i]*pi,
			max_s=max(max_s,s[i]);
		long long l=0,r=max_s,mid;
		while (l<r)
		{
			mid=(l+r+1)/2;
			if (check(mid)) l=mid;
			else r=mid-1;
		}
		cout <<fixed<<setprecision(4)<<l/10000000.0<<endl;
	}
}