[go: up one dir, main page]

login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

a(n) = 1 if A005940(1+n) is deficient, 0 otherwise.
10

%I #51 Aug 03 2022 23:25:18

%S 1,1,1,1,1,0,1,1,1,1,1,0,1,0,1,1,1,1,1,0,1,0,1,0,1,1,1,0,1,0,1,1,1,1,

%T 1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,1,1,0,1,0,1,0,1,1,1,0,1,0,1,1,1,1,1,1,

%U 1,0,1,0,1,1,1,0,1,0,1,0,1,1,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,1,1,0,1,0,1,0,1,0

%N a(n) = 1 if A005940(1+n) is deficient, 0 otherwise.

%C Any node of Doudna-tree (A005940) which is deficient, is marked here with 1, and has also its all leftward descendants (those that are obtained only by prime shifting, A003961) deficient.

%C Any odd node marked with 0 here (that is, nondeficient nodes in A005940) has all its ancestors up at least to the first even ancestor also nondeficient. None is the seen in the illustration below, because the first case is 945, which is the first odd abundant number (A005231), whose even parent is A064989(945) = 120, also an abundant number. (But please see also _Michael De Vlieger_'s illustration of the tree down to level 8 where it is included in the bottom row).

%C Specifically, any odd perfect number, if they exist, should have an abundant mother (which should be also a member of A336930), and also abundant grand-ancestors (on the prime-shift line), up to and including the first even ancestor, at least.

%C On the other hand, when going towards the right, any right child has larger abundancy index than its mother has, and on all rightward leaning edges, except on the rightmost (powers of 2), a nondeficient number is eventually reached. See A336915.

%H Antti Karttunen, <a href="/A336834/b336834.txt">Table of n, a(n) for n = 0..65535</a>

%H Michael De Vlieger, <a href="/A336834/a336834.png">Tree illustrated down to level 7</a> (red circles stand for 1's, white circles for 0's, and the numbers inside circles are those in A005940)

%H Michael De Vlieger, <a href="/A336834/a336834_2.png">Tree illustrated down to level 8</a> (otherwise like in above. Note how the first three odd abundant numbers 945, 1575, 2205 can be seen in the bottom row, inscribed in white circles, being left children of other white circles)

%H Michael De Vlieger, <a href="/A336834/a336834_1.png">Chart of levels 0 <= j <= 16</a>, vertically exaggerated 120X.

%H Michael De Vlieger, <a href="/A336834/a336834_3.png">4096-pixel square raster showing 2^24 terms</a>, where black = 1 and white = 0.

%H <a href="/index/Ch#char_fns">Index entries for characteristic functions</a>

%H Michael De Vlieger, <a href="/A336834/a336834_4.png">Fan-style tree</a> showing levels j = 0..13. Partial annotation of A005940(n) to levels j = 0..7.

%F a(n) = A294934(A005940(1+n)).

%e a(n) tells whether the n-th node in A005940 (here counted with offset 0) is deficient (here marked with 1), or not (marked with 0).

%e The first six levels of the binary tree (compare also to the illustration given at A005940):

%e 1

%e |

%e 1

%e ............../ \..............

%e 1 1

%e ....../ \...... ....../ \......

%e 1 0 1 1

%e / \ / \ / \ / \

%e / \ / \ / \ / \

%e 1 1 1 0 1 0 1 1

%e / \ / \ / \ / \ / \ / \ / \ / \

%e 1 1 1 0 1 0 1 0 1 1 1 0 1 0 1 1

%t Boole[DivisorSigma[1, #] < 2 #] & /@ Nest[Append[#1, Prime[1 + BitLength[#2] - DigitCount[#2, 2, 1]]*#1[[#2 - 2^Floor@ Log2@ #2 + 1]]] & @@ {#, Length@ #} &, {1}, 105]

%t (* Second program: decode chart of levels 0 <= j <= 16 to render 2^16 + 1 terms of a(n) in binary tree form *)

%t Block[{img = Import["https://oeis.org/A336834/a336834_1.png"], s, r, w, d}, s = ImageData[img]; r = (-1 + FirstPosition[s, 1.][[1]])/3; Set[{w, d}, {#1, -1 + #2/r}] & @@ ImageDimensions[img]; Array[Function[{i, j, k}, Array[s[[k, # i + 1]] /. {0. -> 1, 1.0 -> 0} &, j, 0]] @@ {w/#1, #1, #2} & @@ {If[# == 0, 1, 2^(# - 1)], r # + 1} &, d + 1, 0] ]

%t (* Third program: decode raster of 2^24 terms *)

%t ImageData[Import["https://oeis.org/A336834/a336834_3.png"]] /. {0. -> 1, 1. -> 0} // Flatten

%t (* _Michael De Vlieger_, Aug 18 2020 *)

%o (PARI)

%o A005940(n) = { my(p=2, t=1); n--; until(!n\=2, if((n%2), (t*=p), p=nextprime(p+1))); (t) };

%o A294934(n) = ((n+n) > sigma(n));

%o A336834(n) = A294934(A005940(1+n));

%Y Cf. A003961, A005231, A005940, A064989, A294934, A336389, A336835, A336915, A336917, A336930, A337377, A337386.

%Y Cf. A252743, A342000 for similarly constructed sequences, and also arrays A341605/A341606 for a more promising development of the ideas presented here.

%K nonn,tabf

%O 0

%A _Antti Karttunen_ and _Michael De Vlieger_, Aug 07 2020