|
|
A256111
|
|
a(n) = squared distance to the origin of the n-th vertex on a Babylonian Spiral.
|
|
6
|
|
|
0, 1, 5, 13, 26, 50, 65, 85, 116, 100, 97, 85, 85, 90, 128, 205, 293, 409, 481, 586, 730, 845, 890, 841, 833, 745, 514, 244, 65, 17, 106, 338, 698, 1117, 1225, 1193, 1040, 986, 1037, 1060, 850, 477, 197, 85, 80, 232, 530, 757, 650, 522, 225, 16, 50, 333, 797
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
COMMENTS
|
A Babylonian spiral is constructed by starting with a zero vector and progressively concatenating the next longest vector with integral endpoints on a Cartesian grid. (The squares of the lengths of these vectors are A001481.) The direction of the new vector is chosen to create the clockwise spiral which minimizes the change in direction from the previous vector.
. . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . .
. . . 14. . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . .
. . 13. . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . 2 . 3 . . . . . . .
. . . . . . . . . . . 1 . . . . 4 . . . . .
. . . . . . . . . . . o . . . . . . . . . .
. . . . . . . . . . . . . . . . . . 5 . . .
. . 12. . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . 6 . . .
. . . . . . . . . . . . . . . . . . . . . .
. . . . 11. . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . 7 . . . .
. . . . . . . . . . . . . . . . . . . . . .
. . . . . . . 10. . . . . . . . . . . . . .
. . . . . . . . . . . 9 . . . 8 . . . . . .
. . . . . . . . . . . . . . . . . . . . . .
The name is chosen to mislead school students into making an incorrect hypothesis about the Babylonian Spiral's long-term behavior.
|
|
LINKS
|
|
|
EXAMPLE
|
On the above diagram, point 4 is distance sqrt(26) from the origin, so a(4) = 26.
|
|
MATHEMATICA
|
NextVec[{x_, y_}] :=
Block[{n = x^2 + y^2 + 1}, While[SquaresR[2, n] == 0, n++];
TakeSmallestBy[
Union[Flatten[(Transpose[
Transpose[Tuples[{1, -1}, 2]] #] & /@
({{#[[1]], #[[2]]}, {#[[2]], #[[1]]}})) & /@
PowersRepresentations[n, 2, 2], 2]],
Mod[ArcTan[#[[2]], #[[1]]] - ArcTan[y, x], 2 Pi] &, 1][[1]]
]
Norm[#]^2 & /@ Accumulate[NestList[NextVec, {0, 1}, 50]] (* Alex Meiburg, Dec 29 2017 *)
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|