A268752 - OEIS

A268752

Cubefree numbers n such that n^2 + 1 is prime.

1, 2, 4, 6, 10, 14, 20, 26, 36, 66, 74, 84, 90, 94, 110, 116, 124, 126, 130, 134, 146, 150, 156, 170, 180, 204, 206, 210, 230, 236, 260, 284, 300, 306, 314, 326, 340, 350, 386, 396, 406, 420, 430, 436, 444, 466, 470, 474, 490, 556, 570, 634, 636, 644, 646, 654

(list; graph; refs; listen; history; text; internal format)

OFFSET

1,2

COMMENTS

Subsequence of A004709.

Intersection of A004709 and A005574. - Robert Israel, Feb 12 2016

LINKS

K. D. Bajpai, Table of n, a(n) for n = 1..10000

EXAMPLE

a(5) = 10 = 2 * 5 that is cubefree. 10^2 + 1 = 101 which is a prime.

a(6) = 14 = 2 * 7 that is cubefree. 14^2 + 1 = 197 which is a prime.

MAPLE

select(t -> isprime(t^2+1) and max(map(f->f[2], ifactors(t)[2]))<=2, [$1..1000]); # Robert Israel, Feb 12 2016

MATHEMATICA

Select[Range[1000], FreeQ[FactorInteger[#], {_, k_ /; k > 2}] && PrimeQ[#^2 + 1] &]

Select[Range[1000], Max[FactorInteger[#][[;; , 2]]]<3&&PrimeQ[#^2+1]&] (* Harvey P. Dale, May 30 2023 *)

PROG

(PARI) for(n=1, 1000, f = factor(n)[, 2]; if ((#f == 0) || vecmax(f) < 3, if (isprime(n^2+1), print1(n, ", "))));

CROSSREFS

Cf. A000040, A002496, A004709, A005574, A268697.

Sequence in context: A088932 A088954 A000123 * A277277 A241337 A103257

Adjacent sequences: A268749 A268750 A268751 * A268753 A268754 A268755

KEYWORD

nonn

AUTHOR

K. D. Bajpai, Feb 12 2016

STATUS

approved