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.

Previous Showing 21-30 of 49 results. Next

A132263 Product{0<=k<=floor(log_11(n)), floor(n/11^k)}, n>=1.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 44, 46, 48, 50, 52, 54, 56, 58, 60, 62, 64, 99, 102, 105, 108, 111, 114, 117, 120, 123, 126, 129, 176, 180, 184, 188, 192, 196, 200, 204, 208, 212, 216, 275, 280, 285, 290, 295, 300, 305, 310
Offset: 1

Views

Author

Hieronymus Fischer, Aug 20 2007

Keywords

Comments

If n is written in base-11 as n=d(m)d(m-1)d(m-2)...d(2)d(1)d(0) (where d(k) is the digit at position k) then a(n) is also the product d(m)d(m-1)d(m-2)...d(2)d(1)d(0)*d(m)d(m-1)d(m-2)...d(2)d(1)*d(m)d(m-1)d(m-2)...d(2)*...*d(m)d(m-1)d(m-2)*d(m)d(m-1)*d(m).

Examples

			a(50)=floor(50/11^0)*floor(50/11^1)=50*4=200; a(63)=315 since 63=58(base-11) and so a(63)=58*5(base-11)=63*5=315.
		

Crossrefs

For formulas regarding a general parameter p (i.e. terms floor(n/p^k)) see A132264.
For the product of terms floor(n/p^k) for p=2 to p=12 see A098844(p=2), A132027(p=3)-A132033(p=9), A067080(p=10), A132264(p=12).
For the products of terms 1+floor(n/p^k) see A132269-A132272, A132327, A132328.

Formula

Recurrence: a(n)=n*a(floor(n/11)); a(n*11^m)=n^m*11^(m(m+1)/2)*a(n).
a(k*11^m)=k^(m+1)*11^(m(m+1)/2), for 0
Asymptotic behavior: a(n)=O(n^((1+log_11(n))/2)); this follows from the inequalities below.
a(n)<=b(n), where b(n)=n^(1+floor(log_11(n)))/p^((1+floor(log_11(n)))*floor(log_11(n))/2); equality holds for n=k*11^m, 0=0. b(n) can also be written n^(1+floor(log_11(n)))/11^A000217(floor(log_11(n))).
Also: a(n)<=3^((1-log_11(3))/2)*n^((1+log_11(n))/2)=1.346673852...^((1-log_11(3))/2)*11^A000217(log_11(n)), equality holds for n=3*11^m, m>=0.
a(n)>c*b(n), where c=0.4751041275076031053975644472... (see constant A132265).
Also: a(n)>c*(sqrt(2)/2^log_11(sqrt(2)))*n^((1+log_11(n))/2)=0.607848303...*11^00217(log_11(n)).
lim inf a(n)/b(n)=0.4751041275076031053975644472..., for n-->oo.
lim sup a(n)/b(n)=1, for n-->oo.
lim inf a(n)/n^((1+log_p(n))/2)=0.4751041275076031...*sqrt(2)/2^log_11(sqrt(2)), for n-->oo.
lim sup a(n)/n^((1+log_p(n))/2)=sqrt(3)/3^log_11(sqrt(3))=1.346673852..., for n-->oo.
lim inf a(n)/a(n+1)=0.4751041275076031053975644472... for n-->oo (see constant A132265).

A132264 Product{0<=k<=floor(log_12(n)), floor(n/12^k)}, n>=1.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 48, 50, 52, 54, 56, 58, 60, 62, 64, 66, 68, 70, 108, 111, 114, 117, 120, 123, 126, 129, 132, 135, 138, 141, 192, 196, 200, 204, 208, 212, 216, 220, 224, 228, 232, 236, 300, 305, 310, 315
Offset: 1

Author

Hieronymus Fischer, Aug 20 2007

Keywords

Comments

If n is written in base-12 as n=d(m)d(m-1)d(m-2)...d(2)d(1)d(0) (where d(k) is the digit at position k) then a(n) is also the product d(m)d(m-1)d(m-2)...d(2)d(1)d(0)*d(m)d(m-1)d(m-2)...d(2)d(1)*d(m)d(m-1)d(m-2)...d(2)*...*d(m)d(m-1)d(m-2)*d(m)d(m-1)*d(m).

Examples

			a(50)=floor(50/12^0)*floor(50/12^1)=50*4=200.
a(65)=325 since 65=55(base-12) and so a(65)=55*5(base-12)=65*5=325.
		

Crossrefs

For the product of terms floor(n/p^k) for p=2 to p=11 see A098844(p=2), A132027(p=3)-A132033(p=9), A067080(p=10), A132263(p=11).
For the products of terms 1+floor(n/p^k) see A132269-A132272, A132327, A132328.

Formula

The following formulas are given for a general parameter p considering the product of terms floor(n/p^k) for 0<=k<=floor(log_p(n)), where p=12 for this sequence.
Recurrence: a(n)=n*a(floor(n/p)); a(n*p^m)=n^m*p^(m(m+1)/2)*a(n).
a(k*p^m)=k^(m+1)*p^(m(m+1)/2), for 0
Asymptotic behavior: a(n)=O(n^((1+log_p(n))/2)); this follows from the inequalities below.
a(n)<=b(n), where b(n)=n^(1+floor(log_p(n)))/p^((1+floor(log_p(n)))*floor(log_p(n))/2); equality holds for n=k*p^m, 0=0. b(n) can also be written n^(1+floor(log_p(n)))/p^A000217(floor(log_p(n))).
Also: a(n)<=q^((1-log_p(q))/2)*n^((1+log_p(n))/2)=q^((1-log_p(q))/2)*p^A000217(log_p(n)), equality holds for n=q*p^m, m>=0, where q=floor(sqrt(p)+1/2). Also, equality holds for n=(q+1)*p^m, provided p is a A002378-number (in this case we have p=q*(q+1) and so q^((1-log_p(q))/2)=(q+1)^((1-log_p(q+1))/2)).
a(n)>c*b(n), where c=product{k>0, 1-1/(2*p^k)}=0.47735217025489380... (for p=12 see constant A132265).
Also: a(n)>c*(sqrt(2)/2^log_p(sqrt(2)))*n^((1+log_p(n))/2)=0.612870619...*p^A000217(log_p(n)), (p=12).
lim inf a(n)/b(n)=product{k>0, 1-1/(2*p^k)}=0.47735217025489380198334286365820..., for n-->oo (for p=12 see constant A132265).
lim sup a(n)/b(n)=1, for n-->oo.
lim inf a(n)/n^((1+log_p(n))/2)=(sqrt(2)/2^log_p(sqrt(2)))*product{k>0, 1-1/(2*p^k)}=0.612870619..., for n-->oo, (p=12).
lim sup a(n)/n^((1+log_p(n))/2)=sqrt(q)/q^log_p(sqrt(q))=1.358593737..., for n-->oo, (p=12, q=round(sqrt(p))=3).
lim inf a(n)/a(n+1)=product{k>0, 1-1/(2*p^k)}=0.47735217025489380... for n-->oo (for p=12 see constant A132265).

A132269 a(n) = Product_{k>=0} (1 + floor(n/2^k)).

Original entry on oeis.org

1, 2, 6, 8, 30, 36, 56, 64, 270, 300, 396, 432, 728, 784, 960, 1024, 4590, 4860, 5700, 6000, 8316, 8712, 9936, 10368, 18200, 18928, 21168, 21952, 27840, 28800, 31744, 32768, 151470, 156060, 170100, 174960, 210900, 216600, 234000, 240000, 340956, 349272, 374616
Offset: 0

Author

Hieronymus Fischer, Aug 20 2007

Keywords

Comments

If n is written in base 2 as n=d(m)d(m-1)d(m-2)...d(2)d(1)d(0) (where d(k) is the digit at position k) then a(n) is also the product (1+d(m)d(m-1)d(m-2)...d(2)d(1)d(0))*(1+d(m)d(m-1)d(m-2)...d(2)d(1))*(1+d(m)d(m-1)d(m-2)...d(2))*...*(1+d(m)d(m-1)d(m-2))*(1+d(m)d(m-1))*(1+d(m)).
From Gary W. Adamson, Aug 25 2016: (Start)
Given the following production matrix M =
1, 0, 0, 0, 0, ...
2, 0, 0, 0, 0, ...
0, 3, 0, 0, 0, ...
0, 4, 0, 0, 0, ...
0, 0, 5, 0, 0, ...
0, 0, 6, 0, 0, ...
0, 0, 0, 7, 0, ...
...
the sequence is the left-shifted vector as lim_{n->infinity} M^n. (End)

Examples

			a(10) = (1 + floor(10/2^0))*(1 + floor(10/2^1))*(1 + floor(10/2^2))*(1 + floor(10/2^3)) = 11*6*3*2 = 396;
a(17) = 4860 since 17 = 10001_2 and so a(17) = (1+10001_2)*(1+1000_2)*(1+100_2)*(1+10_2)*(1+1) = 18*9*5*3*2 = 4860.
		

Crossrefs

For formulas regarding a general parameter p (i.e., terms 1+floor(n/p^k)) see A132271.
For the product of terms floor(n/p^k) see A098844, A067080, A132027-A132033, A132263, A132264.

Programs

  • Magma
    [1] cat [n le 1 select 2 else (1+n)*Self(Floor(n/2)): n in [1..50]]; // Vincenzo Librandi, Aug 26 2016
  • Maple
    f:= proc(n) option remember; (1+n)*procname(floor(n/2)) end proc:
    f(0):= 1:
    map(f, [$0..100]); # Robert Israel, Aug 26 2016
  • Mathematica
    Table[Product[1 + Floor[2 n/2^k], {k, 2 n}], {n, 0, 42}] (* or *)
    Table[Function[w, Times @@ Map[1 + FromDigits[PadRight[w, #], 2] &, Range@ Length@ w]]@ IntegerDigits[n, 2], {n, 0, 42}] (* Michael De Vlieger, Aug 26 2016 *)

Formula

Recurrence: a(n)=(1+n)*a(floor(n/2)); a(2n)=(1+2n)*a(n); a(n*2^m) = (Product_{k=1..m} (1 + n*2^k))*a(n).
a(2^m-1) = 2^(m*(m+1)/2), a(2^m) = 2^(m*(m+1)/2)*Product_{k=0..m} (1 + 1/2^k), m>=1.
a(n) = A132270(2n) = (1+n)*A132270(n).
Asymptotic behavior: a(n) = O(n^((1+log_2(n))/2)); this follows from the inequalities below.
a(n) <= A098844(n)*Product_{k=0..floor(log_2(n))} (1 + 1/2^k).
a(n) >= A098844(n)/Product_{k=1..floor(log_2(n))} (1 - 1/2^k).
a(n) < c*n^((1+log_2(n))/2) = c*2^A000217(log_2(n)), where c = Product_{k>=0} (1 + 1/2^k) = 4.7684620580627... (see constant A081845).
a(n) > n^((1+log_2(n))/2) = 2^A000217(log_2(n)),
lim sup a(n)/A098844(n) = Product_{k>=0} (1 + 1/2^k) = 4.7684620580627..., for n->oo (see constant A081845).
lim inf a(n)/A098844(n) = 1/Product_{k>=1} (1 - 1/2^k) = 1/0.288788095086602421..., for n->oo (see constant A048651).
lim inf a(n)/n^((1+log_2(n))/2) = 1, for n->oo.
lim sup a(n)/n^((1+log_2(n))/2) = Product_{k>=0} (1 + 1/2^k) = 4.7684620580627..., for n->oo (see constant A081845).
lim inf a(n+1)/a(n) = Product_{k>=0} (1 + 1/2^k) = 4.7684620580627... for n->oo (see constant A081845).
G.f. g(x) satisfies g(x) = (1+2x)*g(x^2) + 2*x^2*(1+x)*g'(x^2). - Robert Israel, Aug 26 2016

A132327 a(n) = Product{k>=0} (1 + floor(n/3^k)).

Original entry on oeis.org

1, 2, 3, 8, 10, 12, 21, 24, 27, 80, 88, 96, 130, 140, 150, 192, 204, 216, 399, 420, 441, 528, 552, 576, 675, 702, 729, 2240, 2320, 2400, 2728, 2816, 2904, 3264, 3360, 3456, 4810, 4940, 5070, 5600, 5740, 5880, 6450, 6600, 6750, 8832, 9024, 9216, 9996, 10200
Offset: 0

Author

Hieronymus Fischer, Aug 20 2007

Keywords

Comments

If n is written in base-3 as n=d(m)d(m-1)d(m-2)...d(2)d(1)d(0) (where d(k) is the digit at position k) then a(n) is also the product (1+d(m)d(m-1)d(m-2)...d(2)d(1)d(0))*(1+d(m)d(m-1)d(m-2)...d(2)d(1))*(1+d(m)d(m-1)d(m-2)...d(2))*...*(1+d(m)d(m-1)d(m-2))*(1+d(m)d(m-1))*(1+d(m)).

Examples

			a(12)=(1+floor(12/3^0))*(1+floor(12/3^1))*(1+floor(12/3^2))=13*5*2=130; a(20)=441 since 20=202(base-3) and so
a(20)=(1+202)*(1+20)*(1+2)(base-3)=21*7*3=441.
		

Crossrefs

Cf. A100220, A132027, A132038, A132264, A132269(for p=2), A132271(for p=10).
For formulas regarding a general parameter p (i.e. terms 1+floor(n/p^k)) see A132271.
For the product of terms floor(n/p^k) see A098844, A067080, A132027-A132033, A132263, A132264.

Programs

  • Mathematica
    Table[Product[1+Floor[n/3^k],{k,0,n}],{n,0,49}] (* James C. McMahon, Mar 07 2025 *)

Formula

Recurrence: a(n)=(1+n)*a(floor(n/3)); a(3n)=(1+3n)*a(n); a(n*3^m)=product{1<=k<=m, 1+n*3^k}*a(n).
a(k*3^m-j)=(k*3^m-j+1)*3^m*p^(m(m-1)/2), for 0=1, a(3^m)=3^(m(m+1)/2)*product{0<=k<=m, 1+1/3^k}, m>=1.
a(n)=A132328(3*n)=(1+n)*A132328(n).
Asymptotic behavior: a(n)=O(n^((1+log_3(n))/2)); this follows from the inequalities below.
a(n)<=A132027(n)*product{0<=k<=floor(log_3(n)), 1+1/3^k}.
a(n)>=A132027(n)/product{1<=k<=floor(log_3(n)), 1-1/3^k}.
a(n)A000217(log_3(n)), where c=product{k>=0, 1+1/p^k}=3.12986803713402307587769821345767... (see constant A132323).
a(n)>n^((1+log_3(n))/2)=3^A000217(log_3(n)).
lim sup a(n)/A132027(n)=2*product{k>0, 1+1/3^k}=3.12986803713402307587769821345767..., for n-->oo (see constant A132323).
lim inf a(n)/A132027(n)=1/product{k>0, 1-1/3^k}=1/0.560126077927948944969792243314140014..., for n-->oo (see constant A100220).
lim inf a(n)/n^((1+log_3(n))/2)=1, for n-->oo.
lim sup a(n)/n^((1+log_3(n))/2)=2*product{k>0, 1+1/3^k}=3.12986803713402307587769821345767..., for n-->oo (see constant A132323).
lim inf a(n+1)/a(n)=2*product{k>0, 1+1/3^k}=3.12986803713402307587769821345767... for n-->oo (see constant A132323).

A054896 a(n) = Sum_{k>0} floor(n/7^k).

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 8, 8, 8, 8, 8, 8, 8, 9, 9, 9, 9, 9, 9, 9, 10, 10, 10, 10, 10, 10, 10, 11, 11, 11, 11, 11, 11, 11, 12, 12, 12, 12, 12, 12, 12, 13, 13, 13, 13, 13, 13
Offset: 0

Author

Henry Bottomley, May 23 2000

Keywords

Comments

Exponent of the highest power of 7 dividing n!.

Examples

			  a(10^0) = 0.
  a(10^1) = 1.
  a(10^3) = 16.
  a(10^3) = 164.
  a(10^4) = 1665.
  a(10^5) = 16662.
  a(10^6) = 166664.
  a(10^7) = 1666661.
  a(10^8) = 16666662.
  a(10^9) = 166666661
		

Crossrefs

Cf. A011371 and A054861 for analogs involving powers of 2 and 3.

Programs

Formula

a(n) = floor(n/7) + floor(n/49) + floor(n/343) + floor(n/2401) + ...
a(n) = (n - A053828(n))/6.
From Hieronymus Fischer, Aug 14 2007: (Start)
a(n) = a(floor(n/7)) + floor(n/7).
a(7*n) = n + a(n).
a(n*7^m) = a(n) + n*(7^m-1)/6.
a(k*7^m) = k*(7^m-1)/6, for 0 <= k < 7, m >= 0.
Asymptotic behavior:
a(n) = n/6 + O(log(n)).
a(n+1) - a(n) = O(log(n)); this follows from the inequalities below.
a(n) <= (n-1)/6; equality holds for powers of 7.
a(n) >= (n-6)/6 - floor(log_7(n)); equality holds for n=7^m-1, m>0. -
lim inf (n/6 - a(n)) = 1/6, for n-->oo.
lim sup (n/6 - log_7(n) - a(n)) = 0, for n-->oo.
lim sup (a(n+1) - a(n) - log_7(n)) = 0, for n-->oo.
G.f.: (1/(1-x))*Sum_{k > 0} x^(7^k)/(1-x^(7^k)). (End)
Partial sums of A214411. - R. J. Mathar, Jul 08 2021
a(n) = A214411(A000142(n)). - Michel Marcus, Oct 07 2024

Extensions

Examples added by Hieronymus Fischer, Jun 06 2012

A132035 Decimal expansion of Product_{k>0} (1-1/7^k).

Original entry on oeis.org

8, 3, 6, 7, 9, 5, 4, 0, 7, 0, 8, 9, 0, 3, 7, 8, 7, 1, 0, 2, 6, 7, 2, 9, 7, 9, 8, 1, 4, 6, 1, 3, 6, 2, 4, 1, 3, 5, 2, 4, 3, 6, 4, 3, 5, 8, 7, 6, 7, 1, 6, 5, 1, 9, 9, 6, 4, 1, 1, 5, 1, 0, 1, 7, 7, 0, 0, 9, 1, 6, 0, 1, 2, 6, 5, 4, 2, 7, 6, 0, 5, 8, 7, 8, 7, 5, 5, 5, 4, 2, 8, 4, 9, 0, 5, 1, 2, 0, 2, 1, 7, 5, 3
Offset: 0

Author

Hieronymus Fischer, Aug 14 2007

Keywords

Examples

			0.8367954070890378710...
		

Programs

  • Mathematica
    digits = 103; NProduct[1-1/7^k, {k, 1, Infinity}, NProductFactors -> 200, WorkingPrecision -> digits+3] // N[#, digits+3]& // RealDigits[#, 10, digits]& // First (* Jean-François Alcover, Feb 18 2014 *)
    RealDigits[QPochhammer[1/7], 10, 100][[1]] (* Amiram Eldar, May 09 2023 *)
  • PARI
    prodinf(k=1, 1 - 1/(7^k)) \\ Amiram Eldar, May 09 2023

Formula

Equals exp(-Sum_{n>0} sigma_1(n)/(n*7^n)) = exp(-Sum_{n>0} A000203(n)/(n*7^n)).
From Amiram Eldar, May 09 2023: (Start)
Equals sqrt(2*Pi/log(7)) * exp(log(7)/24 - Pi^2/(6*log(7))) * Product_{k>=1} (1 - exp(-4*k*Pi^2/log(7))) (McIntosh, 1995).
Equals Sum_{n>=0} (-1)^n/A027875(n). (End)

A054895 a(n) = Sum_{k>0} floor(n/6^k).

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8, 8, 9, 9, 9, 9, 9, 9, 10, 10, 10, 10, 10, 10, 11, 11, 11, 11, 11, 11, 12, 12, 12, 12, 12, 12, 14, 14, 14, 14, 14, 14, 15, 15, 15, 15, 15, 15, 16, 16, 16
Offset: 0

Author

Henry Bottomley, May 23 2000

Keywords

Comments

Different from the highest power of 6 dividing n! (cf. A054861). - Hieronymus Fischer, Aug 14 2007
Partial sums of A122841. - Hieronymus Fischer, Jun 06 2012

Examples

			  a(10^0) = 0.
  a(10^1) = 1.
  a(10^2) = 18.
  a(10^3) = 197.
  a(10^4) = 1997.
  a(10^5) = 19996.
  a(10^6) = 199995.
  a(10^7) = 1999995.
  a(10^8) = 19999994.
  a(10^9) = 199999993.
		

Crossrefs

Cf. A011371 and A054861 for analogs involving powers of 2 and 3.

Programs

  • Haskell
    a054895 n = a054895_list !! n
    a054895_list = scanl (+) 0 a122841_list
    -- Reinhard Zumkeller, Nov 10 2013
    
  • Magma
    function A054895(n)
      if n eq 0 then return n;
      else return A054895(Floor(n/6)) + Floor(n/6);
      end if; return A054895;
    end function;
    [A054895(n): n in [0..100]]; // G. C. Greubel, Feb 09 2023
    
  • Mathematica
    Table[t=0; p=6; While[s=Floor[n/p]; t=t+s; s>0, p *= 6]; t, {n,0,100}]
  • SageMath
    def A054895(n):
        if (n==0): return 0
        else: return A054895(n//6) + (n//6)
    [A054895(n) for n in range(104)] # G. C. Greubel, Feb 09 2023

Formula

a(n) = floor(n/6) + floor(n/36) + floor(n/216) + floor(n/1296) + ...
a(n) = (n - A053827(n))/5.
From Hieronymus Fischer, Aug 14 2007: (Start)
a(n) = a(floor(n/6)) + floor(n/6).
a(6*n) = n + a(n).
a(n*6^m) = n*(6^m-1)/5 + a(n).
a(k*6^m) = k*(6^m-1)/5, for 0 <= k < 6, m >= 0.
Asymptotic behavior:
a(n) = (n/5) + O(log(n)).
a(n+1) - a(n) = O(log(n)); this follows from the inequalities below.
a(n) <= (n-1)/5; equality holds for powers of 6.
a(n) >= ((n-5)/5) - floor(log_6(n)); equality holds for n=6^m-1, m>0.
lim inf (n/5 - a(n)) = 1/5, for n-->oo.
lim sup (n/5 - log_6(n) - a(n)) = 0, for n-->oo.
lim sup (a(n+1) - a(n) - log_6(n)) = 0, for n-->oo.
G.f.: (1/(1-x))*Sum_{k > 0} x^(6^k)/(1-x^(6^k)). (End)

Extensions

An incorrect formula was deleted by N. J. A. Sloane, Nov 18 2008
Examples added by Hieronymus Fischer, Jun 06 2012

A132326 Decimal expansion of Product_{k>=1} (1+1/10^k).

Original entry on oeis.org

1, 1, 1, 2, 2, 3, 4, 5, 6, 9, 1, 3, 7, 0, 5, 0, 6, 3, 2, 1, 2, 6, 0, 7, 8, 0, 6, 7, 0, 9, 4, 4, 0, 5, 8, 0, 3, 7, 4, 7, 5, 0, 7, 4, 6, 7, 5, 7, 7, 5, 9, 2, 8, 3, 5, 7, 8, 7, 9, 5, 8, 2, 3, 7, 0, 3, 3, 2, 5, 3, 4, 6, 9, 4, 8, 8, 1, 4, 1, 1, 0, 4, 3, 7, 6, 4, 7, 2, 2, 2, 2, 6, 4, 2, 1, 3, 5, 2, 3, 5, 5, 6, 4, 7, 4
Offset: 1

Author

Hieronymus Fischer, Aug 20 2007

Keywords

Comments

Half the constant A132325.

Examples

			1.1122345691370506321260780670944...
		

Programs

  • Mathematica
    digits = 105; NProduct[1+1/3^k, {k, 0, Infinity}, NProductFactors -> 100, WorkingPrecision -> digits+3] // N[#, digits+3]& // RealDigits[#, 10, digits]& // First (* Jean-François Alcover, Feb 18 2014 *)
    N[QPochhammer[-1/10,1/10]] (* G. C. Greubel, Dec 01 2015 *)
  • PARI
    prodinf(k=1, 1+1/10^k) \\ Amiram Eldar, May 20 2023

Formula

Equals (1/2)*lim sup_{n->oo} Product_{0<=k<=floor(log_10(n))} (1+1/floor(n/10^k)).
Equals (1/2)*lim sup_{n->oo} A132271(n)/n^((1+log_10(n))/2).
Equals (1/2)*lim sup_{n->oo} A132272(n)/n^((log_10(n)-1)/2).
Equals exp(Sum_{n>0} 10^(-n)*Sum_{k|n} -(-1)^k/k) = exp(Sum_{n>0} A000593(n)/(n*10^n)).
Equals (1/2)*lim sup_{n->oo} A132271(n+1)/A132271(n).
Equals (-1/10; 1/10){infinity}, where (a;q){infinity} is the q-Pochhammer symbol. - G. C. Greubel, Dec 01 2015
Equals (sqrt(2)/2) * exp(log(10)/24 + Pi^2/(12*log(10))) * Product_{k>=1} (1 - exp(-2*(2*k-1)*Pi^2/log(10))) (McIntosh, 1995). - Amiram Eldar, May 20 2023

A067079 Right digitorial function Rd(n) = abcd...*bcd...*cd...*d...*..., for n = abcd... in decimal notation.

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 0, 11, 24, 39, 56, 75, 96, 119, 144, 171, 0, 21, 44, 69, 96, 125, 156, 189, 224, 261, 0, 31, 64, 99, 136, 175, 216, 259, 304, 351, 0, 41, 84, 129, 176, 225, 276, 329, 384, 441, 0, 51, 104, 159, 216, 275, 336, 399, 464, 531, 0, 61, 124, 189
Offset: 0

Author

Amarnath Murthy, Jan 05 2002

Keywords

Comments

Rd(abcd) = abcd*bcd*cd*d, etc.

Examples

			Rd(256) = 256*56*6 = 86016.
		

Crossrefs

Cf. A067080.

Programs

  • Haskell
    import Data.List (tails)
    a067079 = product . map read . init . tails . show :: Integer -> Integer
    -- Reinhard Zumkeller, Apr 23 2013
    
  • Mathematica
    Table[d = IntegerDigits[n]; rd = 1; While[ Length[d] > 0, rd = rd*FromDigits[d]; d = Drop[d, 1]]; rd, {n, 1, 75} ]
    Join[{0},Table[Times@@FromDigits/@NestList[Rest[#]&,IntegerDigits[n], IntegerLength[ n]-1],{n,70}]] (* Harvey P. Dale, Nov 29 2014 *)
  • PARI
    A067079(n)=prod(k=1,logint(n+!n,10)+1,n-n\10^k*10^k) \\ M. F. Hasler, Feb 25 2019

Formula

a(n) = Product_{k=1..length(n)} (n-10^k*floor(n/10^k)). - Vladeta Jovovic, Jan 08 2002
a(n) = 0 iff n == 0 (mod 10). - M. F. Hasler, Feb 25 2019

Extensions

More terms from Robert G. Wilson v, Jan 07 2002

A132020 Decimal expansion of Product_{k>=0} (1 - 1/(2*4^k)).

Original entry on oeis.org

4, 1, 9, 4, 2, 2, 4, 4, 1, 7, 9, 5, 1, 0, 7, 5, 9, 7, 7, 0, 9, 9, 5, 6, 1, 0, 7, 7, 0, 2, 9, 7, 4, 2, 5, 2, 2, 3, 3, 9, 5, 3, 2, 3, 4, 3, 9, 2, 6, 6, 6, 7, 4, 9, 0, 8, 0, 4, 4, 9, 9, 1, 6, 6, 3, 1, 7, 7, 2, 0, 5, 0, 8, 7, 2, 7, 0, 9, 1, 9, 3, 9, 1, 0, 0, 2, 3, 2, 4, 5, 4, 7, 4, 2, 3, 8, 1, 9, 5, 5, 0, 2, 8, 5, 8
Offset: 0

Author

Hieronymus Fischer, Aug 14 2007

Keywords

Comments

This is the limiting probability that a large random symmetric binary matrix is nonsingular (cf. A086812, A048651). In other words, equals Lim_{n->oo} A086812(n)/A006125(n+1).- H. Tracy Hall, Sep 07 2024

Examples

			0.41942244179510759770995610770297425223395323439266674908044991663177...
		

Programs

  • Maple
    evalf(1+sum((-1)^n*2^(n*(n-1)/2)/product(2^k-1, k=1..n), n=1..infinity), 120); # Robert FERREOL, Feb 23 2020
  • Mathematica
    RealDigits[ Product[1 - 1/(2*4^i), {i, 0, 175}], 10, 111][[1]] (* Robert G. Wilson v, May 25 2011 *)
    RealDigits[QPochhammer[1/2, 1/4], 10, 105][[1]] (* Jean-François Alcover, Nov 18 2015 *)
  • PARI
    prodinf(k=0,1-1.>>(2*k+1)) \\ Charles R Greathouse IV, Nov 16 2012

Formula

Equals lim inf_{n->oo} Product_{k=0..floor(log_4(n))} floor(n/4^k)*4^k/n.
Equals lim inf_{n->oo} A132028(n)/n^(1+floor(log_4(n)))*4^((1/2)*(1+floor(log_4(n)))*floor(log_4(n))).
Equals lim inf_{n->oo} A132028(n)/n^(1+floor(log_4(n)))*4^A000217(floor(log_4(n))).
Equals (1/2)*exp(-Sum_{n>0} (4^(-n)*(Sum_{k|n} 1/(k*2^k)))).
Equals lim inf_{n->oo} A132028(n)/A132028(n+1).
Equals Product_{k>0} (1-1/(2^k+1)). - Robert G. Wilson v, May 25 2011
From Robert FERREOL, Feb 23 2020: (Start)
Equals Product_{k>0} (1 + 1/2^k)^(-1) = 2/A081845.
Equals 1 + Sum_{n>=1} (-1)^n*2^(n*(n-1)/2)/((2-1)*(2^2-1)*...*(2^n-1)). (End)
From Peter Bala, Jan 15 2021: (Start)
Constant C = Sum_{n >= 0} 2^n/Product_{k = 1..n} (1 - 4^k).
Faster converging series:
2*C = (1/2)*Sum_{n >= 0} 2^(-n)/Product_{k = 1..n} (1 - 4^k);
(2^4)*C = 7*Sum_{n >= 0} 2^(-3*n)/Product_{k = 1..n} (1 - 4^k);
(2^9)*C = 7*31*Sum_{n >= 0} 2^(-5*n)/Product_{k = 1..n} (1 - 4^k), and so on.
Slower converging series:
C = -Sum_{n >= 0} 2^(3*n)/Product_{k = 1..n} (1 - 4^k);
7*C = Sum_{n >= 0} 2^(5*n)/Product_{k = 1..n} (1 - 4^k);
7*31*C = -Sum_{n >= 0} 2^(7*n)/Product_{k = 1..n} (1 - 4^k), and so on. (End)
Equals Product_{n>=0} (1 - 1/A004171(n)). - Amiram Eldar, May 09 2023

Extensions

Name corrected by Charles R Greathouse IV, Nov 16 2012
Previous Showing 21-30 of 49 results. Next