CSI讲义6-- 有趣的Mod算术运算-欧拉定理

前文再续,书接上一回,我们说到费尔马小定理,这里我们......

欧拉

Mod数为合数时的算术运算

同样的Python代码:

给定任意一个整数n,比如让n=11.
for i in range(1, n):    #i循环从1到n-1
for j in range(1, n):    #j循环从1到n-1
    print ((i * j) % n), #输出 i*j mod n
print                    #只是控制换行

当n=12,我们的程序输出:

1 2 3 4 5 6 7 8 9 10 11
2 4 6 8 10 0 2 4 6 8 10
3 6 9 0 3 6 9 0 3 6 9
4 8 0 4 8 0 4 8 0 4 8
5 10 3 8 1 6 11 4 9 2 7
6 0 6 0 6 0 6 0 6 0 6
7 2 9 4 11 6 1 8 3 10 5
8 4 0 8 4 0 8 4 0 8 4
9 6 3 0 9 6 3 0 9 6 3
10 8 6 4 2 0 10 8 6 4 2
11 10 9 8 7 6 5 4 3 2 1

是不是没有规律?呃,慢着,我看看......第1行怎么那么熟悉?慢着,似乎第5行的结果也很漂亮,因为结果就是1到11到一次重排。哦,原来第7行、第11行,结果都是1到11到一次重排。1、5、7、11都是什么数?不难看出,是与12互素的整数。似乎,规律出来了。我们是不是可以使用上一次课到方法得到:

如果i与n互素,那么i*1 mod n,i*2 mod n,i*3 mod n,...,i*(n-1) mod n[1..n-1]这些数的重排。所以:

    i^(n-1)*(n-1)! ≡ (n-1)! mod n

所以,**i^(n-1) ≡ 1 mod n** ?啊,一个新的“费尔马小定理”?!

慢着,有点不对......因为(n-1)!并不与n互素!!!所以,我们并不可以做两边的“消去”术。似乎有点失望。那为什么费尔马小定理可以那样做?因为n是素数,1到n-1之间到所有数都与n互素。

既然如此,那么,如果n为合数,1到n-1之间的那些数与n互素?它们的相乘是不是还是与n互素?

答案是,1到n-1之间的那些与n互素的数(不就是gcd(i, n)==1 ?)的相乘确实与n互素!

  作业:证明n个与N互素的整数相乘得到的整数与N互素。

现在,如果我们只考虑与n互素的数会怎么样?计算机程序就是帮助我们找规律的利器。为何不尝试用C语言写个程序看看?(虽然我给出的是Python代码。)

 for i in range(1, n): #1到n-1的循环
    if gcd(i, n)==1:   #如果i与n互素
        for j in range(1, n):   #1到n-1的循环
            if gcd(j, n) == 1:   #如果j与n互素
                print ((i * j) % n),  #输出(i * j) mod n
        print                    #内循环结束,输出换行

令n=12,程序输出是:

1 5 7 11
5 1 11 7
7 11 1 5
11 7 5 1

这是什么意思?规律是什么?看多一些数据。令n=21,输出为:

1 2 4 5 8 10 11 13 16 17 19 20
2 4 8 10 16 20 1 5 11 13 17 19
4 8 16 20 11 19 2 10 1 5 13 17
5 10 20 4 19 8 13 2 17 1 11 16
8 16 11 19 1 17 4 20 2 10 5 13
10 20 19 8 17 16 5 4 13 2 1 11
11 1 2 13 4 5 16 17 8 19 20 10
13 5 10 2 20 4 17 1 19 11 16 8
16 11 1 17 2 13 8 19 4 20 10 5
17 13 5 1 10 2 19 11 20 16 8 4
19 17 13 11 5 1 20 16 10 8 4 2
20 19 17 16 13 11 10 8 5 4 2 1

请注意每一行开始的那个数字与n的关系!可以考虑30分钟再往下看。

规律:

一个与n互素的整数i,它分别与所有大于等于1、小于n且与n互素的整数相乘(mod n),所得的整数是所有大于等于1、小于n且与n互素的整数的排列。

记大于等于1、小于n且与n互素的整数的个数为phi(n)。记大于等于1、小于n且与n互素的整数的相乘为Pi(n)。利用在求费尔马小定理时的技巧:

i^phi(n) * Pi(n) ≡ Pi(n) mod n
即:i^phi(n) ≡ 1 mod n

这个公式就是大名鼎鼎的欧拉定理!当n为素数,那么phi(n) == n-1。所以,费尔马小定理只是欧拉定理的一种特殊情形。

欧拉定理

若n,a为正整数,且 gcd(a,n)==1,
则 a^phi(n) ≡ 1 mod n.

注:

phi(n) = phi(n1)*phi(n2), if n == n1 * n2。

如何证明?It is easy,略。(注,如果你想问,容易到底有多容易? 请期待!)

注意,无论是费尔马小定理还是欧拉定理,这里都没有严格证明。请有兴趣的同学自己完成。

2017-06-27整理修订

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 202,607评论 5 476
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 85,047评论 2 379
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 149,496评论 0 335
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,405评论 1 273
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,400评论 5 364
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,479评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,883评论 3 395
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,535评论 0 256
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,743评论 1 295
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,544评论 2 319
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,612评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,309评论 4 318
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,881评论 3 306
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,891评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,136评论 1 259
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 42,783评论 2 349
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,316评论 2 342

推荐阅读更多精彩内容