1. Ha megadunk egy magyarországi települést, akkor megadja, hogy hány településre van Budapesttől, ha egy lépésben légvonalban legfeljebb 30 km távolságra lévő településre utazhatunk.
   
   
2. Ugyanezt a kérdést úgy is oldjuk meg, hogy az légvonalban számolt útnak kell legrövidebbnek lennie, és továbbra is legfeljebb 30 km lehet egy lépés.
   
   
3. Járjuk körbe a magyar megyeszékhelyeket úgy, hogy a körutazás hossza a lehető legkevesebb legyen.
   
   
4. Mind a három részfeladat megoldását jelenítsük meg Google térképen is!

Részfeladatok és megoldásuk

Elmélet és kódok: http://www.movable-type.co.uk/scripts/latlong.html

Clojure

(defn find-distance [lat1 lon1 lat2 lon2]

  (let [radius 6371000

        a (+ (Math/pow (Math/sin (/ (Math/toRadians (- lat2 lat1)) 2)) 2)

             (* (Math/cos (Math/toRadians lat1))

                (Math/cos (Math/toRadians lat2))

                (Math/pow (Math/sin (/ (Math/toRadians (- lon2 lon1)) 2)) 2)))

        c (* 2 (Math/atan2 (Math/sqrt a) (Math/sqrt (- 1 a))))]

    (* radius c)))

Python

import math
# Bemenet varos lat long varos2 lat2 long2 
t = input().strip("\n").split("\t")
lat1 = math.radians(float(t[1]))
long1 = math.radians(float(t[2]))
lat2 = math.radians(float(t[4]))
long2 = math.radians(float(t[5]))
a = math.pow(math.sin((lat1-lat2)/2), 2) + math.cos(lat1)*math.cos(lat2)*math.pow(math.sin((long1-long2)/2), 2)
c = 2 * math.atan2(math.sqrt(a), math.sqrt(1-a))
R = 6371000
d = R*c
print(int(d))

C

#include <math.h>

float rad(float deg) { return deg * M_PI / 180f; }

float fdist(float lat1, float lon1, float lat2, float lon2)

{

    float a =

        sin(rad(lat2 - lat1)) * sin(rad(lat2 - lat1)) +

        cos(rad(lat1)) * cos(rad(lat2)) *

        sin(rad(lon2 - lon1)) * sin(rad(lon2 - lon1));

    float c = 2f * atan2(sqrt(a), sqrt(1.0f - a));

    return 6371000 * c;

}

Kovács Levente megoldása, euklieszi közelítéssel:  megye.szekhely    Korut.java