|
|
A349746
|
|
Numbers k for which k * gcd(sigma(k), u) is equal to sigma(k) * gcd(k, u), where u is obtained by shifting the prime factorization of k two steps toward larger primes [with u = A003961(A003961(k))].
|
|
2
|
|
|
1, 11466, 114660, 411264, 804384, 871416, 4999680, 46332000, 176417280, 378069120, 396168192, 485188704, 709430400, 2004912000, 3921372000, 5600534400, 6128179200, 6956471808, 7556976000, 7746979968, 9904204800, 14092001280, 14182439040, 23423662080, 31998395520
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
Sigma preserves both the 2-adic and 3-adic valuation of the terms of this sequence.
All 65 known 5-multiperfect numbers (A046060) are included in this sequence, as well as the smallest 7-multiperfect number, 141310897947438348259849402738485523264343544818565120000 = A007539(7), and probably the majority of other p-multiperfect numbers as well, where p is a prime > 3. However, any term that is in A349747 is not included in this sequence.
|
|
LINKS
|
|
|
FORMULA
|
|
|
MATHEMATICA
|
f[p_, e_] := NextPrime[p, 2]^e; s[1] = 1; s[n_] := Times @@ f @@@ FactorInteger[n]; q[n_] := n * GCD[(sigma = DivisorSigma[1, n]), (u = s[n])] == sigma * GCD[n, u]; Select[Range[10^6], q] (* Amiram Eldar, Dec 01 2021 *)
|
|
PROG
|
(PARI)
A003961twice(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(1+nextprime(1+f[i, 1]))); factorback(f); };
isA349746(n) = { my(s=sigma(n), u=A003961twice(n)); (n*gcd(s, u) == (s*gcd(n, u))); };
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|