This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A166591 #21 Feb 12 2023 02:08:02
%S A166591 1,5,6,25,8,30,10,125,36,40,14,150,16,50,48,625,20,180,22,200,60,70,
%T A166591 26,750,64,80,216,250,32,240,34,3125,84,100,80,900,40,110,96,1000,44,
%U A166591 300,46,350,288,130,50,3750,100,320,120,400,56,1080,112,1250,132,160
%N A166591 Totally multiplicative sequence with a(p) = p+3 for prime p.
%H A166591 Vaclav Kotesovec, <a href="/A166591/b166591.txt">Table of n, a(n) for n = 1..10000</a>
%F A166591 Multiplicative with a(p^e) = (p+3)^e.
%F A166591 If n = Product p(k)^e(k) then a(n) = Product (p(k)+3)^e(k).
%F A166591 From _Vaclav Kotesovec_, Feb 11 2023: (Start)
%F A166591 Dirichlet g.f.: Product_{p prime} 1 / (1 - p^(1-s) - 3*p^(-s)).
%F A166591 Dirichlet g.f.: zeta(s-1) * (1 + 3/(2^s - 5)) * Product_{p prime, p>2} (1 + 3/(p^s - p - 3)).
%F A166591 Sum_{k=1..n} a(k) has average order 3 * c * zeta(r-1) * n^r / (5*log(5)), where r = log(5)/log(2) = 2.321928094... and c = Product_{p prime, p>2} (1 + 3/(p^r - p - 3)) = 1.68551448153095... (End)
%t A166591 f[p_, e_] := (p + 3)^e; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* _Vaclav Kotesovec_, Feb 11 2023 *)
%o A166591 (PARI) a(n) = my(f=factor(n)); for (i=1, #f~, f[i,1] += 3); factorback(f); \\ _Michel Marcus_, Jun 09 2014
%o A166591 (PARI) for(n=1, 100, print1(direuler(p=2, n, 1/(1-p*X-3*X))[n], ", ")) \\ _Vaclav Kotesovec_, Feb 10 2023
%K A166591 nonn,mult
%O A166591 1,2
%A A166591 _Jaroslav Krizek_, Oct 17 2009
%E A166591 More terms from _Michel Marcus_, Jun 09 2014