首页 > 其他分享 >关于多项式方程所在剩余系的余数循环周期的猜想与推导

关于多项式方程所在剩余系的余数循环周期的猜想与推导

时间:2023-03-14 15:22:06浏览次数:41  
标签:right limits 推导 cdot 多项式 sum alpha 余数 left

前言  

  由于解决这个问题的做法完全是我个人主观完成的,因此可能会存在很严重的错误,如果发现任何的问题与错误请在评论区进行提出指正。还有这个问题在数论中应该有个专业的名词,如果有人知道的话也麻烦在评论区留言告知。

 

正文

  本文尝试解决的问题为:对于 $\forall x \in \mathbb{Z}$ 多项式方程 $f(x) = \sum\limits_{i=0}^{n}{a_i \cdot x^i}$ 在模 $m$ 意义下所对应的剩余系中,寻找一个非最小的余数循环周期 $T$ 使得

$$
\sum\limits_{i=0}^{n}{a_i \cdot (x+T)^i} \equiv \sum\limits_{i=0}^{n}{a_i \cdot x^i} \pmod{m}
$$

  为什么是非最小呢?其实直接取 $T=m$ 就满足上面那个同余方程了。实际上如果要求最小的 $T$ 会涉及到高次同余式的求解,我现在还不会qwq。所以就退而求次,求一个尽可能最小的周期 $T$,同样可以满足 $f(x+T) \equiv f(x) \pmod{m}$。

  比如要找到一个最小的 $T$ 满足高次同余式 $\sum\limits_{i=0}^{n}{a_i \cdot T^i} \equiv 0 \pmod{m}$。正规解法的话我还没学,这里只是记录我想到能得到尽可能小的 $T$ 的做法。就是让每一项的 $T$ 都满足 $a_i \cdot T^i \equiv 0 \pmod{m}$,然后对 $n$ 组方程的 $T$ 取最小公倍数。

  为了方便理解,我们先从数据规模较小的 $n$ 开始模拟,然后再推到更一般的情况。

  当 $n=1$,此时 $f(x) = a_0 + a_1 \cdot x$,对于 $\forall x \in \mathbb{Z}$ 有

$$
\begin{align*}
f(x+T) &= a_0 + a_1 \cdot (x+T) &\pmod{m} \\
&\equiv a_0 + a_1 \cdot x + a_1 \cdot T &\pmod{m} \\
&\equiv f(x) + a_1 \cdot T &\pmod{m}
\end{align*}
$$

  为了使得 $f(x+T) \equiv f(x) \pmod{m}$,很明显我们应该让 $a_1 \cdot T \equiv 0 \pmod{m}$。

  即有 $m \mid a_1 \cdot T$,等价于 $\left. \dfrac{m}{\gcd \left( m, a_1 \right)} \ \middle| \ T \right.$。

  这意味着 $T$ 应该是 $\dfrac{m}{\gcd \left( m, a_1 \right)}$ 的倍数,因此可以令 $T = \dfrac{m}{\gcd \left( m, a_1 \right)}$。

  当 $n=2$,此时 $f(x) = a_0 + a_1 \cdot x + a_2 \cdot x^2$,对于 $\forall x \in \mathbb{Z}$ 有

$$
\begin{align*}
f(x+T) &= a_0 + a_1 \cdot (x+T) + a_2 \cdot (x+T)^2 &\pmod{m} \\
&\equiv a_0 + a_1 \cdot x + a_2 \cdot x^2 + (2 \cdot a_2 \cdot x + a_1) \cdot T + a_2 \cdot T^2 &\pmod{m} \\
&\equiv f(x) + (2 \cdot a_2 \cdot x + a_1) \cdot T + a_2 \cdot T^2 &\pmod{m}
\end{align*}
$$

  此时应该让 $(2 \cdot a_2 \cdot x + a_1) \cdot T + a_2 \cdot T^2 \equiv 0 \pmod{m}$。

  注意到由于是对于 $\forall x \in \mathbb{Z}$ 上面的同余方程都要成立,因此我们在求解 $T$ 时不应该考虑 $x$ 的取值,即可以把 $x$ 从中去除,因此要求解的同余方程就变成了 $(2 \cdot a_2 \cdot x + a_1) \cdot T + a_2 \cdot T^2 \equiv 0 \pmod{m}$。

  令

$$
\begin{cases}
(2 \cdot a_2 \cdot x + a_1) \cdot T &\equiv 0 \pmod{m} \\
a_2 \cdot T^2 &\equiv 0 \pmod{m}
\end{cases}
$$

$$
\begin{cases}
m \mid (2 \cdot a_2 \cdot x + a_1) \cdot T \\
m \mid a_2 \cdot T^2
\end{cases}
$$

$$
\begin{cases}
\dfrac{m}{\gcd \left( m, 2 \cdot a_2 \cdot x + a_1 \right)} &\Bigm| T \\
\dfrac{m}{\gcd \left( m, a_2 \right)} &\Bigm| T^2
\end{cases}
$$

  我们当然可以直接令 $T = \mathop{lcm} \left[ \dfrac{m}{\gcd \left( m, 2 \cdot a_2 \cdot x + a_1 \right)}, \ \dfrac{m}{\gcd \left( m, a_2 \right)} \right]$,但为了使得 $T$ 尽可能的小,对 $\dfrac{m}{\gcd \left( m, a_2 \right)}$ 进行质因数分解,得到
$$
\dfrac{m}{\gcd \left( m, a_2 \right)} = P_{2,1}^{\alpha_{2,1}} \cdot P_{2,2}^{\alpha_{2,2}} \cdots P_{2,k_2}^{\alpha_{2,k_2}}
$$

  可以发现,对于 $\dfrac{m}{\gcd \left( m, a_2 \right)} \Bigm| T^2$ 的情况,当 $T = P_{2,1}^{\left\lceil\frac{\alpha_{2,1}}{2}\right\rceil} \cdot P_{2,2}^{\left\lceil\frac{\alpha_{2,2}}{2}\right\rceil} \cdots P_{2,k_2}^{\left\lceil\frac{\alpha_{2,k_2}}{2}\right\rceil}$,可以取到满足条件的最小的 $T$。

  为了统一起来,对于 $\dfrac{m}{\gcd \left( m, 2 \cdot a_2 \cdot x + a_1 \right)} \Bigm| T$ 的情况,先对 $\dfrac{m}{\gcd \left( m, 2 \cdot a_2 \cdot x + a_1 \right)}$ 进行质因数分解,得到
$$
\dfrac{m}{\gcd \left( m, 2 \cdot a_2 \cdot x + a_1 \right)} = P_{1,1}^{\alpha_{1,1}} \cdot P_{1,2}^{\alpha_{1,2}} \cdots P_{1,k_1}^{\alpha_{1,k_1}}
$$

  此时令 $T = P_{1,1}^{\left\lceil\frac{\alpha_{1,1}}{1}\right\rceil} \cdot P_{1,2}^{\left\lceil\frac{\alpha_{1,2}}{1}\right\rceil} \cdots P_{1,k_1}^{\left\lceil\frac{\alpha_{1,k_1}}{1}\right\rceil}$,就等价于 $T = \dfrac{m}{\gcd \left( m, 2 \cdot a_2 \cdot x + a_1 \right)}$。

  然后对两者取最小公倍数,就有
$$
T = \mathop{lcm} \left[ P_{1,1}^{\left\lceil\frac{\alpha_{1,1}}{1}\right\rceil} \cdot P_{1,2}^{\left\lceil\frac{\alpha_{1,2}}{1}\right\rceil} \cdots P_{1,k_1}^{\left\lceil\frac{\alpha_{1,k_1}}{1}\right\rceil} \ , \ \ P_{2,1}^{\left\lceil\frac{\alpha_{2,1}}{2}\right\rceil} \cdot P_{2,2}^{\left\lceil\frac{\alpha_{2,2}}{2}\right\rceil} \cdots P_{2,k_2}^{\left\lceil\frac{\alpha_{2,k_2}}{2}\right\rceil} \right]
$$

  这时可以考虑更加一般的情况了。

  先对式子 $\sum\limits_{i=0}^{n}{a_i \cdot (x + T)^i}$ 进行等价变形:

\begin{align*}
&\sum\limits_{i=0}^{n}{a_i \cdot (x+T)^i} \\\\
&= a_0 + a_1 \cdot (x+T) + a_2 \cdot (x+T)^2 + \cdots + a_i \cdot (x+T)^i + \cdots + a_n \cdot (x+T)^n \\\\
&= a_0 + a_1 \cdot \sum\limits_{j=0}^{1}{C_{1}^{j} \cdot x^{1-j} \cdot T^j} + a_2 \cdot \sum\limits_{j=0}^{2}{C_{2}^{j} \cdot x^{2-j} \cdot T^j} + \cdots + a_i \cdot \sum\limits_{j=0}^{i}{C_{i}^{j} \cdot x^{i-j} \cdot T^j} + \cdots + a_n \cdot \sum\limits_{j=0}^{n}{C_{n}^{j} \cdot x^{n-j} \cdot T^j} \\\\
&= \left( a_0 + a_1 \cdot x + a_2 \cdot x^2 + \cdots + a_i \cdot x^i + \cdots + a_n \cdot x^n \right) \\\\
& \ \ \ \ + \left( a_1 \cdot \sum\limits_{j=1}^{1}{C_{1}^{j} \cdot x^{1-j} \cdot T^j} + a_2 \cdot \sum\limits_{j=1}^{2}{C_{2}^{j} \cdot x^{2-j} \cdot T^j} + \cdots + a_i \cdot \sum\limits_{j=1}^{i}{C_{i}^{j} \cdot x^{i-j} \cdot T^j} + \cdots + a_n \cdot \sum\limits_{j=1}^{n}{C_{n}^{j} \cdot x^{n-j} \cdot T^j} \right) \\\\
&= \sum\limits_{i=0}^{n}{a_i \cdot x^i} + \sum\limits_{i=1}^{n} \left( {a_i \cdot \sum\limits_{j=1}^{i}{C_{i}^{j} \cdot x^{i-j} \cdot T^j}} \right)
\end{align*}

  即在模 $m$ 意义下要使得
$$
\begin{align*}
& \ \ \ \ \ \sum\limits_{i=0}^{n}{a_i \cdot (x+T)^i} \\
&= \sum\limits_{i=0}^{n}{a_i \cdot x^i} + \sum\limits_{i=1}^{n} \left( {a_i \cdot \sum\limits_{j=1}^{i}{C_{i}^{j} \cdot x^{i-j} \cdot T^j}} \right) \\
&\equiv \sum\limits_{i=0}^{n}{a_i \cdot x^i} \pmod{m}
\end{align*}
$$

  即对于任意的 $\forall x \in \mathbb{Z}$,要求 $T$ 使得

$$
\sum\limits_{i=1}^{n} \left( {a_i \cdot \sum\limits_{j=1}^{i}{C_{i}^{j} \cdot T^j}} \right) \equiv 0 \pmod{m}
$$

  对上式进行等价变形就会得到
$$
\sum\limits_{j=1}^{n}{\left( \sum\limits_{i=j}^{n}{C_{i}^{j} \cdot a_i} \right) \cdot T^j} \equiv 0 \pmod{m}
$$

  每一项在模 $m$ 意义下余 $0$,得到
$$
\begin{cases}
&\left. \dfrac{m}{\gcd \left( m, \ \sum\limits_{i=1}^{n}{C_{i}^{1} \cdot a_i} \right)} \ \middle| \ T \right. \\\\
&\left. \dfrac{m}{\gcd \left( m, \ \sum\limits_{i=2}^{n}{C_{i}^{2} \cdot a_i} \right)} \ \middle| \ T^2 \right. \\\\
& \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \vdots \\
&\left. \dfrac{m}{\gcd \left( m, \ \sum\limits_{i=j}^{n}{C_{i}^{j} \cdot a_i} \right)} \ \middle| \ T^j \right. \\\\
& \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \vdots \\
&\left. \dfrac{m}{\gcd \left( m, \ \sum\limits_{i=n}^{n}{C_{i}^{n} \cdot a_i} \right)} \ \middle| \ T^n \right. \\\\
\end{cases}
$$

  对每一项的除数进行质因数分解,得到

$$
\dfrac{m}{\gcd \left( m, \ \sum\limits_{i=j}^{n}{C_{i}^{j} \cdot a_i} \right)} = \sum\limits_{i=1}^{k_j}{P_{j,i}^{\alpha_{j,i}}}
$$

  那么要求的 $T$ 就是
$$
T = \mathop{lcm}\limits_{1 \leq j \leq n} \left[ {\sum\limits_{i=1}^{k_j}{P_{j,i}^{\left\lceil \frac{\alpha_{j,i}}{j} \right\rceil}}} \right]
$$

  以上。

  感觉很笨是吧,明明有正解(解高次同余式)了,上面的方法不仅复杂,甚至都不是求出余数循环的最小周期 $T$。这篇记录就当作图一乐吧,毕竟很少有机会取思考数学问题。

  当然还有多项式真假分式的情况,这个就太难了没什么思路。

标签:right,limits,推导,cdot,多项式,sum,alpha,余数,left
From: https://www.cnblogs.com/onlyblues/p/17214996.html

相关文章

  • 求商和余数
    需求: 给定两个整数,被除数和除数(都是正数,且不超过int的范围)将两数相除,要求不使用乘法、除法和%运算符。  得到商和余数。分析: 被除数/除数=......
  • 整环上的多变量多项式不是 PID 的证明
    整环上的多变量多项式不是PID.分析与证明:设(R,+,·)是一个整环,需证若n≥2,则R[x1,...,xn]不是PID.为方便起见,使用多重变量简记符号x=x1, ...,xn,以及多......
  • AcWing 199. 余数之和 题解
    做了一下午……题解都看不懂,最后自己比比划划弄懂了。题意:给出\(n,k\),求\(\sum\limits_{i=1}^nk\modi\)。首先取模形式十分不好处理,所以我们可以根据取模运算定义做......
  • 倍数问题(同余定理,对余数的进一步理解)
    题目描述众所周知,小葱同学擅长计算,尤其擅长计算一个数是否是另外一个数的倍数。但小葱只擅长两个数的情况,当有很多个数之后就会比较苦恼。现在小葱给了你n个数,希望你从这......
  • 多项式全家桶
    有时间再修\(\operatorname{FFT}\)P3803【模板】多项式乘法(FFT)点击查看代码#include<bits/stdc++.h>#definecsconst#defineilinline#definefo(i,j,k)for(in......
  • Matlab 多项式的根求解
    ✅作者简介:热爱科研的算法开发者,Python、Matlab项目可交流、沟通、学习。......
  • 基于概率论的MATLAB仿真,内容包括非共轭条件下的后验概率的推导,共轭条件下的非完备集
    1.算法描述1.1先验概率的推导        根据贝叶斯概率论可知,某一事件的后验概率可以根据先验概率来获得,因此,这里首先对事件的先验概率分布进行理论的推导。假设测......
  • 基于概率论的MATLAB仿真,内容包括非共轭条件下的后验概率的推导,共轭条件下的非完备集
    1.算法描述1.1先验概率的推导根据贝叶斯概率论可知,某一事件的后验概率可以根据先验概率来获得,因此,这里首先对事件的先验概率分布进行理论的推导。假设测量的腐蚀数据服从g......
  • 多项式小记
    多项式牛顿迭代对于\(G(f(x))=0\),求解\(f\pmod{x^n}\)$x^{\left\lceil\frac{n}{2}\right\rceil}$意义下的解\(f_{0}\left(x\right)\),要求模\(x^{n}\)意义下的解......
  • pearson总体相关系数到样本相关系数推导过程
    相关系数是根据样本数据计算的度量两个变量之间线性关系强的统计量,若相关关系是根据总体全部数据计算的,成为总体相关系数,记为\(\rho\);若是根据样本数据计算的,则称为样本相......