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.

A006094 Products of 2 successive primes.

Original entry on oeis.org

6, 15, 35, 77, 143, 221, 323, 437, 667, 899, 1147, 1517, 1763, 2021, 2491, 3127, 3599, 4087, 4757, 5183, 5767, 6557, 7387, 8633, 9797, 10403, 11021, 11663, 12317, 14351, 16637, 17947, 19043, 20711, 22499, 23707, 25591, 27221, 28891, 30967, 32399, 34571, 36863
Offset: 1

Views

Author

Keywords

Comments

The Huntley reference would suggest prefixing the sequence with an initial 4 - Enoch Haga. [But that would conflict with the definition! - N. J. A. Sloane, Oct 13 2009]
Sequence appears to coincide with the sequence of numbers n such that the largest prime < sqrt(n) and the smallest prime > sqrt(n) divide n. - Benoit Cloitre, Apr 04 2002
This is true: p(n) < [ sqrt(a(n)) = sqrt(p(n)*p(n+1)) ] < p(n+1) by definition. - Jon Perry, Oct 02 2013
a(n+1) = smallest number such that gcd(a(n), a(n+1)) = prime(n+1). - Alexandre Wajnberg and Ray Chandler, Oct 14 2005
Also the area of rectangles whose side lengths are consecutive primes. E.g., the consecutive primes 7,11 produce a 7 X 11 unit rectangle which has area 77 square units. - Cino Hilliard, Jul 28 2006
a(n) = A001358(A172348(n)); A046301(n) = lcm(a(n), a(n+1)); A065091(n) = gcd(a(n), a(n+1)); A066116(n+2) = a(n+1)*a(n); A109805(n) = a(n+1) - a(n). - Reinhard Zumkeller, Mar 13 2011
See A209329 for the sum of the reciprocals. - M. F. Hasler, Jan 22 2013
A078898(a(n)) = 3. - Reinhard Zumkeller, Apr 06 2015

References

  • H. E. Huntley, The Divine Proportion, A Study in Mathematical Beauty. New York: Dover, 1970. See Chapter 13, Spira Mirabilis, especially Fig. 13-5, page 173.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Subset of the squarefree semiprimes, A006881.
Subsequence of A256617 and A097889.

Programs

  • Haskell
    a006094 n = a006094_list !! (n-1)
    a006094_list = zipWith (*) a000040_list a065091_list
    -- Reinhard Zumkeller, Mar 13 2011
    
  • Haskell
    a006094_list = pr a000040_list
      where pr (n:m:tail) = n*m : pr (m:tail)
            pr _ = []
    -- Jean-François Antoniotti, Jan 08 2020
    
  • Magma
    [NthPrime(n)*NthPrime(n+1): n in [1..41]]; // Bruno Berselli, Feb 24 2011
    
  • Maple
    a:= n-> (p-> p(n)*p(n+1))(ithprime):
    seq(a(n), n=1..43);  # Alois P. Heinz, Jan 02 2021
  • Mathematica
    Table[ Prime[n] Prime[n + 1], {n, 40}] (* Robert G. Wilson v, Jan 22 2004 *)
    Times@@@Partition[Prime[Range[60]], 2, 1] (* Harvey P. Dale, Oct 15 2011 *)
  • MuPAD
    ithprime(i)*ithprime(i+1) $ i = 1..41 // Zerinvary Lajos, Feb 26 2007
    
  • PARI
    g(n) = for(x=1,n,print1(prime(x)*prime(x+1)",")) \\ Cino Hilliard, Jul 28 2006
    
  • PARI
    is(n)=my(p=precprime(sqrtint(n))); p>1 && n%p==0 && isprime(n/p) && nextprime(p+1)==n/p \\ Charles R Greathouse IV, Jun 04 2014
    
  • Python
    from sympy import prime, primerange
    def aupton(nn):
        alst, prevp = [], 2
        for p in primerange(3, prime(nn+1)+1): alst.append(prevp*p); prevp = p
        return alst
    print(aupton(43)) # Michael S. Branicky, Jun 15 2021
    
  • Python
    from sympy import prime, nextprime
    def A006094(n): return (p:=prime(n))*nextprime(p) # Chai Wah Wu, Oct 18 2024

Formula

A209329 = Sum_{n>=2} 1/a(n). - M. F. Hasler, Jan 22 2013
a(n) = A000040(n) * A000040(n+1). - Alois P. Heinz, Jan 02 2021

A066117 Triangle read by rows: T(n,k) = T(n-1,k-1)*T(n,k-1) and T(n,1) = prime(n).

Original entry on oeis.org

2, 3, 6, 5, 15, 90, 7, 35, 525, 47250, 11, 77, 2695, 1414875, 66852843750, 13, 143, 11011, 29674645, 41985913344375, 2806877704512541816406250, 17, 221, 31603, 347980633, 10326201751150285, 433555011900329243987584396875
Offset: 1

Views

Author

Henry Bottomley, Dec 05 2001

Keywords

Comments

As a square array read by descending antidiagonals, A(n, k), n >= 1, k >= 1, gives the encoding defined in A297845 of the polynomial (x+1)^(n-1) * x^(k-1). - Peter Munn, Jul 27 2022

Examples

			T(4,3) = T(3,2)*T(4,2) = 15*35 = 525. Rows start
     2;
    3, 6;
  5, 15, 90;
7, 35, 525, 47250;
...
From _Antti Karttunen_, Sep 18 2016: (Start)
Alternatively, this table can be viewed as a square array. Then the top left 5x4 corner looks as:
    2,       3,        5,         7,         11
    6,      15,       35,        77,        143
   90,     525,     2695,     11011,      31603
47250, 1414875, 29674645, 347980633, 2255916949
(End)
		

Crossrefs

Cf. A000040, A006094 and A066116 (three leftmost diagonal of triangular table = three topmost rows of square array).
Cf. A007188, A267096 (two rightmost diagonals of the triangular table = two leftmost columns of square array).
Cf. also A099884, A255483, A276586, A276588 (other arrays derived from this one).

Programs

Formula

From Antti Karttunen, Sep 19 2016: (Start)
When computed as a square array A(row,col), row >= 1, col >= 1:
A(1,col) = A000040(col), for row > 1, A(row,col) = A(row-1,col)*A(row-1,col+1).
A(row,1) = A007188(row-1), for col > 1, A(row,col) = A003961(A(row,col-1)).
For all row >= 1, col >= 1, A055396(A(row,col)) = col.
(End)
A(1,1) = 2; for n > 1, A(n,k) = A297845(A(n-1,k),6); for k > 1, A(n,k) = A297845(A(n,k-1),3). - Peter Munn, Jul 20 2022
Showing 1-2 of 2 results.