Capital Idea 2 - So Far Away

Rather than query Google Maps to get driving distances between capitals, I found latitude and longitude and use these functions to determine the distance, then ran it against all thousand-some possible distances, so it is ultimately a hash check rather than a calculation to find a distance.

view plainprint?

1. sub haversine {  
2.     my ( $lat1, $lon1, $lat2, $lon2 ) = @_ ;  
3.   
4.     my $theta = $lon1 - $lon2 ;  
5.     my $dist =  
6.         sin( deg2rad( $lat1 ) ) *  
7.         sin( deg2rad( $lat2 ) ) +  
8.         cos( deg2rad( $lat1 ) ) *  
9.         cos( deg2rad( $lat2 ) ) *  
10.         cos( deg2rad( $theta ) ) ;  
11.   
12.     $dist = acos( $dist ) ;  
13.     $dist = rad2deg( $dist ) ;  
14.     $dist = $dist * 60 * 1.1515 ;  
15.     return $dist ;  
16.     }  
17.   
18. sub acos {  
19.     my ( $rad ) = @_ ;  
20.     my $ret = atan2( sqrt( 1 - $rad**2 ), $rad ) ;  
21.     return $ret ;  
22.     }  
23.   
24. sub deg2rad {  
25.     my ( $deg ) = @_ ;  
26.     return ( $deg * $pi / 180 ) ;  
27.     }  
28.   
29. sub rad2deg {  
30.     my ( $rad ) = @_ ;  
31.     return ( $rad * 180 / $pi ) ;  
32.     }  

sub haversine {
    my ( $lat1, $lon1, $lat2, $lon2 ) = @_ ;

    my $theta = $lon1 - $lon2 ;
    my $dist =
        sin( deg2rad( $lat1 ) ) *
        sin( deg2rad( $lat2 ) ) +
        cos( deg2rad( $lat1 ) ) *
        cos( deg2rad( $lat2 ) ) *
        cos( deg2rad( $theta ) ) ;

    $dist = acos( $dist ) ;
    $dist = rad2deg( $dist ) ;
    $dist = $dist * 60 * 1.1515 ;
    return $dist ;
    }

sub acos {
    my ( $rad ) = @_ ;
    my $ret = atan2( sqrt( 1 - $rad**2 ), $rad ) ;
    return $ret ;
    }

sub deg2rad {
    my ( $deg ) = @_ ;
    return ( $deg * $pi / 180 ) ;
    }

sub rad2deg {
    my ( $rad ) = @_ ;
    return ( $rad * 180 / $pi ) ;
    }