login
This site is supported by donations to The OEIS Foundation.

 

Logo

Annual appeal: Please make a donation to keep the OEIS running! Over 6000 articles have referenced us, often saying "we discovered this result with the help of the OEIS".
Other ways to donate

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A166068 a(n) = a(n-1)+ [least square > a(n-1)]. 1
1, 5, 14, 30, 66, 147, 316, 640, 1316, 2685, 5389, 10865, 21890, 43794, 87894, 176103, 352503, 705339, 1410939, 2822283, 5644683, 11290059, 22586380, 45177389, 90362673, 180726709, 361467845, 722962014, 1445926558, 2891903234 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

COMMENTS

This sequence is the base sequence of the map: a(n) = a(n-1)+ [least square > a(n-1)] if a(n) is not divisible by Y, else a(n)=a(n-1)/Y, where Y is a positive integer.

Experimental results shows this map converges to a periodic orbit for all Y.

What is the number and length of periodic orbits for different Y?

What is the trajectory of some input under the map? If Y=2, the map converges to two periodic orbits, {1-5-14-7-16-8-4-2} and {11-27-63-127-271-560-280-140-70-35-71-152-76-38-19-44-22} whose length is L1=8, L2=17.

Two examples of trajectories for initial value 9 resp. 13 under the map for Y=2 are 9-25-61-125-269-558-279-568-284-142-{76-38-19-44-22-11-27-63-127-271-560-280-140-70-35-71-152} and 13-29-65-146-73-154-77-158-79-160-80-40-20-10-{5-14-7-16-8-4-2-1}.

LINKS

Robert Israel, Table of n, a(n) for n = 1..2656

J. C. Lagarias, The 3x+1 problem and its generalizations, Amer. Math. Monthly, 92 (1985), 3-23.

MAPLE

A[1]:= 1:

for n from 1 to 100 do

  A[n+1]:= A[n] + (floor(sqrt(A[n]))+1)^2

od:

seq(A[n], n=1..100); # Robert Israel, Oct 06 2014

PROG

(PARI) lista(n) = {na = 0; for (i=1, n, na += ceil(sqrt(na+1))^2; print1(na, ", "); ); } \\ Michel Marcus, Jun 02 2013

CROSSREFS

Cf. A006370, A048761

Sequence in context: A211804 A283591 A238604 * A070129 A081861 A267167

Adjacent sequences:  A166065 A166066 A166067 * A166069 A166070 A166071

KEYWORD

nonn

AUTHOR

Ctibor O. Zizka, Oct 06 2009

EXTENSIONS

Typo in data corrected by D. S. McNeil, Aug 17 2010

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent | More pages
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy .

Last modified November 18 09:13 EST 2017. Contains 294879 sequences.