A181151 a(n) = prime(n)^3 + prime(n)^2 + 1.

13, 37, 151, 393, 1453, 2367, 5203, 7221, 12697, 25231, 30753, 52023, 70603, 81357, 106033, 151687, 208861, 230703, 305253, 362953, 394347, 499281, 578677, 712891, 922083, 1040503, 1103337, 1236493, 1306911, 1455667, 2064513, 2265253, 2590123, 2704941, 3330151, 3465753, 3894543, 4357317, 4685353, 5207647, 5767381

Offset: 1

Author

Jani Melik, Jan 24 2011

Keywords

• nonn
• easy

Examples

			a(4)=393 because the 4th prime is 7, 7^3=343, 7^2=49, and 343+49+1=393.
		

Links

• Vincenzo Librandi, Table of n, a(n) for n = 1..1000

Crossrefs

Cf. A000040, A001248, A030078, A060800, A135178.

Programs

Magma

[p^3+p^2+1: p in PrimesUpTo(700)]; // Vincenzo Librandi, Jan 26 2011

Maple

A181151 := n -> map (p -> p^(3)+p^(2)+1, ithprime(n)):
seq (A181151(n), n=1..41);

Mathematica

#^3+#^2+1&/@Prime[Range[45]]  (* Harvey P. Dale, Jan 25 2011 *)

Formula

a(n) = A135178(n) + 1. - Elmo R. Oliveira, Jan 27 2023