%I #16 Sep 08 2022 08:45:33
%S 8192,1594323,1220703125,96889010407,34522712143931,302875106592253,
%T 9904578032905937,42052983462257059,504036361936467383,
%U 10260628712958602189,24417546297445042591,243569224216081305397,925103102315013629321,1718264124282290785243
%N a(n) = prime(n)^13.
%H Vincenzo Librandi, <a href="/A138031/b138031.txt">Table of n, a(n) for n = 1..1000</a>
%F From _Amiram Eldar_, Jan 24 2021: (Start)
%F Product_{n>=1} (1 + 1/a(n)) = zeta(13)/zeta(26).
%F Product_{n>=1} (1 - 1/a(n)) = 1/zeta(13) = 1/A013671. (End)
%t Array[Prime[ # ]^13 &, 10]
%t Prime[Range[30]]^13 (* _Vincenzo Librandi_, Mar 27 2014 *)
%o (Magma) [p^13: p in PrimesUpTo(100)]; // _Vincenzo Librandi_, Mar 27 2014
%Y Subsequence of A010801.
%Y Cf. A000040, A013671.
%K nonn,easy
%O 1,1
%A _Vladimir Joseph Stephan Orlovsky_, May 01 2008
%E More terms from _Vincenzo Librandi_, Mar 27 2014