Tag: conic sections

A note on interpretations

Because the problem doesn’t explain how to define “risk”, there are two main ways to interpret the problem.

1. We can assume the laser cut will happen randomly, and the “risk” of a point is related to the likelihood that the cut will hit it. If the point is a geometrical point and the cut is infinitely thin, then all points have zero probability. If assume the point has a finite size, then the problem can make sense, but we still need to worry about how the cuts are distributed. For example:
   • we could start by choosing the angle of the cut uniformly at random and then pick the valid cut at this angle.
   • we could randomly choose a point on the perimeter and then pick the valid cut that passes through that point.
   • we could pick a point uniformly at random on the inside of the area and then choose a valid cut through that point.

   Each method will yield a different answer! This notion is related to Bertrand’s Paradox.

2. Another way to interpret the problem is to avoid probability altogether and to count the number of different admissible cuts that can pass through a point and use this as a measure of how “unsafe” the point is.

I chose interpretation #2 because it’s a well-defined problem that doesn’t require making additional assumptions. I’ll also present the solution to every polygon, because why not!

Even number of sides

Let’s start with an easy case. If the building is a square, or hexagon, or any other polygon with an even number of sides, then it’s clear by symmetry that all cuts will pass through the center of the polygon. See below for an illustration.

We can all agree that the middle of the building is very unsafe, because all cuts go through it. Meanwhile, each other point in the building has precisely one possible cut passing through it. The same holds true if we have infinitely many sides (a circular building).

Odd number of sides

If the building is a polygon with an odd number of sides instead, such as a triangle or pentagon, then things get more interesting. It turns out that not all cuts through the center divide the area into equal halves. To understand why, take a look at the pentagon below with center at $P$ and short radius $|PA’| = 1$. Extend the top and bottom edges to meet at $O$ as shown. A similar figure can be drawn for any number of sides $n$. In any case, we’ll have $\theta = \tfrac{\pi}{n}$.

Clearly $AA’$ and $BB’$ are valid cuts because by symmetry they split the area of the pentagon in half. Suppose we’d like to make $XY$ a valid cut as well. Then it follows that the area $(OAA’)$ must equal the area $(OYX)$. Define the lengths $x = |OX|$ and $y = |OY|$. Then we must have:

\[
xy = |OA’||OA| = \cot(\tfrac{\theta}{2})\left( \cot(\tfrac{\theta}{2}) + \tan(\theta) \right) = \cot^2(\tfrac{\theta}{2})\sec(\theta)
\]

If $P$ is the origin, the Cartesian coordinates of $X$ and $Y$ are

\[
X = \begin{bmatrix}x \,-\, \cot(\tfrac{\theta}{2}) \\ -1 \end{bmatrix}
\qquad
Y = \begin{bmatrix}y\cos(\theta) \,-\, \cot(\tfrac{\theta}{2}) \\ y\sin(\theta) \,-\, 1 \end{bmatrix}
\]

As $X$ slides from $A’$ to $B’$, $Y$ will slide from $A$ to $B$ in such a way that the product $xy$ remains constant. What we’d like to know is the shape of the curve that the line $XY$ traces out as it assumes all allowable configurations. One way to do this is to compute a parametric equation for the line $XY$. To do this, compute the vector $X + t(Y – X)$. Then, eliminate $y$ to obtain the following family of points:

\[
Z(x,t) = \begin{bmatrix}
(1-t)x + \frac{t}{x} \cot^2(\tfrac{\theta}{2}) \,-\, \cot(\tfrac{\theta}{2}) \\ -1 + \frac{t}{x}\cot^2(\tfrac{\theta}{2})\tan(\theta) \end{bmatrix}
\]

Where $t \in [0,1]$ and $x \in [\cot(\tfrac{\theta}{2}) ,\cot(\tfrac{\theta}{2}) +\tan(\theta)]$. We don’t care about all of these points, however. We only want the boundary of the curve. To find it, fix one coordinate and maximize the other. e.g. let $\eta = -1 + \frac{t}{x}\cot^2(\tfrac{\theta}{2})\tan(\theta)$, solve for $x$ in terms of $\eta$ and $t$, substitute this into the first component of $Z(x,t)$, to obtain an expression that only depends on $\eta$ and $t$. This will be a quadratic in $t$ that is maximized when $t=\tfrac{1}{2}$. Therefore, the equation of the locus of points is simply $Z(x,\tfrac{1}{2})$. If the polygon has $n$ sides, there will be $n$ such loci, each rotated about $P$ by an angle $\frac{2k\pi}{n}$, for $k=0,1,\dots,n-1$. In the figure below, I plotted the result for the case $n=3,5,7$ (click to enlarge).

Here are the diagrams repeated, but this time with some cuts drawn in and zoomed in to the middle portion (click to enlarge).

The blue lines are the boundaries to the regions of interest. Inside each region, all the points have the same number of valid cuts through them. For example, in the triangle case:

• the points inside the blue shape belong to three distinct cuts.
• the points on the boundary of the blue shape belong to two distinct cuts.
• the points outside the blue shape each belong to precisely one cut.

In the figure below, I labeled the case $n=5$ with the number of cuts for each portion.

Calculating the area

The bonus question asks about the number of “high-risk” points. In our interpretation, this would be the points in the innermost blue pentagon (it’s not quite a pentagon since its sides are slightly curved). The points in the middle each have five intersecting cuts. In general, in the case of an $n$-sided polygon ($n$ must be odd, of course), the high-risk area will also be (roughly) an $n$-sided polygon located in the middle.

We can calculate the area of this shape, but I wasn’t able to find an elegant way to do it. Roughly, my approach was as follows:

1. Start with parametric equations $x(t)$ and $y(t)$ for each of the curves.
2. Solve for the intersection points between adjacent rotated versions of the curves (these are the vertices of the inner rounded polygon).
3. Convert these intersection points into an interval $t \in [t_1,t_2]$ that allows our original parametric equations to sweep out one side of the inner rounded polygon.
4. Integrate to find the area of the radial sector of the inner polygon using the cross-product formula for area:
   \[
   A = \int_{t_1}^{t_2} |x'(t) y(t) – x(t) y'(t)|\,dt
   \]
5. Multiply the area of the sector by $n$ to get the total area of the inner rounded polygon.
6. Divide by the total area of the polygon ($n \tan(\theta)$) in our case to get the area ratio.

The details aren’t particularly interesting. Thank goodness for Mathematica… The final result for the area ratio is:

Doesn’t appear to be a way to simplify this. For what it’s worth, this function goes to zero very fast, to the tune of $\tfrac{1}{256}\theta^6$. Here is what it looks like graphically.

The area ratio for the case $n=5$ turns out to be $0.000333883\dots$, or about $\frac{1}{3000}$.

Pretty pictures

Congratulations for making it this far! Here is a neat visualization of the regions. I gave each cut thickness and transparency and used evenly spaced angles. The first picture has about 50 thick cuts while the second picture has about 2000 thin cuts. Click the figure to enlarge.