2016-09-17
This slide is based on the article arXiv:1601.06509.
We study the longest cycles consisted by the quadratic residues modulo prime numbers in some way. We give some formulae about the maximum length of the cycles. Especially, the length of the longest cycles for modulo Fermat primes is one.
We find this phenomenon by computation. So let us begin with some examples. Let's consider quadratic residues module $999$ first.
\[ \begin{cases} 454^2\equiv 322\ (\mod 999),\\ 322^2\equiv 787\ (\mod 999),\\ 787^2\equiv 988\ (\mod 999),\\ 988^2\equiv 121\ (\mod 999),\\ 121^2\equiv 655\ (\mod 999),\\ 655^2\equiv 454\ (\mod 999),\\ \end{cases} \qquad \begin{cases} 445^2\equiv 223\ (\mod 999),\\ 223^2\equiv 778\ (\mod 999),\\ 778^2\equiv 889\ (\mod 999),\\ 889^2\equiv 112\ (\mod 999),\\ 112^2\equiv 556\ (\mod 999),\\ 556^2\equiv 445\ (\mod 999).\\ \end{cases} \]
We find that these numbers $\{454,322,787,988,121,655\}$ consist a cycle, also for the numbers $\{445,223,778,889,112,556\}$. The length of the cycle is $6$. In fact it is the maximum number of the elements in the cycles. We called it the length of the
Look another example for modulus $99$, we get
\[ \begin{cases} 22^2\equiv 88\ (\mod 99),\\ 88^2\equiv 22\ (\mod 99),\\ \end{cases} \qquad \begin{cases} 70^2\equiv 49\ (\mod 99),\\ 49^2\equiv 25\ (\mod 99),\\ 25^2\equiv 31\ (\mod 99),\\ 31^2\equiv 70\ (\mod 99).\\ \end{cases} \]
We see that, there exist a smaller cycle. Here we consider the largest cycles since they are more interesting.
Definition(the length of the cycle). Consider the equation $x^2\equiv a\ (\mod m)$. If there exist a series of numbers $\{x_i\}_{i=1}^{k}$ such that \[ \begin{cases} x_1^2 &\equiv x_2\ (\mod m),\\ x_2^2 &\equiv x_3\ (\mod m),\\ &\cdots\\ x_{k-1}^2 &\equiv x_k\ (\mod m),\\ x_k^2 &\equiv x_1\ (\mod m),\\ \end{cases} \] then we call these $k$ numbers consist a cycle modulo $m$. The number $k$ is defined as the length of the cycle.
It infers that. \[ x_i^{2^k}\equiv x_i\ (\mod m),\quad\text{for}\ i=1,2,\ldots,k. \]
There exists at least one largest cycle. We denote the length of it by $L(m)$. For example, $L(99)=4$, $L(999)=6$. In the next section, we try to find out the formula for $L(m)$.
The maximum length of the cycles modulo prime numbers are list as follows:
Obviously, it establishes a map from the set of primes $\mathrm{PRIMES}$ to the set of positive numbers $\mathbb{Z}^{+}$:
\[ \begin{array}{rcl} L:\ \mathrm{PRIMES} &\rightarrow &\mathbb{Z}^{+},\\ p & \mapsto & L(p).\\ \end{array} \]
In fact, E. L. Blanton, Jr., S. P. Hurd and J. S. McCranie [Blanton-Hurd-McCranie] had considered this question. We can find this sequence on the site of
For composites, we have
And when modulo $p^2$, we have
And consider modulo the powers of prime numbers. For examples:
The data listed in these tables are all verified by computations.
$L(p^2)$ is relevant to $p-1$ and $L(p)$. It is rather complicated and we discuss it later. But for few exceptions we have the explicit formula. For example, the Fermat primes.
For the first five Fermat primes \[ F_0=3,\quad F_1=5,\quad F_2=17,\quad F_3=257,\quad F_4=65537, \] we have $L(F_i)=1$, $i=0,1,2,3,4$.
For $p^2=11^2$, we have $L(121)=20$. The largest cycle $\{x_1,x_2,\ldots,x_{20}\}$ is \[ 4, 16, 14, 75, 59, 93, 58, 97, 92, 115, 36, 86, 15, 104, 47, 31, 114, 49, 102, 119. \]
Of course, there are also cycles with length equal to $10$. For example, one of them is \[ 100, 78, 34, 67, 12, 23, 45, 89, 56, 111. \]
For $p^2=17^2$, we have $L(289)=8=\frac{1}{2}(17-1)$. For example, one of the cycles is \[ 256, 222, 154, 18, 35, 69, 137, 273. \]
For $p^2=257^2$, we have $L(66049)=16=\frac{1}{2^4}(257-1)$. For example, one of the cycles is \[ \begin{array}{rrrrrrrr} 65536, & 65022, & 63994, & 61938, & 57826, & 49602, & 33154, & 258,\\ 515, & 1029, & 2057, & 4113, & 8225, & 16449, & 32897, & 65793.\\ \end{array} \]
For $p^2=65537^2$, we have $L(4295098369)=32=\frac{1}{2^{11}}(65537-1)$. For example, one of the cycles is \[ \begin{array}{rrrr} 4294967296, & 4294836222, & 4294574074, & 4294049778, \\ 4293001186, & 4290904002, & 4286709634, & 4278320898, \\ 4261543426, & 4227988482, & 4160878594, & 4026658818, \\ 3758219266, & 3221340162, & 2147581954, & 65538, \\ 131075, & 262149, & 524297, & 1048593, \\ 2097185, & 4194369, & 8388737, & 16777473, \\ 33554945, & 67109889, & 134219777, & 268439553, \\ 536879105, & 1073758209, & 2147516417, & 4295032833 \\ \end{array} \]
Note that \[ 17=2^{2^2}+1, \quad 257=2^{2^3}+1,\quad 65537=2^{2^4}+1. \]
They are the Fermat numbers. For $k=0,1,2,3,4$, $2^{2^k}+1=3,5,17,257,65537$ are all prime numbers. But $2^{32}+1=4294967297$ is not a prime. It equals $641 \times 6700417$.
It is conjectured that there are only 5 Fermat prime numbers. Currently it has been shown that $2^{2^k} + 1$ is composite for $5\leqslant k\leqslant 32$
Prop 4.1. For the Fermat primes, i.e., the prime numbers of the form $p=2^{2^k}+1$, we have $L(p)=1$.
Proof. Suppose $F_k=2^{2^k}+1$ is a prime, $k\geqslant 1$. Then $\varphi(F_k)=2^{2^k}$. To find the largest cycle of the quadratic residue equation of modulo $F_k$. We consider the equation \[ x^{2^m}\equiv x\ (\mod F_k), \] where $x$ is coprime with $F_k$. Then it is equivalent to the equation \begin{equation} x^{2^m-1}\equiv 1\ (\mod F_k).\tag{4.1} \end{equation} By Lemma 3.4, it has a solution. And the number of the solutions is equal to \[ \gcd(2^m-1, 2^{2^k})=1. \] Since $x=1$ is the trivial solution of (4.1), hence there is no other solutions. Hence we have $L(F_k)=1$.
Remark. Since there is a primitive root $\mod F_k$, by Lemma 3.5, for any $1\leqslant m\leqslant 2^k$, $x^{2^m}\equiv 1\ (\mod F_k)$ have exact $2^m$ solutions. Suppose $k\geqslant 1$. Let's consider the equation \[ x^2\equiv -1\ (\mod F_k). \] By Lemma 3.4, it has a solution since $(-1)^{2^{2^k-1}}\equiv 1\ (\mod F_k)$. And there are exactly two solutions of this equation. By Lemma 3.3, these two solutions $\{x_1,x_2\}$ are the primitive roots modulo $F_k$. By definition of primitive root, they have the maximal order. That is, $\mathrm{ord}(x_1)=\mathrm{ord}(x_2)=\varphi(F_k)=2^{2^k}$. Hence, any other number $a$ coprime to $F_k$ (not the primitive root) with order $\mathrm{ord}(a)<2^{2^k}$. And $\mathrm{ord}(a)| 2^{2^k-1}$. Hence, we have \begin{equation}\label{eqn:0} a^{2^{2^k-1}}\equiv 1\ (\mod F_k). \end{equation} It infers that the equation $x^2\equiv a\ (\mod F_k)$ has a solution for any $a$ coprime to $F_k$. (Use Lemma 3.4 again.)
Prop 4.3. For the Fermat primes, i.e., the prime numbers of the form $p=2^{2^k}+1$, we have \[ L(p^2)=L((2^{2^k}+1)^2)=\dfrac{1}{2^{2^k-k-1}}(p-1)=2^{k+1}. \]
Proof. We have verified it for $k=0,1,2,3,4$, the known Fermat primes. Suppose $p=F_k=2^{2^k}+1$ is a prime, $k\geqslant 1$. Similar to the Proposition 1, we consider the equation \[ x^{2^m}\equiv x\ (\mod F_k^2),\quad x\neq F_k, \] and find the nontrivial solution for the maximal $m$. Since $F_k$ is a prime, for $x\neq F_k$, $x$ is coprime with $F_k^2$. Then it is equivalent to \begin{equation} x^{2^m-1}\equiv 1\ (\mod F_k^2).\tag{4.3} \end{equation} By Lemma 3.4, there are $\gcd(2^m-1, \varphi(F_k^2))$ solutions. By Lemma 3.1, $\varphi(F_k^2)=(F_k-1)F_k=2^{2^k}(2^{2^k}+1)$. Thus, if $\gcd(2^m-1, 2^{2^k}+1)>1$, the equation (4.3) has nontrivial solutions. The number of the solutions is equal to $2^{2^k}+1$. Because $F_k=2^{2^k}+1$ is a prime number, \[ \gcd\bigl(2^m-1, 2^{2^k}+1\bigr)>1\ \Leftrightarrow\ (2^{2^k}+1)\,\bigr|\,(2^m-1). \] It infers that $m=h\cdot 2^{k+1}$, here $h\geqslant 1$. And there are $2^{2^k}+1$ solutions include the trivial $x=1$. So, for example, if $h=2$, then \[ x^{2^m}=x^{2^{2^{k+1}}\cdot 2^{2^{k+1}}}=\Bigl(x^{2^{2^{k+1}}}\Bigr)^{2^{2^{k+1}}}\equiv x^{2^{2^{k+1}}}\equiv x\ (\mod F_k^2). \] Hence, $L(F_k^2)=2^{k+1}$ for Fermat prime $F_k$.
Definition. Let $n\in\mathbb{N}$. A primitive root mod $n$ is a residue class $\alpha\in(\mathbb{Z}/n\mathbb{Z})^*$ with maximal order, i.e., $\mathrm{ord}(\alpha)=\varphi(n)$.
Euler's totient function $\varphi(n)$ counts the positive integers up to a given integer $n$ that are relatively prime to $n$
Lemma 3.1. The Euler's totient function $\varphi(m)$ have the following formula for $m=p^n$ (where $p$ is a prime): \[ \varphi(p^n)=(p-1)p^{n-1}. \]
Euler's totient function is a multiplicative function, meaning that if two numbers $m$ and $n$ are coprime, then $\varphi(mn) = \varphi(m) \varphi(n)$. Generally, \[ \varphi(mn)=\varphi(m)\varphi(n)\cdot\frac{d}{\varphi(d)}, \] where $d=\mathrm{gcd}(m,n)$. It infers that \[ \varphi(n^m)=n^{m-1}\varphi(n). \]
Lemma 3.2. (Euler's Criterion, 1748) Let $p$ be an odd prime and $a$ not divisible by $p$. Then $a$ is a quadratic residue modulo $p$ if and only if \[ a^{(p-1)/2}\equiv 1\ (\mod p). \]
The proof uses the fact that the residue classes modulo a prime number are a field.
Lemma 3.3. (
Proof. Let $a$ be a quadratic nonresidue of the Fermat prime $F_n$, and let $e=\mathrm{ord}_{F_n}a$. According to Fermat's little theorem, $a^{F_n-1}\equiv 1\ (\mod F_n)$. So $e | F_n-1=2^{2^n}$. It follows that $e=2^k$ for some nonnegative $k\leqslant 2^n$. On the other hand, by Euler's criterion, \[ a^{(F_n-1)/2}=a^{2^{2^n-1}}\equiv -1\ (\mod F_n). \] On the other hand, if $k < 2^n$, then $2^k | 2^{2^n-1}$ and so $a^{2^{2^n-1}}\equiv 1\ (\mod F_n)$, which is a contradiction. So, $k=2^n$ and $\mathrm{ord}_{F_n}a=2^{2^n}$. Therefore, $a$ is a primitive root modulo $F_n$.
Lemma 3.4.
On one extreme, this means if $\gcd(d,\varphi(n))=1$, then there is a unique $d$-th root $\mod n$ of any $a\in\mathbb{Z}$ coprime to $n$. On the other hand, it means that for any $d|\varphi(n)$, there are exactly $d$ solutions $\mod n$ to the equation $x^d\equiv 1\ (\mod n)$. This in fact characterizes the existence of primitive roots:
Lemma 3.5. Let $n\in\mathbb{N}$. The following are equivalent:
Corollary 3.6. Given any odd prime number $p$. Let $d=\mathrm{gcd}\bigl(2^{p-1}-1,(p-1)p\bigr)$. Then there are exactly $d$ solutions mod $p^2$ to the equation \[ x^d\equiv 1\ (\mod p^2). \]
Since there are $(p-1)/2$ quadratic residues mod $p$, we have $L(p)\leq (p-1)/2$. However, $1$ maps to itself under squaring mod $p$, so we expect $L(p)\leq (p-3)/2$. In the theorem below, we state a condition for this to happen.
Theorem 4.4. (by Dr.Timothy) Let $p$ be a non-Fermat prime. $L(p)=(p-3)/2$ if and only if $q=(p-1)/2$ is prime and $2$ is a generator of $(\mathbb{Z}/q\mathbb{Z})^{*}$.
Proof. Let $d$ be an odd divisor of $p-1$. Then there is a solution to the equations \begin{equation} x^{d}\equiv 1\ (\mod p),\quad x^i\not\equiv 1\ (\mod p),\tag{4.4} \end{equation} $1\leq i < d$. Let $x$ satisfy $(4.4)$. Then such $x$ is a solution to the equation $$ x^{2^m}\equiv x\ (\mod p) $$ for all $m$ such that $\mbox{ord}_d(2)|m$ where $\mbox{ord}_d(2)$ is the multiplicative order of $2$ mod $d$. Furthermore, $$ L(p)=\max_{\substack{d|p-1\\d \equiv 1 \bmod 2}}\mbox{ord}_d(2)\leq\max_{\substack{d|p-1\\d \equiv 1 \bmod 2}}d-1. $$ We need only consider the case where $(p-1)/2$ is odd, since if $d < (p-1)/2$, then $\mbox{ord}_d(2) < (p-3)/2$. Now any odd $d$ dividing $(p-1)$ satisfies $\mbox{ord}_d(2)\leq \mbox{ord}_{(p-1)/2}(2)$, hence $L(p)=\mbox{ord}_{(p-1)/2}(2)$. This is $(p-3)/2$ if and only if $q=(p-1)/2$ is prime and $2$ is a generator of $(\mathbb{Z}/q\mathbb{Z})^{*}$.
Remark 4.5. (by Dr.Timothy) Similar to the case for $L(p)$, we have \[ L(p^2)=\max_{\substack{d|p(p-1)\\d \equiv 1 \bmod 2}}\mbox{ord}_d(2). \] Let $n$ be the largest odd divisor of $p-1$. Then $L(p)=\mbox{ord}_n(2)$ and $L(p^2)=\mbox{ord}_{np}(2)=\mathrm{lcm}(L(p),\mbox{ord}_p(2))$, where $\mathrm{lcm}(a_1,\ldots,a_n)$ stands for the lowest common multiple of $a_1,\ldots,a_n$.
From the above tables, we guess there are formulas for $L(p^n)$ for $n\geqslant 2$.
Prop 4.6. Let $p$ be a prime. Then we have \[ L(p^2)\leqslant(p-1)L(p). \]
Proof. Suppose $p$ is a prime. Similar to the Proposition 1, we consider the equation \[ x^{2^m}\equiv x\ (\mod p^2), \] and find the nontrivial solution for the maximal $m$. Suppose $x$ is coprime with $p$. Then it is equivalent to \begin{equation}\label{eqn:3} x^{2^m-1}\equiv 1\ (\mod p^2). \end{equation} By Lemma 3.4, there are $\gcd(2^m-1, \varphi(p^2))$ solutions. By Lemma 3.1, $\varphi(p^2)=(p-1)p$. By Euler theorem, for every odd prime $p$, we have \[ 2^{p-1}\equiv 1\ (\mod p). \] Thus, $p|(2^{p-1}-1)$ for $p>2$. Hence \[ \begin{split} & \gcd(2^m-1, (p-1)p)>1\\ \Leftrightarrow\ &\gcd(2^{p-1}-1,2^{m}-1)>1\quad\text{or}\quad\gcd(p-1,2^m-1)>1.\\ \end{split} \] Hence for find the largest $m$, we only need to consider $\gcd(2^{p-1}-1,2^{m}-1)>1$ with $m\geqslant p-1$. Which infers that $m=h(p-1)$, $h\geqslant 1$. On the other hand, we can proved that $h\leqslant L(p)$. In fact the following figure describes the relationship between modulo $p$ and modulo $p^2$. Suppose \[ x_1^2\equiv x_2\ (\mod p),\quad x_1^2\equiv x'_2\ (\mod p^2). \] Then there exists some integers $s$ and $t$ such that \[ x_1^2=x_2+sp,\quad x_1^2=x'_2+tp^2. \] Thus, we have \[ x_2-x'_2=p(pt-s). \] It shows that $x'_2$ and $x_2$ lie in the same vertical line in Figure 1. And ${x'}^2_2\equiv x_2^2\equiv x_3\ (\mod p)$. Let ${x'}^2_2\equiv x'_3\ (\mod p^2)$. Then it infers that $ x'_3\equiv x_3\ (\mod p)$. Hence, each $x'_i$ lies in the same line of $x_i$. Therefore, the elements in the biggest cycle all lie in the area bounded by dashed lines below. Because there are no more than $p-1$ vertical lines, the total numbers of the elements in the biggest cycle is less or equal to $(p-1)L(p)$. It means that $L(p^2)\leqslant (p-1)L(p)$.
For $n\geqslant 3$, we have \[ L(p^n)=p^{n-2}\cdot L(p^2). \]
Proof. By using the same idea (as illustrated in Figure 1) in the proof of Proposition 4.6, we know that if $x_i^2\equiv x_j\ (\mod p^2)$ and $x_i^2\equiv x'_j\ (\mod p^3)$, then $x'_j\equiv x_j\ (\mod p^2)$. Thus, following the arguments, we get $L(p^3)\leqslant p\cdot L(p^2)$. Note that the multiplication factor is $p$, not $p-1$. If the largest cycle modulo $p^3$ is less than $pL(p^2)$. That is, the largest cycle doesn't take all of the elements in the area similar in Figure 1. Then, there must be another largest cycle modulo $p^3$. If use curve to describe the cycle, then the curve intersect each vertical line once. They all project onto the "line" $[0,p^2]$. And they are not intersect (i.e., have no common elements). Since $L(p^3)>L(p^2)$, the largest cycle's length must greater than $L(p^2)$. And note that there are $p$ rows($p$ is a prime number), then the number of largest cycles must be $1$. That is $L(p^3)=pL(p^2)$. For the general cases, the idea is the same. By conduction we complete the proof.
Finally, we have
Prop. Suppose $m$ is a composite number, and $m=p_{i_1}^{s_1}p_{i_2}^{s_2}\cdots p_{i_t}^{s_t}$. Here $p_{i_k}$ are all prime numbers. Then we have \[ L(m)=L(p_{i_1}^{s_1}p_{i_2}^{s_2}\cdots p_{i_t}^{s_t})=\mathrm{lcm}\Bigl(L(p_{i_1}^{s_1}), L(p_{i_2}^{s_2}), \ldots, L(p_{i_t}^{s_t})\Bigr). \]
Proof. Let $m=p_{i}^{s_i}p_{j}^{s_j}$. Consider equation \[ x^{2^m}\equiv x\ (\mod p_{i}^{s_i}p_{j}^{s_j}), \] where $(x,p_i)=(x,p_j)=1$. Then it is equivalent to the equation \[ x^{2^m-1}\equiv 1\ (\mod p_{i}^{s_i}p_{j}^{s_j}). \] By Lemma 3.4, it has $\gcd(2^m-1, \varphi(p_{i}^{s_i}p_{j}^{s_j}))$ solutions. \[ \varphi(p_{i}^{s_i}p_{j}^{s_j})=\varphi(p_{i}^{s_i})\varphi(p_{j}^{s_j})=(p_i-1)p_i^{s_i-1}\cdot (p_j-1)p_j^{s_j-1}. \] Hence, we have $L(p_{i}^{s_i}p_{j}^{s_j})=\mathrm{lcm}(L(p_{i}^{s_i}),L(p_{j}^{s_j}))$. By conduction we will complete the proof.
We give a few examples.
\[ \begin{aligned} L(15)&=L(3\cdot 5)=1=\text{lcm}(L(3),L(5))=\text{lcm}(1,1),\\ L(45)&=L(3^2\cdot 5)=2=\text{lcm}(L(3^2),L(5))=\text{lcm}(2,1),\\ L(135)&=L(3^3\cdot 5)=6=\text{lcm}(L(3^3),L(5))=\text{lcm}(6,1),\\ L(75)&=L(3\cdot 5^2)=4=\text{lcm}(L(3),L(5^2))=\text{lcm}(1,4),\\ L(225)&=L(3^2\cdot 5^2)=4=\text{lcm}(L(3^2),L(5^2))=\text{lcm}(2,4),\\ L(675)&=L(3^3\cdot 5^2)=12=\text{lcm}(L(3^3),L(5^2))=\text{lcm}(6,4),\\ L(375)&=L(3\cdot 5^3)=20=\text{lcm}(L(3),L(5^3))=\text{lcm}(1,20),\\ L(1125)&=L(3^2\cdot 5^3)=20=\text{lcm}(L(3^2),L(5^3))=\text{lcm}(2,20),\\ L(3375)&=L(3^3\cdot 5^3)=60=\text{lcm}(L(3^3),L(5^3))=\text{lcm}(6,20),\\ \end{aligned} \]
For $L(p^2)$, the situations are complicated. We list the various formulae here. \noindent(1) The first few primes obey the formula $L(p^2)=(p-1)L(p)$ are: \[ 2,3,5,29,179,293,317,\ldots \] (2) The first few primes obey the formula $L(p^2)=\frac{1}{2}(p-1)L(p)$ are: \[ \begin{aligned} &7,11,13,17,23,47,59,67,71,83,103,107,131,139,\\ &167,173,191,227,239,263,269,347,\ldots\\ \end{aligned} \] (3) The first few primes obey the formula $L(p^2)=\frac{1}{4}(p-1)L(p)$ are: \[ 53,61,97,113,149,193,349,\ldots \] (4) The first few primes obey the formula $L(p^2)=\frac{1}{6}(p-1)L(p)$ are: \[ 19,31,37,43,79,199,211,223,229,277,283,\ldots \] (5) The first few primes obey the formula $L(p^2)=\frac{1}{7}(p-1)L(p)$ are: \[ 197,\ldots \] (6) The first few primes obey the formula $L(p^2)=\frac{1}{8}(p-1)L(p)$ are: \[ 41,89,137,233,281,353,\ldots \] (7) The first few primes obey the formula $L(p^2)=\frac{1}{10}(p-1)L(p)$ are: \[ 311,\ldots \] (8) The first few primes obey the formula $L(p^2)=\frac{1}{12}(p-1)L(p)$ are: \[ 157,181,\ldots \] (9) The first few primes obey the formula $L(p^2)=\frac{1}{18}(p-1)L(p)$ are: \[ 127,271,307,\ldots \] (10) The first few primes obey the formula $L(p^2)=\frac{1}{20}(p-1)L(p)$ are: \[ 101,\ldots \] (11) The first few primes obey the formula $L(p^2)=\frac{1}{24}(p-1)L(p)$ are: \[ 73,313,\ldots \] (12) The first few primes obey the formula $L(p^2)=\frac{1}{40}(p-1)L(p)$ are: \[ 241,\ldots \] (13) The first few primes obey the formula $L(p^2)=\frac{1}{48}(p-1)L(p)$ are: \[ 337,\ldots \] (14) The first few primes obey the formula $L(p^2)=\frac{1}{50}(p-1)L(p)$ are: \[ 151,\ldots \] (15) The first few primes obey the formula $L(p^2)=\frac{1}{54}(p-1)L(p)$ are: \[ 109,163,\ldots \] (16) The first few primes obey the formula $L(p^2)=\frac{1}{110}(p-1)L(p)$ are: \[ 331,\ldots \] (17) The first few primes obey the formula $L(p^2)=\frac{1}{250}(p-1)L(p)$ are: \[ 251,\ldots \]
Although we do not sure whether there are some definite laws about $L(p^2)$, some interesting phenomenons should be noted. If write $L(p^2)=\frac{1}{k}(p-1)L(p)$, then the differences between the adjacent primes in the list of each cases are all divisible by the corresponding $k$.
We would like to express our gratitude to Professor Vilmos Komornik who invite the author to visit the Department of Mathematics in University of Strasbourg. This work is done during the stay of the author in Strasbourg.
We also express our gratitude to Dr. Timothy Foo (a researcher in School of Computing, National University of Singapore, Email: dcsfcht@nus.edu.sg) for sharing his finding. He gives the Theorem 4.4 with the proof and the Remark 4.5.
We thank Professor Vilmos Komornik and Dr. Timothy Foo for discussing this question. And thanks Michel Marcus for giving me the link of