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

A336317
Numbers k such that A122111(k) [conjugated prime factorization of k] is one of Ore's Harmonic numbers (in A001599).
4
1, 6, 40, 126, 154, 204, 1716, 1914, 2772, 8580, 11264, 12090, 12540, 50960, 62790, 64350, 77748, 83200, 104720, 152320, 186116, 193440, 331890, 382720, 432768, 518364, 648788, 684684, 753480, 817344, 895356, 1083852, 1113840, 1619352, 1675044, 1743588, 1759680, 1991340, 2060322, 2360484, 2492028, 2621080, 2932800
OFFSET
1,2
COMMENTS
Numbers k for which A336314(k) = A323173(k).
Sequence A122111(A001599(n)), n >= 1, sorted into ascending order. Positions of zeros in A323174 (corresponding to perfect numbers similarly mapped) is a subsequence.
Note that all terms after 1 seem to be present in A102750. This observation is equal to Ore's conjecture that there are no odd Harmonic numbers larger than one.
Also, all terms after 1 seem to be even, which would imply that apart from its initial 1, A001599 were a subsequence of A102750. However, this is false, as there are terms of A001599 not in A102750, for example 8011798098793361832960 found by David A. Corneth. Note that A122111(8011798098793361832960) = 96922193555635754403846044921625, which is thus an odd term of this sequence.
PROG
(PARI)
isA001599(n) = !((sigma(n, 0)*n)%sigma(n, 1));
isA336317(n) = isA001599(A122111(n)); \\ Program for A122111 given under that entry.
(PARI)
\\ Standalone program:
isA336317(n) = if(1==n, 1, my(f=factor(n), es=Vecrev(f[, 2]), is=concat(apply(primepi, Vecrev(f[, 1])), [0]), pri=0, d=1, s=1, x=1, p, e); for(i=1, #es, pri += es[i]; p = prime(pri); e = 1+is[i]-is[1+i]; d *= e; s *= ((p^e)-1)/(p-1); x *= (p^(e-1))); !((x*d)%s));
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jul 19 2020
STATUS
approved