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-2 of 2 results.

A198244 Primes of the form k^10 + k^9 + k^8 + k^7 + k^6 + k^5 + k^4 + k^3 + k^2 + k + 1 where k is nonprime.

Original entry on oeis.org

11, 10778947368421, 17513875027111, 610851724137931, 614910264406779661, 22390512687494871811, 22793803793211153712637, 79905927161140977116221, 184251916941751188170917, 319465039747605973452001, 1311848376806967295019263, 1918542715220370688851293
Offset: 1

Views

Author

Jonathan Vos Post, Dec 21 2012

Keywords

Comments

Subsequence of A060885.
From Bernard Schott, Nov 01 2019: (Start)
These are the primes associated with the terms k of A308238.
A162861 = A286301 Union {this sequence}.
The numbers of this sequence R_11 = 11111111111_k with k > 1 are Brazilian primes, so belong to A085104. (End)

Examples

			10778947368421 is in the sequence since 10778947368421 = 20^10 + 20^9 + 20^8 + 20^7 + 20^6 + 20^5 + 20^4 + 20^3 + 20^2 + 20 + 1, 20 is not prime, and 10778947368421 is prime.
		

Crossrefs

Similar to A185632 (k^2+ ...), A193366 (k^4+ ...), A194194 (k^6+ ...).

Programs

  • Magma
    [a: n in [0..500] | not IsPrime(n) and IsPrime(a) where a is (n^10+n^9+n^8+n^7+n^6+n^5+n^4+n^3+n^2+n+1)]; // Vincenzo Librandi, Nov 09 2014
    
  • Maple
    f:= proc(n)
    local p,j;
    if isprime(n) then return NULL fi;
    p:= add(n^j,j=0..10);
    if isprime(p) then p else NULL fi
    end proc:
    map(f, [$1..1000]); # Robert Israel, Nov 19 2014
  • PARI
    forcomposite(n=0,10^3,my(t=sum(k=0,10,n^k));if(isprime(t),print1(t,", "))); \\ Joerg Arndt, Nov 10 2014
  • Python
    from sympy import isprime
    A198244_list, m = [], [3628800, -15966720, 28828800, -27442800, 14707440, -4379760, 665808, -42240, 682, 0, 1]
    for n in range(1,10**4):
        for i in range(10):
            m[i+1]+= m[i]
        if not isprime(n) and isprime(m[-1]):
            A198244_list.append(m[-1]) # Chai Wah Wu, Nov 09 2014
    

Formula

{A060885(A018252(n)) which are in A000040}.

Extensions

a(5)-a(6) from Robert G. Wilson v, Dec 21 2012
a(7) from Michael B. Porter, Dec 27 2012
Corrected terms a(6)-a(7) and added terms by Chai Wah Wu, Nov 09 2014

A288939 Nonprime numbers k such that k^6 + k^5 + k^4 + k^3 + k^2 + k + 1 is prime.

Original entry on oeis.org

1, 6, 14, 26, 38, 40, 46, 56, 60, 66, 68, 72, 80, 87, 93, 95, 115, 122, 126, 128, 146, 156, 158, 160, 180, 186, 192, 203, 206, 208, 220, 221, 235, 237, 238, 264, 266, 280, 282, 290, 294, 300, 303, 320, 341, 350, 363, 381, 395, 399, 404, 405, 417, 418, 436, 438, 447, 450
Offset: 1

Views

Author

Bernard Schott, Jun 19 2017

Keywords

Comments

A163268 Union {This sequence} = A100330.
The corresponding prime numbers k^6 + k^5 + k^4 + k^3 + k^2 + k + 1 = 1111111_k are in A194194; all these Brazilian primes belong to A085104 and A285017.

Examples

			6 is in the sequence because 6^6 + 6^5 + 6^4 + 6^3 + 6^2 + 6 + 1 = 1111111_6 = 55987 which is prime.
		

Crossrefs

Programs

  • Maple
    for n from 1 to 200 do s(n):= 1+n+n^2+n^3+n^4+n^5+n^6;
    if not isprime(n) and isprime(s(n)) then print(n,s(n)) else fi; od:
  • Mathematica
    Select[Range@ 450, And[! PrimeQ@ #, PrimeQ@ Total[#^Range[0, 6]]] &] (* Michael De Vlieger, Jun 19 2017 *)
  • PARI
    isok(n) = !isprime(n) && isprime(1+n+n^2+n^3+n^4+n^5+n^6); \\ Michel Marcus, Jun 19 2017
    
  • Python
    from sympy import isprime
    A288939_list = [n for n in range(10**3) if not isprime(n) and isprime(n*(n*(n*(n*(n*(n + 1) + 1) + 1) + 1) + 1) + 1)] # Chai Wah Wu, Jul 13 2017
Showing 1-2 of 2 results.