Pearl necklaces for Maheśvara
Śrīpati’s pearl necklace for Maheśvara
The brāhmaṇa Śrīpati of the Kāśyapa clan was a soothsayer from Rohiṇīkhaṇḍa, which is in the modern Buldhana district of Maharashtra state. Somewhere between 1030 to 1050 CE he composed several works on mathematics, astronomy and divination, one of which is the Gaṇita-tilaka on basic arithmetic and algebra which has come down to us through incomplete manuscripts. In that he posses the following problem (The solution is provided by his Jaina commentator Siṃhatilaka Sūri in 1275 CE):
The brāhmaṇa Śrīpati of the Kāśyapa clan was a soothsayer from Rohiṇīkhaṇḍa, which is in the modern Buldhana district of Maharashtra state. Somewhere between 1030 to 1050 CE he composed several works on mathematics, astronomy and divination, one of which is the Gaṇita-tilaka on basic arithmetic and algebra which has come down to us through incomplete manuscripts. In that he posses the following problem (The solution is provided by his Jaina commentator Siṃhatilaka Sūri in 1275 CE):
viśva-kha-sapta-bhujaṅga-navārkāḥ śaila-turaṅga-samāhata-dehāḥ |
syāt sphuṭa-tāraka-vartula-muktā-bhūṣaṇam atra maheśvara-kaṇṭhe ||
The bodies of all (13), the space (0), the seven (7), the snakes (8), the nine (9) and the suns (12) struck together (means multiplied in mathematical terminology) with the mountains (7) and the horses (7) may now be the clear, sparkling and globular pearls forming an ornament for the neck of Maheśvara.
syāt sphuṭa-tāraka-vartula-muktā-bhūṣaṇam atra maheśvara-kaṇṭhe ||
The bodies of all (13), the space (0), the seven (7), the snakes (8), the nine (9) and the suns (12) struck together (means multiplied in mathematical terminology) with the mountains (7) and the horses (7) may now be the clear, sparkling and globular pearls forming an ornament for the neck of Maheśvara.
The purpose of the problem is two-fold: (1) to make the student familiar with using the Hindu numerical code (the bhūta-saṃkhya) and (2) make the student do some elementary large number multiplication. Thus, the problem is actually a simple multiplication
whose answer is
. Interestingly, Śrīpati offers a clue for the answer in the problem itself: he says the multiplication of the 2 numbers results in a necklace of clear, sparkling, spherical pearls for Maheśvara. This is clearly an allusion to the palindromic structure of the product, with the 0s forming the pearls, since they were written historically as circles and the two flanking 1s form the bindings of the necklace. Some numerical savants are known to exhibit a synesthesia with respect numbers — perhaps such was indeed the situation with Śrīpati for a palindromic number such as this could simultaneously produce in ones mind the vivid image of a pearl necklace.
Moreover, this is not just any pearl necklace but one for Rudra’s neck. It is in this regard we believe he encoded more into that number. As we can see from the above, one of the factors of 1000000001 is 11, which is the characteristic number of the Rudra-s. At the same time, the Rudra-s are also said to be 1000s upon 1000s (Yajurveda: sahasrāṇi sahasraśo ye rudrā adhi bhūmyām |). Thus, this product captures both those aspects. Further, the numbers that yield the product are described in bhūta-saṃkhya (itself eminently amenable synesthetic experience of numbers) as including viśvā (= all); kha (= space); 7, which symbolizes the heavenly realms or vyāhṛti-s, the 8 directional earth-bearing snakes, the 9, which symbolizes the planets, the 12 Āditya-s on one hand and the 7 continental mountain ranges and 7 solar horses on the other. Thus, the two numbers are described by the entities of entire universe pervaded by the 1000s of Rudra-s with their 11-fold essence and their product is seen yielding a necklace for Maheśvara. In this regard, the use of viśva for 13 is curious. In the bhūta-saṃkhya system, viśva represents the viśvedeva-s (all gods). In the gaṇapāṭha database from which Pāṇini constructed his grammar we find viśvadeva as entry 13 in the manojñādi-gaṇa (GP 177.13), thus lending viśva to encode 13.
Maheśvara’s necklace sequence and its factors
Taking the cue from Śrīpati, we can define a general integer sequence
of Maheśvara’s necklaces thusly:
, where
is an integer and
, where
, is also palindromic. This implies that
Thus, the first few terms of our sequence are:
11, 101, 1001, 10001, 100001, 1000001, 10000001, 100000001, 1000000001, 10000000001…
Taking the cue from Śrīpati, we can define a general integer sequence
11, 101, 1001, 10001, 100001, 1000001, 10000001, 100000001, 1000000001, 10000000001…
We can next ask what are the factors of each
. The first few
are factorized and tabulated below:
From the above one may notice a few obvious things. The first 2 terms of the sequence are primes, 11 and 101, but all others are composites. It is also obvious that 3 can never be a factor of
because the sum of the digits of
will never be 3. We also observe that many of the
have a tendency to have a mix of small factors with a very large one. Śrīpati’s original example
is one such:
. We can then ask questions such as: 1) which primes will divide a given
; 2) For which
will a given prime
be a factor of
. 3) Knowing these, we can ask questions, such as, at what further
will we get a
that will be again divisible by 52579, i.e. they will be further Maheśvara’s necklaces of the type specified by Śrīpati.
We notice right away that 11 divides every other term, i.e. whenever
is odd. We also notice that 101 is a factor of
and every 4th term thereafter. Thus, we can formally write that 11 divides every
when
, where
. Similarly, 101 divides every
when
. Thus, it also becomes obvious that 11 and 101 will never be co-factors of the same
. With closer observation we can see that every prime
that divides a subset of
does so at some
, where
: the reason for writing it this way will become clear below. The first few
that divide
for some
are tabulated in the order of how often they do so along with the form
takes when
is divisible by that
:
At first site the order in which the primes which divide
appear beyond 11 and 101 is puzzling — they wildly differ in magnitude and form. However, a closer examination reveals a striking pattern behind this: a prime
appears in the above list as per the multiplicative order of 10 modulo
: consider
where
; when for the first time
, that
is
the multiplicative order of
. Carl Gauss had famously shown in his Disquisitiones Arithmeticae that
is the length of the repeat pattern of the decimal expansion of
. Thus, it is also clear that
divides
for the first time when
. For example, for 7 we get
because 7 divides
for the first time when
to give
. Likewise, the repeat pattern in the decimal expansion of
which is of length 6. From the above we can easily see why a
will divide
first time when
. Thus, the sequence of
that divide
for the first time will be arranged as per the multiplicative order of
:
Now
divides
for the first time when
. Hence, it will not divide
. However, because
divides
, it therefore divides ![10^{m/2}+1 10^{m/2}+1]()
Now what if
divides some factor of
which takes the form
? We can see from polynomial factorization that a polynomial of the form
can frequently have two factors of the form
, namely
or
. For example,
and
. Now, in our case
and
. Those are primes and the first two terms of
; hence, they will not be divided by any other
. Now, less frequently, other numbers of the form
are divisors of another such number for a larger
. For example, 10001 is a divisor of 1000000000001. So, let us assume for a moment that
divides some
which is a factor of
, then
. If
does divides
, then it also divides
. But
; hence,
cannot divide it because it will only divide a number of of the form
when
, i.e. the multiplicative order of
. Thus,
cannot divide any other
where
.
The Maheśvara’s necklaces
are the sequence
; therefore they would be divided by a given
for the first time when
, which is half the multiplicative order of ![10 \mod p \;_{...\blacksquare} 10 \mod p \;_{...\blacksquare}]()
By the procedure we followed above we can see that, after
,
would divide every
where
A corollary to this is only
with even
can be factors of
for only then
would be an integer. Hence, those primes with odd
such as
;
etc will never be factors of any
.
Armed with the above, we can also tell which will be the next
that will have 52579 as a factor as the original example of Śrīpati. For 52579,
Hence, next term would be:
![f[27]=1000000000000000000000000001 f[27]=1000000000000000000000000001]()
We can also see that some
will always come together as factors of
because they have the same
. Thus, 7 and 13 with
or 19 and 52579 with
will always co-occur. Further, if a certain
satisfies the relationship
for a certain
and
also divides that
then
will occur again as a factor of
. For example, consider
. Now,
here; hence
will be divisible by both 7 and 13 as they have
. However, 7 divides 21. Hence, 7 will occur again as the factor of
. Thus, we have:
. Likewise,
will correspond to a
divisible by both 7 and 13. However, as 39 is divisible by 13, we will have 13 occur again as a factor of
. Thus,
. All other
would be square-free.
The families of Maheśvara’s necklaces
We can represent any given
as a clique of its factors. For example, Figure 1 shows the dodecagonal clique formed by the factors of
, which the most composite
for
. We then merge all cliques sharing common nodes for
. The edges and the nodes are then scaled as per their frequency of occurrence across all 50 cliques. The result is a factor graph for
which is shown in Figure 2.
We can see from Figure 2 that there are totally 6 families of
in this range. These families can be described according to their founder member which is then the divisor of the remaining
of that family. The founder member of each family can be described as the
of the form
, where
corresponds to a particular
:
● When
, we get the 11 family. 11, as we saw above, divides every
corresponding to
. Thus, every
corresponding to an odd
is drawn into this family, there making it the largest of them.
● When
, we get the 101 family. 101 draws all
corresponding to
. Thus it becomes the largest of the even
families.
● When
, we get the 1001 family. 1001 being composite is centered on its two factors 73 and 137 and corresponds to the terms where
.
● When
, we get the 100000001 family. This number being composite is centered on its factors 17 and 5882353 and corresponds to the terms where
.
● When
, we get the 10000000000000001 family centered on its factors 353, 449, 641, 1409 and 69857. This encompasses the terms corresponding to
.
● When
, we get the 100000000000000000000000000000001 family centered on its factors 976193, 19841, 6187457 and 834427406578561. This includes the terms corresponding to
.
● When
● When
● When
● When
● When
Thus, we find that the even
terms are split up among the various families that appear as per the powers of 2.
The largest factor of Maheśvara’s necklace
Given the above information, we can cut down the time in which we factorize Maheśvara’s necklaces and gather the set of factors for the first 300 terms of
. We can then ask which is the largest prime
which divides the corresponding
. Figure 3 shows the plot of
against
.
Given the above information, we can cut down the time in which we factorize Maheśvara’s necklaces and gather the set of factors for the first 300 terms of
We see that the general increase of
appears to be linear with
. It is bounded between lines
, where
and
. The upper bounding slope
is easy to understand: as we observed above, some
tend to have factors widely differing in magnitude; thus the large one is closer in magnitude to
. Trivially, first two terms are primes. There after we get
that are minimally composite. These tend to be of a particular form, e.g.:
![f[19]=11 \times 909090909090909091 f[19]=11 \times 909090909090909091]()
![f[31]= 11 \times 909090909090909090909090909091 f[31]= 11 \times 909090909090909090909090909091]()
![f[53]= 11 \times 9090909090909090909090909090909090909090909090909091 f[53]= 11 \times 9090909090909090909090909090909090909090909090909091]()
Thus, these have factors that approach the upper bounding line.
Thus, these have factors that approach the upper bounding line.
The median value of ratio of
to
is approximately 0.5006. This indicates an even distribution with half the number of
being greater than the
and the other half being lesser than
.
We understand the lower bound is less clearly. Is there a way to derive it from theory alone? One can see that for
corresponding to multiples of 15 there is an increased propensity to be close to the lower bound. This is in part expected from the factorization of polynomial of the form
where
is a multiple of 15. For example, we can see that:
As one can see from the above example the polynomial where
tends to have more factors than other polynomial adjacent to it. Thus, these tend to be among highly composite
; hence, they are more likely to have smaller
than their neighbors.
Tailpiece
Finally, let us look at some actual historical pearl ornaments of Maheśvara. The below were likely ātmaliṅga-s of a medieval South Indian ruler or a royal śaiva officiant.
Finally, let us look at some actual historical pearl ornaments of Maheśvara. The below were likely ātmaliṅga-s of a medieval South Indian ruler or a royal śaiva officiant.
In this example we see the first term of Maheśvara necklace sequence. 11 is of course the characteristic number of the Rudra-s. If we wish to add the bottom two pearls which belong to a different register we get 13 which is the next most frequent factor of
.
In this example we get two sets of pearl ornamentation one with 17 and another with 7. Both of these are factors of
, with 7 being the next most frequent and 17 the founder of a distinct family. Of course there are other numbers with other symbolisms in these ornaments.