%I #34 Dec 18 2024 09:21:12
%S 1,3,7,25,33,55,57,69,91,93,2211,4825,12639,28225,32043,68727,89575,
%T 970225,2245557,16322559,22799825,48980427,55037217,60406599,68258725,
%U 325422273,414690595,569173299,794579511,10056372275,10475647197,10759889913,11154517557
%N Odd numbers k for which A003961(k)-2k divides A003961(k)-sigma(k), where A003961 shifts the prime factorization one step towards larger primes, and sigma is the sum of divisors function.
%C Numbers k for which A326057(k) = gcd(A003961(k)-2k, A003961(k)-sigma(k)) is equal to abs(A252748(k)) = |A003961(k)-2k|.
%C The odd terms of A326134 form a subsequence of this sequence. Unlike in A326134, here we don't constrain the value of A252748(k) = A003961(k)-2k, thus allowing also values <= +1. Because of that, the odd terms of A048674 and A348514 are all included here, for example 57 and 68727 that occur in A348514, and 1, 3, 25, 33, 93, 970225, 325422273, 414690595 that occur in A048674.
%C Conjecture (1): This is a subsequence of A319630, in other words, for all terms k, gcd(k, A003961(k)) = 1.
%C Conjecture (2): Apart from 1, there are no common terms with A349169, which would imply that no odd perfect numbers exist.
%C None of the 36 initial terms is Zumkeller, in A083207, because all are deficient (in A005100). See also A337372. - _Antti Karttunen_, Dec 05 2024
%H Amiram Eldar, <a href="/A349753/b349753.txt">Table of n, a(n) for n = 1..36</a> (terms below 10^11)
%H <a href="/index/O#opnseqs">Index entries for sequences where any odd perfect numbers must occur</a>
%H <a href="/index/Pri#prime_indices">Index entries for sequences computed from indices in prime factorization</a>
%H <a href="/index/Si#SIGMAN">Index entries for sequences related to sigma(n)</a>
%t f[p_, e_] := NextPrime[p]^e; s[1] = 1; s[n_] := Times @@ f @@@ FactorInteger[n]; q[n_] := Divisible[(sn = s[n]) - DivisorSigma[1, n], sn - 2*n]; Select[Range[1, 10^6, 2], q] (* _Amiram Eldar_, Dec 04 2021 *)
%o (PARI)
%o A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
%o isA349753(n) = if(!(n%2), 0, my(s = A003961(n), t = (s-(2*n)), u = s-sigma(n)); !(u%t));
%Y Cf. A000203, A003961, A005100, A083207, A252748, A286385, A319630, A326057, A326134, A337372, A349169.
%Y Cf. A048674, A348514.
%Y Subsequence of A378980 (its odd terms).
%K nonn
%O 1,2
%A _Antti Karttunen_, Dec 01 2021