cp's OEIS Frontend

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

%I A317538 #11 Oct 18 2024 18:00:56
%S A317538 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 A317538 73,76,78,81,85,88,91,93,96,100,102,105,108,112,115,117,120,124,126,
%U A317538 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 A317538 Indices m for which A317413(m) = 1.
%F A317538 a(n) = 3*(n-1) + 1 - A073089(n).
%o A317538 (Python)
%o A317538 n, f, i, p, q, base = 1, 1, 0, 0, 1, 2
%o A317538 while i < 10000000:
%o A317538     i, p, q = i+1, p*base, q*base
%o A317538     if i == f:
%o A317538         p, n = p+1, n+1
%o A317538         f = f*n
%o A317538 n, a, j = 0, 0, 0
%o A317538 while p%q > 0:
%o A317538     a, f, p, q = a+1, p//q, q, p%q
%o A317538     if f == 1:
%o A317538         n = n+1
%o A317538         print(n, a-1)
%Y A317538 Cf. A073089, A317413.
%K A317538 nonn
%O A317538 1,2
%A A317538 _A.H.M. Smeets_, Jul 30 2018