分治106Pts MLE求助
查看原帖
分治106Pts MLE求助
616964
Adolfo_North楼主2023/7/27 21:40

RT,P1429 平面最近点对(加强版)的代码直接交了, #7~#50 MLE,#6数组越界RE。求调。

#include<bits/stdc++.h>
using namespace std;
int n;
struct node{
	double x,y;
}a[400001];
bool cmp(node a,node b){
	if(a.x==b.x) return a.y<b.y;
	else return a.x<b.x;
}
bool cmp2(int a2,int b){
	return a[a2].y<a[b].y;
}
double xy(int x,int y){
	double xx=(a[y].x-a[x].x)*(a[y].x-a[x].x);
	double yy=(a[y].y-a[x].y)*(a[y].y-a[x].y);
	return sqrt(xx+yy);
}
double merg(int l,int r){
	vector<int> v;
	double ans=0x7ffffff;
	if(l==r) return ans;
	else if(l+1==r) return xy(l,r);
	int mid=l+r>>1;
	double minm=min(merg(l,mid),merg(mid+1,r));
	for(int i=l;i<=r;i++) if(fabs(a[mid].x-a[i].x)<minm) v.push_back(i);
	sort(v.begin(),v.end(),cmp2);
	int sizev=v.size();
	for(int i=0;i<sizev;i++){
		for(int j=i+1;j<sizev&&fabs(a[v[i]].y-a[v[j]].y)<minm;j++){
			minm=min(minm,xy(v[i],v[j]));
		}
	}
	return minm;
}
int main(){
	ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
	cin>>n;
	for(int i=1;i<=n;i++) cin>>a[i].x>>a[i].y;
	sort(a+1,a+1+n,cmp);
	double ans=merg(1,n);
	printf("%.0f",ans*ans);
	return 0;
}
2023/7/27 21:40
加载中...