牛客网暑期ACM多校训练营(第五场) B Pell方程

题目链接

---

题意：
定义好数$n$满足在区间$[n^2+1,n^2+2n]$存在一个数$x$满足$x|n^4$
给定一个正整数$m$，求不小于$m$的最小好数。
$m\le 10^{1000}$

---

思路：

数论好题~
设$x=n^2+a$
则可得：$n^2+a|n^4$

此时试着构造式子来建立$n^4$和 $n^2+a$的关系，得到：

$$n^4=(n^2+a{)}^2-a^2-2a{n}^2$$

和

$$n^4=(n^2+a)(n^2-a)+a^2$$

因为$(n^2+a)|n^4$，且因为$(n^2+a)|(n^2+a)(n^2-a)$
可推出：

$$(n^2+a)|n^4-(n^2+a)(n^2-a)$$

利用第二个关系式$n^4=(n^2+a)(n^2-a)+a^2$，得：

$$(n^2+a)|a^2$$

等价于：

$$a^2=m(n^2+a)$$

又定义式：$x=n^2+a$

故$1<=a<=2n$，即$1<=a^2<=4n^2<4(n^2+a)$
故
$m=1,2,3$

---

当$m=1$时，$a^2=n^2+a$

即
$$a(a-1)=n^2$$

因$n,a$均为整数，故无解

---

当$m=2$时，$a^2=2n^2+2a$
则：
$$(a-1{)}^2-2n^2=1$$

令 $x=a-1$，$y=n$，则化为：

$$x^2-2y^2=1$$

此为Pell方程的标准形式

利用连分数可以解出最小整数解，解为$(3,2)$
具体原理不再赘述，附上Python代码，输入任意的$d$，则可以输出$x^2-d{y}^2=1$的最小整数解。

def getPell(n):
    j = 1
    while j*j<n:
        j += 1
    if j*j == n:
        print(j,1)
    if j*j > n:
        p = [0 for i in range(0, 1001)]
        q = [0 for i in range(0, 1001)]
        a = [0 for i in range(0, 1001)]
        g = [0 for i in range(0, 1001)]
        h = [0 for i in range(0, 1001)]
        p[1] = q[0] = h[1] = 1
        p[0] = q[1] = g[1] = 0
        a[2] = j - 1
        i = 2
        while 1:
            g[i]=-g[i-1]+a[i]*h[i-1]
            h[i]=(n-g[i]*g[i])/h[i-1]
            a[i+1]=(g[i]+a[2])/h[i]
            p[i]=a[i]*p[i-1]+p[i-2]
            q[i]=a[i]*q[i-1]+q[i-2]
            if(p[i]*p[i]-n*q[i]*q[i]==1):
                print(p[i], q[i])
                break
            i += 1
        
if __name__ == "__main__":
    d = int(input())
    getPell(d)

那应该如何推出其他的解呢，我们可以先看一种通用的解法，即将最小整数解回代：

$$1=3^2-2\ast 2^2=(3-2\sqrt{2})(3+2\sqrt{2})$$

两边同时平方，得：

$$1^2=1=(3-2\sqrt{2}{)}^2(3+2\sqrt{2}{)}^2=(17-12\sqrt{2})(17+12\sqrt{2})$$

故第二组解为：(17, 12)
…
同理易得，对于第$k$组解，有：

$$x_k+\sqrt{2}{y}_k=(3+2\sqrt{2}{)}^k$$

则：

KaTeX parse error: No such environment: eqnarray* at position 7: \begin{̲e̲q̲n̲a̲r̲r̲a̲y̲*̲}̲x_{k+1} + y_{k…

因等式两端整数项对应相等，根号项对应相等，故：

$$\{\begin{array}{llcc}{x}_{k+1}=& 3x_k& +& 4y_k\\ y_{k+1}=& 2x_k& +& 3y_k\end{array}$$

现在考虑如何推出$x,y$单独的通项式：

KaTeX parse error: No such environment: eqnarray* at position 7: \begin{̲e̲q̲n̲a̲r̲r̲a̲y̲*̲}̲x_{k+2} &=&3x…

同理可得：

$$y_{k+2}=6y_{k+1}-y_k$$

---

因为$y=n$，故我们可以令：
$n_0=0,n_1=2$

则：

$$n_k=6n_{k-1}-n_{k-2}$$

故可以推出在$m=2$的情况下所有的$n$

---

$m=3$时，同理可得$n$的通项：
$n_0=0,n_1=6$

$$n_k=14n_{k-1}+n_{k-2}$$

---

故此题只需要递推出最小的不小于给定的$m$值的$n$即可。

代码：

m = (int)(input())
n0, n1 = 0, 2
while n1 < m:
    n0, n1 = n1, 6*n1 - n0
m0, m1 = 0, 6
while m1 < m:
    m0, m1 = m1, 14*m1 - m0
print(min(n1, m1))

---

拓展1：如何更简便地求出Pell方程某一变量的通项公式？

假设已知Pell方程：

$$x^2-d{y}^2=1$$

如果我们已经得到其最小正整数解($x_0,y_0$)
则对于第$k$组解有：

$$x_k+\sqrt{d}{y}_k=(x_0+\sqrt{d}{y}_0{)}^k$$

故可以推出：

KaTeX parse error: No such environment: eqnarray* at position 7: \begin{̲e̲q̲n̲a̲r̲r̲a̲y̲*̲}̲x_{k+1} + \sqr…

由等式两边对应系数相等，得出：

$$\{\begin{array}{llcc}{x}_{k+1}=& x_0{x}_k& +& d{y}_0{y}_k\\ y_{k+1}=& y_0{x}_k& +& x_0{y}_k\end{array}$$

故通项式为：

KaTeX parse error: No such environment: eqnarray* at position 7: \begin{̲e̲q̲n̲a̲r̲r̲a̲y̲*̲}̲x_{k+2} &=&x_…

同理：

$$y_{k+2}=2x_0{y}_{k+1}-y_k$$

故递推式前的系数分别是$(2x_0,-1)$

---

拓展2：对于非标准形式的Pell方程 $x^2-d{y}^2=T$ $(T̸=1)$，我们如何求出其通解呢？

(待更新)