OFFSET
1,2
LINKS
Rok Cestnik, Table of n, a(n) for n = 1..10000
Michael De Vlieger, Log log scatterplot of a(n), n = 1..2^16.
Michael De Vlieger, Plot [b(k) > n] at (x,y) = (n,k), for n = 1..2^12, where k may appear at index b(k) and brackets are Iverson. Hence if b(k) > n, k is banned and shown in black else k is permitted and shown in white. This image is similar to the Example but rotated 90 degrees counterclockwise.
EXAMPLE
a(n) ban 1 2 3 4 5 6 7 ...
1 | | | | | | |
2 x | | | | | |
1 | x | | | | |
3 x x | | | | |
2 x | x | | | |
1 | x x | | | |
4 x x x | | | |
3 x x | x | | |
5 x x x x | | |
1 | | x x x | |
2 x | x x x | |
4 x x x | x | |
6 x x x x x | |
7 x x x x x x |
1 | x | x | x x
3 x x | x | x x
5 x x x x | x x
2 x | x x x x x
8 x x x x x x x
6 x x x x x | x
1 | x x | x x x
.
.
.
MATHEMATICA
c[_] := 0; q[_] := 1; r = k = 1; nn = 120; Do[Set[{a[n], c[k]}, {k, n + k*q[k] + 1}]; q[k]++; If[k > r, r = k]; k = MinimalBy[Range[r], c, 1][[1]]; If[c[k] > n + 1, k = r + 1], {n, nn}]; Array[a, nn] (* Michael De Vlieger, Aug 13 2023 *)
PROG
CROSSREFS
KEYWORD
nonn,look
AUTHOR
Rok Cestnik, Jul 29 2023
STATUS
approved