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

A005382 Primes p such that 2p-1 is also prime.

Original entry on oeis.org

2, 3, 7, 19, 31, 37, 79, 97, 139, 157, 199, 211, 229, 271, 307, 331, 337, 367, 379, 439, 499, 547, 577, 601, 607, 619, 661, 691, 727, 811, 829, 877, 937, 967, 997, 1009, 1069, 1171, 1237, 1279, 1297, 1399, 1429, 1459, 1531, 1609, 1627, 1657, 1759, 1867, 2011
Offset: 1

Views

Author

Keywords

Comments

Sequence gives values of p such Sum_{i=1..p} gcd(p,i) = A018804(p) is prime. - Benoit Cloitre, Jan 25 2002
Let q = 2n-1. For these n (and q), the sum of two cyclotomic polynomials can be written as a product of cyclotomic polynomials and as a cyclotomic polynomial in x^2: Phi(q,x) + Phi(2q,x) = 2 Phi(n,x) Phi(2n,x) = 2 Phi(n,x^2). - T. D. Noe, Nov 04 2003
Primes in A006254. - Zak Seidov, Mar 26 2013
If a(n) is in A168421 then A005383(n) is a twin prime with a Ramanujan prime, A005383(n) - 2. If this sequence has an infinite number of terms in A168421, then the twin prime conjecture can be proved. - John W. Nicholson, Dec 05 2013
Records subsequence of A023509 (n >= 2). - David James Sycamore, May 05 2025

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A010051, A000040, A053685 (subsequence), A006254.
Cf. A023509.

Programs

  • Haskell
    a005382 n = a005382_list !! (n-1)
    a005382_list = filter
       ((== 1) . a010051 . (subtract 1) . (* 2)) a000040_list
    -- Reinhard Zumkeller, Oct 03 2012
    
  • Magma
    [n: n in [0..1000] | IsPrime(n) and IsPrime(2*n-1)]; // Vincenzo Librandi, Nov 18 2010
    
  • Maple
    f := proc(Q) local t1,i,j; t1 := []; for i from 1 to 500 do j := ithprime(i); if isprime(2*j-Q) then t1 := [op(t1),j]; fi; od: t1; end; f(1);
    # second Maple program:
    q:= p-> andmap(isprime, [p, 2*p-1]):
    select(q, [$2..2500])[];  # Alois P. Heinz, Dec 16 2024
  • Mathematica
    Select[Prime[Range[300]], PrimeQ[2#-1]&]
  • PARI
    select(p->isprime(2*p-1),primes(500)) \\ Charles R Greathouse IV, Apr 26 2012
    
  • PARI
    forprime(n=2, 10^3, if(ispseudoprime(2*n-1), print1(n, ", "))) \\ Felix Fröhlich, Jun 15 2014

Formula

a(n) = A129521(n) / A005383(n). - Reinhard Zumkeller, Apr 19 2007
a(n) = (A005383(n) + 1)/2. - Zak Seidov, Nov 04 2010

A047211 Numbers that are congruent to {2, 4} mod 5.

Original entry on oeis.org

2, 4, 7, 9, 12, 14, 17, 19, 22, 24, 27, 29, 32, 34, 37, 39, 42, 44, 47, 49, 52, 54, 57, 59, 62, 64, 67, 69, 72, 74, 77, 79, 82, 84, 87, 89, 92, 94, 97, 99, 102, 104, 107, 109, 112, 114, 117, 119, 122, 124, 127, 129, 132, 134, 137, 139, 142, 144, 147, 149, 152, 154, 157, 159, 162, 164, 167, 169, 172, 174, 177, 179, 182, 184
Offset: 1

Views

Author

N. J. A. Sloane, Dec 11 1999

Keywords

Comments

Conjecture: n such that the characteristic polynomial of M(n) is irreducible over the rationals where M(n) is an n X n matrix with ones on the skew diagonal and below it and the skew line two positions above it and otherwise zeros; see example for one such matrix. Tested up to n=177. - Joerg Arndt, Aug 10 2011

Examples

			The 7 X 7 matrix (dots for zeros):
[....1.1]
[...1.11]
[..1.111]
[.1.1111]
[1.11111]
[.111111]
[1111111]
has the characteristic polynomial x^7 - 5*x^6 - 4*x^5 + 15*x^4 + 5*x^3 - 11*x^2 - x + 1 which is irreducible over the field of rational numbers, and 7 is a term of the sequence. - _Joerg Arndt_, Aug 10 2011
		

Crossrefs

Cf. A053685 (subsequence).

Programs

Formula

a(n) = a(n-1) +a(n-2) -a(n-3).
a(n) = (10*n-3-(-1)^n)/4, (n>=1). [Corrected by Bruno Berselli, Sep 20 2010]
a(n) = 5*floor((n-1)/2) +3 +(-1)^n. - Gary Detlefs, Mar 02 2010
G.f.: x*(2+2*x+x^2)/((1+x)*(1-x)^2). - Paul Barry, Sep 11 2008
a(n) = 5*n-a(n-1)-4 (with a(1)=2). - Vincenzo Librandi, Nov 18 2010
a(n) = floor((5*n-1)/2). - Gary Detlefs, May 14 2011
a(n) = 2*n + floor((n-1)/2). - Arkadiusz Wesolowski, Sep 19 2012
Sum_{n>=1} (-1)^(n+1)/a(n) = sqrt(2+2/sqrt(5))*Pi/10 - sqrt(5)*log(phi)/5, where phi is the golden ratio (A001622). - Amiram Eldar, Dec 07 2021
E.g.f.: 1 + ((10*x - 3)*exp(x) - exp(-x))/4. - David Lovler, Aug 23 2022

Extensions

Conjecture corrected by John M. Campbell, Aug 25 2011

A053635 a(n) = Sum_{d|n} phi(d)*2^(n/d).

Original entry on oeis.org

0, 2, 6, 12, 24, 40, 84, 140, 288, 540, 1080, 2068, 4224, 8216, 16548, 32880, 65856, 131104, 262836, 524324, 1049760, 2097480, 4196412, 8388652, 16782048, 33554600, 67117128, 134218836, 268452240, 536870968, 1073777040, 2147483708, 4295033472, 8589938808
Offset: 0

Views

Author

N. J. A. Sloane, Mar 23 2000

Keywords

Comments

Dirichlet convolution of phi(n) and 2^n. - Richard L. Ollerton, May 06 2021

Crossrefs

Column k=2 of A185651.

Programs

  • Magma
    [0] cat  [&+[EulerPhi(d)*2^(n div d): d in Divisors(n)]: n in [1..40]]; // Vincenzo Librandi, Jul 20 2019
  • Maple
    with(numtheory); A053685:=n->add( phi(n/d)*2^d, d in divisors(n)); # N. J. A. Sloane, Nov 21 2009
  • Mathematica
    a[0] = 0; a[n_] := Sum[EulerPhi[d] 2^(n/d), {d, Divisors[n]}];
    Table[a[n], {n, 0, 31}] (* Jean-François Alcover, Aug 30 2018 *)
  • PARI
    a(n) = if (n, sumdiv(n, d, eulerphi(d)*2^(n/d)), 0); \\ Michel Marcus, Sep 20 2017
    

Formula

a(n) = n * A000031(n).
a(n) = Sum_{k=1..n} 2^gcd(n,k). - Ilya Gutkovskiy, Apr 16 2021
a(n) = Sum_{k=1..n} 2^(n/gcd(n,k))*phi(gcd(n,k))/phi(n/gcd(n,k)). - Richard L. Ollerton, May 06 2021
Showing 1-3 of 3 results.