cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 36 results. Next

A246392 Numbers n such that Phi(10, n) is prime, where Phi is the cyclotomic polynomial.

Original entry on oeis.org

2, 3, 5, 10, 11, 12, 16, 20, 21, 22, 33, 37, 38, 43, 47, 48, 55, 71, 75, 76, 80, 81, 111, 121, 126, 131, 133, 135, 136, 141, 155, 157, 158, 165, 176, 177, 180, 203, 223, 242, 245, 251, 253, 256, 257, 258, 265, 268, 276, 286, 290, 297, 307, 322, 323, 342, 361, 363, 366, 375, 377, 385, 388, 396, 411
Offset: 1

Views

Author

Eric Chen, Nov 13 2014

Keywords

Comments

Numbers n such that (n^5+1)/(n+1) is prime, or numbers n such that A060884(n) is prime.

Crossrefs

Cf. A008864 (1), A006093 (2), A002384 (3), A005574 (4), A049409 (5), A055494 (6), A100330 (7), A000068 (8), A153439 (9), this sequence (10), A162862 (11), A246397 (12), A217070 (13), A006314 (16), A217071 (17), A164989 (18), A217072 (19), A217073 (23), A153440 (27), A217074 (29), A217075 (31), A006313 (32), A097475 (36), A217076 (37), A217077 (41), A217078 (43), A217079 (47), A217080 (53), A217081 (59), A217082 (61), A006315 (64), A217083 (67), A217084 (71), A217085 (73), A217086 (79), A153441 (81), A217087 (83), A217088 (89), A217089 (97), A006316 (128), A153442 (243), A056994 (256), A056995 (512), A057465 (1024), A057002 (2048), A088361 (4096), A088362 (8192), A226528 (16384), A226529 (32768), A226530 (65536).

Programs

  • Magma
    [n: n in [1..500]| IsPrime((n^5+1) div (n+1))]; // Vincenzo Librandi, Nov 14 2014
  • Maple
    A246392:=n->`if`(isprime((n^5+1)/(n+1)),n,NULL): seq(A246392(n), n=1..500); # Wesley Ivan Hurt, Nov 15 2014
  • Mathematica
    Select[Range[700], PrimeQ[(#^5 + 1) / (# + 1)] &] (* Vincenzo Librandi, Nov 14 2014 *)
  • PARI
    for(n=1,10^3,if(isprime(polcyclo(10,n)),print1(n,", "))); \\ Joerg Arndt, Nov 13 2014
    

A217070 Numbers k such that (k^13-1)/(k-1) is prime.

Original entry on oeis.org

2, 3, 5, 7, 34, 37, 43, 59, 72, 94, 98, 110, 133, 149, 151, 159, 190, 207, 219, 221, 251, 260, 264, 267, 282, 286, 291, 319, 355, 363, 373, 382, 397, 398, 402, 406, 408, 412, 436, 442, 486, 489, 507, 542, 544, 552, 553, 582, 585, 592, 603, 610, 614, 634, 643
Offset: 1

Views

Author

Tim Johannes Ohrtmann, Sep 26 2012

Keywords

Crossrefs

Programs

  • Magma
    [n: n in [2..1000] |IsPrime((n^13 - 1) div (n - 1))]; // Vincenzo Librandi, Sep 28 2012
    
  • Mathematica
    Select[Range[2, 1000], PrimeQ[(#^13 - 1)/(# - 1)] &] (* T. D. Noe, Sep 26 2012 *)
  • PARI
    is(n)=isprime(polcyclo(13,n)) \\ Charles R Greathouse IV, Apr 28 2015

A217089 Numbers n such that (n^97-1)/(n-1) is prime.

Original entry on oeis.org

12, 90, 104, 234, 271, 339, 420, 421, 428, 429, 464, 805, 909, 934, 1054, 1114, 1116, 1128, 1144, 1159, 1193, 1364, 1788, 2086, 2215, 2254, 2448, 2461, 2593, 2595, 2771, 2787, 2829, 2859, 2952, 3029, 3075, 3144, 3250, 3265, 3268, 3301, 3701, 3752, 3875, 4026
Offset: 1

Views

Author

Tim Johannes Ohrtmann, Sep 26 2012

Keywords

Crossrefs

Programs

  • Mathematica
    Select[Range[2, 1000], PrimeQ[(#^97 - 1)/(# - 1)] &] (* T. D. Noe, Sep 26 2012 *)
  • PARI
    is(n)=isprime((n^97-1)/(n-1)) \\ Charles R Greathouse IV, Feb 17 2017

Extensions

More terms from T. D. Noe, Sep 26 2012

A085398 Let Cn(x) be the n-th cyclotomic polynomial; a(n) is the least k>1 such that Cn(k) is prime.

Original entry on oeis.org

3, 2, 2, 2, 2, 2, 2, 2, 2, 2, 5, 2, 2, 2, 2, 2, 2, 6, 2, 4, 3, 2, 10, 2, 22, 2, 2, 4, 6, 2, 2, 2, 2, 2, 14, 3, 61, 2, 10, 2, 14, 2, 15, 25, 11, 2, 5, 5, 2, 6, 30, 11, 24, 7, 7, 2, 5, 7, 19, 3, 2, 2, 3, 30, 2, 9, 46, 85, 2, 3, 3, 3, 11, 16, 59, 7, 2, 2, 22, 2, 21, 61, 41, 7, 2, 2, 8, 5, 2, 2
Offset: 1

Views

Author

Don Reble, Jun 28 2003

Keywords

Comments

Conjecture: a(n) is defined for all n. - Eric Chen, Nov 14 2014
Existence of a(n) is implied by Bunyakovsky's conjecture. - Robert Israel, Nov 13 2014

Examples

			a(11) = 5 because C11(k) is composite for k = 2, 3, 4 and prime for k = 5.
a(37) = 61 because C37(k) is composite for k = 2, 3, 4, ..., 60 and prime for k = 61.
		

Crossrefs

Programs

  • Maple
    f:= proc(n) local k;
    for k from 2 do if isprime(numtheory:-cyclotomic(n,k)) then return k fi od
    end proc:
    seq(f(n), n = 1 .. 100); # Robert Israel, Nov 13 2014
  • Mathematica
    Table[k = 2; While[!PrimeQ[Cyclotomic[n, k]], k++]; k, {n, 300}] (* Eric Chen, Nov 14 2014 *)
  • PARI
    a(n) = k=2; while(!isprime(polcyclo(n, k)), k++); k; \\ Michel Marcus, Nov 13 2014

Formula

a(A072226(n)) = 2. - Eric Chen, Nov 14 2014
a(n) = A117544(n) except when n is a prime power, since if n is a prime power, then A117544(n) = 1. - Eric Chen, Nov 14 2014
a(prime(n)) = A066180(n), a(2*prime(n)) = A103795(n), a(2^n) = A056993(n-1), a(3^n) = A153438(n-1), a(2*3^n) = A246120(n-1), a(3*2^n) = A246119(n-1), a(6^n) = A246121(n-1), a(5^n) = A206418(n-1), a(6*A003586(n)) = A205506(n), a(10*A003592(n)) = A181980(n).

A250174 Numbers n such that Phi_14(n) is prime, where Phi is the cyclotomic polynomial.

Original entry on oeis.org

2, 3, 10, 11, 14, 15, 16, 17, 18, 21, 24, 25, 29, 37, 43, 44, 46, 49, 52, 54, 61, 66, 72, 73, 78, 84, 86, 87, 99, 101, 106, 114, 115, 128, 133, 135, 136, 143, 145, 148, 164, 169, 170, 173, 200, 219, 224, 226, 228, 231, 234, 240, 248, 255, 262, 275, 281, 282, 298, 301
Offset: 1

Views

Author

Eric Chen, Dec 24 2014

Keywords

Comments

n = 9069 * 2^64163 + 1 is an example of a rather large member of this sequence. The generated 115914 decimal digit prime is proved by the N-1 method (because n is prime and n*(n-1) is fully factored and this provides for an exactly 33.33...% factorization for Phi_14(n) - 1). - Serge Batalov, Mar 13 2015

Examples

			2 is in the sequence because 2^6-2^5+2^4-2^3+2^2-2+1 = 43 which is prime.
		

Crossrefs

See A250177 for cross-references, A100330 (Phi_7(n) = n^6 + n^5 + n^4 + n^3 + n^2 + n + 1 primes; these two sequences can also be considered an extension of each other into negative n values), A250177 (Phi_21(n) primes).

Programs

  • Mathematica
    a250174[n_] := Select[Range[n], PrimeQ@Cyclotomic[14, #] &]; a250174[256]
  • PARI
    isok(n) = isprime(polcyclo(14, n)); \\ Michel Marcus, Mar 13 2015

A088550 Primes of the form n^6 + n^5 + n^4 + n^3 + n^2 + n + 1.

Original entry on oeis.org

7, 127, 1093, 19531, 55987, 5229043, 8108731, 25646167, 321272407, 917087137, 3092313043, 4201025641, 9684836827, 31401724537, 47446779661, 52379047267, 83925549247, 100343116693, 141276239497, 153436090543, 265462278481
Offset: 1

Views

Author

Cino Hilliard, Nov 17 2003

Keywords

Comments

These numbers, starting with 127, are repunit primes 1111111_n in a base n >= 2, so except 7, they are all Brazilian primes belonging to A085104. In fact, 7 = 111_2 is also Brazilian by this other way. (See Links "Les nombres brésiliens", § V.4 -§ V.5.) A088550 is generated by the bases n present in A100330. - Bernard Schott, Dec 20 2012

Examples

			a(3) = 1093 = 3^6 + 3^5 + 3^4 + 3^3 + 3^2 + 3 + 1 is prime.
		

Crossrefs

Programs

  • Magma
    [a: n in [0..100] | IsPrime(a) where a is 1+n+n^2+n^3+n^4+n^5+n^6] ; // Vincenzo Librandi, Jul 14 2012
  • Maple
    A088550 := proc(n)
        numtheory[cyclotomic](7,A100330(n)) ;
    end proc:
    seq(A088550(n),n=1..30) ;
  • Mathematica
    Select[Table[n^6 + n^5 + n^4 + n^3 + n^2 + n + 1, {n, 100}], PrimeQ] (* Alonso del Arte, Feb 07 2014 *)
    Select[Table[Total[n^Range[0,6]],{n,100}],PrimeQ] (* Harvey P. Dale, Aug 13 2024 *)
  • PARI
    polypn(n,p) = { for(x=1,n, if(p%2,y=2,y=1); for(m=1,p, y=y+x^m; ); if(isprime(y),print1(y",")); ) }
    

A128164 Least k > 2 such that (n^k - 1)/(n-1) is prime, or 0 if no such prime exists.

Original entry on oeis.org

3, 3, 0, 3, 3, 5, 3, 0, 19, 17, 3, 5, 3, 3, 0, 3, 25667, 19, 3, 3, 5, 5, 3, 0, 7, 3, 5, 5, 5, 7, 0, 3, 13, 313, 0, 13, 3, 349, 5, 3, 1319, 5, 5, 19, 7, 127, 19, 0, 3, 4229, 103, 11, 3, 17, 7, 3, 41, 3, 7, 7, 3, 5, 0, 19, 3, 19, 5, 3, 29, 3, 7, 5, 5, 3, 41, 3, 3, 5, 3, 0, 23, 5, 17, 5, 11, 7, 61, 3, 3
Offset: 2

Views

Author

Alexander Adamchuk, Feb 20 2007

Keywords

Comments

a(n) = A084740(n) for all n except n = p-1, where p is an odd prime, for which A084740(n) = 2.
All nonzero terms are odd primes.
a(n) = 0 for n = {4,9,16,25,32,36,49,64,81,100,121,125,144,...}, which are the perfect powers with exceptions of the form n^(p^m) where p>2 and (n^(p^(m+1))-1)/(n^(p^m)-1) are prime and m>=1 (in which case a(n^(p^m))=p). - Max Alekseyev, Jan 24 2009
a(n) = 3 for n in A002384, i.e., for n such that n^2 + n + 1 is prime.
a(152) > 20000. - Eric Chen, Jun 01 2015
a(n) is the least number k such that (n^k - 1)/(n-1) is a Brazilian prime, or 0 if no such Brazilian prime exists. - Bernard Schott, Apr 23 2017
These corresponding Brazilian primes are in A285642. - Bernard Schott, Aug 10 2017
a(152) = 270217, see the top PRP link. - Eric Chen, Jun 04 2018
a(184) = 16703, a(200) = 17807, a(210) = 19819, a(306) = 26407, a(311) = 36497, a(326) = 26713, a(331) = 25033; a(185) > 66337, a(269) > 63659, a(281) > 63421, and there are 48 unknown a(n) for n <= 1024. - Eric Chen, Jun 04 2018
Six more terms found: a(522)=20183, a(570)=12907, a(684)=22573, a(731)=15427, a(820)=12043, a(996)=14629. - Michael Stocker, Apr 09 2020

Examples

			a(7) = 5 because (7^5 - 1)/6 = 2801 = 11111_7 is prime and (7^k - 1)/6 = 1, 8, 57, 400 for k = 1, 2, 3, 4. - _Bernard Schott_, Apr 23 2017
		

Crossrefs

Cf. A002384, A049409, A100330, A162862, A217070-A217089. (numbers b such that (b^p-1)/(b-1) is prime for prime p = 3 to 97)
A126589 gives locations of zeros.

Programs

  • Mathematica
    Table[Function[m, If[m > 0, k = 3; While[! PrimeQ[(m^k - 1)/(m - 1)], k++]; k, 0]]@ If[Set[e, GCD @@ #[[All, -1]]] > 1, {#, IntegerExponent[n, #]} &@ Power[n, 1/e] /. {{k_, m_} /; Or[Not[PrimePowerQ@ m], Prime@ m, FactorInteger[m][[1, 1]] == 2] :> 0, {k_, m_} /; m > 1 :> n}, n] &@ FactorInteger@ n, {n, 2, 17}] (* Michael De Vlieger, Apr 24 2017 *)
  • PARI
    a052409(n) = my(k=ispower(n)); if(k, k, n>1)
    a052410(n) = if (ispower(n, , &r), r, n)
    is(n) = issquare(n) || (ispower(n) && !ispseudoprime((n^a052410(a052409(n))-1)/(n-1)))
    a(n) = if(is(n), 0, forprime(p=3, 2^16, if(ispseudoprime((n^p-1)/(n-1)), return(p)))) \\ Eric Chen, Jun 01 2015, corrected by Eric Chen, Jun 04 2018, after Charles R Greathouse IV in A052409 and Michel Marcus in A052410

Extensions

a(18) = 25667 found by Henri Lifchitz, Sep 26 2007

A217071 Numbers k such that (k^17-1)/(k-1) is prime.

Original entry on oeis.org

2, 11, 20, 21, 28, 31, 55, 57, 62, 84, 87, 97, 107, 109, 129, 147, 149, 157, 160, 170, 181, 189, 191, 207, 241, 247, 251, 274, 295, 297, 315, 327, 335, 349, 351, 355, 364, 365, 368, 379, 383, 410, 419, 423, 431, 436, 438, 466, 472, 506, 513, 527, 557, 571, 597
Offset: 1

Views

Author

Tim Johannes Ohrtmann, Sep 26 2012

Keywords

Crossrefs

Programs

  • Magma
    [n: n in [2..1000] |IsPrime((n^17 - 1) div (n - 1))]; // Vincenzo Librandi, Sep 28 2012
    
  • Mathematica
    Select[Range[2, 1000], PrimeQ[(#^17 - 1)/(# - 1)] &] (* T. D. Noe, Sep 26 2012 *)
  • PARI
    is(n)=isprime((n^17-1)/(n-1)) \\ Charles R Greathouse IV, Feb 17 2017

Extensions

More terms from T. D. Noe, Sep 26 2012

A217072 Numbers k such that (k^19-1)/(k-1) is prime.

Original entry on oeis.org

2, 10, 11, 12, 14, 19, 24, 40, 45, 46, 48, 65, 66, 67, 75, 85, 90, 103, 105, 117, 119, 137, 147, 164, 167, 179, 181, 205, 220, 235, 242, 253, 254, 263, 268, 277, 303, 315, 332, 337, 366, 369, 370, 389, 399, 404, 424, 431, 446, 449, 480, 481, 506, 509, 521, 523
Offset: 1

Views

Author

Tim Johannes Ohrtmann, Sep 26 2012

Keywords

Crossrefs

Programs

  • Magma
    [n: n in [2..1000] |IsPrime((n^19 - 1) div (n - 1))]; // Vincenzo Librandi, Sep 28 2012
    
  • Mathematica
    Select[Range[2, 1000], PrimeQ[(#^19 - 1)/(# - 1)] &] (* T. D. Noe, Sep 26 2012 *)
  • PARI
    is(n)=isprime((n^19-1)/(n-1)) \\ Charles R Greathouse IV, Feb 17 2017

Extensions

More terms from T. D. Noe, Sep 26 2012

A246397 Numbers n such that Phi(12, n) is prime, where Phi is the cyclotomic polynomial.

Original entry on oeis.org

2, 3, 4, 5, 9, 10, 12, 13, 17, 25, 27, 30, 31, 36, 38, 39, 43, 48, 52, 55, 56, 61, 62, 65, 83, 92, 94, 99, 100, 104, 105, 109, 114, 118, 126, 131, 166, 168, 169, 172, 183, 185, 190, 194, 196, 198, 209, 224, 225, 229, 231, 239, 244, 257, 260, 261, 263, 269, 270, 272, 278, 291, 296, 299, 300, 302, 308, 311
Offset: 1

Views

Author

Eric Chen, Nov 13 2014

Keywords

Comments

Numbers n such that n^4-n^2+1 is prime, or numbers n such that A060886(n) is prime.

Crossrefs

Cf. A008864 (1), A006093 (2), A002384 (3), A005574 (4), A049409 (5), A055494 (6), A100330 (7), A000068 (8), A153439 (9), A246392 (10), A162862 (11), this sequence (12), A217070 (13), A006314 (16), A217071 (17), A164989 (18), A217072 (19), A217073 (23), A153440 (27), A217074 (29), A217075 (31), A006313 (32), A097475 (36), A217076 (37), A217077 (41), A217078 (43), A217079 (47), A217080 (53), A217081 (59), A217082 (61), A006315 (64), A217083 (67), A217084 (71), A217085 (73), A217086 (79), A153441 (81), A217087 (83), A217088 (89), A217089 (97), A006316 (128), A153442 (243), A056994 (256), A056995 (512), A057465 (1024), A057002 (2048), A088361 (4096), A088362 (8192), A226528 (16384), A226529 (32768), A226530 (65536).

Programs

  • Maple
    A246397:=n->`if`(isprime(n^4-n^2+1),n,NULL): seq(A246397(n),n=1..300); # Wesley Ivan Hurt, Nov 14 2014
  • Mathematica
    Select[Range[350], PrimeQ[Cyclotomic[12, #]] &] (* Vincenzo Librandi, Jan 17 2015 *)
  • PARI
    for(n=1,10^3,if(isprime(polcyclo(12,n)),print1(n,", "))); \\ Joerg Arndt, Nov 13 2014
Showing 1-10 of 36 results. Next