A067080 If n = ab...def in decimal notation then the left digitorial function Ld(n) = ab...def*ab...de*ab...d*...*ab*a.
1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 40, 42, 44, 46, 48, 50, 52, 54, 56, 58, 90, 93, 96, 99, 102, 105, 108, 111, 114, 117, 160, 164, 168, 172, 176, 180, 184, 188, 192, 196, 250, 255, 260, 265, 270, 275, 280, 285, 290, 295, 360, 366, 372
Offset: 1
Examples
Ld(256) = 256*25*2 =12800. a(31)=floor(31/10^0)*floor(31/10^1)=31*3=93; a(42)=168 since 42=42(base-10) and so a(42)=42*4(base-10)=42*4=168.
Links
- Hieronymus Fischer, Table of n, a(n) for n = 1..10000
Crossrefs
Programs
-
Haskell
a067080 n = if n <= 9 then n else n * a067080 (n `div` 10) -- Reinhard Zumkeller, Nov 29 2012
-
Mathematica
Table[d = IntegerDigits[n]; rd = 1; While[ Length[d] > 0, rd = rd*FromDigits[d]; d = Drop[d, -1]]; rd, {n, 1, 75} ] Table[Times@@NestList[Quotient[#,10]&,n,IntegerLength[n]-1],{n,70}] (* Harvey P. Dale, Dec 16 2013 *)
-
PARI
a(n)=my(t=n);while(n\=10,t*=n); t \\ Charles R Greathouse IV, Nov 20 2012
Formula
a(n) = Product_{k=1..length(n)} floor(n/10^(k-1)). - Vladeta Jovovic, Jan 08 2002
From Hieronymus Fischer, Aug 13 2007: (Start)
a(n) = product{0<=k<=floor(log_10(n)), floor(n/10^k)}, n>=1.
Recurrence:
a(n) = n*a(floor(n/10));
a(n*10^m) = n^m*10^(m(m+1)/2)*a(n).
a(k*10^m) = k^(m+1)*10^(m(m+1)/2), for 0
a(n) <= b(n), where b(n)=n^(1+floor(log_10(n)))/10^(1/2*(1+floor(log_10(n)))*floor(log_10(n))); equality holds for n=k*10^m, m>=0, 1<=k<10. Here b(n) can also be written n^(1+floor(log_10(n)))/10^A000217(floor(log_10(n))).
Also: a(n) <= 3^((1-log_10(3))/2)*n^((1+log_10(n))/2)=1.332718...*10^A000217(log_10(n)), equality for n=3*10^m, m>=0.
a(n) > c*b(n), where c=0.472362443816572... (see constant A132026).
Also: a(n) > c*2^((1-log_10(2))/2)*n^((1+log_10(n))/2) = 0.601839...*10^A000217(log_10(n)).
lim inf a(n)/b(n) = 0.472362443816572..., for n-->oo.
lim sup a(n)/b(n) = 1, for n-->oo.
lim inf a(n)/n^((1+log_10(n))/2) = 0.472362443816572...*sqrt(2)/2^log_10(sqrt(2)), for n-->oo.
lim sup a(n)/n^((1+log_10(n))/2) = sqrt(3)/3^log_10(sqrt(3)), for n-->oo.
lim inf a(n)/a(n+1) = 0.472362443816572... for n-->oo (see constant A132026).
a(n) = O(n^((1+log_10(n))/2)). (End)
Extensions
More terms from Robert G. Wilson v, Jan 07 2002
A098844 a(1)=1, a(n) = n*a(floor(n/2)).
1, 2, 3, 8, 10, 18, 21, 64, 72, 100, 110, 216, 234, 294, 315, 1024, 1088, 1296, 1368, 2000, 2100, 2420, 2530, 5184, 5400, 6084, 6318, 8232, 8526, 9450, 9765, 32768, 33792, 36992, 38080, 46656, 47952, 51984, 53352, 80000, 82000, 88200, 90300
Offset: 1
Keywords
Examples
a(10) = floor(10/2^0)*floor(10/2^1)*floor(10/2^2)*floor(10/2^3) = 10*5*2*1 = 100; a(17) = 1088 since 17 = 10001(base 2) and so a(17) = 10001*1000*100*10*1(base 2) = 17*8*4*2*1 = 1088.
Links
- Hieronymus Fischer, Table of n, a(n) for n = 1..1000
Crossrefs
Programs
-
Mathematica
lst={};Do[p=n;s=1;While[p>1,p=IntegerPart[p/2];s*=p;];AppendTo[lst,s],{n,1,6!,2}];lst (* Vladimir Joseph Stephan Orlovsky, Jul 28 2009 *)
-
PARI
a(n)=if(n<2,1,n*a(floor(n/2)))
-
Python
from math import prod def A098844(n): return n*prod(n//2**k for k in range(1,n.bit_length()-1)) # Chai Wah Wu, Jun 07 2022
Formula
a(2^n) = 2^(n*(n+1)/2) = A006125(n+1).
From Hieronymus Fischer, Aug 13 2007: (Start)
a(n) = Product_{k=0..floor(log_2(n))} floor(n/2^k), n>=1.
Recurrence:
a(n*2^m) = n^m*2^(m(m+1)/2)*a(n).
a(n) <= n^((1+log_2(n))/2) = 2^A000217(log_2(n)); equality iff n is a power of 2.
a(n) >= c(n)*(n+1)^((1 + log_2(n+1))/2) for n != 2,
where c(n) = Product_{k=1..floor(log_2(n))} (1 - 1/2^k); equality holds iff n+1 is a power of 2.
a(n) > c*(n+1)^((1 + log_2(n+1))/2)
where c = 0.288788095086602421... (see constant A048651).
lim inf a(n)/n^((1+log_2(n))/2)=0.288788095086602421... for n-->oo.
lim sup a(n)/n^((1+log_2(n))/2) = 1 for n-->oo.
lim inf a(n)/a(n+1) = 0.288788095086602421... for n-->oo (see constant A048651).
a(n) = O(n^((1+log_2(n))/2)). (End)
Extensions
Formula section edited by Hieronymus Fischer, Jun 13 2012
A132027 a(n) = Product_{k=0..floor(log_3(n))} floor(n/3^k), n>=1.
1, 2, 3, 4, 5, 12, 14, 16, 27, 30, 33, 48, 52, 56, 75, 80, 85, 216, 228, 240, 294, 308, 322, 384, 400, 416, 729, 756, 783, 900, 930, 960, 1089, 1122, 1155, 1728, 1776, 1824, 2028, 2080, 2132, 2352, 2408, 2464, 3375, 3450, 3525, 3840, 3920, 4000, 4335
Offset: 1
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 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(11)=floor(11/3^0)*floor(11/3^1)*floor(11/3^2)=11*3*1=33; a(13)=52 since 13=111(base-3) and so a(13)=111*11*1(base-3)=13*4*1=52.
Links
- Ivan Neretin, Table of n, a(n) for n = 1..10000
Crossrefs
Programs
-
Mathematica
Table[(f = If[# < 3, #, #*f[Quotient[#, 3]]] &)[n], {n, 51}] (* Ivan Neretin, May 29 2016 *)
Formula
Recurrence: a(n)=n*a(floor(n/3)); a(n*3^m)=n^m*3^(m(m+1)/2)*a(n).
a(k*3^m)=k^(m+1)*3^(m(m+1)/2), for k=1 or 2.
a(n)<=b(n), where b(n)=n^(1+floor(log_3(n)))/3^(1/2*(1+floor(log_3(n)))*floor(log_3(n))); equality holds if n is a power of 3 or two times a power of 3.
Also: a(n)<=2^((1-log_3(2))/2)*n^((1+log_3(n))/2)=1.1364507...*3^A000217(log_3(n)), equality for n=2*3^m, m>=0.
a(n)>c*b(n), where c=0.3826631966790330232889550... (see constant A132019).
Also: a(n)>c*2^((1-log_3(2))/2)*n^((1+log_3(n))/2)=0.434877...*3^A000217(log_3(n)).
lim inf a(n)/b(n)=0.3826631966790330232889550..., for n-->oo.
lim sup a(n)/b(n)=1, for n-->oo.
lim inf a(n)/n^((1+log_3(n))/2)=0.3826631966790330232889550...*sqrt(2)/2^log_3(sqrt(2)), for n-->oo.
lim sup a(n)/n^((1+log_3(n))/2)=sqrt(2)/2^log_3(sqrt(2)), for n-->oo.
lim inf a(n)/a(n+1)=0.3826631966790330232889550... for n-->oo (see constant A132019).
a(n)=O(n^((1+log_3(n))/2)).
A132033 Product{0<=k<=floor(log_9(n)), floor(n/9^k)}, n>=1.
1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 36, 38, 40, 42, 44, 46, 48, 50, 52, 81, 84, 87, 90, 93, 96, 99, 102, 105, 144, 148, 152, 156, 160, 164, 168, 172, 176, 225, 230, 235, 240, 245, 250, 255, 260, 265, 324, 330, 336, 342, 348, 354, 360, 366, 372
Offset: 1
Comments
If n is written in base-9 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(85)=floor(85/9^0)*floor(85/9^1)*floor(85/9^2)=85*9*1=765; a(88)=792 since 88=107(base-9) and so a(88)=107*10*1(base-9)=88*9*1=792.
Crossrefs
Programs
-
Mathematica
Table[Product[Floor[n/9^k],{k,0,Floor[Log[9,n]]}],{n,62}] (* James C. McMahon, Mar 03 2025 *)
Formula
Recurrence: a(n)=n*a(floor(n/9)); a(n*9^m)=n^m*9^(m(m+1)/2)*a(n).
a(k*9^m)=k^(m+1)*9^(m(m+1)/2), for 0
Asymptotic behavior: a(n)=O(n^((1+log_9(n))/2)); this follows from the inequalities below.
a(n)<=b(n), where b(n)=n^(1+floor(log_9(n)))/9^((1+floor(log_9(n)))*floor(log_9(n))/2); equality holds for n=k*9^m, 0=0. b(n) can also be written n^(1+floor(log_9(n)))/9^A000217(floor(log_9(n))).
Also: a(n)<=3^(1/4)*n^((1+log_9(n))/2)=1.316074013...*9^A000217(log_9(n)), equality holds for n=3*9^m, m>=0.
a(n)>c*b(n), where c=0.4689451783670236932832800... (see constant A132024).
Also: a(n)>c*2^((1-log_9(2))/2)*n^((1+log_9(n))/2)=0.4689451783670...*1.267747616...*9^A000217(log_9(n)).
lim inf a(n)/b(n)=0.4689451783670236932832800..., for n-->oo.
lim sup a(n)/b(n)=1, for n-->oo.
lim inf a(n)/n^((1+log_9(n))/2)=0.4689451783670236932832800...*sqrt(2)/2^log_9(sqrt(2)), for n-->oo.
lim sup a(n)/n^((1+log_9(n))/2)=3^(1/4)=1.316074013..., for n-->oo.
lim inf a(n)/a(n+1)=0.4689451783670236932832800... for n-->oo (see constant A132025).
A132328 a(n) = Product_{k>0} (1+floor(n/3^k)).
1, 1, 1, 2, 2, 2, 3, 3, 3, 8, 8, 8, 10, 10, 10, 12, 12, 12, 21, 21, 21, 24, 24, 24, 27, 27, 27, 80, 80, 80, 88, 88, 88, 96, 96, 96, 130, 130, 130, 140, 140, 140, 150, 150, 150, 192, 192, 192, 204, 204, 204, 216, 216, 216, 399, 399, 399, 420, 420, 420, 441, 441, 441, 528
Offset: 0
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))*(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^1))*(1+floor(12/3^2))=5*2=10; a(19)=21 since 19=201(base-3) and so a(19)=(1+20)*(1+2)(base-3)=7*3=21.
Links
- Robert Israel, Table of n, a(n) for n = 0..10000
Crossrefs
Programs
-
Maple
f:= proc(n) option remember; local t; t:= floor(n/3); (1+t)*procname(t) end proc: f(0):= 1: f(1):= 1: f(2):= 1: map(f, [$0..100]); # Robert Israel, Oct 20 2020
-
Mathematica
(* Using definition *) Table[Product[1 + Floor[n/3^k], {k, IntegerLength[n, 3] - 1}], {n, 0, 100}] (* Using recurrence -- faster *) a[0] = 1; a[n_] := a[n] = (1 + #)*a[#] & [Floor[n/3]]; Table[a[n], {n, 0, 100}] (* Paolo Xausa, Sep 23 2024 *)
Formula
Recurrence: a(n)=(1+floor(n/3))*a(floor(n/3)); a(3n)=(1+n)*a(n); a(n*3^m)=product{0<=k
a(k*3^m-j)=k^m*3^(m(m-1)/2), for 0=1. a(3^m)=p^(m(m-1)/2)*product{0<=k
Asymptotic behavior: a(n)=O(n^((log_3(n)-1)/p)); this follows from the inequalities below.
a(n)<=A132027(n)/(n+1)*product{0<=k<=floor(log_3(n)), 1+1/3^k}.
a(n)>=A132027(n)/((n+1)*product{0
a(n)A000217(log_3(n))/(n+1), where c=product{k>0, 1+1/3^k}=3.12986803713402307587769821345767... (see constant A132323).
a(n)>n^((1+log_3(n))/2)/(n+1)=3^A000217(log_3(n))/(n+1).
lim sup n*a(n)/A132027(n)=2*product{k>0, 1+1/3^k}=3.12986803713402307587769821345767..., for n-->oo (see constant A132323).
lim inf n*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).
A132263 Product{0<=k<=floor(log_11(n)), floor(n/11^k)}, n>=1.
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
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
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.
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
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
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)).
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
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.
Links
- Robert Israel, Table of n, a(n) for n = 0..10000
Crossrefs
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.
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
A132326 Decimal expansion of Product_{k>=1} (1+1/10^k).
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
Comments
Half the constant A132325.
Examples
1.1122345691370506321260780670944...
Links
- G. C. Greubel, Table of n, a(n) for n = 1..1200
- Richard J. McIntosh, Some Asymptotic Formulae for q-Hypergeometric Series, Journal of the London Mathematical Society, Vol. 51, No. 1 (1995), pp. 120-136; alternative link.
Crossrefs
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/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
A179051 Number of partitions of n into powers of 10 (cf. A011557).
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10
Offset: 0
Keywords
Examples
a(19) = #{10 + 9x1, 19x1} = 2; a(20) = #{10 + 10, 10 + 10x1, 20x1} = 3; a(21) = #{10 + 10 + 1, 10 + 11x1, 21x1} = 3.
Links
Crossrefs
Programs
-
Haskell
a179051 = p 1 where p _ 0 = 1 p k m = if m < k then 0 else p k (m - k) + p (k * 10) m -- Reinhard Zumkeller, Feb 05 2012
-
Mathematica
terms = 10001; CoefficientList[Product[1/(1 - x^(10^k)) + O[x]^terms, {k, 0, Log[10, terms] // Ceiling}], x] (* Jean-François Alcover, Dec 12 2021, after Ilya Gutkovskiy *)
Comments