32
$\begingroup$

Choose n points randomly from a circle, how to calculate the probability that all the points are in one semicircle? Any hint is appreciated.

$\endgroup$
  • 2
    $\begingroup$ Have I taken too much of a simplistic view on the problem by thinking the probability is $\left\(\dfrac{1}{2}\right\)^n$? $\endgroup$ – Noble. Mar 9 '13 at 1:10
  • 2
    $\begingroup$ @Noble: Yes, you have -- that's the probability that the points are all in one particular semicircle. $\endgroup$ – joriki Mar 9 '13 at 1:10
  • $\begingroup$ @joriki I thought as much, it's a much more interesting problem then! $\endgroup$ – Noble. Mar 9 '13 at 1:11
  • 1
    $\begingroup$ Hint: Start with a point randomly on the circle and draw a diameter from that point. All you got to do now is ensure that rest of the $n-1$ points lie on the same side of the diameter (i.e., on a semi-circle). You can place the $n-1$ points using a coin toss. $\endgroup$ – jay-sun Mar 9 '13 at 1:12
  • 2
    $\begingroup$ @jay-sun I dont think this is the correct way, three points could still be in one semicircle even if the last two are on two different sides of the diameter joining the first point and the center. $\endgroup$ – NECing Mar 9 '13 at 1:14
26
$\begingroup$

A variation on @joriki's answer (and edited with help from @joriki):

Suppose that point $i$ has angle $0$ (angle is arbitrary in this problem) -- essentially this is the event that point $i$ is the "first" or "leading" point in the semicircle. Then we want the event that all of the points are in the same semicircle -- i.e., that the remaining points end up all in the upper halfplane.

That's a coin-flip for each remaining point, so you end up with $1/2^{n-1}$. There's $n$ points, and the event that any point $i$ is the "leading" point is disjoint from the event that any other point $j$ is, so the final probability is $n/2^{n-1}$ (i.e. we can just add them up).

A sanity check for this answer is to notice that if you have either one or two points, then the probability must be 1, which is true in both cases.

$\endgroup$
  • 2
    $\begingroup$ I don't understand this answer. (Strange, since you think it's a variation on mine. :-) I presume the "if" in "if the remaining points end up all in the upper halfplane" is intended to mean "if and only if"? If so, why is that? And why is the final probability simply the number of points times this one probability you calculated? $\endgroup$ – joriki Mar 9 '13 at 1:46
  • 1
    $\begingroup$ @joriki Basically I'm breaking this down into conditional probabilities. The angle around the circle is just an arbitrary assignment, so conditionally pick $i$. All of these conditional probabilities are going to be identical, and there's $n$ of them, so whatever that probability is, multiply it by $n$. That's the easy part. (cont'd...) $\endgroup$ – John Moeller Mar 9 '13 at 1:48
  • 9
    $\begingroup$ I think I see now -- I find it rather confusingly formulated, but if I understand correctly, you mean something like this: The probability of the remaining $n-1$ points being in the semicircle clockwise of a given point is $1/2^{n-1}$. These $n$ events (one for each given point) are disjoint, and exactly one of them has to occur for the points to lie in a semicircle; thus the desired probability is their sum. That's a nice argument :-) $\endgroup$ – joriki Mar 9 '13 at 2:02
  • 4
    $\begingroup$ There's a slight variation of this answer that uses the inherent symmetry: instead of picking $n$ points at random, pick $n$ random diameters of the circle and pick the $n$ points by randomly picking one of the 2 poles of each diameter. By essentially the same argument, you have the probability given by $(2n)/2^n=n/2^{n-1}$. $\endgroup$ – sai Mar 9 '13 at 2:16
  • 1
    $\begingroup$ @Arrow: The events mentioned in the sentence before that: "the remaining $n−1$ points being in the semicircle clockwise of a given point". That defines one event per point, so $n$ events in all, and they're disjoint. $\endgroup$ – joriki Apr 5 '16 at 15:08
7
$\begingroup$

See

https://mathoverflow.net/questions/33112/estimate-probability-0-is-in-the-convex-hull-of-n-random-points

for the general problem (when the points have any distribution that is invariant w.r.t. rotation about the origin) and

https://mathoverflow.net/questions/2014/if-you-break-a-stick-at-two-points-chosen-uniformly-the-probability-the-three-re/2016#2016

for a nice application.

As a curiosity, this answer can be expressed as a product of sines:

Prove that $\prod_{k=1}^{n-1}\sin\frac{k \pi}{n} = \frac{n}{2^{n-1}}$

$\endgroup$
  • $\begingroup$ "when the points have any distribution that is invariant w.r.t. rotation about the origin", actually the theorem proved by J. G. Wendel pointed out in your first link is more general. We only need a distribution symmetric with respect to $0$ and the result apply to any dimension. $\endgroup$ – Gilles Bonnet May 9 '14 at 19:06
  • $\begingroup$ Also, I think it would be nice if you quote directly in your answer the theorem of Wendel (just in order people looking at the present question don't miss the fact this is a specific case of a much more general "classical" result) $\endgroup$ – Gilles Bonnet May 9 '14 at 19:31
5
$\begingroup$

Here's another way to do this:

Divide the circle into $2k$ equal sectors. There are $2k$ contiguous stretches of $k$ sectors each that form a semicircle, and $2k$ slightly shorter contiguous stretches of $k-1$ sectors that almost form a semicircle. The number of the semicircles containing all the points minus the number of slightly shorter stretches containing all the points is $1$ if the points are contained in at least one of the semicircles and $0$ otherwise; that is, it's the indicator variable for the points all being contained in at least one of the semicircles. The probability of an event is the expected value of its indicator variable, which in this case is

$$2k\left(\frac k{2k}\right)^n-2k\left(\frac{k-1}{2k}\right)^n=\frac k{2^{n-1}}\left(1-\left(1-\frac1k\right)^n\right)\;.$$

The limit $k\to\infty$ yields the desired probability:

$$ \lim_{k\to\infty}\frac k{2^{n-1}}\left(1-\left(1-\frac1k\right)^n\right)=\lim_{k\to\infty}\frac k{2^{n-1}}\cdot\frac nk=\frac n{2^{n-1}}\;. $$

$\endgroup$
  • $\begingroup$ Why is $ \lim_{k\to\infty}k2^{-n}\left(1-\left(1-\frac2k\right)^n\right)=\lim_{k\to\infty}k2^{-n}\left(\frac{2n}k\right)$ true? $\endgroup$ – sai Mar 9 '13 at 5:09
  • 2
    $\begingroup$ @sai: Apply the binomial theorem to the $n$-th power -- the first term cancels the $1$, the second term yields $2n/k$ and the remaining terms have more than one inverse power of $k$ and thus go to zero as $k\to\infty$. Note that $n$ is fixed; it's not a question of taking $n$ and $k$ to infinity simultaneously; we're just adding a finite number of terms, so the standard rules for adding convergent sequences apply. $\endgroup$ – joriki Mar 9 '13 at 8:23
  • $\begingroup$ sorry for pinging you, but can we use this method for a quadrant (1/4 circle)? $\endgroup$ – user2838619 Jun 30 '16 at 16:54
  • 1
    $\begingroup$ In the first equation mentioned in the answer, shouldn't we get 1 - 1/k, instead of 1 - 2/k on the RHS? $\endgroup$ – user10 Aug 8 '16 at 18:32
  • $\begingroup$ It seems to me that there are two errors (that happen to given the correct limit) in the first equation. I think should be $$2k\left(\frac k{2k}\right)^n-2k\left(\frac{k-1}{2k}\right)^n= \frac{k}{2^{n-1}}\left(1-\left(1-\frac1k\right)^n\right) \to \frac{n}{2^{n-1}}$$ $\endgroup$ – leonbloy Jan 6 '17 at 1:47
1
$\begingroup$

Bull, 1948, Mathematical Gazette, Vol 32 No 299 (Dec), pp87-88 solves this problem in the context of the broken stick problem (he uses polytopes and relative volumes in his argument). Rushton, 1949, Mathematical Gazette, Vol 33 No 306 (May), pp286-288 points out that the problem can be re-stated in terms of placing points at random on the circumference of a circle. Ruston's answer is the clearest I have seen. Place $n$ points randomly on the circumference. Label them $X_1, X_2, ..., X_n$. Open up the circle at $X_n$ and produce a straight line. Label the line $OX_n$ (where $O$ is the part of the circle previously immediately adjacent to $X_n$). There are $n$ line segments: $OX_1, X_1X_2, ..., X_{n-1}X_n$. Each segment is equally likely to be longer than half the length of $OX_n$ (and thus correspond to greater than a semi-circle of the orginal circle). The probability that the first segment fulfils this condition is the probability that the remaining $n-1$ points lie upon the second half of the line $OX_n$. That is $(\frac{1}{2})^{(n-1)}$. The probability that there is one segment (note there can be at most one) greater than half the length of the circumference is the sum of the probabilities that each particular segment could be so (because these are mutually exclusive): $n(\frac{1}{2})^{(n-1)}$. So, the favorable probability is $1 -n(\frac{1}{2})^{(n-1)}$.

$\endgroup$
1
$\begingroup$

Another simpler approach,

1) Randomly pick $1$ out of $n$ points and call it $A$ : $\binom n1$ ways

2) Starting from $A$, mark another point $B$ on circumference, such that $length(AB) = \frac12(Cirumference)$ [so that $AB$ and $BA$ are two semi-circles]

3) Now out of remaining $(n-1)$ points, each point can lie on either $AB$ or $BA$ with probability $\frac12$

4) For ALL the remaining $(n-1)$ points to lie on EITHER $AB$ OR $BA$ (i.e., all $(n-1)$ lie on same semi-circle), the joint probability is $\frac12*\frac12 ...(n-1) times$ $=$ $(\frac12)^{(n-1)}$

Since #1 above (randomly picking $A$) is an independent event, $\therefore$ $(\frac12)^{(n-1)}$ (expression in #4) will add $\binom n1$ times

$\implies$ Required probability is $\binom n1(\frac12)^{(n-1)}$ $=$ $n(\frac12)^{(n-1)}$

$\endgroup$
  • $\begingroup$ @joriki, Perhaps this would be simpler to understand. I have tried to put the same idea as John Moeller in a different way. $\endgroup$ – Rahul May 25 '16 at 16:05
  • $\begingroup$ logic has flaws at step 4, as it should be $(1/2)^{n-2}$ $\endgroup$ – fizis Dec 2 '16 at 9:23
  • $\begingroup$ @fizis You are right, he is missing a factor of 2 here. $\endgroup$ – James Apr 22 '18 at 16:15

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.