|
%I #55 Aug 24 2020 22:29:06
%S 1,3,7,15,31,127,1023,8191,34335,57855,131071,524287,2147483647
%N Numbers k for which A000265(1+A000265(sigma(k))) is equal to A000265(1+k).
%C Numbers k such that A336698(k) [= A000265(1+A161942(k))] is equal to A000265(1+k).
%C Numbers k such that A337194(k) = 2^e * A000265(1+k), for some e >= 1, where that e = A337195(k).
%C Any odd perfect number would trivially satisfy this condition.
%C Also, all hypothetical quasiperfect numbers, numbers k that satisfy sigma(k) = 2k+1, would be members.
%C Question: Is A066175 a subsequence of this sequence?
%C From _Antti Karttunen_, Aug 23 2020: (Start)
%C Numbers k such that (1+k) = 2^e * A336698(k), for some e >= 0.
%C Thus numbers k such that for some e >= 0, (1+k) = 2^(e-A337195(k)) * A337194(k), or equally, that A337194(k) = 2^(A337195(k)-e) * (1+k).
%C Conjecture: There are no even terms. This is equivalent to claim that there are no k such that A336698(k) = 1+k: If we assume that k is even, then in above equations we set e=0, and the requirement will then become that A337194(k) = 2^A337195(k)*(1+k), thus 1+k = A336698(k) = A000265(1+A000265(sigma(k))).
%C (End)
%H Paolo Cattaneo, <a href="http://www.bdim.eu/item?fmt=pdf&id=BUMI_1951_3_6_1_59_0">Sui numeri quasiperfetti</a>, Bollettino dell’Unione Matematica Italiana, Serie 3, Vol.6(1951), n.1, p. 59-62.
%H P. Hagis and G. L. Cohen, <a href="http://dx.doi.org/10.1017/S1446788700018401">Some Results Concerning Quasiperfect Numbers</a>, J. Austral. Math. Soc. Ser. A 33, 275-286, 1982.
%H V. Siva Rama Prasad and C. Sunitha, <a href="http://nntdm.net/papers/nntdm-23/NNTDM-23-3-073-078.pdf">On quasiperfect numbers</a>, Notes on Number Theory and Discrete Mathematics, Vol. 23, 2017, No. 3, 73-78.
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/QuasiperfectNumber.html">Quasiperfect Number</a>
%H <a href="/index/O#opnseqs">Index entries for sequences where odd perfect numbers must occur, if they exist at all</a>
%H <a href="/index/Si#SIGMAN">Index entries for sequences related to sigma(n)</a>
%t Block[{f}, f[n_] := n/2^IntegerExponent[n, 2]; Select[Range[2^20], f[1 + f[DivisorSigma[1, #]]] == f[1 + #] &] ] (* _Michael De Vlieger_, Aug 22 2020 *)
%o (PARI)
%o A000265(n) = (n>>valuation(n,2));
%o isA336701(n) = (A000265(1+A000265(sigma(n))) == A000265(1+n));
%Y Subsequence of A336700.
%Y Cf. A000203, A000265, A066175, A161942, A336698, A336699, A336937, A337194, A337195.
%Y Cf. A000668 (a subsequence).
%Y See also comments in A326042, A332223.
%K nonn,more
%O 1,2
%A _Antti Karttunen_, Aug 02 2020
| 