login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A125508 Lowest number having a particular shape of factor decomposition binary tree. 2
1, 2, 4, 8, 16, 20, 32, 40, 64, 72, 88, 128, 160, 176, 200, 220, 256, 272, 288, 320, 336, 360, 400, 420, 460, 480, 512, 540, 544, 640, 704, 864, 880, 920, 1024, 1056, 1152, 1184, 1200, 1280, 1344, 1440, 1600, 1640, 1680, 1800, 1840, 1920, 2048, 2176, 2368 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

COMMENTS

For a natural number, n, make it the root node of a binary tree. The left child node (L) is the largest divisor of n which is greater than 1 but less than or equal to the square root of n, if this exists. The right child node is n/L, if the left node exists. Thus if n is a prime it is a leaf node; otherwise if it is composite then it is the product of its two children. If n = 1 then we have an empty tree. A term in this sequence a(n) is such that no number (x) exists 1<=x<a(n) with an isomorphic binary tree.

LINKS

Alois P. Heinz, Table of n, a(n) for n = 1..1000

Gordon Hamilton, Integral Fission, Video for grade 7 teachers

EXAMPLE

a(8) = 40, make this the root node.

40 is first decomposed into 5 and 8; 5 is a leaf node.

8 is decomposed into 2 and 4; 2 is a leaf node.

4 is decomposed into 2 and 2; both are leaf nodes.

There is no number x (where 1 <= x < 40) which has a factor decomposition binary tree isomorphic to this.

PROG

(PARI) fissr(n) = {if (isprime(n), n, my(d = divisors(n)); my(nddt = #d\2 + 1); if (#d % 2, [fissr(d[nddt]), n, fissr(d[nddt])], [fissr(d[nddt-1]), n, fissr(d[nddt])]); ); }

fiss(n) = if (n==1, [], if (type(fv = fissr(n))== "t_INT", [fv], fv));

empty(v) = {my(nv = vector(#v)); for (i=1, #v, if (type(v[i]) == "t_INT", nv[i] = 0, nv[i] = empty(v[i])); ); nv; }

eqvec(va, vb) = {if (type(va) != type(vb), return (0)); if (type(va) == "t_INT", return (va == vb)); if (#va != #vb, return (0)); for (i=1, #va, if (!eqvec(va[i], vb[i]), return (0)); ); return (1); }

isalready(erep, vrep) = {if (! #vrep, return (0)); for (j=1, #vrep, if (eqvec(erep, vrep[j]), return (1); ); ); return (0); }

addrep(erep, vrep) = {nvrep = vector(#vrep+1, i, if (i <= #vrep, vrep[i], 0)); nvrep[#vrep+1] = erep; nvrep; }

lista(nn) = {vrep = []; for (n=1, nn, rep = fiss(n); erep = empty(rep); if (! isalready(erep, vrep), print1(n, ", "); vrep = addrep(erep, vrep); ); ); } \\ Michel Marcus, May 25 2014

CROSSREFS

Sequence in context: A241582 A082003 A093107 * A067945 A166156 A089473

Adjacent sequences:  A125505 A125506 A125507 * A125509 A125510 A125511

KEYWORD

nonn

AUTHOR

Paul Richards, Jan 18 2007

STATUS

approved

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

License Agreements, Terms of Use, Privacy Policy .

Last modified May 24 21:25 EDT 2017. Contains 287008 sequences.