多项式除法的C语言实现

问题的引出

在该博客中我将着手解决这样一个问题，在求解模素数$p$的同余式：
$$f(x)=a_n{x}^n+\cdots +a_{1}x+a_0\equiv 0\left(modp\right)$$
时，其中$a_n$恒不为零. 通常需要借助费马小定理和多项式的欧几里得除法将上述一个复杂的问题转化为一个次数更低的模素数同余式。

具体来说，由费马小定理可知，多项式$x^p-x\left(modp\right)$对任何整数取值为零.进一步，利用多项式欧几里得除法可知存在整系数多项式$q(x)$,$r(x)$使得：
$$f(x)=q(x)(x^p-x)+r(x),deg(r(x))<deg(g(x)).$$
从而容易将上述同余式：
$$f(x)=a_n{x}^n+\cdots +a_{1}x+a_0\equiv 0\left(modp\right)$$
降次为同余式：
$$r(x)\equiv 0\left(modp\right)$$
不难证明上述两个模$p$同余式等价.这时我便将初始的问题转化为一个次数更低的模$p$同余式，进而方便原复杂问题的进一步求解。在这个化简过程中，起到核心作用的便是多项式的欧几里得除法，所以为了最终成功编程实现模素数$p$任意次同余式的求解，我需要首先编程实现整系数的多项式欧几里得除法.

多项式欧几里得除法

引理（多项式欧几里得除法）：
设
$$f(x)=a_n{x}^n+\cdots +a_{1}x+a_0$$
为$n$次整系数多项式，
$$g(x)=x^n+\cdots +b_{1}x+b_0$$
为$m\ge 1$次首一整系数多项式，则存在整系数多项式 q ( x ) q(x)