大衍求一术

大衍求一术，又名求一术，是中国数学史中通常被用来泛指南宋数学家秦九韶发明的求解中国剩余定理的历史算法（不是中国剩余定理的现代算法）。秦九韶原来的《数书九章》求解一次同余式组的算法的总称叫做大衍总数术，或大衍术、而其中一个计算乘率的子程序，才是大衍求一术。

大衍术是秦九韶最得意的创作，特放在《数书九章》之首。欧洲直到18世纪德国数学家高斯，才有相类的结果。秦九韶大衍术领先世界五百余年。

历史

一次同余式组问题，最早见于《孙子算经》卷下第二十六问： Template:Quote

宋周密《志雅堂杂钞》述鬼谷算，又名隔墙算。

杨辉《继古摘奇算法》记载剪管术解孙子物不知数的算法：

Template:Quote Template:Quote 明严恭《通源算法》：^[1]Template:Quote

明朝数学家程大位有《孙子歌》如下：^[2] Template:Quote

清代学者张敦仁在《求一算术》中提出大衍求一术源自《孙子算经》物不知数，后世学者，多从其说。但近年李俨、钱宝琮等学者提出大衍求一术很可能源自西汉《三统历》中计算上元积年的“通其率”近似法，其后《古四分历》和《乾象历》沿用此法。^[3]^[4]

事实上，秦九韶在《数书九章序》中就写道：“独大衍法不载《九章》，未有能推之者，历家推演法颇用之”，“历家虽用，其用不知，小试经世，姑推所为，述大衍第一。”

大衍总数术

秦九韶大衍总数术原载《数书九章》第一卷上大衍类《蓍卦发微》： Template:Quote

程序如下：^[5]^[6]^[7]。

1。将问数（有关问题中的数字）[

m_{1}

，

m_{2}

……

m_{n}

]，分为四大类

元数：整数
收数：带小数的有理数。
通数：分数
复数：10的倍数。

2.如果问数[

m_{1}

，

m_{2}

……

m_{n}

]是收数，则以10的倍数乘之化为元数列，按元数计算。

3：如果问数[

m_{1}

，

m_{2}

……

m_{n}

]包含分数，则化为元数列，按元数计算。

4。当问数[

m_{1}

，

m_{2}

……

m_{n}

]是整数

4.1:将问数成对以其最大公约数约化，约奇数不约偶数；如二数都是偶数，约化其中之一。如二数中一个数是10的倍数，另一数的个位是5，则约化偶数不约化奇数。最后得到一组

定母：[ $m'_{1}$ ， $m'_{2}$ …… $m'_{n}$ ]。

定母之中避免过多1，得到一个“1”即可。各定母为相应问数的因子；定母两两互为质数。定母的乘积称为衍母，是问数的最小公倍数。

衍母 v=∏

m'_{i}

; i=1…n;

4.3：将各定母分别除衍母，得到一组衍数：[

M_{1}

，

M_{2}

……

M_{n}

]

其中衍数

M_{i} =

v/

m'_{i}

4.4 用定母

m'_{i}

约化衍数

M_{i}

(i=1……n)；余数

d_{i}

称为奇数。

4.4 当奇数

d_{i}

=1时，乘率

x_{i}

=1。

当奇数 $d_{i}$ ≠1时，从定母 $m'_{i}$ 和奇数 $d_{i}$ ，用大衍求一术计算乘率 $x_{i}$ ；（i=1……n)。

4.5 衍數乘以乘率，稱為用數。

4.6 各餘數乘以相應的用數，得出各總。

4.7 各總的總和，稱為總數。

4.8 把總數除以衍母所得的餘數，便是所求數。

例一：《孙子算经》“物不知数”，

問数：3，5，7；因两两互素，也即定母。
定母：3，5，7
衍母=3*5*7=105
衍数：35，21，15
奇数：2，1，1
乘率：2，1，1
用數：70，21，15
餘數：2，3，2
各總：140，63，30
總數=140+63+30=233
所求數=23

例二

問数：108，57，75，40
定母：27，19，25，8
衍母 =27*19*25*8=102600
衍数：3800，5400，4104，12825
奇数：20，4，4，1
乘率：23，5，19，1