A324316 -id:A324316 - cp's OEIS frontend

A324975 Rank of the n-th Carmichael number.

6, 10, 12, 8, 8, 10, 6, 6, 8, 18, 52, 12, 12, 18, 98, 164, 22, 6, 50, 8, 96, 34, 52, 46, 52, 6, 6, 156, 20, 46, 36, 32, 16, 8, 304, 36, 20, 36, 10, 316, 76, 468, 8, 30, 24, 1580, 84, 54, 8, 12, 250, 28, 92, 36, 20, 418, 456, 928, 188, 16, 8, 276, 284, 56, 144

Offset: 1

Bernd C. Kellner and Jonathan Sondow, Mar 24 2019

nonn

See A324974 for definition and explanation of rank of a special polygonal number, hence of rank of a Carmichael number A002997 by Kellner and Sondow 2019.

The ranks of the primary Carmichael numbers A324316 form the subsequence A324976.

			If m = A002997(1) = 561 = 3*11*17, then p = 17, so a(1) = 2+2*((561/17)-1)/(17-1) = 6.

Amiram Eldar, Table of n, a(n) for n = 1..10000
Bernd C. Kellner and Jonathan Sondow, On Carmichael and polygonal numbers, Bernoulli polynomials, and sums of base-p digits, Integers 21 (2021), #A52, 21 pp.; arXiv:1902.10672 [math.NT], 2019.
Bernd C. Kellner, On primary Carmichael numbers, Integers 22 (2022), #A38, 39 pp.; arXiv:1902.11283 [math.NT], 2019.
Wikipedia, Polygonal number

Subsequence of A324974.
A324976 is a subsequence.
Cf. also A002997, A324316, A324972, A324973, A324977.

T = Cases[Range[1, 10000000, 2], n_ /; Mod[n, CarmichaelLambda[n]] == 1 && ! PrimeQ[n]];
GPF[n_] := Last[Select[Divisors[n], PrimeQ]];
Table[2 + 2*(T[[i]]/GPF[T[[i]]] - 1)/(GPF[T[[i]]] - 1), {i, Length[T]}]

a(n) = 2+2*((m/p)-1)/(p-1), where m = A002997(n) and p is its greatest prime factor. (See Formula in A324974.) Hence a(n) is even, by Carmichael's theorem that p-1 divides (m/p)-1, for any prime factor p of a Carmichael number m.

A324976 Rank of the n-th primary Carmichael number.

Original entry on oeis.org

12, 8, 18, 12, 52, 52, 20, 32, 16, 54, 8, 36, 124, 34, 12, 72, 96, 26, 28, 76, 98, 1804, 108, 124, 18, 72, 172, 120, 10, 104, 32, 244, 130, 376, 18, 92, 780, 36, 172, 92, 284, 24, 198, 12, 244, 64, 234, 340, 100, 284, 24, 124, 44, 518, 364, 16, 82, 148, 8, 206

Offset: 1

Bernd C. Kellner and Jonathan Sondow, Mar 24 2019

See A324974 for definition and explanation of rank of a special polygonal number, hence of rank of a primary Carmichael number A324316 by Kellner and Sondow 2019.

			If m = A324316(1) = 1729 = 7*13*19, then p = 19, so a(1) = 2+2*((1729/19)-1)/(19-1) = 12.

Amiram Eldar, Table of n, a(n) for n = 1..10000
Bernd C. Kellner and Jonathan Sondow, On Carmichael and polygonal numbers, Bernoulli polynomials, and sums of base-p digits, Integers 21 (2021), #A52, 21 pp.; arXiv:1902.10672 [math.NT], 2019.
Bernd C. Kellner, On primary Carmichael numbers, Integers 22 (2022), #A38, 39 pp.; arXiv:1902.11283 [math.NT], 2019.
Wikipedia, Polygonal number

Subsequence of A324975 (rank of the n-th Carmichael number A002997) and of A324974 (rank of the n-th special polygonal number A324973).

Cf. also A324316, A324972.

SD[n_, p_] := If[n < 1 || p < 2, 0, Plus @@ IntegerDigits[n, p]];
LP[n_] := Transpose[FactorInteger[n]][[1]];
TestCP[n_] := (n > 1) && SquareFreeQ[n] && VectorQ[LP[n], SD[n, #] == # &];
T = Select[Range[1, 10^7, 2], TestCP[#] &];
GPF[n_] := Last[Select[Divisors[n], PrimeQ]];
Table[2 + 2*(T[[i]]/GPF[T[[i]]] - 1)/(GPF[T[[i]]] - 1), {i, Length[T]}]

a(n) = 2+2*((m/p)-1)/(p-1), where m = A324316(n) and p is its greatest prime factor. Hence a(n) is even; see Formula in A324975.

More terms from Amiram Eldar, Mar 27 2019

A324857 Numbers m > 1 such that there exists a prime divisor p of m with s_p(m) = p.

Original entry on oeis.org

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

Offset: 1

Jonathan Sondow, Mar 17 2019

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

			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.

Robert Israel, Table of n, a(n) for n = 1..10000
Bernd C. Kellner, On primary Carmichael numbers, Integers 22 (2022), Article #A38, 39 pp.; arXiv:1902.11283 [math.NT], 2019.

A324456 is the union of A324857 and A324858.

Includes A083558.

S:= (p,m) -> convert(convert(m,base,p),`+`):
filter:= proc(m) ormap(p -> S(p,m) = p, numtheory:-factorset(m)) end proc:
select(filter, [$2..500]); # Robert Israel, Mar 20 2019

s[n_, b_] := If[n < 1 || b < 2, 0, Plus @@ IntegerDigits[n, b]];
f[n_] := AnyTrue[Divisors[n], PrimeQ[#] && s[n, #] == # &];
Select[Range[400], f[#] &]n (* simplified by Bernd C. Kellner, Apr 02 2019 *)

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

A324974 Rank of the n-th special polygonal number A324973(n).

Original entry on oeis.org

3, 3, 3, 5, 3, 3, 6, 3, 6, 3, 11, 5, 3, 3, 8, 10, 5, 6, 12, 3, 15, 9, 3, 5, 3, 8, 3, 8, 19, 14, 5, 7, 3, 6, 6, 36, 21, 66, 22, 3, 10, 5, 6, 3, 3, 50, 10, 20, 5, 14, 11, 51, 3, 10, 21, 6, 13, 5, 16, 25, 3, 3, 6, 6, 12, 14, 10, 68, 5, 28, 3, 11, 29, 3, 56, 6, 19

Offset: 1

Bernd C. Kellner and Jonathan Sondow, Mar 24 2019

nonn

While two polygonal numbers of different ranks can be equal (e.g., P(6,n) = P(3,2n-1)), that cannot occur for special polygonal numbers, since for fixed p the value of P(r,p) is strictly increasing with r. Thus the rank of a special polygonal number is well-defined.

The Carmichael numbers A002997 and primary Carmichael numbers A324316 are special polygonal numbers (see Kellner and Sondow 2019). Their ranks form the subsequences A324975 and A324976.

			If m = A324973(4) = 70 = 2*5*7, then p = 7, so a(4) = 2+2*((70/7)-1)/(7-1) = 5.

Amiram Eldar, Table of n, a(n) for n = 1..10000
Bernd C. Kellner and Jonathan Sondow, On Carmichael and polygonal numbers, Bernoulli polynomials, and sums of base-p digits, Integers 21 (2021), #A52, 21 pp.; arXiv preprint, arXiv:1902.10672 [math.NT], 2019-2021.
Bernd C. Kellner, On primary Carmichael numbers, Integers 22 (2022), #A38, 39 pp.; arXiv preprint, arXiv:1902.11283 [math.NT], 2019-2022.
Wikipedia, Polygonal number.

A324975 and A324976 are subsequences.

Cf. A002997, A324316, A324972, A324973, A324977.

GPF[n_] := Last[Select[Divisors[n], PrimeQ]];
T = Select[Flatten[Table[{p, (p^2*(r - 2) - p*(r - 4))/2}, {p, 3, 150}, {r, 3, 100}], 1], SquareFreeQ[Last[#]] && First[#] == GPF[Last[#]] &];
TT = Take[Union[Table[Last[T[[i]]], {i, Length[T]}]], 47];
Table[2 + 2*(t/GPF[t] - 1)/(GPF[t] - 1), {t, TT}]

a(n) = 2 + 2*((m/p)-1)/(p-1), where m = A324973(n) and p is its greatest prime factor. (Proof: solve m = P(r,p) = (p^2*(r-2) - p*(r-4))/2 for r.)

Several missing terms inserted by and more terms from Jinyuan Wang, Feb 18 2021

cp's OEIS Frontend

A324975 Rank of the n-th Carmichael number.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A324976 Rank of the n-th primary Carmichael number.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A324857 Numbers m > 1 such that there exists a prime divisor p of m with s_p(m) = p.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

A324974 Rank of the n-th special polygonal number A324973(n).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

Extensions