|
%I #7 May 21 2016 22:34:41
%S 1,17,26,289,44,442,62,4913,676,748,98,7514,116,1054,1144,83521,152,
%T 11492,170,12716,1612,1666,206,127738,1936,1972,17576,17918,260,19448,
%U 278,1419857,2548,2584,2728,195364,332,2890,3016,216172,368,27404,386
%N Totally multiplicative sequence with a(p) = 9p-1 for prime p.
%H G. C. Greubel, <a href="/A166658/b166658.txt">Table of n, a(n) for n = 1..10000</a>
%F Multiplicative with a(p^e) = (9p-1)^e. If n = Product p(k)^e(k) then a(n) = Product (9*p(k)-1)^e(k).
%t a[1] = 1; a[n_] := (fi = FactorInteger[n]; Times @@ ((9*fi[[All, 1]] - 1)^fi[[All, 2]])); Table[a[n], {n, 1, 100}] (* _G. C. Greubel_, May 21 2016 *)
%K nonn,mult
%O 1,2
%A _Jaroslav Krizek_, Oct 18 2009
|
