login
Positions of 0's in A342585.
5

%I #17 Oct 12 2021 21:54:37

%S 1,4,8,14,20,28,37,46,57,69,82,95,110,125,142,159,177,196,216,238,260,

%T 285,310,335,362,390,418,448,478,511,544,578,613,648,685,722,761,800,

%U 842,884,927,971,1018,1065,1112,1161,1210,1259,1309,1361,1413,1467,1521

%N Positions of 0's in A342585.

%C This sequence is infinite.

%H Rémy Sigrist, <a href="/A343880/b343880.txt">Table of n, a(n) for n = 1..10000</a>

%H Rémy Sigrist, <a href="/A343880/a343880.gp.txt">PARI program for A343880</a>

%H Rémy Sigrist, <a href="/A343880/a343880.png">Colored scatterplot of (n, a(n+1)-a(n)) for n = 1..9999</a> (where the color is function of a(n+1)-a(n))

%e A342585(8) = 0, so 8 belongs to the sequence.

%t Position[Block[{c, k, m}, c[0] = 1; {0}~Join~Reap[Do[k = 0; While[IntegerQ[c[k]], Set[m, c[k]]; Sow[m]; If[IntegerQ@ c[m], c[m]++, c[m] = 1]; k++]; Sow[0]; c[0]++, 52]][[-1, -1]]], 0][[All, 1]] (* _Michael De Vlieger_, Oct 12 2021 *)

%o (PARI) See Links section.

%o (Python)

%o from collections import Counter

%o def aupton(terms):

%o num, A342585lst, inventory, alst, idx = 0, [0], Counter([0]), [1], 1

%o while len(alst) < terms:

%o idx += 1

%o c = inventory[num]

%o if c == 0:

%o num = 0

%o alst.append(idx)

%o else:

%o num += 1

%o A342585lst.append(c)

%o inventory.update([c])

%o return alst

%o print(aupton(53)) # _Michael S. Branicky_, Oct 12 2021

%Y Cf. A342585.

%K nonn

%O 1,2

%A _Rémy Sigrist_, May 02 2021