%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