OFFSET
1,2
COMMENTS
From Antti Karttunen, Jul 23 2022: (Start)
Numbers k such that k is a multiple of A350073(k).
For any square s in this sequence, A349162(s) = 1, i.e. sigma(s) divides A003961(s), and also A286385(s). Question: Is 1 the only square in this sequence? (see the conjecture in A350072).
If both x and y are terms and gcd(x, y) = 1, then x*y is also present.
After 2, the only primes present are Mersenne primes, A000668.
(End)
LINKS
MATHEMATICA
f[p_, e_] := NextPrime[p]^e; s[1] = 1; s[n_] := Times @@ f @@@ FactorInteger[n]; oddpart[n_] := n/2^IntegerExponent[n, 2]; q[n_] := oddpart[(sigma = DivisorSigma[1, n])] == GCD[sigma, s[n]]; Select[Range[1500], q] (* Amiram Eldar, Dec 04 2021 *)
PROG
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Dec 03 2021
STATUS
approved