|
|
A324857
|
|
Numbers m > 1 such that there exists a prime divisor p of m with s_p(m) = p.
|
|
4
|
|
|
6, 10, 12, 15, 18, 20, 21, 24, 33, 34, 36, 39, 40, 45, 48, 57, 63, 65, 66, 68, 72, 80, 85, 87, 91, 93, 96, 99, 105, 111, 117, 130, 132, 133, 135, 136, 144, 145, 160, 165, 171, 175, 185, 189, 192, 205, 217, 225, 231, 249, 255, 258, 259, 260, 261, 264, 265, 272, 273, 279, 285, 288, 297, 301, 305, 320, 325, 327, 333, 341, 351, 384, 385
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
The function s_p(m) gives the sum of the base-p digits of m.
m must have at least 2 prime factors, since s_p(p^k) = 1 < p.
The sequence contains the primary Carmichael numbers A324316.
The main entry for this sequence is A324456 = numbers m > 1 such that there exists a divisor d > 1 of m with s_d(m) = d. It appears that d is usually prime: compare the sparser sequence A324858 = numbers m > 1 such that there exists a composite divisor c of m with s_c(m) = c. However, d is usually composite for higher values of m.
The sequence contains the 3-Carmichael numbers A087788, but not all Carmichael numbers A002997. This is a nontrivial fact. The smallest Carmichael number that is not a member is 173085121 = 11*31*53*61*157. For further properties of the terms see A324456 and Kellner 2019. - Bernd C. Kellner, Apr 02 2019
|
|
LINKS
|
|
|
EXAMPLE
|
s_p(m) = 1 < p for m = 2, 3, 4, 5 with prime p dividing m, but if m = 6 and p = 2 then s_p(m) = s_2(2 + 2^2) = 1 + 1 = 2 = p, so a(1) = 6.
|
|
MAPLE
|
S:= (p, m) -> convert(convert(m, base, p), `+`):
filter:= proc(m) ormap(p -> S(p, m) = p, numtheory:-factorset(m)) end proc:
|
|
MATHEMATICA
|
s[n_, b_] := If[n < 1 || b < 2, 0, Plus @@ IntegerDigits[n, b]];
f[n_] := AnyTrue[Divisors[n], PrimeQ[#] && s[n, #] == # &];
|
|
PROG
|
(PARI) isok(n) = {if (n>1, my(vp=factor(n)[, 1]); for (k=1, #vp, if (sumdigits(n, vp[k]) == vp[k], return (1)))); } \\ Michel Marcus, Mar 19 2019
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,base
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|