|
%I #24 Oct 31 2022 11:47:57
%S 0,1,17,33,289,305,561,817,4913,4929,5185,5441,9537,9793,13889,17985,
%T 83521,83537,83793,84049,88145,88401,92497,96593,162129,162385,166481,
%U 170577,236113,240209,305745,371281,1419857,1419873,1420129,1420385
%N a(0)=1, a(1)=1, a(n) = 17*a(n/2) for n=2,4,6,..., a(n) = 16*a((n-1)/2) + a((n+1)/2) for n=3,5,7,....
%C A 17-divide version of A084230.
%C The Harborth : f(2^k)=3^k suggests that a family of sequences of the form: f(2^k)=Prime[n]^k There does indeed seem to be an infinite family of such functions.
%H G. C. Greubel, <a href="/A116523/b116523.txt">Table of n, a(n) for n = 0..2500</a>
%H H. Harborth, <a href="http://dx.doi.org/10.1090/S0002-9939-1977-0429714-1">Number of Odd Binomial Coefficients</a>, Proc. Amer. Math. Soc. 62, 19-22, 1977.
%H Hsien-Kuei Hwang, Svante Janson, and Tsung-Hsi Tsai, <a href="https://arxiv.org/abs/2210.10968">Identities and periodic oscillations of divide-and-conquer recurrences splitting at half</a>, arXiv:2210.10968 [cs.DS], 2022, pp. 27, 33.
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/Stolarsky-HarborthConstant.html">Stolarsky-Harborth Constant</a>
%F a(n) = Sum_{k=0..n-1} 16^wt(k), where wt = A000120. - _Mike Warburton_, Mar 22 2019
%p a:=proc(n) if n=0 then 0 elif n=1 then 1 elif n mod 2 = 0 then 17*a(n/2) else 16*a((n-1)/2)+a((n+1)/2) fi end: seq(a(n),n=0..38);
%t b[0] := 0; b[1] := 1; b[n_?EvenQ] := b[n] = 17*b[n/2]; b[n_?OddQ] := b[n] = 16*b[(n - 1)/2] + b[(n + 1)/2]; a = Table[b[n], {n, 1, 25}]
%Y Cf. A000120, A006046, A077465, A084230.
%K nonn
%O 0,3
%A _Roger L. Bagula_, Mar 15 2006
%E Edited by _N. J. A. Sloane_, Apr 16 2005
|
