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.

A300406 Primes of the form 13*2^n + 1.

Original entry on oeis.org

53, 3329, 13313, 13631489, 3489660929, 62864142619960717084721153, 5100145160001678120616578906356228963083163798627028041729, 6779255729241169695101387251026410519979286814120235842117075415451380965612384558178346467329, 1735489466685739441945955136262761093114697424414780375581971306355553527196770446893656695635969
Offset: 1

Views

Author

Martin Renner, Mar 05 2018

Keywords

Comments

For the corresponding exponents n see A032356.

Crossrefs

Programs

  • GAP
    Filtered(List([1..500],n->13*2^n + 1),IsPrime); # Muniru A Asiru, Mar 06 2018
    
  • Magma
    [a: n in [1..400] | IsPrime(a) where a is 13*2^n + 1]; // Vincenzo Librandi, Mar 06 2018
    
  • Maple
    a:=(n,k)->`if`(isprime(k*2^n+1), k*2^n+1, NULL):
    seq(a(n,13), n=1..316);
  • Mathematica
    Select[Table[13 2^n + 1, {n, 400}], PrimeQ] (* Vincenzo Librandi, Mar 06 2018 *)
  • PARI
    lista(nn) = {for(k=1, nn, if(ispseudoprime(p=13*2^k+1), print1(p, ", ")));} \\ Altug Alkan, Mar 29 2018

Formula

a(n) = A168596(A032356(n)). - Michel Marcus, Mar 29 2018

A291049 Primes of the form 2^r * 17^s + 1.

Original entry on oeis.org

2, 3, 5, 17, 137, 257, 65537, 157217, 295937, 557057, 1336337, 96550277, 1212153857, 2281701377, 5473632257, 395469930497, 1401249857537, 2637646790657, 4964982194177, 28572702478337, 1271035441709057, 38280596832649217, 1872540629620228097, 6634884445436379137
Offset: 1

Views

Author

Muniru A Asiru, Sep 15 2017

Keywords

Comments

Primes of the forms a^r * b^s + 1 where (a, b) = (2,1), (2,3), (2,5), (2,7), (2,11) and (2,13) are A092506, A005109, A077497, A077498, A077499 and A173236.
Fermat prime exponents r are 0, 1, 2, 4, 8, 16.
For n > 2, all terms are congruent to 5 (mod 6).
Also, these are prime numbers p for which (p*34^p)/(p-1) is an integer.

Examples

			With n = 1, a(1) = 2^0 * 17^0 + 1 = 2.
With n = 5, a(5) = 2^3 * 17^1 + 1 = 137.
list of (r,s): (0,0), (1,0), (2,0), (4,0), (3,1), (8,0), (16,0), (5,3), (10,2), (15,1), (4,4), (2,6).
		

Crossrefs

Cf. Sequences of primes of form 2^n * q^u + 1: A092506 (q=1), A005109 (q=3), A077497 (q=5), A077498 (q=7), A077499 (q=11), A173236 (q=13).

Programs

  • GAP
    K:=26*10^7+1;; # to get all terms <= K.
    A:=Filtered(Filtered([1,3..K],i-> i mod 6=5),IsPrime);;    I:=[17];;
    B:=List(A,i->Elements(Factors(i-1)));;
    C:=List([0..Length(I)],j->List(Combinations(I,j),i->Concatenation([2],i)));;
    A291049:=Concatenation([2,3],List(Set(Flat(List([1..Length(C)],i->List([1..Length(C[i])],j->Positions(B,C[i][j]))))),i->A[i]));
    
  • Maple
    N:= 10^20: # to get all terms <= N+1
    S:= NULL:
    for r from 0 to ilog2(N) do
      for s from 0 to floor(log[17](N/2^r)) do
        p:= 2^r*17^s +1;
        if isprime(p) then
         S:= S, p
        fi
    od od:
    sort([S]); # Robert Israel, Sep 26 2017
  • Mathematica
    With[{nn = 10^19, q = 17}, Select[Sort@ Flatten@ Table[2^i*q^j + 1, {i, 0, Log[2, nn]}, {j, 0, Log[q, nn/2^i]}], PrimeQ]] (* Michael De Vlieger, Sep 18 2017, after Robert G. Wilson v at A005109 *)
  • PARI
    lista(nn) = my(t, v=List([])); for(r=0, logint(nn, 2), t=2^r; for(s=0, logint(nn\t, 17), if(isprime(t+1), listput(v, t+1)); t*=17)); Vec(vecsort(v)) \\ Jinyuan Wang, Jun 26 2022
Showing 1-2 of 2 results.