A317538 Indices m for which A317413(m) = 1.

1, 4, 6, 9, 12, 16, 19, 21, 24, 28, 30, 33, 37, 40, 43, 45, 48, 52, 54, 57, 60, 64, 67, 69, 73, 76, 78, 81, 85, 88, 91, 93, 96, 100, 102, 105, 108, 112, 115, 117, 120, 124, 126, 129, 133, 136, 139, 141, 145, 148, 150, 153, 156, 160, 163, 165, 169, 172, 174, 177, 181, 184, 187, 189, 192, 196, 198, 201, 204

Offset

1,2

Links

Table of n, a(n) for n=1..69.

Formula

a(n) = 3*(n-1) + 1 - A073089(n).

Prog

(Python)
n, f, i, p, q, base = 1, 1, 0, 0, 1, 2
while i < 10000000:
....i, p, q = i+1, p*base, q*base
....if i == f:
........p, n = p+1, n+1
........f = f*n
n, a, j = 0, 0, 0
while p%q > 0:
....a, f, p, q = a+1, p//q, q, p%q
....if f == 1:
........n = n+1
........print(n, a-1)

Crossrefs

Cf. A073089, A317413.

Sequence in context: A144922 A298864 A310667 * A310668 A185595 A310669

Adjacent sequences: A317535 A317536 A317537 * A317539 A317540 A317541

Keyword

nonn

Author

A.H.M. Smeets, Jul 30 2018

Status

approved