P3132 [USACO16JAN] Angry Cows G 题解

老实说，我一开始真止步外层二分了，里层可以二分的性质确实没想到。

思维上的问题还是蛮大的。

来一发很有道理的二分套二分做法。

首先不难发现如果 $R_1$ 是合法的，则 $R_2(R_1<R_2)$ 一定合法，因此可以二分答案。

判定函数怎么写？其实也是一个二分，我们先尝试一个点 $x$,然后分类讨论：

• 若 $x$ 左右都可以覆盖，返回答案合法。
• 若 $x$ 左边无法覆盖，为了覆盖，考虑减小 $x$ 的值。
• 若 $x$ 右边无法覆盖，为了覆盖，考虑增大 $x$ 的值。

二分终止未能找到答案则返回答案不合法。

由此不难写出代码。

#include<bits/stdc++.h>
#define LL long long
#define LD long double
using namespace std;
const LL N=3e5+5;
const LD eps=1e-6;
LL n;
LD a[N],f[N],g[N],ans,l,r;
int check(LD x,LD y)
{
	int t=0;
	while(a[t+1]+eps<x)t++;
	LD r=y,pos=x;
	for(int i=t;i>=1;i--)
	{
		if(pos>r+a[i])return -1;
		while(1<=i&&pos<r+a[i])i--;
		i++,pos=a[i],r-=1.000;
	}
	r=y,pos=x;
	for(int i=t+1;i<=n;i++)
	{
		if(a[i]>r+pos)return 1;
		while(i<=n&&a[i]<r+pos)i++;
		i--,pos=a[i],r-=1.000;
	}
	return 0;
	
}
bool pd(LD x)
{
	LD l=a[1],r=a[n];
	while(l+eps<=r)
	{
		LD mid=(l+r)/2;
		int t=check(mid,x);
		if(t==0)return true;
		if(t==-1)r=mid;
		if(t==1)l=mid;
	}
	return false;
}
int main()
{
	scanf("%lld",&n);
	for(int i=1;i<=n;i++)
	{
		scanf("%Lf",&a[i]);
	}
	sort(a+1,a+n+1);
	l=0.000,r=a[n]-a[1];
	while(l+eps<=r)
	{
		LD mid=(l+r)/2;
		if(pd(mid))
		{
			ans=mid;
			r=mid;
		}
		else l=mid;
	}
	printf("%.1Lf",ans);
	return 0;
}