login
Indices m for which A317413(m) = 1.
3

%I #11 Oct 18 2024 18:00:56

%S 1,4,6,9,12,16,19,21,24,28,30,33,37,40,43,45,48,52,54,57,60,64,67,69,

%T 73,76,78,81,85,88,91,93,96,100,102,105,108,112,115,117,120,124,126,

%U 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

%N Indices m for which A317413(m) = 1.

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

%o (Python)

%o n, f, i, p, q, base = 1, 1, 0, 0, 1, 2

%o while i < 10000000:

%o i, p, q = i+1, p*base, q*base

%o if i == f:

%o p, n = p+1, n+1

%o f = f*n

%o n, a, j = 0, 0, 0

%o while p%q > 0:

%o a, f, p, q = a+1, p//q, q, p%q

%o if f == 1:

%o n = n+1

%o print(n, a-1)

%Y Cf. A073089, A317413.

%K nonn

%O 1,2

%A _A.H.M. Smeets_, Jul 30 2018