The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.



(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A063511 a(n) = a(floor(square root(n))) * 2. 2


%S 1,2,2,4,4,4,4,4,4,4,4,4,4,4,4,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,

%T 8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,

%U 8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8

%N a(n) = a(floor(square root(n))) * 2.

%C From _Kevin Ryde_, May 11 2020: (Start)

%C The sqrt steps in the definition are equivalent to A211667 but here factors of 2 instead of counting, so a(n) = 2^A211667(n). A211667 is a double logarithm and the effect of power 2^ is to turn the second into a rounding. So a(n) is the bit length of n (see A070939) increased to the next power of 2 if not already a power of 2. Each n = 2^(2^k) is a new high a(n) = 2^(k+1), since such an n is bit length 2^k+1.

%C In a microcomputer, it's common for machine words to be power-of-2 sizes such as 16, 32, 64, 128 bits. a(n) can be thought of as the word size needed to contain integer n. Some algorithms by their nature expect power-of-2 sizes, for example Schönhage and Strassen's big integer multiplication.

%C This sequence differs from A334789 (2^log*(n)) for n>=256. For example a(256)=16 whereas A334789(256)=8. The respective exponent sequences are A211667 (for here) and A001069 (for A334789) which likewise differ for n>=256.

%C (End)

%H Kevin Ryde, <a href="/A063511/b063511.txt">Table of n, a(n) for n = 1..8192</a>

%H Martin Fürer, <a href="http://web.archive.org/web/1id_/http://www.cse.psu.edu/~furer/Papers/mult.pdf">Faster integer multiplication</a>, Proceedings of the 39th Annual ACM Symposium on Theory of Computing, June 11-13 2007. And <a href="http://dx.doi.org/10.1137/070711761">in SIAM Journal of Computing</a>, volume 30, number 3, 2009, pages 979-1005. (See size "n" calculation at the start of Algorithm Integer-Multiplication.)

%H <a href="/index/Di#divseq">Index to divisibility sequences</a>

%F a(n) = 2^A211667(n) = 2^ceiling(log_2(log_2(n+1))). - _Kevin Ryde_, May 11 2020

%o (PARI) a(n) = if(n==1,1, 2<<logint(logint(n,2),2)); \\ _Kevin Ryde_, May 11 2020

%Y Cf. A001146 (indices of new highs), A334789.

%K easy,nonn,changed

%O 1,2

%A _Reinhard Zumkeller_, Jul 30 2001

%E Formula and code by Charles R Greathouse IV moved to A334789 where they apply. - _Kevin Ryde_, May 11 2020

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified June 4 01:54 EDT 2020. Contains 334812 sequences. (Running on oeis4.)