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

A060723 a(n) is the denominator of r(n) where r(n) is the sequence of rational numbers defined by the recursion: r(0) = 0, r(1) = 1 and for n>1 r(n) = r(n-1) + r(n-2)/2. From this definition it is clear that a(n) is always a power of 2 (see A060755).

Original entry on oeis.org

1, 1, 1, 2, 1, 4, 4, 8, 1, 16, 16, 32, 8, 64, 64, 128, 8, 256, 256, 512, 128, 1024, 1024, 2048, 256, 4096, 4096, 8192, 2048, 16384, 16384, 32768, 1024, 65536, 65536, 131072, 32768, 262144, 262144, 524288, 65536, 1048576, 1048576, 2097152
Offset: 0

Views

Author

Avi Peretz (njk(AT)netvision.net.il), Apr 21 2001

Keywords

Comments

It can be proved that r(n) is an integer (i.e. a(n) = 1) if and only if n is one of 0, 1, 2, 4, 8.

Examples

			The sequence r(n) begins 0, 1, 1, 3/2, 2, 11/4, 15/4, 41/8, 7, 153/16, 209/16, 571/32, 363/16, 2023/64, 2749/64, 7521/128, 5135/64, ...
		

Crossrefs

Cf. A060755, A305491 (numerators).

Programs

  • Mathematica
    Denominator[RecurrenceTable[{a[0]==0,a[1]==1,a[n]==a[n-1]+a[n-2]/2},a,{n,50}]] (* Harvey P. Dale, Mar 07 2016 *)
    Table[Denominator[Simplify[((1/2(1 + Sqrt[3]))^x - (1/2(Sqrt[3] - 1))^x Cos[Pi x])/ Sqrt[3]]], {x, 0, 43}]  (* Peter Luschny, Jun 02 2018 *)

Formula

r(n) = (((1/2)*(sqrt(3) + 1))^n - ((1/2)*(sqrt(3) - 1))^n * cos(Pi*n))/sqrt(3). - Peter Luschny, Jun 02 2018

Extensions

More terms from Vladeta Jovovic, Apr 24 2001

A060796 Upper central divisor of n-th primorial.

Original entry on oeis.org

2, 3, 6, 15, 55, 182, 715, 3135, 15015, 81345, 448630, 2733549, 17490603, 114388729, 785147363, 5708795638, 43850489690, 342503171205, 2803419704514, 23622001517543, 201817933409378, 1793779635410490, 16342166369958702, 154171363634898185, 1518410187442699518, 15259831781575946565
Offset: 1

Views

Author

Labos Elemer, Apr 27 2001

Keywords

Comments

Also: Write product of first n primes as x*y with x < y and x maximal; sequence gives value of y. Indeed, p(n)# = primorial(n) = A002110(n) is never a square for n >= 1; all exponents in the prime factorization are 1. Therefore primorial(n) has N = 2^n distinct divisors. Since this is an even number, the N divisors can be grouped in N/2 pairs {d(k), d(N+1-k)} with product equal to p(n)#. One of the two is always smaller and one is larger than sqrt(p(n)#). This sequence gives the (2^(n-1)+1)-th divisor, which is the smallest one larger than sqrt(p(n)#). - M. F. Hasler, Sep 20 2011

Examples

			n = 8, q(8) = 2*3*5*7*11*13*17*19 = 9699690. Its 128th and 129th divisors are {3094, 3135}: a(8) = 3135, and 3094 < A000196(9699690) = 3114 < 3135. [Corrected by _M. F. Hasler_, Sep 20 2011]
		

Crossrefs

Programs

  • Mathematica
    k = 1; Do[k *= Prime[n]; l = Divisors[k]; x = Length[l]; Print[l[[x/2 + 1]]], {n, 1, 24}] (* Ryan Propper, Jul 25 2005 *)
  • PARI
    A060796(n) = divisors(prod(k=1,n,prime(k)))[2^(n-1)+1] \\ Requires stack size > 2^(n+5). - M. F. Hasler, Sep 20 2011

Formula

a(n) = A033677(A002110(n)).
a(n) = A002110(n)/A060795(n). - M. F. Hasler, Mar 21 2022

Extensions

More terms from Ryan Propper, Jul 25 2005
a(24)-a(37) in b-file calculated from A182987 by M. F. Hasler, Sep 20 2011
a(38) from David A. Corneth, Mar 21 2022
a(39)-a(70) in b-file from Max Alekseyev, Apr 20 2022

A060795 Write product of first n primes as x*y with x

Original entry on oeis.org

1, 2, 5, 14, 42, 165, 714, 3094, 14858, 79534, 447051, 2714690, 17395070, 114371070, 783152070, 5708587335, 43848093003, 342444658094, 2803119896185, 23619540863730, 201813981102615, 1793779293633437, 16342050964565645, 154170926013430326, 1518409177581024365
Offset: 1

Views

Author

Labos Elemer, Apr 27 2001

Keywords

Comments

Or, lower central divisor of n-th primorial.
Subsequence of A005117 (squarefree numbers). - Michel Marcus, Feb 22 2016

Examples

			n = 8: q(8) = 2*3*5*7*11*13*17*19 = 9699690. Its 128th and 129th divisors are {3094, 3135}: a(8) = 3094 and 3094 < A000196(9699690) = 3114 < 3135. [Corrected by _Colin Barker_, Oct 22 2010]
2*3*5*7 = 210 = 14*15 with difference of 1, so a(4) = 14.
		

Crossrefs

Programs

  • Maple
    F:= proc(n) local P,N,M;
         P:= {seq(ithprime(i),i=1..n)};
         N:= floor(sqrt(convert(P,`*`)));
         M:= map(convert, combinat:-powerset(P),`*`);
         max(select(`<=`,M,N))
    end proc:
    map(F, [$1..20]); # Robert Israel, Feb 22 2016
  • Mathematica
    a[n_] := (m = Times @@ Prime[Range[n]] ; dd = Divisors[m]; dd[[Length[dd]/2 // Floor]]); Table[Print[an = a[n]]; an, {n, 1, 25}] (* Jean-François Alcover, Oct 15 2016 *)
  • PARI
    a(n) = my(m=prod(i=1, n, prime(i))); divisors(m)[numdiv(m)\2]; \\ Michel Marcus, Feb 22 2016

Formula

a(n) = A060775(A002110(n)). - Labos Elemer, Apr 27 2001
a(n) = A002110(n)/A060796(n). - M. F. Hasler, Mar 21 2022

Extensions

More terms from Ed Pegg Jr, May 28 2001
a(16)-a(23) computed by Jud McCranie, Apr 15 2000
a(24) and a(25) from Robert Israel, Feb 22 2016
a(25) corrected by Jean-François Alcover, Oct 15 2016
a(26)-a(33) in b-file from Amiram Eldar, Apr 09 2020
Up to a(38) using b-file of A060796, by M. F. Hasler, Mar 21 2022
a(39)-a(70) in b-file from Max Alekseyev, Apr 20 2022

A060797 Integer part of square root of n-th primorial, A002110(n).

Original entry on oeis.org

1, 1, 2, 5, 14, 48, 173, 714, 3114, 14936, 80434, 447839, 2724103, 17442771, 114379899, 784149081, 5708691485, 43849291330, 342473913399, 2803269796341, 23620771158594, 201815957246321, 1793779464521955, 16342108667160301, 154171144824008979
Offset: 0

Views

Author

Labos Elemer, Apr 27 2001

Keywords

Comments

Integer part of square root of product of n first primes.

Examples

			n=8, q(8) = 2*3*5*7*11*13*17*19 = 9699690, a(8)=3114. This is between the 128th and 129th divisors of the 8th primorial: 3094 < A000196(9699690)=3114 < 3135.
(In general, x=A002110(n) always has 2^n divisors, and A000196(x) always lies between the k-th and (k+1)-th divisors of x, where k=ceiling(tau(x)/2) and tau(x) is the number of divisors of x.) - _M. F. Hasler_, Sep 02 2012
		

Crossrefs

Programs

  • Mathematica
    a = {}; Do[b = 1; Do[b = b Prime[x], {x, 1, n}]; AppendTo[a, Floor[b^(1/2)]], {n, 1, 100}]; a (* Artur Jasinski *)
    Join[{1},Floor[Sqrt[#]]&/@FoldList[Times,Prime[Range[30]]]] (* Harvey P. Dale, Nov 22 2023 *)
  • PARI
    A060797(n)=sqrtint(prod(k=1, n, prime(k))) \\ M. F. Hasler, Sep 02 2012

Formula

a(n) = A000196(A002110(n)) = floor(sqrt(A002110(n))).

Extensions

a(23) correction by Hans Havermann, Dec 02 2010
Extended to a(0)=1=sqrt(A002110(0)) by M. F. Hasler, Sep 02 2012

A305491 a(n) = numerator(r(n)) where r(n) = (((1/2)*(sqrt(3) + 1))^n - ((1/2)*(sqrt(3) - 1))^n * cos(Pi*n))/sqrt(3).

Original entry on oeis.org

0, 1, 1, 3, 2, 11, 15, 41, 7, 153, 209, 571, 195, 2131, 2911, 7953, 679, 29681, 40545, 110771, 37829, 413403, 564719, 1542841, 263445, 5757961, 7865521, 21489003, 7338631, 80198051, 109552575, 299303201, 12776743, 1117014753, 1525870529, 4168755811, 1423656585
Offset: 0

Views

Author

Peter Luschny, Jun 02 2018

Keywords

Comments

Let f(x, y) = ((y+1)^x - (y-1)^x * cos(Pi*x))/(y * 2^x). Then f(n, sqrt(3)) are the rational numbers a(n)/A060723(n) and f(n, sqrt(5)) the Fibonacci numbers A000045(n).
From Paul Curtz, Dec 05 2018: (Start)
The binomial inverse of the rational sequence r(n) starts 0, 1, -1, 3/2, -2, 11/4, -15/4, 41/8, -7, 153/16, -209/16, ... and is up to signs equal to r(n). The difference table starts:
0, 1, 1, 3/2, 2, 11/4, 15/4, 41/8, ...
1, 0, 1/2, 1/2, 3/4, 1, 11/8, 15/8, ...
-1, 1/2, 0, 1/4, 1/4, 3/8, 1/2, 11/16, ...
3/2, -1/2, 1/4, 0, 1/8, 1/8, 3/16, 1/4, ...
...
Let s(n) = 2*r(n+1) - r(n) then s(n) = 1, 2, 5/2, 7/2, 19/4, 13/2, ... = A173299(n)/A173300(n) for n >= 1. (End)

Crossrefs

Cf. A060723 (denominators), A060755, A000045, A305492.

Programs

  • Mathematica
    Table[Numerator[Simplify[((1/2 (Sqrt[3] + 1))^x - (1/2 (Sqrt[3] - 1))^x Cos[Pi  x])/Sqrt[3]]], {x, 0, 36}]

Formula

A recurrence for r(n) is given in A060723.
Showing 1-5 of 5 results.