login
A350742
Sisyphus with squares: a(1) = 1; thereafter, if a(n-1) is even, a(n) = a(n-1)/2; otherwise a(n) = a(n-1) + k^2, where k^2 is the smallest square (starting from 1) that has not yet been added.
2
1, 2, 1, 5, 14, 7, 23, 48, 24, 12, 6, 3, 39, 88, 44, 22, 11, 75, 156, 78, 39, 139, 260, 130, 65, 209, 378, 189, 385, 610, 305, 561, 850, 425, 749, 1110, 555, 955, 1396, 698, 349, 833, 1362, 681, 1257, 1882, 941, 1617, 2346, 1173, 1957, 2798, 1399, 2299, 3260, 1630
OFFSET
1,2
COMMENTS
Allan C. Wechsler's initial remarks: (Start)
This variation onA350877was proposed byKeith F. Lynchon 2022 January 28. Because the former sequence uses prime increments, and because squares grow much faster than primes, Lynch felt that my conjecture about the former sequence (that all integers eventually appear) would not be true of this sequence. I responded that because squares nevertheless grow much more slowly than exponentials, we should still expect the same basic argument, that descents to the vicinity of 1 should happen on the order of once per "era", to apply, and I still thought that all values would appear, albeit perhaps very slowly.
But later, Gareth McCaughan did the "naive" probabilistic calculation more carefully, and I am now convinced that he and Lynch are correct and that the sequence eventually permanently exceeds any given value, most often without ever equaling that value exactly.
Primes > 2 are all odd, so in the former sequence every ascents will be immediately followed by at least one descent. But squares alternate parity, so in this sequence, ascents (after the first) will always occur in pairs, followed by at least one descent. (End)
LINKS
Michael De Vlieger,Log-log scatterplot of a(n)for n = 1..2^16, annotating the first terms, showing records in red and terms preceded by an even number in blue.
Hans Havermann,A350742 minimain the first 150 billion terms.
EXAMPLE
Given a(10) = 12, the even rule gives a(11) = 6 and a(12) = 3, and then the odd rule must govern. Because the odd rule has already acted 5 times before a(10), we must add 6^2 = 36, so a(13) = 39. Now the odd rule must act again, giving a(14) = a(13) + 7^2 = 39 + 49 = 88.
MATHEMATICA
j = 1; q = 1; {j}~Join~Reap[Do[If[EvenQ[j], k = j/2, k = j + q^2; q++]; Sow[k]; j = k, {i, 55}]][[-1, -1]] (*Michael De Vlieger,Jan 28 2022 *)
PROG
(PARI) k=0; for (n=1, 56, print1 (v=if (n==1, 1, v%2, v+k++^2, v/2) "," )) \\Rémy Sigrist,Jan 29 2022
CROSSREFS
A variation of the Sisyphus sequenceA350877.
KEYWORD
nonn,easy
AUTHOR
Allan C. Wechsler,Jan 28 2022
STATUS
approved