%I #21 Jan 27 2023 11:22:58
%S 13,37,151,393,1453,2367,5203,7221,12697,25231,30753,52023,70603,
%T 81357,106033,151687,208861,230703,305253,362953,394347,499281,578677,
%U 712891,922083,1040503,1103337,1236493,1306911,1455667,2064513,2265253,2590123,2704941,3330151,3465753,3894543,4357317,4685353,5207647,5767381
%N a(n) = prime(n)^3 + prime(n)^2 + 1.
%H Vincenzo Librandi, <a href="/A181151/b181151.txt">Table of n, a(n) for n = 1..1000</a>
%F a(n) = A135178(n) + 1. - _Elmo R. Oliveira_, Jan 27 2023
%e a(4)=393 because the 4th prime is 7, 7^3=343, 7^2=49, and 343+49+1=393.
%p A181151 := n -> map (p -> p^(3)+p^(2)+1, ithprime(n)):
%p seq (A181151(n), n=1..41);
%t #^3+#^2+1&/@Prime[Range[45]] (* _Harvey P. Dale_, Jan 25 2011 *)
%o (Magma) [p^3+p^2+1: p in PrimesUpTo(700)]; // _Vincenzo Librandi_, Jan 26 2011
%Y Cf. A000040, A001248, A030078, A060800, A135178.
%K nonn,easy
%O 1,1
%A _Jani Melik_, Jan 24 2011
|