[BZOJ3707][坐标旋转][随机化]圈地

BZOJ3707

分析：
考虑暴力做法：$n^3$枚举选法，然后验证
考虑优化：发现有些枚举是没有必要的：比如枚举到x最小的和x最大的（但有可能答案在其中），所以我们考虑随机化贪心然后取最小值，即每次把坐标轴随机旋转一个角度
我们可以以x为第一关键字排序，然后将点分块，对于每个块内暴力询问，这样在数据随机的情况下表现的比较好
注意在n很小的时候一定要把所有点丢进一个块里，等于是分段暴力（不然你连样例都过不了）

Code：

#include<bits/stdc++.h>
#define db double
using namespace std;
inline int read(){
	int res=0,f=1;char ch=getchar();
	while(!isdigit(ch)) {if(ch=='-') f=-f;ch=getchar();}
	while(isdigit(ch)) {res=(res<<1)+(res<<3)+(ch^48);ch=getchar();}
	return res*f;
}
const int N=1e3+5;
struct point{
	db x,y;
	point(){}
	point(db _x,db _y): x(_x),y(_y) {}
	friend inline point operator - (const point &a,const point &b) {return point(a.x-b.x,a.y-b.y);}
	friend inline db operator * (const point &a,const point &b) {return a.x*b.y-a.y*b.x;}
}p[N];
inline bool operator < (const point &a,const point &b) {if(a.x==b.x) return a.y<b.y;return a.x<b.x;}
int sqr,n,lim;
db ans=1000000000;
inline point rotate(point a,db th) {return point(a.x*cos(th)-a.y*sin(th),a.y*cos(th)+a.x*sin(th));}
inline void calc(int a,int b){
	for(int i=a;i<=b-2;i++)
		for(int j=i+1;j<=b-1;j++)
			for(int k=j+1;k<=b;k++)
				ans=min(ans,abs((p[i]-p[k])*(p[j]-p[k]))/2);
}
inline void work(){
	db t=rand()/10000;
	for(int i=1;i<=n;i++) p[i]=rotate(p[i],t);
	sort(p+1,p+n+1);
	for(int i=1;i<=lim;i++){
		int l=sqr*(i-1)+1,r=min(sqr*i,n);
		calc(l,r);
	}
}
int main(){
	srand(time(0));
	n=read();sqr=sqrt(n)+10;
	lim=n/sqr+(n%sqr!=0);
	for(int i=1;i<=n;i++) p[i].x=read(),p[i].y=read();
	for(int i=1;i<=1000;i++) work();
	printf("%.2lf",ans);
	return 0;
}