1. 程式人生 > >BZOJ 2458 BeiJing2011 最小三角形 計算幾何+分治

BZOJ 2458 BeiJing2011 最小三角形 計算幾何+分治

題目大意:給定平面上的一個點集,求這個點集所能組成的周長最小的三角形

與平面最近點對一個道理- - 這個題也是分治做法

做法如下:

1.記錄全域性答案ans

2.將所有點按照x值排序

3.定義Solve(l,r)為處理[l,r]區間內的最小三角形

4.對於每層Solve(l,r),將當前區間分成左右兩部分,分別遞迴處理

5.兩側的最小三角形都以處理完畢,現在我們要處理的就是兩區間之間的點構成的三角形

6.將本層中與點mid的橫座標之差不超過ans/2的點拎出來,按照縱座標排序 (其實這步可以直接遞迴退出時歸併排序)

7.維護雙指標,保證內部縱座標之差不超過ans/2,對於每個點在範圍內暴力查詢最小三角形即可

時間複雜度期望O(nlogn)...吧?

為何隨便寫了一發就進第一篇了- -

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define M 200200
#define INF 1e10
#define EPS 1e-7
using namespace std;
struct Point{
	int x,y;
	Point() {}
	Point(int _,int __):
		x(_),y(__) {}
	friend istream& operator >> (istream &_,Point &p)
	{
		scanf("%d%d",&p.x,&p.y);
		return _;
	}
	friend double Distance(const Point &p1,const Point &p2)
	{
		return sqrt( (double)(p1.x-p2.x)*(p1.x-p2.x) + (double)(p1.y-p2.y)*(p1.y-p2.y) );
	}
}points[M],_points[M];
bool Compare1(const Point &p1,const Point &p2)
{
	return p1.x < p2.x;
}
bool Compare2(const Point &p1,const Point &p2)
{
	return p1.y < p2.y;
}
int n;
double ans=INF;
void Solve(int l,int r)
{
	if(l==r) return ;
	if(l+1==r)
	{
		if( Compare2(points[r],points[l]) )
			swap(points[l],points[r]);
		return ;
	}
	int i,j,k,tail,mid=l+r>>1,temp=points[mid].x;
	Solve(l,mid);Solve(mid+1,r);
	int l1=l,l2=mid+1;
	for(i=l;i<=r;i++)
	{
		if( l2>r || l1<=mid && Compare2(points[l1],points[l2]) )
			_points[i]=points[l1++];
		else
			_points[i]=points[l2++];
	}
	memcpy(points+l,_points+l,sizeof(Point)*(r-l+1) );
	for(i=l,tail=l;i<=r;i++)
		if(abs(points[i].x-temp)<ans/2)
		{
			for(;points[i].y-points[tail].y>ans/2;tail++);
			for(j=tail;j<i-1;j++)
				if(abs(points[j].x-temp)<ans/2)
					if(points[i].y-points[j].y<ans/2)
						for(k=j+1;k<i;k++)
							ans=min(ans,Distance(points[i],points[j])+Distance(points[i],points[k])+Distance(points[j],points[k]) );
		}
}
int main()
{
	
	//freopen("2458.in","r",stdin);
	//freopen("2458.out","w",stdout);
	
	int i;
	cin>>n;
	for(i=1;i<=n;i++)
		cin>>points[i];
	sort(points+1,points+n+1,Compare1);
	Solve(1,n);
	printf("%.6lf\n",ans);
	return 0;
}