【poj3714】 Raid

原创 2016年08月30日 10:28:47

http://poj.org/problem?id=3714 (题目链接)

现在才搞平面最近点对。。感觉有点尴尬

题意:给出平面上两组点,每组n个,求两组点之间最短距离

Solution1
  平面最近点对,分治即可。
  将点按横坐标排序,然后每次二分成左边和右边分别计算最小距离,再计算中间的最小距离,这里需要把中间符合条件的点按照纵坐标排序,然后当当前枚举的两点的纵坐标之差大于答案时break,否则会TLE。
  
代码:

// poj3714
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
#include<vector>
#define inf 2147483640
#define LL long long
#define free(a) freopen(a".in","r",stdin);freopen(a".out","w",stdout);
using namespace std;
inline LL getint() {
    LL x=0,f=1;char ch=getchar();
    while (ch>'9' || ch<'0') {if (ch=='-') f=-1;ch=getchar();}
    while (ch>='0' && ch<='9') {x=x*10+ch-'0';ch=getchar();}
    return x*f;
}

const int maxn=1000010;
struct point {double x,y;int flag;}p[maxn];
int n,tmp[maxn];

bool cmpx(point a,point b) {
    return a.x==b.x ? a.y<b.y : a.x<b.x;
}
bool cmpy(int a,int b) {
    return p[a].y==p[b].y ? p[a].x<p[b].x : p[a].y<p[b].y;
}
double dis(point a,point b) {
    return sqrt((double)(a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
double solve(int l,int r) {
    double res=1e60;
    if (l==r) return res;
    if (l+1==r) {
        if (p[l].flag==p[r].flag) return res;
        return dis(p[l],p[r]);
    }
    int mid=(l+r)>>1;
    res=solve(l,mid);
    res=min(res,solve(mid+1,r));
    int num=0;
    for (int i=l;i<=r;i++)
        if (fabs(p[i].x-p[mid].x)<=res) tmp[++num]=i;
    sort(tmp+1,tmp+num+1,cmpy);
    for (int i=1;i<=num;i++)
        for (int j=i+1;j<=num;j++) {
            if (fabs(p[tmp[i]].y-p[tmp[j]].y)>=res) break; //剪枝
            if (p[tmp[i]].flag!=p[tmp[j]].flag) res=min(res,dis(p[tmp[i]],p[tmp[j]]));
        }
    return res;
}
int main() {
    int T;
    scanf("%d",&T);
    while (T--) {
        scanf("%d",&n);
        for (int i=1;i<=n;i++) scanf("%lf%lf",&p[i].x,&p[i].y),p[i].flag=0;
        for (int i=1;i<=n;i++) scanf("%lf%lf",&p[i+n].x,&p[i+n].y),p[i+n].flag=1;
        n<<=1;
        sort(p+1,p+1+n,cmpx);
        printf("%.3f\n",solve(1,n));
    }
    return 0;
}

Solution2
  hzwer上惊现平面最近点对的随机化算法(貌似是随机分块),于是我就蒯了过来,虽然并不知道为什么可以这样写,但是好像很厉害的样子。
  上网搜了下,发现期望复杂度是O(n)的。度娘链接
  然而= =:
  这里写图片描述
  这里写图片描述
  比分治还跑的慢,坑比东西。

代码:

// poj3714
#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<cmath>
#define LL long long
#define inf 2147483640
#define Pi acos(-1.0)
#define free(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout);
using namespace std;
inline LL getint() {
    int f,x=0;char ch=getchar();
    while (ch<='0' || ch>'9') {if (ch=='-') f=-1;else f=1;ch=getchar();}
    while (ch>='0' && ch<='9') {x=x*10+ch-'0';ch=getchar();}
    return x*f;
}

const int maxn=1000010;
struct point {double x,y;int flag;}p[maxn];
int n,block,m;

bool cmp(point a,point b) {
    return a.x==b.x ? a.y<b.y : a.x<b.x;
}
point rotate(point a,double x) {
    return (point){(double)a.x*cos(x)-(double)a.y*sin(x),(double)a.y*cos(x)+(double)a.x*sin(x),a.flag};
}
double dis(point a,point b) {
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
int main() {
    int T;scanf("%d",&T);
    while (T--) {
        scanf("%d",&n);
        for (int i=1;i<=n;i++) scanf("%lf%lf",&p[i].x,&p[i].y),p[i].flag=0;
        for (int i=1;i<=n;i++) scanf("%lf%lf",&p[i+n].x,&p[i+n].y),p[i+n].flag=1;
        n<<=1;
        block=(int)sqrt(n);
        m=n/block+(n%block!=0);
        double t=rand()/10000;
        for (int i=1;i<=n;i++) p[i]=rotate(p[i],t);
        sort(p+1,p+1+n,cmp);
        double ans=1e60;
        for (int i=1;i<=m;i++) {
            int t1=block*(i-1),t2=min(block*i,n);
            for (int j=t1;j<=t2;j++)
                for (int k=t1+1;k<=t2;k++) if (p[j].flag!=p[k].flag) ans=min(ans,dis(p[j],p[k]));
        }
        printf("%.3f\n",ans);
    }
    return 0;
}
版权声明:本文为博主原创文章,未经博主允许不得转载。

POJ 3714 Raid (最近点对)感觉数据怪怪的

传送门题目大意就是要你从两组点钟选出距离最小的两点,输出他们的距离。 一道经典最小点对问题+模板我看了晚上很多人的代码,这题交的时候感觉怪怪的 比如G++wa , C++ac这类情况我的代码也是用...

poj 3714 Raid (类型不同点间的最近点对)

http://poj.org/problem?id=3714   #include "stdio.h" #include "cmath" #include "iostream" #inclu...
  • Risir
  • Risir
  • 2011年08月15日 12:08
  • 179

【分治】【二分】POJ 3714 Raid

After successive failures in the battles against the Union, the Empire retreated to its last strongh...

POJ 3714 Raid(变种最近点对问题:分治+剪枝)

Raid Time Limit: 5000MS   Memory Limit: 65536K Total Submissions: 11728   Accepted...

POJ3714 Raid

最近点对
  • ivy_uu
  • ivy_uu
  • 2016年07月05日 21:59
  • 166

POJ 3714 Raid 最近点对

求最近点对,只不过这两个点需要属于不同的集合,那么就给两个集合的点分别标记一个id号,在计算时,两个集合合并起来,并排序,递归求解,只不过,求两点距离时,如果id号是同一集合的,直接返回一个很大的数就...

poj 3714 Raid【(暴力+剪枝) || (分治法+剪枝)】

题目:  http://poj.org/problem?id=3714 http://acm.hust.edu.cn/vjudge/contest/view.action?cid=27048#pr...

【解题报告】 POJ 3714 Raid -- 点对最小距离 + 分治法

最近点对_分治算法               O(nlgn)思路:对所有点先按x不减排序,               二分x,得到点集S1,点集S2,通过递归求得S1,S2的最小点对距离d1,d2...

POJ 3714 Raid(最近点对)

DescriptionAfter successive failures in the battles against the Union, the Empire retreated to its l...

POJ 3714 Raid(平面最近点对,不同类型点之间)

Raid Time Limit: 5000MS   Memory Limit: 65536K Total Submissions: 7289   Accepted:...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:【poj3714】 Raid
举报原因:
原因补充:

(最多只允许输入30个字)