A001912 -id:A001912 - cp's OEIS frontend

A002496 Primes of the form k^2 + 1.

2, 5, 17, 37, 101, 197, 257, 401, 577, 677, 1297, 1601, 2917, 3137, 4357, 5477, 7057, 8101, 8837, 12101, 13457, 14401, 15377, 15877, 16901, 17957, 21317, 22501, 24337, 25601, 28901, 30977, 32401, 33857, 41617, 42437, 44101, 50177

Offset: 1

N. J. A. Sloane

It is conjectured that this sequence is infinite, but this has never been proved.
An equivalent description: primes of form P = (p1*p2*...*pm)^k + 1 where p1..pm are primes and k > 1, since then k must be even for P to be prime.
Also prime = p(n) if A054269(n) = 1, i.e., quotient-cycle-length = 1 in continued fraction expansion of sqrt(p). - Labos Elemer, Feb 21 2001
Also primes p such that phi(p) is a square.
Also primes of form x*y + z, where x, y and z are three successive numbers. - Giovanni Teofilatto, Jun 05 2004
It is a result that goes back to Mirsky that the set of primes p for which p-1 is squarefree has density A, where A = A005596 denotes the Artin constant. More precisely, Sum_{p <= x} mu(p-1)^2 = A*x/log x + o(x/log x) as x tends to infinity. Conjecture: Sum_{p <= x, mu(p-1)=1} 1 = (A/2)*x/log x + o(x/log x) and Sum_{p <= x, mu(p-1)=-1} 1 = (A/2)*x/log x + o(x/log x). - Pieter Moree (moree(AT)mpim-bonn.mpg.de), Nov 03 2003
Also primes of the form x^y + 1, where x > 0, y > 1. Primes of the form x^y - 1 (x > 0, y > 1) are the Mersenne primes listed in A000668(n) = {3, 7, 31, 127, 8191, 131071, 524287, 2147483647, ...}. - Alexander Adamchuk, Mar 04 2007
With the exception of the first two terms {2,5}, the continued fraction (1 + sqrt(p))/2 has period 3. - Artur Jasinski, Feb 03 2010
With the exception of the first term {2}, congruent to 1 (mod 4). - Artur Jasinski, Mar 22 2011
With the exception of the first two terms, congruent to 1 or 17 (mod 20). - Robert Israel, Oct 14 2014
From Bernard Schott, Mar 22 2019: (Start)
These primes are the primitive terms which generate the sequence of integers with only one prime factor and whose Euler's totient is a square: A054755. So this sequence is a subsequence of A054755 and of A039770. Additionally, the terms of this sequence also have a square cototient, so this sequence is a subsequence of A063752 and A054754.
If p prime = n^2 + 1, phi(p) = n^2 and cototient(p) = 1^2.
Except for 3, the four Fermat primes in A019434 {5, 17, 257, 65537}, belong to this sequence; with F_k = 2^(2^k) + 1, phi(F_k) = (2^(2^(k-1)))^2.
See the file "Subfamilies and subsequences" (& I) in A039770 for more details, proofs with data, comments, formulas and examples. (End)
In this sequence, primes ending with 7 seem to appear twice as often as primes ending with 1. This is because those with 7 come from integers ending with 4 or 6, while those with 1 come only from integers ending with 0 (see De Koninck & Mercier reference). - Bernard Schott, Nov 29 2020
The set of odd primes p for which every elliptic curve of the form y^2 = x^3 + d*x has order p-1 over GF(p) for those d with (d,p)=1 and d a fourth power modulo p. - Gary Walsh, Sep 01 2021 [edited, Gary Walsh, Apr 26 2025]

Jean-Marie De Koninck & Armel Mercier, 1001 Problèmes en Théorie Classique des Nombres, Problème 211 pp. 34 and 169, Ellipses, Paris, 2004.
Leonhard Euler, De numeris primis valde magnis (E283), reprinted in: Opera Omnia. Teubner, Leipzig, 1911, Series (1), Vol. 3, p. 22.
G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, 5th ed., Oxford Univ. Press, 1979, th. 17.
Hugh L. Montgomery, Ten Lectures on the Interface Between Analytic Number Theory and Harmonic Analysis, Amer. Math. Soc., 1996, p. 208.
C. Stanley Ogilvy, Tomorrow's Math. 2nd ed., Oxford Univ. Press, 1972, p. 116.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, page 118.
David Wells, The Penguin Dictionary of Curious and Interesting Numbers (Rev. ed. 1997), p. 134.

T. D. Noe, Table of n, a(n) for n = 1..10000
Tewodros Amdeberhan, Luis A. Medina and Victor H. Moll, Arithmetical properties of a sequence arising from an arctangent sum, J. Numb. Theory, Vol. 128, No. 6 (2008), pp. 1807-1846, eq. (1.10).
William D. Banks, John B. Friedlander, Carl Pomerance and Igor E. Shparlinski, Multiplicative structure of values of the Euler function, in High Primes and Misdemeanours: Lectures in Honour of the Sixtieth Birthday of Hugh Cowie Williams (A. Van der Poorten, ed.), Fields Inst. Comm. 41 (2004), pp. 29-47.
Paul T. Bateman and Roger A. Horn, A heuristic asymptotic formula concerning the distribution of prime numbers, Mathematics of Computation, Vol. 16, No. 79 (1962), pp. 363-367.
Joel E. Cohen, Conjectures about Primes and Cyclic Numbers, arXiv:2508.08335 [math.NT], 2025. See p. 11.
Frank Ellermann, Primes of the form (m^2)+1 up to 10^6.
Dan Ismailescu and Yunkyu James Lee, Polynomially growing integer sequences all whose terms are composite, arXiv:2501.04851 [math.NT], 2025. See p. 9.
Leon Mirsky, The number of representations of an integer as the sum of a prime and a k-free integer, Amer. Math. Monthly, Vol. 56, No. 1 (1949), pp. 17-19.
Maxie D. Schmidt, New Congruences and Finite Difference Equations for Generalized Factorial Functions, arXiv:1701.04741 [math.CO], 2017.
Apoloniusz Tyszka, On sets X subset of N for which we know an algorithm that computes a threshold number t(X) in N such that X is infinite if and only if X contains an element greater than t(X), 2019.
Apoloniusz Tyszka and Sławomir Kurpaska, Open problems that concern computable sets X, subset of N, and cannot be formally stated as they refer to current knowledge about X, (2020).
Eric Weisstein's World of Mathematics, Landau's Problems.
Eric Weisstein's World of Mathematics, Near-Square Prime.
Wikipedia, Bateman-Horn Conjecture.
Marek Wolf, Search for primes of the form m^2+1, arXiv:0803.1456 [math.NT], 2008-2010.

Cf. A083844 (number of these primes < 10^n), A199401 (growth constant).
Cf. A001912, A005574, A054964, A062325, A088179, A090693, A141293, A172168.
Cf. A000668 (Mersenne primes), A019434 (Fermat primes).
Subsequence of A039770.
Cf. A010051, subsequence of A002522.
Cf. A237040 (an analog for n^3 + 1).
Cf. A010051, A000290; subsequence of A028916.
Subsequence of A039770, A054754, A054755, A063752.
Primes of form n^2+b^4, b fixed: A243451 (b=2), A256775 (b=3), A256776 (b=4), A256777 (b=5), A256834 (b=6), A256835 (b=7), A256836 (b=8), A256837 (b=9), A256838 (b=10), A256839 (b=11), A256840 (b=12), A256841 (b=13).
Cf. A030430 (primes ending with 1), A030432 (primes ending with 7).

a002496 n = a002496_list !! (n-1)
a002496_list = filter ((== 1) . a010051') a002522_list
-- Reinhard Zumkeller, May 06 2013

[p: p in PrimesUpTo(100000)| IsSquare(p-1)]; // Vincenzo Librandi, Apr 09 2011

select(isprime, [2, seq(4*i^2+1, i= 1..1000)]); # Robert Israel, Oct 14 2014

Select[Range[100]^2+1, PrimeQ]
Join[{2},Select[Range[2,300,2]^2+1,PrimeQ]] (* Harvey P. Dale, Dec 18 2018 *)

isA002496(n) = isprime(n) && issquare(n-1) \\ Michael B. Porter, Mar 21 2010

is_A002496(n)=issquare(n-1)&&isprime(n) \\ For "random" numbers in the range 10^10 and beyond, at least 5 times faster than the above. - M. F. Hasler, Oct 14 2014

# Python 3.2 or higher required
from itertools import accumulate
from sympy import isprime
A002496_list = [n+1 for n in accumulate(range(10**5),lambda x,y:x+2*y-1) if isprime(n+1)] # Chai Wah Wu, Sep 23 2014

# Python 2.4 or higher required
from sympy import isprime
A002496_list = list(filter(isprime, (n*n+1 for n in range(10**5)))) # David Radcliffe, Jun 26 2016

There are O(sqrt(n)/log(n)) terms of this sequence up to n. But this is just an upper bound. See the Bateman-Horn or Wolf papers, for example, for the conjectured for what is believed to be the correct density.
a(n) = 1 + A005574(n)^2. - R. J. Mathar, Jul 31 2015
Sum_{n>=1} 1/a(n) = A172168. - Amiram Eldar, Nov 14 2020
a(n+1) = 4*A001912(n)^2 + 1. - Hal M. Switkay, Apr 03 2022

Formula, reference, and comment from Charles R Greathouse IV, Aug 24 2009

Edited by M. F. Hasler, Oct 14 2014

A005574 Numbers k such that k^2 + 1 is prime.

Original entry on oeis.org

1, 2, 4, 6, 10, 14, 16, 20, 24, 26, 36, 40, 54, 56, 66, 74, 84, 90, 94, 110, 116, 120, 124, 126, 130, 134, 146, 150, 156, 160, 170, 176, 180, 184, 204, 206, 210, 224, 230, 236, 240, 250, 256, 260, 264, 270, 280, 284, 300, 306, 314, 326, 340, 350, 384, 386, 396

Offset: 1

N. J. A. Sloane

Hardy and Littlewood conjectured that the asymptotic number of elements in this sequence not exceeding n is approximately c*sqrt(n)/log(n) for some constant c. - Stefan Steinerberger, Apr 06 2006
Also, nonnegative integers such that a(n)+i is a Gaussian prime. - Maciej Ireneusz Wilczynski, May 30 2011
Apparently Goldbach conjectured that any a > 1 from this sequence can be written as a=b+c where b and c are in this sequence (Lemmermeyer link below). - Jeppe Stig Nielsen, Oct 14 2015
No term > 2 can be both in this sequence and in A001105 because of the Aurifeuillean factorization (2*k^2)^2 + 1 = (2*k^2 - 2*k + 1) * (2*k^2 + 2*k + 1). - Jeppe Stig Nielsen, Aug 04 2019

Harvey Dubner, "Generalized Fermat primes", J. Recreational Math., 18 (1985): 279-280.
R. K. Guy, "Unsolved Problems in Number Theory", 3rd edition, A2.
G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, 5th ed., Oxford Univ. Press, 1979, p. 15, Thm. 17.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

T. D. Noe, Table of n, a(n) for n = 1..10000
H. Dubner, Generalized Fermat primes, J. Recreational Math. 18.4 (1985-1986), 279. (Annotated scanned copy)
F. Ellermann, Primes of the form (m^2)+1 up to 10^6.
L. Euler, Lettre CXLIX (to Goldbach), 1752.
L. Euler, De numeris primis valde magnis, Novi Commentarii academiae scientiarum Petropolitanae 9 (1764), pp. 99-153. See pp. 123-125.
R. K. Guy, Letter to N. J. A. Sloane, 1988-04-12 (annotated scanned copy).
F. Lemmermeyer, Primes of the form a^2+1, Math Overflow question (2010).
Eric Weisstein's World of Mathematics, Landau's Problems.
Eric Weisstein's World of Mathematics, Power.
Eric Weisstein's World of Mathematics, Near-Square Prime.
Marek Wolf, Search for primes of the form m^2+1, arXiv:0803.1456 [math.NT], 2008-2010.

Cf. A002522, A001912, A002496, A062325, A090693, A000068, A006314, A006313, A006315, A006316, A056994, A056995, A057465, A057002, A088361, A088362, A226528, A226529, A226530, A251597, A253854, A244150, A243959, A321323.
Other sequences of the type "Numbers k such that k^2 + i is prime": this sequence (i=1), A067201 (i=2), A049422 (i=3), A007591 (i=4), A078402 (i=5), A114269 (i=6), A114270 (i=7), A114271 (i=8), A114272 (i=9), A114273 (i=10), A114274 (i=11), A114275 (i=12).
Cf. A010051, A259645, A295405 (characteristic function).

a005574 n = a005574_list !! (n-1)
a005574_list = filter ((== 1) . a010051' . (+ 1) . (^ 2)) [0..]
-- Reinhard Zumkeller, Jul 03 2015

[n: n in [0..400] | IsPrime(n^2+1)]; // Vincenzo Librandi, Nov 18 2010

Select[Range[350], PrimeQ[ #^2 + 1] &] (* Stefan Steinerberger, Apr 06 2006 *)
Join[{1},2Flatten[Position[PrimeQ[Table[x^2+1,{x,2,1000,2}]],True]]]  (* Fred Patrick Doty, Aug 18 2017 *)

isA005574(n) = isprime(n^2+1) \\ Michael B. Porter, Mar 20 2010

for(n=1, 1e3, if(isprime(n^2 + 1), print1(n, ", "))) \\ Altug Alkan, Oct 14 2015

from sympy import isprime; [print(n, end = ', ') for n in range(1, 400) if isprime(n*n+1)] # Ya-Ping Lu, Apr 23 2025

a(n) = A090693(n) - 1.

a(n) = 2*A001912(n-1) for n > 1. - Jeppe Stig Nielsen, Aug 04 2019

A002384 Numbers m such that m^2 + m + 1 is prime.

Original entry on oeis.org

1, 2, 3, 5, 6, 8, 12, 14, 15, 17, 20, 21, 24, 27, 33, 38, 41, 50, 54, 57, 59, 62, 66, 69, 71, 75, 77, 78, 80, 89, 90, 99, 101, 105, 110, 111, 117, 119, 131, 138, 141, 143, 147, 150, 153, 155, 161, 162, 164, 167, 168, 173, 176, 188, 189, 192, 194, 203, 206, 209, 215

Offset: 1

N. J. A. Sloane

A002383 lists the corresponding primes. - Bernard Schott, Dec 22 2012
This is also the list of bases where 111 represents a prime number. - Christian N. K. Anderson, Mar 28 2013
If d>1 divides m^2 + m + 1, then m + k*d is not in the sequence, for all k>=1. - Gionata Neri, Mar 04 2017

A. J. C. Cunningham, Binomial Factorisations, Vols. 1-9, Hodgson, London, 1923-1929; see Vol. 1, pp. 245-259.
D. H. Lehmer, Guide to Tables in the Theory of Numbers. Bulletin No. 105, National Research Council, Washington, DC, 1941, p. 46.
L. Poletti, Le serie dei numeri primi appartenente alle due forme quadratiche (A) n^2+n+1 e (B) n^2+n-1 per l'intervallo compreso entro 121 milioni, e cioè per tutti i valori di n fino a 11000, Atti della Reale Accademia Nazionale dei Lincei, Memorie della Classe di Scienze Fisiche, Matematiche e Naturali, s. 6, v. 3 (1929), pages 193-218.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Zak Seidov, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)
A. J. C. Cunningham, Binomial Factorisations, Vols. 1-9, Hodgson, London, 1923-1929. [Annotated scans of a few pages from Volumes 1 and 2]
Pantelis A. Damianou, On prime values of cyclotomic polynomials, arXiv:1101.1152 [math.NT], Jan 06 2011.
Oliver Knill, Goldbach for Gaussian, Hurwitz, Octavian and Eisenstein primes, arXiv preprint arXiv:1606.05958 [math.NT], 2016.
Oliver Knill, Some experiments in number theory, arXiv preprint arXiv:1606.05971 [math.NT], 2016.

Cf. A002383, A049407, A049408, A075723, A088503, A110284.

[ n: n in [1..300] | IsPrime(n^2+n+1)] // Vincenzo Librandi, Nov 21 2010

Select[Range@ 216, PrimeQ[#^2 + # + 1] &] (* Michael De Vlieger, Mar 06 2017 *)

is(n)=isprime(n^2+n+1) \\ Charles R Greathouse IV, Jan 21 2014

a(n) = (A088503(n) - 1)/2. - Ray Chandler

Extended by Ray Chandler, Sep 07 2005

A002407 Cuban primes: primes which are the difference of two consecutive cubes.

Original entry on oeis.org

7, 19, 37, 61, 127, 271, 331, 397, 547, 631, 919, 1657, 1801, 1951, 2269, 2437, 2791, 3169, 3571, 4219, 4447, 5167, 5419, 6211, 7057, 7351, 8269, 9241, 10267, 11719, 12097, 13267, 13669, 16651, 19441, 19927, 22447, 23497, 24571, 25117, 26227, 27361, 33391

Offset: 1

N. J. A. Sloane

Primes of the form p = (x^3 - y^3)/(x - y) where x=y+1. See A007645 for generalization. I first saw the name "cuban prime" in Cunningham (1923). Values of x are in A002504 and y are in A111251. - N. J. A. Sloane, Jan 29 2013
Prime hex numbers (cf. A003215).
Equivalently, primes of the form p=1+3k(k+1) (and then k=floor(sqrt(p/3))). Also: primes p such that n^2(p+n) is a cube for some n>0. - M. F. Hasler, Nov 28 2007
Primes p such that 4p = 1+3s^2 for some integer s (A121259). - Michael Somos, Sep 15 2005
This sequence is believed to be infinite. - N. J. A. Sloane, May 07 2020

			a(1) = 7 = 1+3k(k+1) (with k=1) is the smallest prime of this form.
a(10^5) = 1792617147127 since this is the 100000th prime of this form.

Allan Joseph Champneys Cunningham, On quasi-Mersennian numbers, Mess. Math., 41 (1912), 119-146.
Allan Joseph Champneys Cunningham, Binomial Factorisations, Vols. 1-9, Hodgson, London, 1923-1929; see Vol. 1, pp. 245-259.
J.-M. De Koninck & A. Mercier, 1001 Problèmes en Théorie Classique des Nombres, Problem 241 pp. 39; 179, Ellipses Paris 2004.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

David A. Corneth, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)
A. J. C. Cunningham, On quasi-Mersennian numbers, Mess. Math., 41 (1912), 119-146. [Annotated scan of page 144 only]
A. J. C. Cunningham, Binomial Factorisations, Vols. 1-9, Hodgson, London, 1923-1929. [Annotated scans of a few pages from Volumes 1 and 2]
R. K. Guy, Letter to N. J. A. Sloane, 1987.
G. L. Honaker, Jr., Prime curio for 127.
Michael Penn, Nearly cubic primes., YouTube video, 2023.
Project Euler, Problem 131: Prime cube partnership.
Eric Weisstein's World of Mathematics, Cuban Prime
Wikipedia, Cuban prime.

Cf. A000217, A002504, A003215, A111251, A113478, A145203.

Cf. A002648 (with x=y+2), A003627, A007645, A201477, A334520.

[a: n in [0..100] | IsPrime(a) where a is (3*n^2+3*n+1)]; // Vincenzo Librandi, Jan 20 2020

lst={}; Do[If[PrimeQ[p=(n+1)^3-n^3], AppendTo[lst, p]], {n, 10^2}]; lst (* Vladimir Joseph Stephan Orlovsky, Aug 21 2008 *)
Select[Table[3x^2+3x+1,{x,100}],PrimeQ] (* or *) Select[Last[#]- First[#]&/@ Partition[Range[100]^3,2,1],PrimeQ] (* Harvey P. Dale, Mar 10 2012 *)
Select[Differences[Range[100]^3],PrimeQ] (* Harvey P. Dale, Jan 19 2020 *)

{a(n)= local(m, c); if(n<1, 0, c=0; m=1; while( cMichael Somos, Sep 15 2005 */

A002407(n,k=1)=until(isprime(3*k*k+++1) && !n--,);3*k*k--+1
list_A2407(Nmax)=for(k=1,sqrt(Nmax/3),isprime(t=3*k*(k+1)+1) && print1(t",")) \\ M. F. Hasler, Nov 28 2007

from sympy import isprime
def aupto(limit):
    alst, k, d = [], 1, 7
    while d <= limit:
        if isprime(d): alst.append(d)
        k += 1; d = 1+3*k*(k+1)
    return alst
print(aupto(34000)) # Michael S. Branicky, Jul 19 2021

a(n) = 6*A000217(A111251(n)) + 1. - Christopher Hohl, Jul 01 2019
From Rémi Guillaume, Nov 07 2023: (Start)
a(n) = A003215(A111251(n)).
a(n) = (3*(2*A002504(n) - 1)^2 + 1)/4.
a(n) = (3*A121259(n)^2 + 1)/4.
a(n) = prime(A145203(n)). (End)

More terms from James Sellers, Aug 08 2000

Entry revised by N. J. A. Sloane, Jan 29 2013

A002645 Quartan primes: primes of the form x^4 + y^4, x > 0, y > 0.

Original entry on oeis.org

2, 17, 97, 257, 337, 641, 881, 1297, 2417, 2657, 3697, 4177, 4721, 6577, 10657, 12401, 14657, 14897, 15937, 16561, 28817, 38561, 39041, 49297, 54721, 65537, 65617, 66161, 66977, 80177, 83537, 83777, 89041, 105601, 107377, 119617, 121937

Offset: 1

N. J. A. Sloane

The largest known quartan prime is currently the largest known generalized Fermat prime: The 1353265-digit 145310^262144 + 1 = (145310^65536)^4 + 1^4, found by Ricky L Hubbard. - Jens Kruse Andersen, Mar 20 2011

Primes of the form (a^2 + b^2)/2 such that |a^2 - b^2| is a square. - Thomas Ordowski, Feb 22 2017

			a(1) =   2 = 1^4 + 1^4.
a(2) =  17 = 1^4 + 2^4.
a(3) =  97 = 2^4 + 3^4.
a(4) = 257 = 1^4 + 4^4.

A. J. C. Cunningham, Binomial Factorisations, Vols. 1-9, Hodgson, London, 1923-1929; see Vol. 1, pp. 245-259.
N. D. Elkies, Primes of the form a^4 + b^4, Mathematical Buds, Ed. H. D. Ruderman Vol. 3 Chap. 3 pp. 22-8 Mu Alpha Theta 1984.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Zak Seidov, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)
A. J. C. Cunningham, High quartan factorisations and primes, Messenger of Mathematics 36 (1907), pp. 145-174.
A. J. C. Cunningham, Binomial Factorisations, Vols. 1-9, Hodgson, London, 1923-1929. [Annotated scans of a few pages from Volumes 1 and 2]
Ernest G. Hibbs, Component Interactions of the Prime Numbers, Ph. D. Thesis, Capitol Technology Univ. (2022), see p. 33.

Subsequence of A002313 and of A028916.
Intersection of A004831 and A000040.
Cf. A002646, A000583, A003336, A000290, A256852.

a002645 n = a002645_list !! (n-1)
a002645_list = 2 : (map a000040 $ filter ((> 1) . a256852) [1..])
-- Reinhard Zumkeller, Apr 11 2015

nn = 100000; Sort[Reap[Do[n = a^4 + b^4; If[n <= nn && PrimeQ[n], Sow[n]], {a, nn^(1/4)}, {b, a}]][[2, 1]]]
With[{nn=20},Select[Union[Flatten[Table[x^4+y^4,{x,nn},{y,nn}]]],PrimeQ[ #] && #<=nn^4+1&]] (* Harvey P. Dale, Aug 10 2021 *)

upto(lim)=my(v=List(2),t);forstep(x=1,lim^.25,2,forstep(y=2,(lim-x^4)^.25,2,if(isprime(t=x^4+y^4),listput(v,t))));vecsort(Vec(v)) \\ Charles R Greathouse IV, Jul 05 2011

list(lim)=my(v=List([2]),x4,t); for(x=1,sqrtnint(lim\=1,4), x4=x^4; forstep(y=1+x%2,min(sqrtnint(lim-x4,4), x-1),2, if(isprime(t=x4+y^4), listput(v,t)))); Set(v) \\ Charles R Greathouse IV, Aug 20 2017

A000040 INTERSECTION A003336. - Jonathan Vos Post, Sep 23 2006

A256852(A049084(a(n))) > 1 for n > 1. - Reinhard Zumkeller, Apr 11 2015

More terms from Victoria A Sapko (vsapko(AT)canes.gsw.edu), Nov 07 2002

A006686 Octavan primes: primes of the form p = x^8 + y^8.

Original entry on oeis.org

2, 257, 65537, 2070241, 100006561, 435746497, 815730977, 832507937, 1475795617, 2579667841, 4338014017, 5110698017, 6975822977, 16983628577, 17995718017, 25605764801, 32575757441, 37822859617, 37839636577, 54875880097, 54876264161, 103910985281, 110081078977

Offset: 1

N. J. A. Sloane

The largest known octavan prime is currently the largest known generalized Fermat prime: The 1353265-digit 145310^262144+1 = (145310^32768)^8+1^8, found by Ricky L Hubbard. - Jens Kruse Andersen, Mar 20 2011

			65537 = 1^8 + 4^8.

A. J. C. Cunningham, High quartan factorisations and primes, Messenger of Mathematics, 36, 11 (1907), pp. 145-174.
A. J. C. Cunningham, Binomial Factorisations, Vols. 1-9, Hodgson, London, 1923-1929; see Vol. 1, pp. 245-259.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Charles R Greathouse IV, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)
A. J. C. Cunningham, Binomial Factorisations, Vols. 1-9, Hodgson, London, 1923-1929. [Annotated scans of a few pages from Volumes 1 and 2]

Intersection of A003380 and A000040. Subsequence of A291206.

lst={}; Do[If[PrimeQ[a^8+b^8], AppendTo[lst, a^8+b^8]], {a, 100}, {b, a, 100}]; Sort[lst] (* T. D. Noe *)
Union[Select[Total/@(Tuples[Range[30],2]^8),PrimeQ]] (* Harvey P. Dale, Apr 06 2013 *)

list(lim)=my(v=List([2]),x8,t); for(x=1,sqrtnint(lim\=1,8), x8=x^8; forstep(y=1+x%2,min(sqrtnint(lim-x8,8), x-1),2, if(isprime(t=x8+y^8), listput(v,t)))); Set(v) \\ Charles R Greathouse IV, Aug 20 2017

Corrected and extended by Jud McCranie, Jan 04 2001

A199307 Primes of the form 4n^3 + 1.

Original entry on oeis.org

5, 109, 257, 1373, 2917, 4001, 27437, 62501, 157217, 202613, 237277, 296353, 470597, 629857, 665501, 1492993, 1556069, 1898209, 2456501, 2634013, 3217429, 3322337, 4244833, 5038849, 5180117, 6572129, 10512289, 11453153, 12706093

Offset: 1

N. J. A. Sloane, Nov 05 2011

Dirichlet's theorem on primes in arithmetic progressions tells us, for example, that there are infinitely many primes of the form 4n+1. For primes represented by polynomials of degree greater than 1, the Bateman-Horn paper gives a conjecture on the density.

Vincenzo Librandi, Table of n, a(n) for n = 1..1000
P. Bateman and R. A. Horn, A heuristic asymptotic formula concerning the distribution of prime numbers, Mathematics of Computation, 16 (1962), 363-367.
Wikipedia, Bateman-Horn Conjecture

Cf. A115104, A001912, A002496, A005574.

[ a: n in [0..200] | IsPrime(a) where a is 4*n^3+1 ]; // Vincenzo Librandi, Nov 08 2011

Select[Table[4n^3+1,{n,1100}],PrimeQ] (* Vincenzo Librandi, Aug 01 2012 *)

for(n=1,1e3,if(isprime(t=4*n^3+1),print1(t", "))) \\ Charles R Greathouse IV, Nov 21 2011

A094550 Numbers n such that there are integers a < b with a+(a+1)+...+(n-1) = (n+1)+(n+2)+...+b.

Original entry on oeis.org

4, 6, 9, 11, 14, 15, 16, 17, 19, 21, 22, 23, 24, 25, 26, 29, 30, 31, 32, 34, 35, 36, 38, 39, 40, 41, 43, 44, 46, 48, 49, 50, 51, 52, 53, 54, 56, 57, 59, 61, 64, 66, 68, 69, 70, 71, 72, 74, 76, 77, 79, 81, 82, 83, 84, 86, 87, 89, 91, 93, 94, 95, 96, 97, 98, 99, 100, 101, 104

Offset: 1

T. D. Noe, May 10 2004

Liljestrom shows that n is in this sequence if and only if 4n^2+1 is composite.
Complement of A001912.
From Hermann Stamm-Wilbrandt, Sep 16 2014: (Start)
For n > 1, A047209 is a subset of this sequence [ 4*n^2+1 is divisible by 5 if n is (1 or 4) mod 5].
A092464 is a subset of this sequence [4*n^2+1 is divisible by 13 if n is (4 or 9) mod 13].
The above are for divisibility by 5, 13; notation (1,4,5), (4,9,13). Divisibility by p for a and p-a; notation (a,p-a,p). These are the next tuples: (2,15,17), (6,23,29), (3,34,37), (16,25,41), ... . The corresponding sequences are a subset of this sequence [ 2,15,19,32,36,49,... for (2,15,17) ]. These sequences have no entries in the OEIS yet. For any prime of the form 4*k+1 there is exactly one of these tuples/sequences [solution to 4*a^2+1=0 (mod p)].
For n>1, A000290 (squares) is a subset of this sequence (4,9,16,25,...) [ 4*(m^2)^2+1 is divisible by  m^2+(m+1)^2, tuple (m^2, (m+1)^2, m^2+(m+1)^2) ].
(End)

			6 is in this sequence because 1+2+3+4+5 = 7+8.

Vincenzo Librandi, Table of n, a(n) for n = 1..1000
R. J. Liljestrom and Richard Zucker, Numerical Fulcrums (PowerPoint Format)

Cf. A094551, A094552, A094553.

Cf. A001912.

[n: n in [1..100] |not IsPrime(4*n^2 + 1)]; // Vincenzo Librandi, Sep 27 2012

lst={}; Do[i1=n-1; i2=n+1; s1=i1; s2=i2; While[i1>1 && s1!=s2, If[s1T. D. Noe, Nov 12 2010 *)

cp's OEIS Frontend

A140126 Partial sums of A001912.

Views

Author

Keywords

Examples

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A002496 Primes of the form k^2 + 1.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

Python

Python

Formula

Extensions

A005574 Numbers k such that k^2 + 1 is prime.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Magma

Mathematica

PARI

PARI

Python

Formula

A002384 Numbers m such that m^2 + m + 1 is prime.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

Extensions

A002407 Cuban primes: primes which are the difference of two consecutive cubes.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Mathematica

PARI

PARI

Python

Formula

Extensions

A002645 Quartan primes: primes of the form x^4 + y^4, x > 0, y > 0.

Views