%I#7 2018年1月22日03:07:29
%S 1,2,3,4,5,6,7,8,9,10,11,12,13,15,16,17,18,19,20,22,23,24,25,27,29,30,
%电话31,32,33,36,37,40,41,43,44,45,47,48,49,50,53,54,55,59,60,61,62,64,66,
%U 67,71,72,73,75,79,80,81,83,88,89,90,91,93,96,97,99100单位
%N Matula-Goebel有根树的数量,这样根的每个分支都有相同数量的叶子。
%e树木顺序开始:
%e 1 o
%e 2(o)
%e 3((o))
%e 4(oo)
%e 5((o))
%e 6(o(o))
%e 7((oo))
%e 8(ooo)
%e 9((o)(o))
%e 10(o(o))
%e 11(((o)))
%e 12(oo(o))
%e 13((o(o)))
%e 15((o)(o))
%e 16(oooo)
%e 17(((oo))
%e 18(o(o)(o))
%e 19((ooo))
%e 20(oo((o)))
%e 22(o((o)))
%e 23((o)(o))
%e 24(ooo(o))
%e 25((o)(o))
%e 27((o)(o)(o))
%e 29((o(o)))
%e 30(o(o)(o))
%t nn=2000;
%t素数MS[n_]:=如果[n===1,{},平坦[Cases[FactorInteger[n],{p_,k_}:>表[PrimePi[p],{k}]]];
%t叶数[n_]:=如果[n===1,1,其中[{m=primeMS[n]},如果[Length[m]===1,叶数[First[m]],总[leafcount/@m]]];
%t选择[Range[nn],SameQ@@leafcount/@primeMS[#]&]
%Y参见A000081、A007097、A061775、A111299、A214577、A276625、A290760、A290822、A291442、A298533、A29853。
%K非n
%O 1,2号机组
%A _Gus Wiseman_,2018年1月20日