Hilbert Curve Intersections
时间限制:C/C++/Rust/Pascal 1秒,其他语言2秒
空间限制:C/C++/Rust/Pascal 9.765625 M,其他语言19.53125 M
64bit IO Format: %lld

题目描述


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). H n is defined recursively in terms of H n-1, for n = 2, 3, ... by four transformations:

1.Halve all coordinates in H n-1.
2.Add a copy rotated 90 degrees counterclockwise about the point (0, 1/2).
3.Add the reflection across the line x = 1/2.
4.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.

The coordinates of vertices of Hn are odd multiples of 1/2 n+1. The coordinates of horizontal segment endpoints will always be multiples of 1/2 n. Hence the specified horizontal segment can only cross vertical segments in H n.

输入描述:

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 H n 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.

输出描述:

The output is one integer per line for each data set: the number of intersections of H n 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.
示例1

输入

复制
3 2 7 7
4 0 16 1
30 1 1073741823 1
0

输出

复制
3
16
1073741822