%I #8 Aug 02 2018 11:44: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
|