- 在线时间
- 17 小时
- 最后登录
- 2016-8-29
- 注册时间
- 2009-1-19
- 听众数
- 3
- 收听数
- 0
- 能力
- 0 分
- 体力
- 423 点
- 威望
- 0 点
- 阅读权限
- 30
- 积分
- 178
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 100
- 主题
- 20
- 精华
- 0
- 分享
- 0
- 好友
- 8
升级   39% TA的每日心情 | 开心 2016-8-29 17:02 |
|---|
签到天数: 18 天 [LV.4]偶尔看看III
 |
中国剩余定理新解/ H% b X" ~+ l! q6 ~
为什么要编写本文?因为,中国剩余定理属于算术问题,算术的特征是准确性和唯一性。编写的目的是为了让人们确认它的准确性和唯一性,探索最简单的计算方法,便于人们推广应用。 h1 b! D9 k, J ?% T( U' P
本文的主要内容:1、原题,2、基本原理,3、计算方法比较,4、中国剩余定理的扩展,5、趣谈中国剩余定理。
! S1 r7 b" y0 ?5 P! `- E一、原题5 Q& \ ^' o7 a3 \& |& U( m; t
中国剩余定理,又叫孙子定理,原文是:6 p; m4 P% q7 W4 C$ @7 U
“《孙子算经》中的题目:有物不知其数,三个一数余二,五个一数余三,七个一数又余二,问该物总数几何? 8 R. a" m0 ^$ F" X6 a; v$ q/ I# j
《孙子算经》中的解法:三三数之,取数七十,与余数二相乘;五五数之,取数二十一,与余数三相乘;七七数之,取数十五,与余数二相乘。将诸乘积相加,然后减去一百零五的倍数。”对于这种解法(以下简称前人的计算方法)。
3 ?8 g2 D, o; s5 Q7 G那么,前人计算方法的原理是什么?只有从计算原理,才知道该计算方法是否正确,所求之数是否存在,存在于什么范围之内?计算方法最简单、最方便的是什么?
# o4 X1 W- X+ V( J中国剩余定理,是从一种现象开始,寻找解决现象的方法。本人认为应该从基础原理开始,便于人们掌握和应用:
6 D$ ?9 o0 T/ R+ t' M& E9 w9 T二、基本原理
) V# P( D* N9 a) |" f; I+ @3 \1、在整数中,用不同的器具,如2,3,5,7,11,13,17,…,N(素数)量一个固定的数,每一个量器都有一个固定的余数。
% k- m2 s0 Q: r8 j如某数为53,用2,3,5,7进行衡量,有53/2余1,53/3余2,53/5余3,53/7余4,即,53除以2,3,5,7都有一个固定的余数,也只有一个固定的余数。
+ N# [6 b7 r+ J, t& H( z" C6 h2、用2,3,5,7,11,13,17,…,N作为器具,每一个量器取一个固定的余数,在2*3*5*7*11*13*17*…*N之内(也就是在除数的最小公倍数之内),满足这些余数条件的数只有一个(具体分析见后)。8 ?( H! n/ _) x1 r8 O; n
如用2,3,5,7量M,当M/2余1,M/3余2,M/5余4,M/7余5,得M=89,89必然存在于2*3*5*7=210之内,在每210个连续数之内满足M/2余1,M/3余2,M/5余4,M/7余5的数只有一个数,可以用89+210N表示。
( o: B" s5 X& l" c: V4 q% K必然存在的理由依据:用2,3,5,7,11,13,17,…,N作为器具时。自然数除以2可以取2种余数,余0,余1;除以3可以取3种余数,余0,余1,余2;除以5可以取5种余数,余0,余1,余2,余3,余4;…,除以N可以取N种余数,余0,余1,余2,余3,余4,…,余N-1。不同的余数排列组合共有2*3*5*7*11*13*17*…*N个,正好对应2*3*5*7*11*13*17*…*N之内的数,每一组不同的余数对应这之内的一个数。为什么存在一一对应呢?请看下面的具体计算。3 E" p4 P+ _* B6 ^$ }
3、除以A余C的数,为等差数列C+AN中的数,也只有C+AN中的数满足除以A余C。, z* i ?- x9 c k, ]6 s V1 t8 h
如除以7余3的数,为等差数列3+7N中的数,也只有3+7N中的数满足除以7余3。
* E6 Y( N& L2 w0 d! u, d4 L' e. _4、除以A,B,…,C都余K,称为同余项。满足除以A,B,…,C都余K的数为等差数列K+(A*B*…*C)N,也只有K+(A*B*…*C)N中的数满足分别除以A,B,…,C都余K。
# b; z ]* z+ t8 ~5 F& I+ g4 Q三、计算方法比较
; G* T/ T+ X% l计算方法,指求除以素数2,3,5,7,11,13,17,…,N各取一个固定的余数,在2*3*5*7*11*13*17*…*N之中,寻找满足这些条件的这一个数的方法。( \: i" Q4 `' l V
上面虽然说的是除以素因子2,3,5,7,11,…,N的依次排列,换为不同的素数A,B,C,D,E,…,Z这种说法仍然成立;排列顺序不一定依次,可以按自己的须要,认为怎样计算方便就怎样进行排列。
; X' i& X5 ~8 ~2 O. F Y: c8 U解决这一问题,前人有两种方法,本人有一种计算方法,下面以同一个题进行比较:
: ~8 {+ s3 j1 U方法一、/ q/ r" z, O0 J) c4 g
前人的计算方法为什么成立?其解题思路为:
" E+ W1 g/ t. \; R6 ~$ w6 C- a令某数为M,令素数为A,B,C,D,…,Z,已知M/A余a,M/B余b,M/C余c,M/D余d,…,M/Z余z。求M=?
7 ?& a. j4 r% ^4 l+ J9 d, V6 L因为A,B,C,D,…,Z为不同的素数,故,B*C*D*…*Z不可能被A整除,有(B*C*D*…*Z)N中取A个连续项,这A个连续项分别除以A的余数必然存在0,1,2,3,…,A-1,每一个余数都不缺少。所以,从这A个连续项中必然能寻找到除以A余1的数。再用除以A余1的这个数*a,其积必然除以A余a,这个除以A余a的数,为能够被素数B*C*D*…*Z整除的数,为第一个数;/ K. T2 ^' }7 C9 s' e0 t
再按同样的道理,从A*C*D*…*Z的倍数中寻找除以B余b的数,该数具备被素数A,C,D,…,Z整除的特性,为第二个数;% [4 c8 P3 p! T9 M
因为,第一个数除以A余a,第二个数能被素数A,C,D,…,Z整除,即能被A整除,所以,第一个数+第二个数之和,仍然保持除以A余a;
: @; W& N7 S2 P6 S t同理,第二个数除以B余b,因第一个数能被B整除,所以,第二个数+第一个数之和,仍然保持除以B余b。即,第一个数+第二个数之和,为满足除以A余a,除以B余b,并且,能被C,D,…,Z整除。- w, Q( @' ~ V. `. C
按同样的方法,从A*B*D*…*Z的倍数中寻找除以C余c的数,该数具备能被A,B,D,…,Z整除的特性,为第三个数;0 v8 }" L8 g' `8 X" |( w5 P
因第一个数和第二个数,都能被C整除,故第三个数+(第一个数+第二个数),仍然保持除以C余c;又因第三个数能被A、B整除,所以,(第一个数+第二个数)+第三个数之和,仍然保持除以A余a,除以B余b。即第一个数+第二个数+第三个数之和,为除以A余a,除以B余b,除以C余c的数,而且,能被素数D,…,Z整除的数;; E! b2 @# m- Y5 E) n& f
依此类推,按上面的方法寻找到除以各素数余数的数之总和,为满足除以各素数余数的条件的数。总和再减去能被这几个素数共同整除的数(A*B*C*D*…*Z)N后,其差仍然保持除以各素数余数的条件的数。由此构成前人对中国剩余定理的解法。: T6 G( G3 X5 _ R, w4 `: F4 L# R
例:某数为M,M/3余2,M/5余4,M/7余3,M/11余6,求M=? `' b* S' G5 f0 z! n5 K- J _
解:
: Q9 B% C' Q8 [ u4 _1、5*7*11=385,对385N取3项:385,770,1155,因385/3余1,有385*2=770,为除以3余2的数;" j% ~3 @! R' B/ o( ^. T! D) d
2、3*7*11=231,对231N取5项:231,462,693,924,1155,因231/5余1,有231*4=924,为除以5余4的数;( e3 T: Q( r Q6 A0 j) Q5 O
3、3*5*11=165,对165N取7项:165,330,495,660,825,990,1155,有330/7余1,得330*3=990,为除以7余3的数;
0 c/ @$ i- a4 H+ Y! u4、3*5*7=105,对105N取11项:105,210,315,420,525,630,735,840,945,1050,1155,因210/11余1,有210*6=1260,为除以11余5的数;' f, G, f+ g4 X8 s
5、770+924+990+1260=3944,因3*5*7*11=1155,3944/1155余479,即479为满足这些条件的数。
; Q. v* o6 A/ M, @其实,前人的计算方法,完全可以简化一个步骤,直接寻找满足余数条件的数,省略一步乘以余数得寻找之数:8 @7 G' p) ?# G; F( c
1、5*7*11=385,对385N取3项:385,770,1155,因770/3余2;
6 T! S6 o* Q6 N+ T2、3*7*11=231,对231N取5项:231,462,693,924,1155,因924/5余4;
* L& q. d& d7 c3、3*5*11=165,对165N取7项:165,330,495,660,825,990,1155,因990/7余3;0 g% H b) I# W! g$ i8 {4 g
4、3*5*7=105,对105N取11项:105,210,315,420,525,630,735,840,945,1050,1155,因105/11余6。, \% n. t5 y! {
5、770+924+990+105=2789,因3*5*7*11=1155,2789/1155余479,即479为满足这些条件的数。4 J+ N. A/ q" P7 t$ Q _; y6 x
方法二、" _7 C' Y% y# i: p9 t7 q2 i4 [
除以3余2为等差数列2+3N有:2,5,8,11,14,…,479,…,1154,要取5*7*11=385项,# U- z" W: C( x; A( J# |5 o
除以5余4为等差数列4+5N有:4,9,14,19,…,479,…,1154,要取3*7*11=231项。) w* ~% g0 s. d0 q
除以7余3为等差数列3+7N有:3,10,17,24,…,479,…,1151,要取3*5*11=165项,; X8 m/ z* r( H
除以11余6为等差数列6+11N有:5,16,27,38,…,479,…,1149,要取3*5*7=105项。
( ]# A0 X) u' w7 s4 T因4个等差数列在1155内,只有479都同时存在,所以,只有479为满足这些条件的数。
) l" `4 c( a7 q7 `方法二对于计算较大的剩余数,很不适用,不可取。
, G% z* {. y8 V5 V# @$ l0 w! e本人的方法:
) Y: q: i( f- O& c( x该计算方法,可以直观地说明剩余数在范围内存在的必然性和唯一性。我们边计算,边说明。7 \+ W" V! Y+ f, o, m
(1),初步计算:+ p' t; ]- a( h: U
对于素数11来说,1,2,3,4,5,6,7,8,9,10,11,在这11个连续数中,只有6除以11余6,满足除以11余6为等差数列6+11N;
& t/ K/ `" C( L. o0 Y# D素数7,对6+11N取7项,6,17,28,39,50,61,72,这7个数除以7余1,2,3,4,5,6,0都存在,这就是存在的必然性;只有17/7余3,这就是唯一性。因11*7=77,即在77之内满足除以11余6,除以7余3的数只有17,满足这两个条件的数为17+77N,表明每77个连续数之内只有一个数满足这两个条件。
( ~; O/ q% {- T8 x& L# t素数5,对17+77N取5项:17,94,171,248,325,这5个数除以5余1,2,3,4,0都存在,这就是存在的必然性;只有94除以5余4,这就是唯一性。因77*5=385,即在385之内满足除以11余6,除以7余3,除以5余4的数只有94,满足这三个条件的数为94+385N,表明每385个连续数中只有一个数满足这三个条件。
' {) b& N. O+ ]' C素数3,对94+385N取3项:94,479,864,这3个数除以3余1,2,0都存在,这就是存在的必然性;只有479除以3余2,这就是唯一性。因385*3=1155,满足这4个条件的数为479+1155N,表明每1155个连续数中只有一个数满足这四个条件,也必然有一个数满足这四个条件。. B; a) w, T* ?1 g1 e
(二),简化计算(该方法实用于计算大数):9 j1 b9 p; w; e' K: |% g- g8 h
1,素数11,除以11余6为等差数列6+11N;
! m7 d) g, D7 y# t$ _ W2,素数7,将等差数列6+11N的首项和公差同除以7的余数为6和4,组成新的等差数列6+4N,取7项之内必然有除以7余3的数,第一项6,第二项(6+4)-7=3,当出现与余数相同的数时,就不需要再往下计算了。代入原等差数列6+11*(2-1)=17,因11*7=77,得新的等差数列17+77N;
9 Z6 ^; ]6 |* w6 a$ r3,素数5,将17+77N同除以5变为2+2N,取5项之内:2,4,得第二项满足条件,代入有17+77=94,因77*5=385,得新的等差数列94+385N;
" H* X8 r9 L, ~, n5 k6 t) e4,素数3,将94+385N同除以3变为1+1N,取3项之内:1,2,得第二项满足条件,代入有94+385=497,因11*7*3*5=1155,即497+1155N数列的数都满足这4个条件。8 n, p4 I+ ]1 j! z- T' D% L: A3 S
本人的这种方法,是在3+5+7+11=26个数之内,寻找在1155之内满足除以这4个素数余数条件的数,这种方法也叫滚雪球的方法,容易理解、操作都比前人的简单、方便。8 G0 t4 e; C3 u% t* L8 c: R
) g* g$ p, o' L/ j方法二、同余的解法:) K9 D6 L( n. {7 ?( S/ f
例M/5余3,M/11余4,M/19余3,M/29余4,M/37余3,求M=?
" ^, b6 S8 `% x8 q0 e4 P& IM除以5,19,37同余3,满足除以5,19,37都余3为等差数列3+(5*19*37)N,即3+3515N;7 N& y5 I1 q2 K( J5 a* Q! c
M除以11,29同余4,满足除以11和29都余4为等差数列4+(11*29)N′,即4+319 N′。: \ I* P4 Z. L( i" c& ?
这里又出现三种解法:
( \9 l4 C1 t. a5 W' J% o6 C1、将3+3515N取319项,4+319 N′取3515项,共同有之数为此题的解;' _; M: Q& x2 ]
2、将3+3515N取319项,寻找除以319余4的数,为此题的解; H2 a; G# y5 {: ^
方法1和2比较烦锁(略)。# S. g3 L. \. C* r1 A
3、将3+3515N同除以29得等差数列3+6N取29项之内:3,9,15,21,27,4(当≥29时,减去29再算,下类似),得第6项除以29余4,代入原等差数列3+3515(6-1)=17578,因3515*29=101935,得等差数列17578+101935N。1 t# d& N' g3 v9 T; F* u
将等差数列17578+101935N同除以11得等差数列0+9N,取11项内:0,9,7,5,3,1,10,8,6,4,得第10项除以11余4,代入原等差数列17578+101935*(10-1)=934993。又因5*11*19*29*37=1121285,即在1121285之内只有M为934993时,满足M/5余3,M/11余4,M/19余3,M/29余4,M/37余3。等差数列934993+1121285N的数,都满足这些条件。
' X# b# H I9 J9 f四、中国剩余定理的扩展
& z" D7 G; [% @& {4 z6 Z4 C' W: \中国剩余定理在初始阶段,都是除以素数的余数。后来,不可避免地扩展到除以合数的余数问题。当涉及除以合数的余数后,就出现了有的题有解,有的题无解。这到底是为什么?0 y+ g) ? z, O4 k& y+ [1 Z
前面说的是除以不同素数的余数,那么,除以任意数的余数是否存在必然性和唯一性呢?
) B) Q! p5 A* g, h; Y) i* \这里的任意数,不包括0和1,任意数分为单合数、多合数、混合数、素数:
* s; R. t" [# \单合数,指由单个素因子组成的合数,如8,16等由单个素因子2组成;9,27,81等由单个素因子3组成。 ( g* ~0 c5 S6 t1 o! ~: n& Z
多合数指两个以上不同素因子组成的合数,如30=2*3*5。5 o1 L( V+ }- n, v, q
混合数指由单合数与多合数组成的合数,如840=8*3*5*7。
- i% p3 _7 }3 E4 n% f6 b例一,M/16余3,M/105余4,M/11余5,M/3余1,求M=?
, Z, d6 ]" ~) Q% l) M I1 Y; S首先审题:16为单个素因子2,105=3*5*7,还有11,3。这里只有一个重复素因子3,对于单个素因子3已经明确M/3余1;在多合数中M/105余4,用该余数4/3余1,没有矛盾,该题为可解之题。
) B3 D2 y0 W% c% d/ n9 H! R/ ^$ c抛开重复素数3后有:16*105*11=18480,表明在18480内该题有解,有唯一的解:- G5 r9 t' X6 c& s: d' b
满足除以105余4为等差数列4+105N,将4+105N化为4+9N取16项内:4,13,6,15,8,1,10,3,得第8项满足除以16余3,代入原等差数列4+105(8-1)=739,因105*16=1680,即在1680内同时满足除以105余4,除以16余3,只有739,得等差数列739+1680N;
, X- l6 q4 b" |' V7 F Q1 w+ c将该等差数列739+1680N化为2+8N,取11项:2,10,7,4,1,9,6,3,0,8,5,满足除以11余5为11项,代入的等差数列739+1680(11-1)=17539,因1680*11=18480,满足这些条件的数只有17539+18480N等差数列的数。, q1 t; O+ O0 f, M8 |3 M
例二,M/135余37,M/63余28,M/105余7,求M=?
y& l0 n8 c) ]& J0 o这里的135=3*3*3 *5,63=3*3*7,105=3*5*7。单合数有3的平方,3的立方;素数有3,5,7。计算最小公倍数时,单合数取素数N次方中最高的,即3*3*3=27,表明该题有除数27,5,7。有27*5*7=945,那么,满足这些条件的最小数在它们的公倍数945之内。
$ c I! a7 h L& [/ F/ _# r审题:37/3余1,28/3余1,7/3余1;37/5余2,7/5余2;28/7余0,7/7余0。其余数除以素因子的余数,都没有矛盾,表明该题有解。1 d* o# V9 z* H: r* ^- X' M; H
因M/135包括M/27和M/5,所以,该题变为:M/135余37,M/7余0。! J4 p, a+ g7 ^2 r( I: ]4 i
M/135余37为等差数列37+135N,化为2+2N取7项内:2,4,6,1,3,5,0,得第7项除以7余1,代入原等差数列37+135*(7-1)=847,即847+945N数列的数都满足这些条件。
% e4 ?( t0 C- }8 ?说明:) S/ m% }0 ?- B* g" L
1、 单个素因子组成的合数是不能拆分的,其除数取素数N次方中最高的;' v% n. m: c" I+ A/ p1 o0 y h7 X: s
2、 多个素因子组成的合数是可以拆分的,拆分方法:为总余数除以各素因子;- [, v% r& @: D; Y+ L5 g
3、 重复素因子,应先检验是否有矛盾,无矛盾是可解之题,有矛盾是无解之题;; b. U. n9 _; X! e# H+ t% N& f$ j
4、 重复素因子,重复合数,应当抛开重复后,再进行计算。
: m) \5 k% T$ I) A再举一个重复数的抛开,如M/105余19,M/165余139,求M=?
& {* N' g I5 @, N$ U! U8 {审题:105=3*5*7,165=3*5*11,两个合数中都有素因子3和5,因总余数19/3余1,139/3余1;19/5余4,139/5余4,余数除以素因子的余数无矛盾,该题有解。
" R% q" O6 D5 o4 F5 r. ]8 I选择抛开105中的3和5,剩余7,有19/7余5,原题变为M/7余5,M/165余139。
, @7 ^0 {/ \+ g: x9 h/ r由M/165余139为等差数列139+165N,先化为6+4N取7项内:6,3,0,4,1,5,第6项满足除以7余5,代入原等差数列139+165*(6-1)=964,因165*7=1155,得满足该题的解为964+1155N。4 J( l, _% \( J( r" T. y* D3 l8 {( L
中国剩余定理的结论:
/ }) S9 }1 X: j% s- s7 q令任意固定整数为M,当M/A余a,M/B余b,M/C余c,M/D余d,…,M/Z余z时,这里的A,B,C,D,…,Z为除数,除数为任意自然数(不包括0和1)时;余数a,b,c,d,z为自然整数时。
5 }- c$ {5 d" y 1、当题正确时,在这些除数的最小公倍数内有解,有唯一的解,每一个最小公倍数内都有唯一的解;当题错误时,在整个自然数范围内都无解。 2 h( y" t) |& h5 S+ r& b, q
2、当M在两个或两个以上的除数的最小公倍数内时,这两个或两个以上的除数和余数同时定位M在最小公倍数内的具体位置,也就是M的大小。8 H$ D+ K6 ^/ |+ R6 J* a3 D0 w
五、趣谈中国剩余定理6 G, B% [* p& U) v
我们以57为例,有57/2余1,57/3余0,57/5余2,57/7余1,57/11余2,57/13余5,57/17余6,57/19余0,57/23余11,57/29余28,57/31余26,57/37余20,57/41余16,57/43余14,57/47余10。: n6 {6 A- D6 i- \
这里有15个素数作为57的参照,那么,所有素数的乘积大于57许多;分别乘积,大于57的也有许多,素数如3*5*7=105,2*3*11=66,2*47=94都大于57。所有素数的乘积与部分素数的乘积,对于57的定位如何呢?
. Z& w3 t7 s) v& ?4 z0 z" }7 G8 @(一)部分素数
- i2 l: @" ?: w) n9 C& B: S1、按3*5*7=105,为:
: ]3 S5 c+ I1 s/ P) J4 t满足除以7余1为等差数列1+7N,取5项有1,8,15,22,29,满足除以5余2为22,因7*5=35,得等差数列22+35N,
" r( B' T& Q' k3 }对等差数列22+35N取3项为,22,57,92,满足除以3余0为57,因35*3=105,得等差数列57+105N满足除以3,5,7的余数条件;
; B: ~* \0 M5 M2、按2*3*11=66,为:4 T, `: M( M! C% ~' X% o, e0 t
满足除以11余2为等差数列2+11N,取3项有2,13,24,满足除以3余0为24,因11*3=33,得等差数列24+33N,取2项24,57,有57满足除以2余1,因33*2=66,得等差数列57+66N满足除以2,3,11的余数条件;
) n4 g% }) Z2 B; z/ U j; U1 G3、按2*47=94为:满足除以47余10为等差数列10+47N,取2项10,57,有57满足除以2余1,因2*47=94,得等差数列57+94N满足除以2和47的余数条件。
3 a* }$ j) Y5 `2 ?% U9 M* a(二)、全部素数
, l3 g/ R. x {前面计算了除以素因子3,5,7的余数后,得等差数列57+105N,- j; x6 g5 e( x! i5 P& r9 V I
再计算11时,将57+105N取11项,57,162,267,…,只有第一项满足除以11余2,因105*11=1155,得等差数列57+1155N,再计算13时,将该等差数列取13项,57,1212,2367,…,只有第一项满足除以13余5;……。! k$ {; t( z: I2 I
也就是说,在614889782588491410之内的任意一个数,锁定了它除以这些素因子的余数;反过来除以这些素因子的余数,也锁定了在614889782588491410之内的这个数。不论你先计算除以哪一个素因子的余数,后计算除以哪一个素因子的余数,结果都是一样的。( k0 m6 y: z! {- P. r6 ^$ k
从这里也可以看出,当几个素因子的乘积大于这个数时,从除以这几个素因子的余数,就可以准确地计算出这个数。
2 ]: H& Y8 V& @5 C9 Z: c; B3 A四川省三台县工商局:王志成。+ s1 O. A8 T1 Y% w& d+ s7 ^" R
|
zan
|