The quotient triangle
Consider a positive integer . Then for all do the floor operation . Say , we get , a sequence of quotients of the division . If we do this for all we get the quotient triangular array whose top few elements are show below.
It is obvious that the first column is the sequence of itself. The second column is the sequence of positive integers, each repeated twice; the third column is the sequence of positive integers, each repeated thrice; so on.
It we linearize this triangular array we get the sequence . Since each row of the triangular array adds elements to the sequence, it grows as the sum of numbers from 1:n. We see that the successive maxima are attained at the following terms of : 1, 2, 4, 7, 11, 16, 22, 29, 37, 46, 56… We can plot by placing the first term at 0, i.e., its x-coordinate will be 0, while its y-coordinate will be ; for and so on. Based on the sum of integers 1:n we can show that the successive maxima attained by the sequence would be bounded in this plot by the parabola (Figure 1):
This sequence also features cycles of decaying values (Figure 1) between successive maxima, starting from one maximum and going to 1 before jumping back to the next maximum and decaying again. Since the quotients are successively determined by , each cycle of the sequence of length is a discrete form of the rectangular hyperbola (Figure 2):
Thus, we may term sequence the parabola-hyperbola sequence. Since any will produce quotient 1, the fraction of 1s in will converge to ; similarly, with produce quotient 2; hence, the fraction of 2s in will converge to . Thus, generally .
We then define a summatory sequence such that each term is the sum of the nth row of the quotient triangular array :
This sequence goes thus: 1, 3, 5, 8, 10, 14, 16, 20, 23, 27… and is shown graphically in figure 3.
Each term of this sequence is the sum of the terms in a given hyperbolic cycle (Figure 2) of sequence (Figure 1). Thus, it can be approximated in continuous form by the area under the hyperbola corresponding to each cycle . Hence, we can create a continuous approximation for by the integral function:
This is shown by the red curve in Figure 3. We notice that this continuous approximation falls short of actual discrete sequence. This correction factor can be empirically determined using the method of least squares to be a linear term. Thus, with this correction the below function provides a good fit for the average behavior of :
The empirical determination of this correction term prompted us to try to determine it from first principles. The function is a continuous approximation of the area in one cycle of the sequence . But in reality our cycles are approximately the discrete harmonic series as opposed to a smooth hyperbola. From Euler’s work we know that the difference between the discrete form and the continuous integral converges to that mysterious number Euler’s constant for a unit hyperbola. To make use of this we tried some experiments with the hyperbola and realized that the best way to capture the area of the cycle was to use the symmetry of the hyperbola about the line (Figure 2). Here, we can count all in a cycle vertically till (see vertical line in Figure 2). The due to symmetry we can make the count again from the other end horizontally till . This way we would have completely covered all the discrete counts except that we would have counted the square twice. Hence, need to subtract from our sum. We can get the two symmetric discrete sums now by the addition of to the integral of the continuous form. Thus, we can write the approximate area for a cycle and hence the function approximating as:
Thus, we get our correction term to be . It gave us great pleasure to have figured this out from scratch without any study of mathematical literature in this regard. We noticed that this already good fit can be made even better by addition of a further constant term . We do not know if this is really a constant or is some subtle term currently beyond our reach. Thus, we may write the final function approximating summatory sequence as:
This function gives a mean difference of with for (Figure 3: light blue dashed line). While this captures the average behavior of , examination of the specific behavior of shows that it exhibits saltations that are greater than usual at certain values. To better understand this we created the difference sequence (Figure 4):
The sequence shows several interesting features which we consider in detail below:
1) The lowest value it ever attains is 2. A closer examination of the indices at which reveals that they correspond to being a prime in the parent sequence. Thus, there are as many minima in as there are primes. This can be explained thus: Since, a prime is completely divisible only by , these two values will generate quotients of 1 and respectively to add to the quotient sum. The remaining quotients will be the same as previous number as none of the other between 1 and will perfectly divide the prime. Thus, the quotient sum will minimally differ from that of the previous number by 2.
2) Jumps above the median value of 6 (for this range; Figure 4: violet line) have propensity to increase with increasing . Analysis of these jumps revealed that they have a significantly higher propensity to occur at . Figure 5 shows a box-plot of , which indicates that the median value of this sequence for is significantly higher than median value of the sequence for any other as well as the median value of the overall sequence. This can be explained by considering the following corresponds to in the parent sequence. A number of the form will undergo complete divisions by at least 1, 2, 3, 6 out of the first 6 . However, the number before it will not be divisible by 2, 3 or their multiples. Since it cannot undergo complete division by them, its quotients will be less than the corresponding ones of by 1. Thus there will be a jump in corresponding to in . As e.g. consider 12:
Thus, causes a spike in due to completion of divisions by 2,3,4,6.
3) Notably, the successive record values of for are attained at: 1, 3, 5, 11, 23, 35, 47, 59, 119, 179, 239, 359, 719, 839, 1259, 1679, 2519, 5039, 7559, 10079, 15119, 20159, 25199, 27719 (Figure 4). These values of correspond to 2, 4, 6, 12, 24, 36, 48, 60, 120, 180, 240, 360, 720, 840, 1260, 1680, 2520, 5040, 7560, 10080, 15120, 20160, 25200, 27720 in . Remarkably, these latter numbers are the highly composite numbers first defined by Ramanujan in 1915 CE. They are numbers with more divisors than any integer smaller than them. We can see why these highly composite numbers spawn maxima in — they will be maximally completely divided and thus yield the largest quotient sum for any number up to them. Thus, presents them as antipodes to the primes. Further, Ramanujan defined a second more inclusive set of numbers, the largely composite numbers, which are numbers with more or the same number of divisors than any integer smaller than them. This more inclusive set spawns high-points in , which are not maxima but still stand out (Figure 4).
As an aside one may note that the above at which the maxima occur are often primes: 15 out of the 24. Most of these are also in particular Sophie Germain or safe primes (primes of form where if is a prime and is also a prime. Then the former is Sophie Germain prime and the latter a safe prime).
4) The successive record values of are: 2, 3, 4, 6, 8, 9, 10, 12, 16, 18, 20, 24, 30, 32, 36, 40, 48, 60, 64, 72, 80, 84, 90, 96… It turns out this sequence (another recorded by Ramanujan) is the number of divisors of the corresponding highly composite numbers which spawn the maximum. Thus, the successive maximum values of are attained at one less than a highly composite number and is equal to the number of divisors of the HCN.
5) From the the above discussion of the specific cases of the maxima and minima it should be apparent that:
where is the number of divisors including 1 and itself for .
The remainder triangle
The above triangular array can be converted to a second triangular array by the below transformation:
This turns out to be the triangular array of remainders. This can be alternatively obtained thus: Consider a positive integer . Then for all , . Say , we get , a sequence of remainders of the division . If we do this for all we get the remainder triangular array , which is the same as the array obtained by the above-stated transformation of . Shown below are its initial terms.
We observe that each column in can be described thus: 1st, all 0s; 2nd a 2-term cycle sequence of 0,1; 3rd a 3-term cycle of sequence 0,1,2. In general terms the kth column is a term cycle of form
If we linearize we get the remainders sequence , which goes 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 2, 1, 0, 0, 0, 0, 2, 1… As the linearization of a triangular array also grows as , where is the row number of . The sequence is plotted in Figure 6 with the x-coordinate of the first term as 0 and y=coordinate as .
It shows cycles of length . Each pair of cycles shows the same maximum value. Examination of the sequence shows that successively greater maxima are all bounded by the parabola:
Just as we did with , we can similarly define the summatory sequence (we are using this notation after we discovered post facto that has been used for related remainder sequences in the mathematical literature, where each term is the sum of the numbers in a row of . This sequence goes thus: 0, 0, 1, 1, 4, 3, 8, 8, 12, 13, 22, 17, 28, 31, 36, 36, 51, 47, 64, 61… It is plotted in Figure 7.
We determined that the sequence can be described by a continuous function, a parabola of the form (shown as red line in Figure 7), where we empirically determined .
We next sought to obtain the exact value of the constant of the parabola from first principles. For this consider the plot a single cycle of length of sequence (Figure 8). We observe that as it is a series of right triangles.
The largest of these triangles is an isosceles triangle with . The next triangle has , the next has . Thus, the general expression for the heights and bases of these triangles is where and the area of the triangles would be . Thus at limit the area under the triangles in a single cycle, which will be the value of the nth term of for large , will be given by the sum (assuming ):
We realized, much to our pleasure on discovering it, that this sum is convergent and remarkably can be expressed in terms of the Riemann function as:
Thus we can express the parabola fitting the sequence precisely as:
Thus, again we stumble upon one of those deep links we seen in mathematics, in this case between remainders of division and the number . Now, to study the exact behavior of we define the difference sequence :
This difference sequence is plotted in Figure 9.
The sequence shows the following notable features:
1) takes positive values more often than negative values (3.03:1) in the first 30000 terms. This corresponds to the ratio of numbers with the sum of their divisors smaller than themselves being less than them by at least 2 (known as deficient numbers) to those whose sum of divisors smaller than them is equal or greater than them. This will become clear from the discussion that follows further down.
2) Successive maxima, i.e. the record value of the sequence up to certain occur only at even values of . Closer examination revealed that the maxima are associated with being a prime in the parent sequence ; then will be a maximum. If we take a number then the numbers will be even. will be divisible by 3. Thus, these will not be primes. That leaves and ; the latter is the same as for another . Thus, will not be divisible by 2 and 3 and could be prime . Thus, all primes can be expressed as . Hence, in the sequence , will correspond and to . Thus, at , the successive maxima occur at or corresponding to every that related to a prime in the sequence. That successive maxima will occur at primes can be explained easily. A prime will not be divisible by any number other than 1 and itself so all other remainders will be non-zero. The number before it will have at least 2 or at least 2, 3 and 6 as divisors that will yield remainder 0. Thus, a prime will result in remainders adding to a larger number relative to the remainder sum of its prior number, thereby causing a maximum to attained. Thus, there will as many successive maxima in as primes.
3) only when , which corresponds to in .
4) Negative values of show a more complex behavior. They are seen where the in is a number which might be perfect numbers, abundant numbers and highly abundant numbers. We discuss these in detail below.
5) The successive minima, i.e. record negative values assumed by for a given are fewer than the successive maxima. While we see a succession of 3245 maxima, corresponding the same number of primes in the first 30000 terms of , we only see a succession of 60 minima. We did some experimentation to discover where these minima occur. Let us define an arithmetic function as below for all :
. Where is a divisor of from
Whenever is greater than of all integers lesser than it, then is considered part of the sequence. It goes: 2, 6, 12, 24, 36, 48, 60, 72, 84, 96… and may be called “remainder anti-primes”. When belongs to this sequence in it spawns a minimum at . Thus, for the successive minima occur only at a of the form . Remainder anti-primes are a subset of a class of numbers, highly abundant numbers, defined by the mathematician S. Sivasankaranarayana Pillai. These are numbers whose divisors smaller than themselves sum up to a value, which is bigger than the sum of the divisors for each of the numbers smaller than that number. More generally, whenever of is a highly abundant number it results in strongly negative values in .
6) Definition of also helps us understand the value assumed by :
Thus, for a prime since the divisors are only , we have . Thus successive maxima in occur at and attain the value . Thus, after 1 they are either of the form or and will always be odd. Further, this implies that the successive maxima are bounded by the line (Figure 9). The successive minima usually take the form . The middle of these, the only even values of minima, are very rare: there being on 3 in the first 30000 terms, namely . The successive minima increase in magnitude with increasing but can attain much greater magnitudes than the maxima in their vicinity. Their increase with is not linear unlike the maxima as the for this specialized subset of highly abundant numbers grows rapidly. We empirically found a curve of the form to approximately fit them (Figure 9).
7) When it corresponds to , i.e. in of 6, 28, 496, 8128 in our range of 1:30000. These numbers are known after the Pythagorean sage Nicomachus as perfect numbers because their divisors other than themselves add up exactly to them (one can see why this is so from the above, ):
One notices that they have a curious form with respect to the powers of 2:
More generally, if is a prime and is also a prime (such primes are called Mersenne primes) then is a perfect number. The search for such perfect numbers goes on this date computationally and checking online while writing this article shows that to date 50 such numbers have been found. All known perfect numbers are of a form ending in 6 or 28. Yet it has not be formally proven that no odd perfect numbers exist or if there are an infinite number of perfect numbers!
8) Abundant numbers are those numbers whose divisors smaller than themselves add up to a number greater than themselves. e.g. . When is an abundant number in it results in a negative value of . The odd abundant numbers like 945, 1575, 2205 in result in the only even indices ( ) at which the sequence is negative.
9) is unknown to us. This would imply that the sum of the divisors of a number which are smaller than it can never be 1 more than the number. Similarly , appears to be unknown, i.e. corresponding to deficient numbers with a deficiency of 3. We do not know if these have been proven or disproven. are also rare. We get only 3 such numbers from 1:30000 which correspond the below of and observe that:
These are numbers whose sum of divisors smaller than them fall short of them by 2. Similarly, we have only 5 numbers resulting in : 20, 104 464, 650, 1952. These are numbers whose sum of divisors smaller than them exceed them by 2.
The and sequences (Figures 4 and 9) are striking in being simple means of illustrating the fundamental asymmetry in the number world between primes and their anti-numbers. reveals the highly composite numbers, and more weakly the largely composite numbers of Ramanujan as the antipodes of primes. suggests that the remainder anti-primes, a subset of the highly abundant numbers of Pillai, are the anti-numbers of primes. From the vantage point of these sequences the primes appear to define a more regular pattern in terms of minima and maxima respectively while these anti-numbers define more mysterious patterns.
There are some questions for which we do not have ready answers are (though they might have been answered/studied by mathematicians):
1) Can you give an exact form for the constant correction term used to recapitulating ?
2) Can you derive a function that bounds the maxima of ?
3) Can you derive the exact form of the bounding function of the minima of ?
We carried out all this quite independently of the mathematical literature and were pleased to see the organic emergence of the two classes of anti-prime numbers including those previously discovered by Ramanujan. We showed an earlier variant of this material to a person with extraordinary mathematical talent. He brought to our attention a formal proof for the relationship between perfect numbers and the remainder sequence, which was given to him by a poorly known genius Suryanarayana from Vishakhapatnam who studied such issues extensively in the last century. However, such a proof might have been first published by the mathematician Lucas in the late 1800s. Searching the internet, we came across a simple presentation of that proof in an excellent paper by Spivey in the Mathematics Magazine in 2005. Hence, we are not providing that here. Otherwise our account here provides a summary of the main interesting results concerning these sequences we found as part of our investigation.