hiho一下第125周 GeoHash一·编码解码

最新推荐文章于 2024-03-25 11:33:00 发布

yuhong_liu

最新推荐文章于 2024-03-25 11:33:00 发布

阅读量726

点赞数

分类专栏： hihocoder

本文链接：https://blog.csdn.net/viphong/article/details/53321109

版权

hihocoder 专栏收录该内容

19 篇文章 0 订阅

订阅专栏

题目1 : GeoHash一·编码解码

时间限制: 10000ms

单点时限: 1000ms

内存限制: 256MB

描述

小Hi：上一次我们讲到了在一个城市里，利用四叉树来查找一个坐标附近的点。假如我们把范围扩大到整个地球呢？

小Ho：扩大到整个地球，那坐标怎么办？

小Hi：坐标的话，我们就用经纬度好了。纬度从-90°到90°，经度从-180°到180°。

小Ho：四叉树里面我们用的都是整数坐标，如果变成实数坐标感觉就不那么容易了。

小Hi：没错，而且在扩展到全球的情况下，可以预见点数也会变得非常巨大。因此用四叉树来存储这么多点显然不太现实。

小Ho：那有什么好一点的方法么？

小Hi：当然有了，这一次我们采用编码的方式来解决。

提示：geohash

输入

第1行：2个整数N,M。1≤N,M≤100。

第2..N+1行：每行2个实数x,y，第i+1行表示第i个点的纬度和经度。-90≤x≤90，-180≤y≤180

第N+2..N+M+1行：每行1个字符串，表示需要解码的geohash代码

输出

第1..N行：每行1个字符串，第i行表示以第i个点的geohash编码，长度为10。

第N+1..N+M行：每行2个实数x,y，第N+i行表示以第i个geohash编码所表示的区域中心点坐标，保留6位小数。

样例输入

3 2
36.255833 117.103367
27.293056 112.688922
34.477861 110.082600
wx0csn0ng82y
ww0k7k0et784

样例输出

ww7q2b0jd1
wsb5k299d4
wqnk0ux8n7
39.672792 113.730620
34.519663 112.995297

#include<bits/stdc++.h>
using namespace std;


char Base32[]=
{
    '0','1','2','3','4','5','6','7','8','9','b','c','d','e','f','g','h','j','k','m','n','p','q','r','s','t','u','v','w','x','y','z'
};

string encode(double latt,double lonn,int precision)
{
    double lat[]= {-90,90};
    double lon[]= {-180,180};
    int len =precision*5;
    string geohash;
    int bits=0;
    for(int i=1; i<=len; i++)
    {
        if (i&1)
        {
            double mid=(lon[0]+lon[1])/2;
            if (lonn>mid)
                bits=bits<<1|1,lon[0]=mid;
            else
                bits<<=1,lon[1]=mid;
        }
        else
        {
            double mid=(lat[0]+lat[1])/2;
            if (latt>mid)
                bits=bits<<1|1,lat[0]=mid;
            else
                bits<<=1,lat[1]=mid;
        }
        if (i%5==0)
            geohash+=Base32[bits],bits=0;
    }
    return geohash;

}
int indexofBase32(char s)
{
    return lower_bound(Base32,Base32+32,s)-Base32;
}
void decode(char  s[],double &x,double&y)
{
    bool odd=1;
    double lat[]= {-90,90};
    double lon[]= {-180,180};
    double mid;
    int len=strlen(s);
    for(int i=0; i<len; i++)
    {

        int bits=indexofBase32(s[i]);
        for(int j=4; j>=0; j--)
        {
                int bit=(bits>>j)&1;
            if (odd&1)
            {
                mid=(lon[0]+lon[1])/2;
                lon[1-bit]=mid;
            }
            else
            {
                mid=(lat[0]+lat[1])/2;
                lat[1-bit]=mid;
            }
            odd^=1;
        }
    }
    x=(lat[0]+lat[1] )/2;
    y=(lon[0]+lon[1] )/2;
}
int main()
{
    double x,y;
    int n,m;
    cin>>n>>m;
    for(int i=1; i<=n; i++)
    {
        scanf("%lf%lf",&x,&y);
        string ret=encode(x,y,10);
        printf("%s\n",ret.c_str());
    }
    char ss[15];
    for(int i=1; i<=n; i++)
    {
        scanf("%s",ss);
        decode(ss,x,y);
        printf("%.6lf %.6lf\n",x,y);
    }


}