%I #8 May 22 2016 12:57:01
%S 1,4,7,16,13,28,19,64,49,52,31,112,37,76,91,256,49,196,55,208,133,124,
%T 67,448,169,148,343,304,85,364,91,1024,217,196,247,784,109,220,259,
%U 832,121,532,127,496,637,268,139,1792,361,676
%N Totally multiplicative sequence with a(p) = 3p-2 for prime p.
%H G. C. Greubel, <a href="/A166669/b166669.txt">Table of n, a(n) for n = 1..10000</a>
%F Multiplicative with a(p^e) = (3p-2)^e. If n = Product p(k)^e(k) then a(n) = Product (3*p(k)-2)^e(k).
%t a[1] = 1; a[n_] := (fi = FactorInteger[n]; Times @@ ((3*fi[[All, 1]] - 2)^fi[[All, 2]])); Table[a[n], {n, 1, 100}] (* _G. C. Greubel_, May 22 2016 *)
%K nonn,mult
%O 1,2
%A _Jaroslav Krizek_, Oct 18 2009
|