Pyrknot Planet Visibility Problem
An alien on planet Pyrknot (a perfect sphere of radius R) observes six planets that appear nightly in uniformly random locations across the sky. The alien says:
"If at a given moment there exists somewhere on Pyrknot that all six planets are visible, then from my friend's position on the surface they have an α probability of also seeing all of the planets."
A tower is then considered. It aggregates visibility: the tower can see every planet that is visible from at least one surface point within distance r of the tower's base.
Task: In the limit where r≪R, the new observation probability is α+β⋅(r/R). Find exact values of α and β.
Pyrknot is a sphere of radius R. A planet at direction d∈S2 (unit vector from the center) is visible from surface point p∈S2 if it is above the local horizon, i.e.:
p⋅d>0
The set of surface points from which planet i is visible is the open hemisphere:
H(di)={p∈S2:p⋅di>0}
The region where all 6 planets are simultaneously visible is:
V=⋂i=16H(di)
The condition "there exists somewhere on Pyrknot that all six planets are visible" is exactly V=∅.
Geometrically, V=∅ iff the six direction vectors d1,…,d6 all lie in some common open hemisphere — since a point p sees all 6 iff p is in the hemisphere "opposite" to all di.
The problem asks for a conditional probability:
α=P(friend sees all 6∣∃p∈S2:all 6 visible from p)
By Bayes' theorem (since "friend sees all 6" ⇒ "someone sees all 6"):
α=P(V=∅)P(friend sees all 6)
The six planets are placed independently and uniformly on S2. For a fixed friend position q, each planet independently satisfies:
P(q⋅di>0)=21
by symmetry (half the sphere is above any horizon). By independence:
P(friend sees all 6)=(21)6=641
We need the probability that 6 uniform random points on S2 all lie in some common open hemisphere. This is a classical result.
Wendel / Cover–Efron Formula: For n i.i.d. uniform points on Sd−1 (a sphere in Rd), the probability they all lie in some open hemisphere is:
p(n,d)=2n−11∑k=0d−1(kn−1)
For our sphere S2 (so d=3) with n=6 planets:
p(6,3)=251[(05)+(15)+(25)]=321(1+5+10)=3216=21
So P(V=∅)=21.
α=1/21/64=321
The tower at base position q aggregates visibility from all surface points within distance r of the base. Letting ε=r/R (the angular radius of the disk), the tower sees planet i iff:
∃p∈D(q,ε):p⋅di>0
where D(q,ε)={p∈S2:∠(q,p)<ε} is the spherical cap of angular radius ε centered at q.
The maximum of p⋅di over p∈D(q,ε) is achieved at the point in the cap closest to di. If θ=∠(q,di):
maxp∈D(q,ε)p⋅di=cos(max(0,θ−ε))
This is positive iff max(0,θ−ε)<π/2, i.e., θ<π/2+ε, i.e.:
q⋅di>−sin(ε)
So: planet i is visible from the tower iff ∠(q,di)<π/2+ε.
Since di is uniform on S2, the probability that a uniform point lies within angle π/2+ε of the fixed point q is:
P(∠(q,d)<2π+ε)=21−cos(π/2+ε)=21+sin(ε)≈21+2ε
Note on 3D vs 2D: In a 2D cross-sectional model (circular sky), the analogous formula gives 21+πε. This is incorrect for our 3D sphere. The correct coefficient is 21, not π1, because we integrate over the spherical measure sinθdθ, not the circular measure dθ.
Since the 6 planets are independent:
P(tower sees all 6)=∏i=16P(q⋅di>−sinε)=(21+sinε)6
Expanding to first order in ε=r/R:
(21+2ε)6=641(1+ε)6≈641(1+6ε)=641+646ε=641+323⋅Rr
One must verify that the event ¬A (no common viewpoint exists) does not contribute at first order. When one planet is in the thin band q⋅d1∈(−ε,0) and the remaining five are strictly above the equator of q, a small tilt of the hemisphere toward d1 generically includes all six. The probability of the tilt failing (some dj dropping below the horizon) requires two planets near the equator simultaneously, contributing O(ε2). Therefore, the correction from ¬A does not affect the first-order term.
P(tower sees all 6∣V=∅)≈21641+323⋅Rr=321+163⋅Rr
β=163
α=321,β=163
| Issue | 2D Approach | Correct 3D Approach |
|---|
| Missing conditional | Read off single-planet fraction as α | Divide by P(V=∅)=1/2 |
| Sky geometry | Circular sky → coefficient 1/π | Spherical sky → coefficient 1/2 |
| Number of planets | Computed for 1 planet | Raise to the 6th power |
| α result | 1/2 (wrong) | 1/32 |
| β result | 1/π (wrong) | 3/16 |
The denominator P(V=∅)=1/2 comes from the Wendel formula, a non-trivial combinatorial geometry result. Without it, the conditional probability structure of the problem is missed entirely.