傻逼数学题(math)

傻逼数学题

题目描述

由于乱码,复制不下来,只好截图了

输入格式

第一行一个正整数n

接下来n行每行两个整数x,y,表示第i个点的坐标

输出格式

一行一个数表示最小距离和,保留6位小数

样例输入

4

1 1

2 3

3 3

3 4

样例输出

3.414214

样例解释

数据约定

对于30%的数据,n<=400.

对于100%的数据,3<=n<=200000.

反思

这一题其实算是个原题了,就是昨天的最近点对,只不过我们需要寻找的变成了3个,思路都是差不多的,但是我就是觉得太复杂,所以没有去实现.所以在考场上我就只打了一个n三方的暴力,也就只拿了30pts,所以我们有任务千万别拖欠,一朝拖欠,就要付出惨痛的代价......

思路

就是原题啦,我们只需要在最近点对的基础改几个地方

1.边界条件(l+2==r)因为有3个

2.在越过中间寻找时,注意3重循环

代码

 #include<bits/stdc++.h>
 #define ll long long
 #define FOR(i,a,b) for(register int i=a;i<=b;i++)
 #define ROF(i,a,b) for(register int i=a;i>=b;i--)
 using namespace std;
 int n;
 const ll N=;
 int tmpt[];
 struct ss
 {
     double x;double y;
 }pot[];
 int scan()
 {
     int as=,f=;
     char c=getchar();
     while(c>''||c<''){if(c=='-') f=-;c=getchar();}
     while(c>=''&&c<=''){as=(as<<)+(as<<)+c-'';c=getchar();}
     return as*f;
 }
 bool cmp(ss i,ss j)
 {
     return i.x<j.x;
 }
 bool cmp2(int i,int j)
 {
     return pot[i].y<pot[j].y;
 }
 double dis(int i,int j)
 {
     double f=sqrt(double((pot[i].x-pot[j].x)*(pot[i].x-pot[j].x)+(pot[i].y-pot[j].y)*(pot[i].y-pot[j].y)));
     return f;
 }
 double par(int l,int r)
 {
 //    cout<<"y"<<endl;
     double d=N;
     if(l==r||l+==r) return d;
     if(l+==r) return dis(l,r)+dis(l+,r)+dis(l+,l);
     int mid=(l+r)>>;//除以2
     double d1=par(l,mid);
     double d2=par(mid+,r);
     d=min(d1,d2);
     int k=;
     FOR(i,l,r)
     {
         if(abs(pot[i].x-pot[mid].x)<=d)
         {
             tmpt[++k]=i;//保存位置
         }
     }
     sort(tmpt+,tmpt+k+,cmp2);
     FOR(i,,k)
     {
         FOR(j,i+,k)
         {
             if(pot[tmpt[j]].y-pot[tmpt[i]].y>=d/) break;
             FOR(j1,j+,k)
             {
                 if(pot[tmpt[j1]].y-pot[tmpt[i]].y>=d/) break;
                 double d3=dis(tmpt[i],tmpt[j])+dis(tmpt[i],tmpt[j1])+dis(tmpt[j1],tmpt[j]);
                 d=min(d,d3);
             }
         }
     }
     return d;
 }
 int main()
 {
 //    freopen("copy.in","r",stdin);
 //    freopen("copy.out","w",stdout);
     n=scan();
     FOR(i,,n)
     {
         scanf("%lf%lf",&pot[i].x,&pot[i].y);
     }
     sort(pot+,pot+n+,cmp);
     double h=par(,n);
     printf("%.6lf",h);//求1~n的最近点对
     return ;
 }