BZOJ 2458 BeiJing2011 最小三角形 計算幾何+分治
阿新 • • 發佈:2019-02-03
題目大意:給定平面上的一個點集,求這個點集所能組成的周長最小的三角形
與平面最近點對一個道理- - 這個題也是分治做法
做法如下:
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; }