login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 56th year, we are closing in on 350,000 sequences, and we’ve crossed 9,700 citations (which often say “discovered thanks to the OEIS”).

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A100002 Start with a sequence of 1's, then replace every other 1 with a 2; then replace every third of the remaining 1's with a 3 and every third of the remaining 2's with a 3; then replace every fourth remaining 1, 2 or 3 with a 4; and so on. The limiting sequence is shown here. 6
1, 2, 1, 2, 3, 3, 1, 2, 4, 4, 3, 4, 1, 2, 5, 5, 3, 5, 1, 2, 4, 5, 3, 4, 6, 6, 1, 2, 6, 3, 7, 7, 6, 4, 7, 7, 5, 6, 1, 2, 5, 3, 8, 8, 7, 4, 8, 8, 1, 2, 6, 7, 3, 6, 5, 8, 4, 8, 5, 6, 9, 9, 1, 2, 9, 3, 10, 10, 9, 4, 10, 10, 7, 8, 9, 5, 7, 10, 1, 2, 9, 7, 3, 4, 9, 6, 11, 11, 10, 11 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

COMMENTS

The position of the 1's is given by A000960. - T. D. Noe, Oct 26 2004

LINKS

T. D. Noe, Table of n, a(n) for n=1..10000

T. D. Noe, Plot of first 5000 terms

A post on sci.math.research newsgroup.

FORMULA

a(1, j)=1 for all j>=1; a(n, j)=a(n-1, j) except when #{i<=j s.t. a(n-1, i)=a(n-1, j)} is multiple of n, in which case a(n, j)=n; a(j) is the limit of the (stationary) a(n, j) when n tends to infinity.

It appears that the maximal value among the first n terms grows like sqrt(4n/3).

Note that the first occurrence of n is bounded by A000960; that is, A100287(n) <= A000960(n), with equality only for n=1. - T. D. Noe, Nov 12 2004

EXAMPLE

Here are the first 6 stages in the construction:

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1...

1 2 1 2 1 2 1 2 1 2 1 2 1 2 1 2 1 2 1 2 1 2 1 2 1 2 1 2 1 2...

1 2 1 2 3 3 1 2 1 2 3 3 1 2 1 2 3 3 1 2 1 2 3 3 1 2 1 2 3 3...

1 2 1 2 3 3 1 2 4 4 3 4 1 2 1 2 3 3 1 2 4 4 3 4 1 2 1 2 3 3...

1 2 1 2 3 3 1 2 4 4 3 4 1 2 5 5 3 5 1 2 4 5 3 4 1 2 1 2 3 3...

1 2 1 2 3 3 1 2 4 4 3 4 1 2 5 5 3 5 1 2 4 5 3 4 6 6 1 2 6 3...

...

MATHEMATICA

nn=100; t=Table[1, {nn}]; done=False; k=1; While[ !done, k++; cnt=Table[0, {k-1}]; Do[If[t[[i]]<k, cnt[[t[[i]]]]++; If[Mod[cnt[[t[[i]]]], k]==0, t[[i]]=k]], {i, nn}]; done=(Max[cnt]<k)]; t (* T. D. Noe *)

a[n_] := Fold[Function[{b1, b2}, Fold[Function[{a1, a2}, ReplacePart[a1, Pick[Position[a1, a2], Take[Flatten[Array[{Array[0 &, b2 - 1], 1} &, Length[a1]]], Length[Position[a1, a2]]], 1] -> b2]], b1, Range[b2 - 1]]], Array[1 &, n], Range[2, 2 Sqrt[n/Pi] + 1]]; a[100] (* Birkas Gyorgy, Feb 06 2011 *)

PROG

/* C */ #define MAXVAL 2048 /* Large enough... */ unsigned int counts[MAXVAL][MAXVAL]; /* Initialized at all 0 */ unsigned int seq_value (void) /* Successive calls return values in the sequence, in order. */ { unsigned int value; unsigned int i; value = 1; for ( i=2; i<MAXVAL; i++ ) if ( ++counts[i][value] >= i ) { counts[i][value] = 0; value = i; } return value; }

CROSSREFS

Cf. A100287 (first occurrence of n).

Sequence in context: A071766 A007305 A112531 * A348330 A328471 A227909

Adjacent sequences:  A099999 A100000 A100001 * A100003 A100004 A100005

KEYWORD

easy,nice,nonn

AUTHOR

David A. Madore, Oct 25 2004

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
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified November 30 19:02 EST 2021. Contains 349424 sequences. (Running on oeis4.)