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”).

A364962
Odd numbers k such that A005941(k) is either k itself or its descendant in Doudna-tree, A005940.
2
1, 3, 5, 11, 19, 23, 31, 37, 41, 43, 47, 53, 59, 61, 71, 73, 79, 83, 85, 89, 97, 101, 103, 107, 109, 113, 127
OFFSET
1,2
COMMENTS
Questions: Is 85 the only composite in this sequence? (See also A364565). Are there any more terms after 127, or is the sequence finite?
EXAMPLE
85 = 5*17 is a term, because A005941(85) = 133 = 7*19 = A003961(85), thus 133 is a left hand side child of 85 in the tree depicted in A005940, and therefore 85 is included in this sequence. (See also the last example in A364959).
PROG
(PARI)
A005941(n) = { my(f=factor(n), p, p2=1, res=0); for(i=1, #f~, p = 1 << (primepi(f[i, 1])-1); res += (p * p2 * (2^(f[i, 2])-1)); p2 <<= f[i, 2]); (1+res) }; \\ (After David A. Corneth's program for A156552)
A064989(n) = {my(f); f = factor(n); if((n>1 && f[1, 1]==2), f[1, 2] = 0); for (i=1, #f~, f[i, 1] = precprime(f[i, 1]-1)); factorback(f)};
A252463(n) = if(!(n%2), n/2, A064989(n));
isA364962(n) = if(!(n%2), 0, my(k=A005941(n)); while(k>n, k = A252463(k)); (k==n));
CROSSREFS
KEYWORD
nonn,hard,more
AUTHOR
Antti Karttunen, Aug 14 2023
STATUS
approved