|
%I #8 May 22 2016 12:56:41
%S 1,17,25,289,41,425,57,4913,625,697,89,7225,105,969,1025,83521,137,
%T 10625,153,11849,1425,1513,185,122825,1681,1785,15625,16473,233,17425,
%U 249,1419857,2225,2329,2337,180625,297,2601,2625,201433,329,24225,345
%N Totally multiplicative sequence with a(p) = 8p+1 for prime p.
%H G. C. Greubel, <a href="/A166666/b166666.txt">Table of n, a(n) for n = 1..10000</a>
%F Multiplicative with a(p^e) = (8p+1)^e. If n = Product p(k)^e(k) then a(n) = Product (8*p(k)+1)^e(k).
%t a[1] = 1; a[n_] := (fi = FactorInteger[n]; Times @@ ((8*fi[[All, 1]] + 1)^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
|
