Efficient Algorithm for Exponentiation of Imaginary Number



==Quickest way to find 'i' without calculator==

'''Hassan Raza Khan'''<span id="fnc-1"></span>[[#fn-1|<sup>1</sup>]]

<span id="fn-1"></span>
<span style="text-align: center; font-size: 75%;">([[#fnc-1|<sup>1</sup>]]) Hassan Raza Khan is a high-school student at The City School DHA Campus, Lahore.</span>

==1 Introduction==

Say you want to find the value of  i) <math display="inline">i^{8237918273}</math>; ii) <math display="inline">i^{823718932}</math>; iii) <math display="inline">i^{1291393767}</math>; iv) <math display="inline">i^{812731274}</math>;

and there's like 20 seconds left in the exam, and this paper does not allow calculators. You start solving it something like this:

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> i^{8237918273} = \sqrt{-1} \cdot \sqrt{-1} \cdot \sqrt{-1} \cdot \sqrt{-1} \cdot \sqrt{-1} \cdot \sqrt{-1} \cdot \sqrt{-1} \cdot \sqrt{-1} \ldots  </math>
|}
|}

and before you know it, the 20 seconds are over. You couldn't even answer any of the question.

The student next to you, the class topper, recognizes that we can find the value of such a large number of <math display="inline">i</math> using polar coordinates, by employing Euler's formula. He starts writing:

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> e^{i\theta } = \cos (\theta ) + i \sin (\theta{).} </math>
|}
|}

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> i = \cos \left(\frac{\pi }{2}\right)+ i \sin \left(\frac{\pi }{2}\right) </math>
|}
|}

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> i^n = \left(\cos \left(\frac{\pi }{2}\right)+ i \sin \left(\frac{\pi }{2}\right)\right)^n </math>
|}
|}

He now starts using De Moivre's Theorem <math display="inline">n = 8237918273</math>.

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> i^{8237918273} = \cos \left(8237918273 \times \frac{\pi }{2}\right)+ i \sin \left(8237918273 \times \frac{\pi }{2}\right) </math>
|}
|}

By now, the proctor comes and snatches his paper. He could not even solve one of the 4 parts. None of the students was able to solve any part within less than 20 seconds&#8211; except a boy named Draco, who solved all 4 of these parts within less than 20 seconds. Fun fact: he didn't even do any calculations on the paper. He only used the last two digits and performed The Algorithm in his mind.

==2 The Algorithm==

The algorithm involves three steps:

<ol>

<li>'''Real/Imaginary number determination'''

If the last two digits are even, the value remains unchanged (either 1 or -1). If odd, subtract 1 (resulting in either i or -i).     </li>

<li>'''Division by 2'''

Divide the resultant number by 2.     </li>

<li>'''Sign determination'''

If the remaining number after division by 2 is odd, the sign is minus; if it's even, the sign is plus. </li>

</ol>

==3 Proving that only the last two digits are necessary==

To demonstrate that only the last two digits are necessary for determining the value of <math display="inline">i^n</math>, where <math display="inline">n</math> is a large exponent, we can employ the periodicity of powers of <math display="inline">i</math>.

Consider the powers of <math display="inline">i</math> when raised to successive positive integer exponents:

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> \begin{align}     i^1 & = i \\     i^2 & = -1 \\     i^3 & = -i \\     i^4 & = 1 \\     i^5 & = i \\     i^6 & = -1 \\     i^7 & = -i \\     i^8 & = 1 \\     \vdots  \end{align} </math>
|}
|}

We observe that the powers of <math display="inline">i</math> repeat in cycles of four: <math display="inline">i, -1, -i, 1</math>. This periodicity implies that the value of <math display="inline">i^n</math> depends only on the remainder when <math display="inline">n</math> is divided by 4.

Now, let's consider the exponent <math display="inline">n = 8237918273</math>. We want to find <math display="inline">i^n</math>:

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> i^{8237918273} = i^{4k + r}, </math>
|}
|}

Consider <math display="inline">i^{8237918273}</math>. We can express this as <math display="inline">i^{4k + r}</math>, where <math display="inline">k</math> is an integer and <math display="inline">r</math> is the remainder when <math display="inline">8237918273</math> is divided by <math display="inline">4</math>. In this case, <math display="inline">r = 73</math>.

Now, let's break down <math display="inline">i^{8237918273}</math> further:

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>i^{8237918273} = i^{8237918200+73}.</math>
|}
|}

Since <math display="inline">8237918200</math> has two zeros and is divisible by 100, and since 100 is divisible by 4, we can conclude that <math display="inline">8237918200</math> is divisible by 4.

Therefore, <math display="inline">i^{4k+r} = i^{8237918200+73} = i^{8237918200} \cdot i^{73}.</math>

Now, <math display="inline">i^{8237918200}</math> simplifies to <math display="inline">1</math> since any power of <math display="inline">i</math> with an exponent divisible by 4 is equal to 1. Therefore, we have:

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>1 \cdot i^{73} = i^{73}.</math>
|}
|}

So, essentially, <math display="inline">i^{8237918273}</math> is congruent to <math display="inline">i^{73}</math>. This result demonstrates a general pattern: for any large exponent <math display="inline">n</math>, <math display="inline">i^n</math> is congruent to <math display="inline">i^{r}</math>, where <math display="inline">r</math> is the remainder when <math display="inline">n</math> is divided by 4. Thus, only the last two digits of the exponent (<math display="inline">73</math> in this case) are crucial for determining the value of <math display="inline">i^n</math>.

This observation allows for a more efficient approach when faced with large exponentiation of <math display="inline">i</math> without the need to calculate the entire exponent. The periodic nature of <math display="inline">i</math> simplifies the computation, making it feasible to focus solely on the last two digits of the exponent.

==4 Demonstrating the algorithm==

What Draco essentially did was use the above three-step algorithm like so:

<ol>

<li><math display="inline">i^{8237918273}</math>;

<ol>

</li>
<li>'''Real/Imaginary number determination'''

The last two digits are 73, an odd number. So it must be an imaginary number (either <math display="inline">i</math> or <math display="inline">-i</math>). Since the number is odd, we subtract it by 1 to make it parallel to its even counterpart, resulting in 72.         </li>
<li>'''Division by 2'''

Dividing 72 by 2 results in 36, an even number.         </li>
<li>'''Sign determination'''

Since the resultant number after division by 2 is even, we can establish that the sign of this number would be positive.

Since the number is both imaginary and positive, it must be <math display="inline">i</math>.     </li>

</ol>

<li><math display="inline">i^{823718932}</math>;

<ol>

</li>
<li>'''Real/Imaginary number determination'''

The last two digits are 32, an even number. So it must be a real number (either 1 or -1). Since the number is even, we do not subtract it, so it remains 32.         </li>
<li>'''Division by 2'''

Dividing 32 by 2 results in 16, an even number.         </li>
<li>'''Sign determination'''

Since the resultant number after division by 2 is even, we can establish that the sign of this number would be positive.

Since the number is both real and positive, it must be 1.     </li>

</ol>

<li><math display="inline">i^{1291393767}</math>;

<ol>

</li>
<li>'''Real/Imaginary number determination'''

The last two digits are 67, an odd number. So it must be an imaginary number (either <math display="inline">i</math> or <math display="inline">-i</math>). Since the number is odd, we subtract it by 1 to make it parallel to its even counterpart, resulting in 66.         </li>
<li>'''Division by 2'''

Dividing 66 by 2 results in 33, an odd number.         </li>
<li>'''Sign determination'''

Since the resultant number after division by 2 is odd, we can establish that the sign of this number would be negative.

Since the number is both imaginary and negative, it must be <math display="inline">-i</math>.     </li>

</ol>

<li><math display="inline">i^{812731274}</math>;

<ol>

</li>
<li>'''Real/Imaginary number determination'''

The last two digits are 74, an even number. So it must be a real number (either 1 or -1). Since the number is even, we do not subtract it, so it remains 74.         </li>
<li>'''Division by 2'''

Dividing 74 by 2 results in 37, an odd number.         </li>
<li>'''Sign determination'''

Since the resultant number after division by 2 is odd, we can establish that the sign of this number would be negative.

Since the number is both real and negative, it must be <math display="inline">-1</math>.     </li>

</ol>

</ol>

In Table 1, we present a matrix illustrating the relationship between the final and initial numbers corresponding to different values of the complex unit <math display="inline">i</math>. The table categorizes the initial and final numbers as either even or odd and displays the resulting values of <math display="inline">i</math> for each combination. This matrix provides a concise reference for understanding the periodicity of <math display="inline">i</math> and its cyclic behavior based on the evenness or oddness of the exponents. Memorizing this table can make you the Draco in the same circumstance.


{|  class="floating_tableSCP wikitable" style="text-align: left; margin: 1em auto;min-width:50%;"
|+ style="font-size: 75%;" |Table. 1 Matrix of Initial and Final Numbers Corresponding to Values of <math>i</math>
|-
| rowspan='2' colspan='2' style="border-right: 2px solid;" | 
| colspan='2' style="text-align: center;border-right: 2px solid;border-left: 2px solid;" | Initial Number
|-
| colspan='1' style="text-align: center;border-right: 2px solid;border-left: 2px solid;" | Even
| colspan='1' style="text-align: center;border-right: 2px solid;border-right: 2px solid;" | Odd
|- style="border-top: 2px solid;"
| rowspan='2' colspan='1' style="border-left: 2px solid;border-right: 2px solid;border-bottom: 2px solid;" | Final Number
| style="border-right: 2px solid;" | Even 
| colspan='1' style="text-align: center;border-right: 2px solid;border-left: 2px solid;" | 1
| colspan='1' style="text-align: center;border-right: 2px solid;border-right: 2px solid;" | i
|- style="border-bottom: 2px solid;"
| style="border-right: 2px solid;" | Odd  
| colspan='1' style="text-align: center;border-right: 2px solid;border-left: 2px solid;" | -1
| colspan='1' style="text-align: center;border-right: 2px solid;border-right: 2px solid;" | -i

|}

==5 Purpose of Odd/Even Analysis==

The divisional remainder and its congruent complex numeber counterpart can be illustrated as below:

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>1 \equiv i \pmod{4} </math>
|-
| style="text-align: center;" | <math> 2 \equiv -1 \pmod{4} </math>
|-
| style="text-align: center;" | <math> 3 \equiv -i \pmod{4} </math>
|-
| style="text-align: center;" | <math> 0 \equiv 1 \pmod{4} </math>
|}
|}

The odd/even analysis serves as a strategic approach to efficiently determine the divisional remainder of the numbers in the context of complex exponentiation. By focusing solely on the last two digits of the exponent, we can employ a systematic method to evaluate whether the number is divisible by 4.

For instance, consider an exponent with the last two digits being 32. If we divide 32 by 2, we obtain 16, an even number. This indicates that the original exponent is divisible by 2 twice, meaning it is divisible by 4. Therefore, the number is congruent to 0 modulo 4.

On the contrary, let's take another example with the last two digits being 34. If we divide 34 by 2, we get 17, an odd number. In this case, 34 can only be divided by 2 once while remaining a whole number. Consequently, it is congruent to 2 modulo 4.

The same approach is applied to odd numbers, where they are subtracted by 1 in order to align them with their even counterparts. This adjustment ensures that odd numbers maintain their distinctive properties in the modulo 4 arithmetic.

It's worth noting that one could also add 1 instead of subtracting 1 from the odd number. However, by adopting the subtraction approach, the signs become universally positive for even final numbers and negative for odd final numbers. This consistency simplifies the association of odd numbers with negative results and even numbers with positive results, aiding in a more straightforward interpretation of the outcomes.

In summary, the odd/even analysis not only streamlines the computation of complex exponentiation but also provides a clear and consistent method for determining the divisional remainder based on the last two digits of the exponent.

==Acknowledgements==

I would like to express my gratitude for the support and encouragement I received from my friends, family, and teachers throughout the course of this research paper.

Lastly, I owe a special debt of gratitude to Dr Javed Hussain and Dr Zarqa Bano from IBA Sukkur University for their expert guidance. As a high school student new to research papers, their assistance was invaluable. They not only helped me with the paper but also instilled in me the importance of intellectual curiosity and hard work.

===BIBLIOGRAPHY===

Khan, H. R. (2023, October 27). Efficient Algorithm for Identifying Repeating Patterns Modulo 4 and its Application in Complex Numbers. https://doi.org/10.31219/osf.io/8zc7e

Euler, L. (1748). ''Introductio in analysin infinitorum''. Apud Marcum-Michaelem Bousquet & Socios.

De Moivre, A. (1730). ''Miscellanea Analytica de Seriebus et Quadraturis'' [Analytical Miscellany on Series and Integration].

Zhou, X. (2017). ''Number Theory - Modular Arithmetic: Math for Gifted Students''. CreateSpace Independent Publishing Platform.

Halmos, P. (1980). The heart of mathematics. ''American Mathematical Monthly, 87'', 519–524.

Nahin, P. (1998). ''An Imaginary Tale: The Story of <math>\sqrt{-1}</math>''. New Jersey: Princeton University Press.