%I #16 Sep 08 2022 08:46:15
%S 4,9,125,5047,39916811,6227020813,355687428096017,121645100408832019,
%T 25852016738884976640023,8841761993739701954543616000029,
%U 8222838654177922817725562880000031,13763753091226345046315979581580902400000037
%N a(n) = prime(n)! + prime(n).
%F a(n) = A039716(n) + A000040(n).
%F a(n) = A005095(A000040(n)).
%t Array[Prime[#]! + Prime[#]&, 20]
%t #!+#&/@Prime[Range[20]] (* _Harvey P. Dale_, Aug 04 2016 *)
%o (Magma) [Factorial(p)+p: p in PrimesUpTo(40)];
%o (PARI) lista(nn) = forprime(p=2, nn, print1(p! + p, ", ")); \\ _Altug Alkan_, Jan 23 2016
%Y Cf. A000040, A005095, A039716, A114574.
%K nonn,easy
%O 1,1
%A _Vincenzo Librandi_, Jan 22 2016
|