>>835
R: 地球半径
P: 起点の位置ベクトル
Q: 終点の位置ベクトル (平面の交線とか外積でも何でも使って求める)

弧長 = R*arccos(P・Q/(|P||Q|))

もう来んなよ