|
|
A130135
|
|
Numbers k such that prime(k^3) - k is prime.
|
|
2
|
|
|
2, 4, 12, 30, 40, 42, 62, 66, 70, 72, 84, 88, 98, 112, 130, 150, 180, 248, 264, 298, 336, 342, 354, 408, 438, 450, 468, 482, 486, 546, 550, 558, 602, 606, 682, 686, 700, 710, 720, 732, 744, 770, 774, 804, 840, 848, 854, 948, 968, 976, 994, 996, 1048, 1056
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
LINKS
|
|
|
EXAMPLE
|
a(1)=2 because prime(2^3) - 2 = 19 - 2 = 17 (a prime),
a(2)=4 because prime(4^3) - 4 = 311 - 4 = 307 (a prime),
a(3)=12 because prime(12^3) - 12 = 14741 (a prime), ...
|
|
MAPLE
|
ts_pra_3:=proc(n) local i, ans; ans := [ ]: for i from 1 to n do if (isprime(ithprime(i^(3))-i)) then ans := [ op(ans), i ]: fi: od; RETURN(ans) end: ts_pra_3(200);
|
|
MATHEMATICA
|
Select[Range[1100], PrimeQ[Prime[#^3]-#]&] (* Harvey P. Dale, Mar 24 2023 *)
|
|
PROG
|
(PARI) is(k) = isprime(prime(k^3)-k); \\ Jinyuan Wang, Apr 10 2020
|
|
CROSSREFS
|
Cf. A064269 (prime(k) - k is prime), A105962 (prime(k^2) - k is prime).
|
|
KEYWORD
|
nonn,less
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|