Hilbert Curve Intersections
Time Limit: 1 Second Memory Limit: 32768 KB
Starting curve | Transformation 1 | Transformation 2 | Transformation 3 | Transformation 4 | |
H1 | |||||
H2 | |||||
H3 | |||||
H4 |
David Hilbert proved the existence of a very counter-intuitive curve that fills space. The construction of the Hilbert curve is based on a sequence of curves, H1, H2, H3, H4, ... composed of horizontal and vertical segments. Each curve lies in the unit square [0, 1] - [0, 1]. H1 contains just three segments, connecting the points (1/4, 3/4) to (1/4, 1/4) to (3/4, 1/4) to (3/4, 3/4). Hn is defined recursively in terms of Hn-1, for n = 2, 3, ... by four transformations:
Halve all coordinates in Hn-1.
Add a copy rotated 90 degrees counterclockwise about the point (0, 1/2).
Add the reflection across the line x = 1/2.
Let m = (1/2)^(n+1). Add segments connecting endpoints (1/2 - m, 1/2 - m) to
(1/2 + m, 1/2 - m), (m, 1/2 - m) to (m, 1/2 + m), and (1 - m, 1/2 - m) to (1
- m, 1/2 + m).
Your job is to count the number of intersections of horizontal line segments
with these curves. For example, consider Figures 1 and 2, which illustrate the
first two example input data sets below.
Figure 1
Segment from (2/8, 7/8) to (7/8, 7/8) crossing H3 three times.
Figure 2
Segment from (0/16, 1/16) to (16/16, 1/16) crossing H4 sixteen times.
The coordinates of vertices of Hn are odd multiples of ?n+1. The coordinates of horizontal segment endpoints will always be multiples of ?n. Hence the specified horizontal segment can only cross vertical segments in Hn.
Input
Input consists of one to 100 data sets, one per line, followed by a final line
containing only 0. Each data set consists of four integers separated by blanks
in the form "n x1 x2 y"
which represents Hn and the segment from (x1/2^n, y/2^n) to (x2/2^n, y/2^n),
where 0 < n < 31, x1 < x2, and each of x1, x2, and y lie in the range
0 to 2^n, inclusive.
Output
The output is one integer per line for each data set: the number of intersections of Hn with the segment. Caution: A brute force solution that computes each intersection individually will not finish within the one minute time limit. As you can see below, there may be more than one billion intersections for any data set.
Sample Input
3 2 7 7 4 0 16 1 30 1 1073741823 1 0
Sample Output
3 16 1073741822Submit
Source: Mid-Central USA 2002