A237627 -id:A237627 - cp's OEIS frontend

A053698 a(n) = n^3 + n^2 + n + 1.

1, 4, 15, 40, 85, 156, 259, 400, 585, 820, 1111, 1464, 1885, 2380, 2955, 3616, 4369, 5220, 6175, 7240, 8421, 9724, 11155, 12720, 14425, 16276, 18279, 20440, 22765, 25260, 27931, 30784, 33825, 37060, 40495, 44136, 47989, 52060, 56355, 60880

Offset: 0

Henry Bottomley, Mar 23 2000

a(n) = 1111 in base n.

n^3 + n^2 + n + 1 = (n^2 + 1)*(n + 1), therefore a(n) is never prime. - Alonso del Arte, Apr 22 2014

			a(2) = 15 because 2^3 + 2^2 + 2 + 1 = 8 + 4 + 2 + 1 = 15.
a(3) = 40 because 3^3 + 3^2 + 3 + 1 = 27 + 9 + 3 + 1 = 40.
a(4) = 85 because 4^3 + 4^2 + 4 + 1 = 64 + 16 + 4 + 1 = 85.
From _Bruno Berselli_, Jan 02 2017: (Start)
The terms of the sequence are provided by the row sums of the following triangle (see the seventh formula in the previous section):
.   1;
.   3,   1;
.   9,   5,   1;
.  19,  13,   7,   1;
.  33,  25,  17,   9,   1;
.  51,  41,  31,  21,  11,   1;
.  73,  61,  49,  37,  25,  13,  1;
.  99,  85,  71,  57,  43,  29, 15,  1;
. 129, 113,  97,  81,  65,  49, 33, 17,  1;
. 163, 145, 127, 109,  91,  73, 55, 37, 19,  1;
. 201, 181, 161, 141, 121, 101, 81, 61, 41, 21, 1;
...
Columns from the first to the fifth, respectively: A058331, A001844, A056220 (after -1), A059993, A161532. Also, eighth column is A161549.
(End)

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1).

Cf. A237627 (subset of semiprimes).
Cf. A056106 (first differences).
Cf. A062158, A027444, A104878, A055129.

[n^3+n^2+n+1: n in [0..50]]; // Vincenzo Librandi, May 01 2011

A053698:=n->n^3 + n^2 + n + 1; seq(A053698(n), n=0..50); # Wesley Ivan Hurt, Apr 22 2014

Table[n^3 + n^2 + n + 1, {n, 0, 39}] (* Alonso del Arte, Apr 22 2014 *)
FromDigits["1111", Range[0, 50]] (* Paolo Xausa, May 11 2024 *)

Vec((1 + 5*x^2) / (1 - x)^4 + O(x^50)) \\ Colin Barker, Jan 02 2017

def a(n): return (n**3+n**2+n+1) # Torlach Rush, May 08 2024

For n >= 2, a(n) = (n^4-1)/(n-1) = A024002(n)/A024000(n) = A002522(n)*(n+1) = A002061(n+1) + A000578(n).
G.f.: (1+5*x^2) / (1-x)^4. - Colin Barker, Jan 06 2012
a(n) = -A062158(-n). - Bruno Berselli, Jan 26 2016
a(n) = Sum_{i=0..n} 2*n*(n-i)+1. - Bruno Berselli, Jan 02 2017
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4) for n > 3. - Colin Barker, Jan 02 2017
a(n) = A104878(n+3,n) = A055129(4,n) for n > 0. - Mathew Englander, Jan 06 2021
E.g.f.: exp(x)*(x^3+4*x^2+3*x+1). - Nikolaos Pantelidis, Feb 06 2023

A241554 Semiprimes generated by the polynomial 2 * n^2 + 29.

Original entry on oeis.org

1711, 1829, 2077, 2479, 3071, 3901, 5029, 6527, 6757, 7471, 7967, 8479, 10397, 10981, 11581, 14141, 15167, 15517, 15871, 16591, 16957, 17701, 18079, 18847, 19631, 20837, 22927, 23791, 25567, 26941, 27877, 28829, 29797, 30287, 31279, 31781, 32287, 35941, 38117

Offset: 1

K. D. Bajpai, Apr 25 2014

nonn

2 * n^2 + 29 is a well-known Legendre prime-producing polynomial which generates 29 distinct primes for n = 0, 1, ..., 28. For n = 29, it yields the first semiprime, 1711 = 29 * 59.

The number n = 185 is the least positive integer for which 2*n^2 + 29 = 68479 = 31 * 47 * 47 is not squarefree.

			2 * 30^2 + 29 = 1829 = 31 * 59, which is a semiprime and is a term.
2 * 35^2 + 29 = 2479 = 37 * 67, which is a semiprime and is a term.

K. D. Bajpai, Table of n, a(n) for n = 1..10000
R. A. Mollin, Prime-Producing Quadratics, The American Mathematical Monthly, Vol. 104, No. 6 (Jun. - Jul., 1997), pp. 529-544.

Cf. A007641 (for primes).

Cf. A001358, A072381, A082919, A145292, A228183, A237627, A353004, A353388.

with(numtheory):A241554:= proc() local k; k:=2*x^2+29;if bigomega(k)=2 then RETURN (k); fi; end: seq(A241554(), x=0..500);

A241554 = {}; Do[k = 2 * n^2 + 29; If[PrimeOmega[k] == 2, AppendTo[A241554, k]], {n,200}]; A241554

s=[]; for(n=1, 200, t=2*n^2+29; if(bigomega(t)==2, s=concat(s, t))); s \\ Colin Barker, Apr 26 2014

A241659 Primes p such that p^3 + 2 is semiprime.

Original entry on oeis.org

2, 11, 13, 17, 19, 23, 31, 41, 53, 59, 89, 101, 131, 137, 149, 193, 211, 223, 227, 229, 233, 239, 251, 271, 293, 317, 331, 359, 401, 449, 461, 557, 563, 571, 593, 599, 619, 641, 659, 677, 691, 719, 739, 751, 809, 821, 853, 929, 971, 991, 1009, 1013, 1039, 1051

Offset: 1

K. D. Bajpai, Apr 26 2014

nonn

			11 is prime and appears in the sequence because 11^3 + 2 = 1333 = 31 * 43, which is a semiprime.
17 is prime and appears in the sequence because 17^3 + 2 = 4915 =  5 * 983, which is a semiprime.
37 is prime but does not appear in the sequence because 37^3 + 2 = 50655 =  3 * 5 * 11 * 983, which is not a semiprime.

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

Cf. A001358, A063637, A063638, A228183, A237627, A241483, A241493.

with(numtheory): KD:= proc() local a, b, k; k:=ithprime(n); a:=bigomega(k^3+2); if a=2 then RETURN (k); fi; end: seq(KD(), n=1..500);

A241659 = {}; Do[t = Prime[n]; If[PrimeOmega[t^3 + 2] == 2, AppendTo[A241659, t]], {n, 500}]; A241659
(*For the b-file*) c = 0; Do[t = Prime[n]; If[PrimeOmega[t^3 + 2] == 2, c++; Print[c, "  ", t]], {n, 1,6*10^4}];
Select[Prime[Range[200]],PrimeOmega[#^3+2]==2&] (* Harvey P. Dale, Feb 05 2025 *)

s=[]; forprime(p=2, 1200, if(bigomega(p^3+2)==2, s=concat(s, p))); s \\ Colin Barker, Apr 27 2014

A176070 Numbers of the form k^3+k^2+k+1 that are the product of two distinct primes.

Original entry on oeis.org

15, 85, 259, 1111, 4369, 47989, 65641, 291919, 2016379, 2214031, 3397651, 3820909, 5864581, 9305311, 13881841, 15687751, 16843009, 19756171, 22030681, 28746559, 62256349, 64160401, 74264821, 79692331, 101412319, 117889591, 172189309, 185518471, 191435329, 270004099, 328985791

Offset: 1

Vladimir Joseph Stephan Orlovsky, Apr 07 2010

nonn

As k^3 + k^2 + k + 1 = (k + 1) * (k^2 + 1) and k <= 1 does not give a term, k + 1 and k^2 + 1 must be prime so k must be even. - David A. Corneth, May 30 2023

			15 is in the sequence as 15 = 3*5 = 2^3+2^2+2+1; 15 is a product of two distinct primes and of the form k^3 + k^2 + k + 1.

Cf. A002496, A006093, A006881, A053698, A070689, A174969, A176069, A237627 (semiprimes of that form).

f[n_]:=Last/@FactorInteger[n]=={1,1};Select[Array[ #^3+#^2+#+1&,7! ],f[ # ]&]

upto(n) = {my(res = List(), u = sqrtnint(n, 3) + 1); forprime(p = 3, u, c = (p-1)^2 + 1; if(isprime(c), listput(res, c*p))); res} \\ David A. Corneth, May 30 2023

a(n) = (A070689(n + 1) + 1) * (A070689(n + 1)^2 + 1). - David A. Corneth, May 30 2023

Name corrected by and more terms from David A. Corneth, May 30 2023

A241716 Primes p such that p^3 - 2 is semiprime.

Original entry on oeis.org

2, 3, 5, 7, 11, 13, 17, 41, 43, 47, 61, 79, 89, 101, 107, 139, 157, 181, 199, 239, 271, 307, 311, 331, 337, 347, 349, 379, 397, 409, 421, 431, 479, 487, 499, 521, 523, 541, 571, 607, 613, 641, 643, 661, 673, 701, 719, 761, 769, 811, 823, 829, 839, 877, 881, 883

Offset: 1

K. D. Bajpai, Apr 27 2014

nonn

			11 is prime and appears in the sequence because 11^3 - 2 = 1329 = 3 * 443, which is a semiprime.
17 is prime and appears in the sequence because 17^3 - 2 = 4911 = 3 * 1637, which is a semiprime.
23 is prime but does not appear in the sequence because 23^3 - 2 = 12165 =  3 * 5 * 811, which is not a semiprime.

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

Cf. A001358, A063637, A063638, A072381, A082919, A145292, A228183, A237627, A241483, A241493, A241659.

with(numtheory):A241716:= proc() local k; k:=ithprime(x); if bigomega(k^3-2)=2 then RETURN (k); fi; end: seq(A241716(), x=1..500);

A241716 = {}; Do[t = Prime[n]; If[PrimeOmega[t^3 - 2] == 2, AppendTo[A241716, t]], {n, 500}]; A241716
Select[Prime[Range[200]],PrimeOmega[#^3-2]==2&] (* Harvey P. Dale, Dec 09 2018 *)

A241732 Primes p such that p^3 + 2 and p^3 - 2 are semiprime.

Original entry on oeis.org

2, 11, 13, 17, 41, 89, 101, 239, 271, 331, 571, 641, 719, 1051, 1231, 1321, 1549, 1559, 1721, 1741, 1831, 1993, 1999, 2029, 2311, 2459, 2749, 2837, 2861, 2939, 3389, 3467, 3671, 4049, 4111, 4273, 4787, 4919, 4969, 5657, 5689, 5861, 6221, 6679, 6691, 6829, 7109

Offset: 1

K. D. Bajpai, Apr 27 2014

nonn

			11 is prime and appears in the sequence because 11^3 + 2 = 1333 = 31 * 43 and 11^3 - 2 = 1329 = 3 * 443, both are semiprime.
41 is prime and appears in the sequence because 41^3 + 2 = 68923 = 157 * 439 and 41^3 - 2 = 68919 = 3 * 22973, both are semiprime.

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

Cf. A001358, A063637, A063638, A072381, A082919, A145292, A228183, A237627, A241483, A241493, A241659.

with(numtheory): KD:= proc() local k; k:=ithprime(n); if bigomega(k^3+2)=2 and bigomega(k^3-2)=2 then k; fi; end: seq(KD(), n=1..2000);

A241732 = {}; Do[t = Prime[n]; If[PrimeOmega[t^3 + 2] == 2 && PrimeOmega[t^3 - 2] == 2, AppendTo[A241732, t]], {n, 500}]; A241732
Select[Prime[Range[1000]],PrimeOmega[#^3+2]==PrimeOmega[#^3-2]==2&] (* Harvey P. Dale, Jun 20 2019 *)

A241607 Semiprimes generated by the polynomial (1/4)(n^5 - 133n^4 + 6729n^3 - 158379n^2 + 1720294*n - 6823316).

Original entry on oeis.org

5141923, 6084557, 11403823, 13201987, 17488411, 20017609, 33239291, 37446979, 42070423, 47139347, 72512623, 88747907, 118408673, 129881707, 169708339, 184952323, 201267887, 278376073, 324881567, 406044923, 436421497, 538566199, 616639427, 658920007, 750410069

Offset: 1

K. D. Bajpai, Apr 26 2014

nonn

(1/4)*(n^5 - 133*n^4 + 6729*n^3 - 158379*n^2 + 1720294*n - 6823316) is a well known prime producing polynomial found by Shyam Sunder Gupta, which generates 57 distinct primes for n = 0,1,...,55,56.

For n = 57, this polynomial yields the first semiprime: 5141923 = 821 * 6263.

			For n=57: (1/4)*(n^5 - 133*n^4 + 6729*n^3 - 158379*n^2 + 1720294*n - 6823316) = 5141923 = 821 * 6263, which is a semiprime and is included in the sequence.
For n=58: (1/4)*(n^5 - 133*n^4 + 6729*n^3 - 158379*n^2 + 1720294*n - 6823316) = 6084557 = 131 * 46447, which is a semiprime and is included in the sequence.

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

Cf. A007641 (for primes).

Cf. A001358, A072381, A082919, A121887, A145292, A228183, A237627.

with(numtheory): KD:= proc() local a,b,k; k:=(1/4)*(n^5 - 133*n^4 + 6729*n^3 - 158379*n^2 + 1720294*n - 6823316); a:=bigomega(k); if a=2 then RETURN (k); fi; end: seq(KD(), n=0..200);

A241607 = {}; Do[k= (1/4) * (n^5 - 133 * n^4 + 6729 * n^3 - 158379 * n^2 + 1720294 * n - 6823316); If[PrimeOmega[k] ==2, AppendTo[A241607, k]], {n,200}]; A241607
(*For the b-file:*) n=0;Do[t=((1/4) * (k^5 - 133 * k^4 + 6729 * k^3 - 158379 * k^2 + 1720294 * k - 6823316));If[PrimeOmega[t]==2, n++; Print[n," ",t]], {k,10^6}]

s=[]; for(n=1, 200, t=(1/4)*(n^5-133*n^4+6729*n^3-158379*n^2+1720294*n-6823316); if(bigomega(t)==2, s=concat(s, t))); s \\ Colin Barker, Apr 26 2014

cp's OEIS Frontend

A053698 a(n) = n^3 + n^2 + n + 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Python

Formula

A241554 Semiprimes generated by the polynomial 2 * n^2 + 29.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

A241659 Primes p such that p^3 + 2 is semiprime.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

A176070 Numbers of the form k^3+k^2+k+1 that are the product of two distinct primes.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A241716 Primes p such that p^3 - 2 is semiprime.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

A241732 Primes p such that p^3 + 2 and p^3 - 2 are semiprime.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

A241607 Semiprimes generated by the polynomial (1/4)*(n^5 - 133*n^4 + 6729*n^3 - 158379*n^2 + 1720294*n - 6823316).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

A241607 Semiprimes generated by the polynomial (1/4)(n^5 - 133n^4 + 6729n^3 - 158379n^2 + 1720294*n - 6823316).