|
%I #14 Feb 28 2017 22:35:26
%S 0,0,1,0,0,1,1,1,3,0,0,1,0,0,1,2,2,4,1,1,2,1,1,2,2,2,5,0,0,1,0,0,1,1,
%T 1,3,0,0,1,0,0,1,5,3,4,2,2,3,2,2,3,3,3,6,1,1,2,1,1,2,2,2,4,1,1,2,1,1,
%U 2,3,3,6,2,2,3,2,2,3,3,3,7
%N a(n) = log_3(A002426(n)/numerator(A002426(n)/3^n)).
%C (-1)^a(n) = A186039(n).
%H G. C. Greubel, <a href="/A186038/b186038.txt">Table of n, a(n) for n = 0..1000</a>
%p A186038 := proc(n)
%p local m,ml ;
%p m := A002426(n) ;
%p ml := %/3^n ;
%p m/numer(ml) ;
%p ilog[3](%) ;
%p end proc: # _R. J. Mathar_, Feb 13 2015
%t b[n_] := If[n < 0, 0, 3^n Hypergeometric2F1[1/2, -n, 1, 4/3]];
%t a[n_] := Log[3, b[n]/Numerator[b[n]/3^n]]; Table[a[n], {n,0,50}] (* _G. C. Greubel_, Feb 27 2017 *)
%K nonn,easy
%O 0,9
%A _Paul Barry_, Feb 11 2011
|
