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.

A215658 Primes p such that the smallest positive integer k for which p# + k is square satisfies p# + k = k^2, where p# = 2*3*5*7*11*...*p is a primorial.

Original entry on oeis.org

2, 3, 5, 7, 17
Offset: 1

Views

Author

Jonathan Sondow, Sep 02 2012

Keywords

Comments

The corresponding values of k are 2, 3, 6, 15, 715 = A215659.
The equation p# + k = k^2 has an integer solution k if and only if 1 + 4*p# is a square.
Conjecture: Not the same sequence as A192579, which is finite.
When p is in this sequence, p# = k(k-1) is in A161620, the intersection of A002110 and A002378. - Jeppe Stig Nielsen, Mar 27 2018

Examples

			The smallest square > 17# = 510510 is 715^2 = 17# + 715, so 17 is a member.
		

Crossrefs

Programs

  • Mathematica
    t = {}; pm = 1; Do[pm = pm*p; s = Floor[Sqrt[pm]]; If[pm == s*(s+1), AppendTo[t, p]], {p, Prime[Range[100]]}]; t (* T. D. Noe, Sep 05 2012 *)
  • PARI
    for (n=1, 10, if (ceil(sqrt(prod(i=1, n, prime(i))))^2 - prod(i=1, n, prime(i)) - ceil(sqrt(prod(i=1, n, prime(i)))) == 0, print(prime(n)));); \\ Michel Marcus, Sep 05 2012
    
  • Python
    from sympy import primorial, integer_nthroot, prime
    A215658_list = [prime(i) for i in range(1,10**2) if integer_nthroot(4*primorial(i)+1,2)[1]] # Chai Wah Wu, Apr 01 2021

Formula

A145781(n) = A216144(n) if and only if prime(n) is a member.
a(n)# = A215659(n)*(A215659(n)-1).

A145781 Least m >= 0 which when added to primorial(n) yields a square.

Original entry on oeis.org

2, 3, 6, 15, 91, 246, 715, 3535, 21099, 95995, 175470, 4468006, 31516774, 192339970, 212951314, 5138843466, 76699112491, 103728576730, 3051100701874, 14417674958635, 245230361204214, 2296196521511806, 10940476546738414
Offset: 1

Views

Author

Keywords

Examples

			a(1)= 4 - 2 = 2, a(2)= 9 - 2*3 = 3, a(3)= 36 - 2*3*5 = 6, a(4)= 225 - 2*3*5*7 = 15, and a(5)= 2401 - 2*3*5*7*11 = 91.
		

Crossrefs

Programs

  • Mathematica
    (IntegerPart[Sqrt[#]]+1)^2-#&/@Rest[FoldList[Times,1,Prime[Range[30]]]]  (* Harvey P. Dale, Jan 14 2011 *)
  • PARI
    j=[];for (n=1,30,j=concat(j,(ceil(sqrt(prod(i=1,n,prime(i))))^2 -prod(i=1,n,prime(i)))));j

Formula

a(n) = A216144(n)^2 - A002110(n). - Jonathan Sondow, Sep 02 2012

Extensions

Definition shortened and clarified by Jonathan Sondow, Sep 02 2012
Showing 1-2 of 2 results.