The number associated with my diurnal age today (28262) is what might be termed a "super sphenic number" as I'll explain in this post. Firstly however, its factorisation:$$28262=2 \times 13 \times 1087$$If we reverse its digits, we get the number 26282 and this number is also sphenic:$$26282 = 2 \times 17 \times 773$$Let's now concatenate the factors of 28262 in ascending order. This gives us the number 2131087 which is also sphenic:$$2131087 = 7 \times 167 \times 1823$$28262 has a sum of digits of 20 and if we add this to the original number we get the palindromic number 28282 which is also sphenic:$$28262+20=28282 = 2 \times 79 \times 179$$The number has a product of digits of 384 and if we subtract this from the original number we get 27878 which is sphenic:$$28262 - 384 = 27878 = 2 \times 53 \times 263$$If we consider only the internal digits of 28262, we get the number 826 which is also sphenic:$$2\, 826 \, 2 \rightarrow826=2 \times 7 \times 59$$When a sphenic number is considered as a sphenic brick then it has an associated number in the form of the brick's surface area. In the case of 28262, this associated surface area is 32662 square units and this number too is sphenic:$$32662 = 2 \times 7 \times 2333$$28262 has a sum of proper divisors that is also sphenic:$$ \sigma(28262) - 28262 =17434 = 2 \times 23 \times 379 $$The number has a sum of prime factors (1102) that is sphenic:$$2 +13+1087=1102 = 2 \times 19 \times 29$$28262 has a totient of 13032 which is not sphenic but its cototient (number - totient) of 15230 is:$$28262 - 13032 = 15230 = 2 \times 5 \times 1523$$So we can see that 28262 may well be termed a super sphenic number because of the above associations.
Wednesday, 19 August 2026
Monday, 13 April 2026
Divisors and Antidivisors: A Fresh Perspective
I've written about antidivisors before in posts titled Anti-Divisors and More on Anti-Divisors. Here is a fresh take. My diurnal age today is 28134 and the antidivisors of this number are 4, 12, 36, 108, 2084, 6252 and 18756. If we concatenate these numbers from left to right we get a rather large number (412361082084625218756) with the property that:$$ \frac{412361082084625218756}{28134} = 14657037111133334$$Numbers with this property form OEIS A249764:
Up to 40000, these numbers are relatively rare (permalink):
15, 30, 105, 120, 150, 222, 375, 585, 1500, 1695, 1755, 1800, 2700, 3449, 3750, 3840, 4891, 6720, 7680, 12000, 13583, 14400, 15000, 18750, 19200, 20940, 28134, 30000, 34800, 35625
The OEIS lists an analog of this, namely OEIS A069872:
Comparatively, members of this sequence are rather more numerous. There are 181 numbers in the range up to 40000 (permalink):
1, 2, 4, 5, 6, 8, 10, 15, 16, 20, 24, 25, 30, 32, 40, 50, 60, 64, 80, 90, 96, 100, 104, 120, 124, 125, 128, 150, 160, 200, 240, 250, 255, 256, 288, 320, 360, 375, 380, 384, 400, 425, 464, 480, 495, 500, 512, 600, 618, 625, 640, 750, 795, 800, 864, 875, 960, 1000, 1024, 1110, 1230, 1250, 1280, 1300, 1390, 1400, 1408, 1440, 1469, 1500, 1525, 1536, 1600, 1632, 1920, 2000, 2048, 2050, 2250, 2400, 2500, 2556, 2560, 2910, 2944, 2952, 3000, 3040, 3125, 3200, 3330, 3360, 3625, 3750, 3825, 3840, 4000, 4096, 4304, 4625, 4650, 4992, 5000, 5120, 5250, 5280, 5300, 5568, 5760, 6000, 6144, 6150, 6250, 6400, 6528, 6750, 7168, 7200, 7560, 7680, 8000, 8192, 9000, 9330, 9375, 9600, 10000, 10240, 10500, 10752, 11875, 12000, 12500, 12800, 13410, 13600, 14000, 15000, 15360, 15625, 15680, 16000, 16384, 16640, 17500, 17920, 18432, 18750, 19710, 20000, 20200, 20250, 20480, 20610, 21600, 21760, 22752, 23040, 23375, 24000, 24120, 24576, 25000, 25600, 25984, 27000, 28125, 30720, 31250, 32000, 32768, 33930, 34480, 34800, 35000, 36000, 37500, 37830, 38400, 39000, 40000
A variation on this is OEIS A240265:
These are less numerous. In the range up to 40000, they are:
4, 15, 16, 255, 375, 495, 795, 1469, 3825, 9375, 28125
All these numbers are also in OEIS A069872 and are marked in red in that sequence's list of members. It can be seen that 28125 is rather special and this number marked my diurnal age only 9 days ago. However, at the time, I missed its membership in these two sequences.
Monday, 11 August 2025
Prime Sum of Divisors
Surprisingly, for me at least, there are not many numbers whose sum of divisors is prime. However, the number associated with my diurnal age today, \( \textbf{27889}\), is one such number. The numbers up to 40000 are 2, 4, 9, 16, 25, 64, 289, 729, 1681, 2401, 3481, 4096, 5041, 7921, 10201, 15625, 17161, 27889, 28561 and 29929. These are the initial terms of OEIS A023194. Here is the breakdown (with 2 excluded - permalink):
number square root sigma4 2 79 3 1316 4 3125 5 3164 8 127289 17 307729 27 10931681 41 17232401 49 28013481 59 35414096 64 81915041 71 51137921 89 801110201 101 1030315625 125 1953117161 131 1729327889 167 2805728561 169 3094129929 173 30103
The members of this sequence are all square numbers with the exception of 2. Some are powers of 2 as well such as 16, 64 and 4096. The OEIS comments contain the following observations (that I'm not sure I understand):
All terms except the first are squares. Why?
Answer from Gabe Cunningham (gcasey@mit.edu): "From the fact that the sigma (the sum-of-divisors function) is multiplicative, we can derive that the \( \sigma(n)\) is even except when \(n\) is a square or twice a square.
"If \(n = 2 \times (2 \times k + 1)^2\), that is, \(n\) is twice an odd square, then:$$ \sigma(n) = 3 \times \sigma((2 \times k + 1)^2)$$If \(n = 2 \times (2 \times k)^2 \), that is, \(n\) is twice an even square, then \( \sigma(n)\) is only prime if \(n\) is a power of 2; otherwise we have:$$ \sigma(n) = \sigma(8 \times 2^m) \times \sigma(k/2^m)$$ for some positive integer \(m\).
"So the only possible candidates for values of n other than squares such that \( \sigma(n)\) is prime are odd powers of 2. But:$$ \begin{align} \sigma(2^{(2 \times m + 1)}) &= 2^{(2 \times m + 2)} - 1 \\ &= (2^{m + 1} + 1) \times (2^{m + 1 }- 1) \end{align}$$which is only prime when \(m = 0\), that is, when \(n = 2\). So 2 is the only nonsquare \(n\) such that \( \sigma(n)\) is prime."
I looked to Gemini to explain what the above means by asking the following question;
What does it mean to say that "the sigma function is multiplicative" where sigma is the sum of divisors function? The fuller quote is as follows:
... I just included the quote above ...
The quote is taken from the comments to OEIS A023194 (https://oeis.org/A023194). Can you break this down for me as I simply don't understand what the author is on about?
Here was the response (link):
Of course. Let's break down this number theory argument step by step. It's a great example of how a single property of a function can lead to a powerful conclusion.
What does "Multiplicative" Mean?
First, let's tackle the main concept. In number theory, a function
f(n)is called multiplicative if, for any two coprime integersaandb, the following equation holds true:
f(a * b) = f(a) * f(b)Coprime means that the two numbers share no common factors other than 1. For example, 4 and 9 are coprime, but 6 and 9 are not (they share a common factor of 3).
The sigma function, written as σ(n), calculates the sum of all the positive divisors of a number
n.Let's test the multiplicative property with an example, n = 12.
We can split 12 into the coprime factors 3 and 4.
σ(12): The divisors of 12 are 1, 2, 3, 4, 6, and 12. Their sum is
1 + 2 + 3 + 4 + 6 + 12 = 28.σ(3): The divisors of 3 are 1 and 3. Their sum is
1 + 3 = 4.σ(4): The divisors of 4 are 1, 2, and 4. Their sum is
1 + 2 + 4 = 7.Now let's check if the property holds:
σ(3) * σ(4) = 4 * 7 = 28
It works!
σ(12) = σ(3) * σ(4). This property is incredibly useful because it allows us to find σ(n) by first finding the prime factorization ofn. Ifn = p₁^a₁ * p₂^a₂ * ..., thenσ(n) = σ(p₁^a₁) * σ(p₂^a₂) * ....Deconstructing the Argument
Now, let's walk through the author's logic for finding the non-square numbers
nwhereσ(n)is prime. The goal is to eliminate all possibilities until only one remains.Part 1: When is σ(n) an odd number?
The author starts with the fact that
σ(n)is almost always even. Forσ(n)to be a prime number (other than 2), it must be an odd prime. So, when isσ(n)odd?
σ(n)is odd only if all the factors in its prime decomposition,σ(p₁^a₁), σ(p₂^a₂), ..., are also odd.
For powers of 2 (p=2):
σ(2^a) = 1 + 2 + 4 + ... + 2^a = 2^(a+1) - 1. This sum is always odd. (e.g., σ(4) = 7, σ(8) = 15).6 For powers of an odd prime (p is odd):
σ(p^a) = 1 + p + p² + ... + p^a. This sum of numbers is odd only if there is an odd number of terms. The number of terms isa+1. Fora+1to be odd,amust be even.So, for
σ(n)to be odd, any odd prime factor innmust have an even exponent. This meansnmust be of the form:
n = 2^k * (an odd number that is a perfect square)This leads to two cases:
k is even:
nis a perfect square (e.g.,2² * 3² = 6²).k is odd:
nis twice a perfect square (e.g.,2³ * 5² = 2 * (2*5)² = 2 * 10²).The author states:
σ(n)is odd only whennis a square or twice a square. This is the conclusion. Since we are looking for non-square values ofn, this immediately tells us we only need to check numbers that are twice a square.Part 2: Analyzing
nas "Twice a Square"We've narrowed our search to
n = 2k². The author now splits this into two sub-cases based on whetherkis odd or even.
Case A: n is twice an odd square. (n = 2 * m², where m is odd).
Since 2 and m² are coprime, we can use the multiplicative property:
σ(n) = σ(2 * m²) = σ(2) * σ(m²) = 3 * σ(m²)
For this result (3 * σ(m²)) to be a prime number, it must be equal to 3. This is only possible if σ(m²) = 1. The only number whose sum of divisors is 1 is the number 1 itself. So m² = 1, which means m = 1.
This gives us our first candidate: n = 2 * 1² = 2.
Let's check: σ(2) = 1 + 2 = 3, which is a prime number. And n=2 is not a square. So, 2 is a solution.
Case B: n is twice an even square. (n = 2 * m², where m is even).
Let the even number m = 2^j * s, where s is an odd number and j ≥ 1.
Then n = 2 * (2^j * s)² = 2 * 2^(2j) * s² = 2^(2j+1) * s².
This n is composed of an odd power of 2 (2^(2j+1)) and an odd square (s²). These two parts are coprime.
σ(n) = σ(2^(2j+1)) * σ(s²)
As we established, both σ(power of 2) and σ(odd square) are odd numbers greater than 1 (assuming s > 1). When you multiply two integers that are both greater than 1, the result is a composite (non-prime) number.
Therefore, for σ(n) to have any chance of being prime, the s² part must not exist. This means s must be 1.
Part 3: The Final Elimination
The previous step forces us to conclude that the only possible candidates left are when
s=1. This meansnmust be an odd power of 2.
n = 2^(2m+1)form ≥ 0.Let's find the sigma of this number:
σ(n) = σ(2^(2m+1)) = 2^((2m+1)+1) - 1 = 2^(2m+2) - 1
The author then uses a standard algebraic factorization for a difference of squares (a² - b² = (a-b)(a+b)):
2^(2m+2) - 1 = (2^(m+1))² - 1² = (2^(m+1) - 1) * (2^(m+1) + 1)
For this product to be a prime number, one of its factors must be equal to 1.
The second factor,
(2^(m+1) + 1), is clearly larger than 1.Therefore, the first factor must be 1:
2^(m+1) - 1 = 1
2^(m+1) = 2
m+1 = 1
m = 0
Plugging m=0 back into our formula for n:
n = 2^(2*0 + 1) = 2¹ = 2.
This confirms that the only possible candidate that survived this entire process of elimination is the number 2.
Summary
The author brilliantly uses the multiplicative nature of
σ(n)to prove that 2 is the only non-square number whose sum of divisors is a prime number.
I'm not sure I'm any the wiser but we definitely all such numbers are squares!
Anyway, by contrast, numbers for which the sum of the \( \textbf{proper}\) divisors is prime are far more common. There are 3771 of them in the range up to 40000.