oyjpArt ACM/ICPC算法程序设计空间
// I am new in programming, welcome to my blog
I am oyjpart(alpc12, 四城)

【转】质数初步

<转>

模p运算

给定一个正整数p,任意一个整数n,一定存在等式

          n = kp + r

其中k、r是整数,且 0 ≤ r < p,称呼k为n除以p的商,r为n除以p的余数。

对于正整数p和整数a,b,定义如下运算:

  • 取模运算:a mod p 表示a除以p的余数。
  • 模p加法:(a + b) mod p ,其结果是a+b算术和除以p的余数,也就是说,(a+b) = kp +r,则 (a+b) mod p = r。
  • 模p减法:(a-b) mod p ,其结果是a-b算术差除以p的余数。
  • 模p乘法:(a × b) mod p,其结果是 a × b算术乘法除以p的余数。

可以发现,模p运算和普通的四则运算有很多类似的规律,如:
规律 公式
结合率 ((a+b) mod p + c)mod p = (a + (b+c) mod p) mod p
((a*b) mod p * c)mod p = (a * (b*c) mod p) mod p
交换率 (a + b) mod p = (b+a) mod p
(a × b) mod p = (b × a) mod p
分配率 ((a +b)mod p × c) mod p = ((a × c) mod p + (b × c) mod p) mod p

简单的证明其中第一个公式:

 ((a+b) mod p + c) mod p = (a + (b+c) mod p) mod p
假设
a = k1 p + r1
b = k2 p + r2
c = k3 p + r3

a+b = (k1 + k2) p + (r1 + r2)
如果(r1 + r2) >= p ,则
(a+b) mod p = (r1 + r2) -p
否则
(a+b) mod p = (r1 + r2)
再和c进行模p和运算,得到
结果为 r1 + r2 + r3的算术和除以p的余数。
对右侧进行计算可以得到同样的结果,得证。

模p相等

如果两个数a、b满足a mod p = b mod p,则称他们模p相等,记做

 a ≡ b mod p
可以证明,此时a、b满足 a = kp + b,其中k是某个整数。

对于模p相等和模p乘法来说,有一个和四则运算中迥然不同得规则。在四则运算中,如果c是一个非0整数,则

       ac = bc 可以得出  a =b

但是在模p运算中,这种关系不存在,例如:

 (3 x 3) mod 9 = 0
(6 x 3) mod 9 = 0
但是
3 mod 9 = 3
6 mod 9 =6

定理(消去律):如果gcd(c,p) = 1 ,则 ac ≡ bc mod p 可以推出 a ≡ b mod p

 证明:
因为ac ≡ bc mod p
所以ac = bc + kp,也就是c(a-b) = kp
因为c和p没有除1以外的公因子,因此上式要成立必须满足下面两个条件中的一个
1) c能整除k
2) a = b
如果2不成立,则c|kp
因为c和p没有公因子,因此显然c|k,所以k = ck'
因此c(a-b)kp可以表示为c(a-b) =ck'p
因此a-b = k'p,得出a ≡ b mod p
如果a = b,则a ≡ b mod p 显然成立
得证

欧拉函数

欧拉函数是数论中很重要的一个函数,欧拉函数是指:对于一个正整数n,小于n且和n互质的正整数的个数,记做:φ(n),其中φ(1)被定义为1,但是并没有任何实质的意义。

定义小于n且和n互质的数构成的集合为Zn,称呼这个集合为n的完全余数集合。

显然,对于素数p,φ(p)= p -1.对于两个素数p、q,他们的乘积n = pq 满足φ(n) =(p-1)(q-1)

        证明:对于质数p,q,满足φ(n) =(p-1)(q-1)
考虑n的完全余数集Zn = { 1,2,....,pq -1}
而不和n互质的集合由下面三个集合的并构成:
1) 能够被p整除的集合{p,2p,3p,....,(q-1)p} 共计q-1个
2) 能够被q整除的集合{q,2q,3q,....,(p-1)q} 共计p-1个
3) {0}
很显然,1、2集合中没有共同的元素,因此Zn中元素个数 = pq - (p-1 + q- 1 + 1) = (p-1)(q-1)

欧拉定理

对于互质的整数a和n,有aφ(n) ≡ 1 mod n

        证明:
首先证明下面这个命题:
对于集合Zn={x1,x2,...,xφ(n)},考虑集合
S = {ax1 mod n,ax2mod n,...,axφ(n)mod n}
则S = Zn

1) 由于a,n互质,xi也与n互质,则axi也一定于p互质,因此
任意xi,axi mod n 必然是Zn的一个元素
2) 对于Zn中两个元素xi和xj,如果xi ≠ xj
则axi mod n ≠ axi mod n,这个由a、p互质和消去律可以得出。
所以,很明显,S=Zn

既然这样,那么
(ax1 × ax2×...×axφ(n))mod n
= (ax1 mod n × ax2mod n × ... × axφ(n)mod n)mod n
= (x1 × x2 × ... × xφ(n))mod n
考虑上面等式左边和右边
左边等于(aφ(n) × (x1 × x2 × ... × xφ(n))mod n) mod n
右边等于x1 × x2 × ... × xφ(n))mod n
而x1 × x2 × ... × xφ(n))mod n和p互质
根据消去律,可以从等式两边约去,就得到:
aφ(n) ≡ 1 mod n
推论:对于互质的数a、n,满足aφ(n)+1 ≡ a mod n

费马定理


a是不能被质数p整除的正整数,则有ap-1 ≡ 1 mod p


证明这个定理非常简单,由于φ(p) = p-1,代入欧拉定理即可证明。


同样有推论:对于不能被质数p整除的正整数a,有ap ≡ a mod p

发表于 2008-03-21 20:55 oyjpart 阅读(687) 评论(0)  编辑 收藏 引用 所属分类: ACM/ICPC
 

标题  
姓名  
主页
验证码 *
内容(提交失败后,可以通过“恢复上次提交”恢复刚刚提交的内容)  
  登录  使用高级评论  新用户注册  返回页首  恢复上次提交      
[使用Ctrl+Enter键可以直接提交]
相关文章:
  • SRM406 PTS500 FoldThePaper
  • 点错题,看错题,做错题!可以去死了!!!!!!!!!
  • 线性筛素数方法
  • 基本参数搜索
  • 湖大邀请赛
  • 1周又过去了
  • 杭州邀请赛 ACM在何方?
  • SRM401 550 PTS ParticleCollision
  • 中南赛A题 Accumulation Degree
  • PKU2504 Rounding Box
博客园      BlogJava   博客生活   IT博客网   C++博客   PHP博客   博客园社区
管理博客   教师博客   天文博客   汽车博客   足球博客   股票博客   电子技术博客
 
<2008年3月>
日一二三四五六
2425262728291
2345678
9101112131415
16171819202122
23242526272829
303112345

 公告


Contact me via POJ mail : Alpc12 and email(MSN also) : yescrystalblue@sina.com

 导航

  • C++博客
  • 首页
  • 发新文章
  • 联系
  • 聚合
  • 管理

 统计

  • 随笔: 172
  • 文章: 5
  • 评论: 492
  • 引用: 0

常用链接

  • 我的随笔
  • 我的评论
  • 我参与的随笔

留言簿(28)

  • 给我留言
  • 查看公开留言
  • 查看私人留言

随笔分类(134)

  • ACM/ICPC(88) (rss)
  • 程序设计(JAVA)(6) (rss)
  • 心情书写 艺术体验(39) (rss)
  • 英语学习(1) (rss)

随笔档案(172)

  • 2008年7月 (1)
  • 2008年6月 (11)
  • 2008年5月 (11)
  • 2008年4月 (3)
  • 2008年3月 (7)
  • 2008年2月 (4)
  • 2008年1月 (8)
  • 2007年12月 (2)
  • 2007年11月 (4)
  • 2007年10月 (2)
  • 2007年9月 (6)
  • 2007年8月 (4)
  • 2007年7月 (5)
  • 2007年6月 (8)
  • 2007年5月 (9)
  • 2007年4月 (10)
  • 2007年3月 (8)
  • 2007年2月 (10)
  • 2007年1月 (6)
  • 2006年12月 (4)
  • 2006年11月 (10)
  • 2006年10月 (10)
  • 2006年9月 (9)
  • 2006年8月 (20)

文章档案(5)

  • 2006年8月 (5)

相册

  • ACM/ICPC
  • MY PICTURES
  • YR
  • 贴图

My friends

  • ZZN's blog
  • !ASP' blog!
  • alpc01's blog
  • alpc02
  • alpc16's blog
  • alpc50's blog
  • alpc55's blog
  • alpc60's blog
  • alpc62's blog
  • byron's blog
  • Cockerel's Treasures
  • coldwings/alpc42
  • drizzlecrj
  • encu_zp's blog
  • felicia's blog
  • javaman's blog
  • owen's blog
  • qywyh's blog
  • richardxx's blog
  • roba's blog
  • Warrior's blog
  • wywcgs's blog
  • yiyiyi4321's blog

最新随笔

  • 1. 期待爱
  • 2. 《the negotiator》
  • 3. 《In.Bruges》
  • 4. 《The Oxford Murders》
  • 5. SRM406 PTS500 FoldThePaper
  • 6. 无题
  • 7. 向Roba学习
  • 8. 生活
  • 9. 点错题,看错题,做错题!可以去死了!!!!!!!!!
  • 10. 线性筛素数方法
  • 11. 基本参数搜索
  • 12. 这样的生活
  • 13. 湖大邀请赛
  • 14. [转]女人常说的32句谎话及详解...
  • 15. Tomcat JSP Web 开发中的乱码问题小结
  • 16. 【转】一页纸多一点的博士论文拿到诺奖
  • 17. 1周又过去了
  • 18. 对死者致以最大的尊敬 对生者致以最深的关怀
  • 19. 杭州邀请赛 ACM在何方?
  • 20. SRM401 550 PTS ParticleCollision
  • 21. 中南赛A题 Accumulation Degree
  • 22. PKU2504 Rounding Box
  • 23. Classical Movies
  • 24. Petr疯了?
  • 25. 08年中南赛--失意后的反思
  • 26. 有信息的搜索--传教士与食人魔的故事
  • 27. 手机终于修好了
  • 28. 【转】一个可以气死日本人的中国大学生!!
  • 29. 校赛裁判 : Judge ACs!
  • 30. 【转】质数初步
  • 31. 【转】正则表达式30分钟入门教程
  • 32. 【转】强大的poj分类
  • 33. 蔡蕾(male)的20个问答
  • 34. recent competitions
  • 35. 累哦
  • 36. Mario Forever
  • 37. SRM389, SRM390, Qual1
  • 38. 被TC说犯规了...
  • 39. JAVA常用设计模式2 MVC [ Model-View-Controller ]
  • 40. SRM388

搜索

  •  

积分与排名

  • 积分 - 158844
  • 排名 - 3

最新随笔

  • 1. 期待爱
  • 2. 《the negotiator》
  • 3. 《In.Bruges》
  • 4. 《The Oxford Murders》
  • 5. SRM406 PTS500 FoldThePaper
  • 6. 无题
  • 7. 向Roba学习
  • 8. 生活
  • 9. 点错题,看错题,做错题!可以去死了!!!!!!!!!
  • 10. 线性筛素数方法
  • 11. 基本参数搜索
  • 12. 这样的生活
  • 13. 湖大邀请赛
  • 14. [转]女人常说的32句谎话及详解...
  • 15. Tomcat JSP Web 开发中的乱码问题小结
  • 16. 【转】一页纸多一点的博士论文拿到诺奖
  • 17. 1周又过去了
  • 18. 对死者致以最大的尊敬 对生者致以最深的关怀
  • 19. 杭州邀请赛 ACM在何方?
  • 20. SRM401 550 PTS ParticleCollision
  • 21. 中南赛A题 Accumulation Degree
  • 22. PKU2504 Rounding Box
  • 23. Classical Movies
  • 24. Petr疯了?
  • 25. 08年中南赛--失意后的反思
  • 26. 有信息的搜索--传教士与食人魔的故事
  • 27. 手机终于修好了
  • 28. 【转】一个可以气死日本人的中国大学生!!
  • 29. 校赛裁判 : Judge ACs!
  • 30. 【转】质数初步
  • 31. 【转】正则表达式30分钟入门教程
  • 32. 【转】强大的poj分类
  • 33. 蔡蕾(male)的20个问答
  • 34. recent competitions
  • 35. 累哦
  • 36. Mario Forever
  • 37. SRM389, SRM390, Qual1
  • 38. 被TC说犯规了...
  • 39. JAVA常用设计模式2 MVC [ Model-View-Controller ]
  • 40. SRM388

最新评论

  • 1. re: 期待爱
  • 推荐歌曲啊
    Of course OK.~~
  • --oyjpart
  • 2. re: 期待爱
  • Are you OK??
  • --alpc62
  • 3. re: PKU3121 Sum of Different Primes
  • 恩?
  • --oyjpart
  • 4. re: PKU3121 Sum of Different Primes
  • 自己都没交成功
  • --ssadwll
  • 5. re: 闲来切题 呵呵
  • Contact me via POJ mail : alpc12
    email(MSN also) : yescrystalblue@sina.com
  • --oyjpart
  • 6. re: 闲来切题 呵呵
  • 能留个联系方式吗? 最近也在为ACM努力着。 有不懂的,希望你能帮我!

    谢谢了。。。
  • --lihao102
  • 7. re: 闲来切题 呵呵
  • 评论内容较长,点击标题查看
  • --oyjpart
  • 8. re: 闲来切题 呵呵
  • 1724 roads 双调路径能说进更详细些吗?

       我不是很明白,有讲这方面的资料吗?

               谢谢了。。。。。。
  • --lihao102
  • 9. re: PKU1465 Multiple
  • shayne ward的歌有几首都不错^^
  • --sdfond
  • 10. re: 线性筛素数方法
  • 呃 ...
  • --oyjpart

阅读排行榜

  • 1. KMP算法浅析(3833)
  • 2. 奶牛题比赛(2399)
  • 3. PKU 1456 Supermarket 并查集加速实例(2264)
  • 4. 闲来切题 呵呵(2098)
  • 5. PKU 1328 Radar Installation (1934)

评论排行榜

  • 1. 基本参数搜索(13)
  • 2. 一诀成都,金牌!(12)
  • 3. 闲来切题 呵呵(11)
  • 4. 肖申克的救赎(11)
  • 5. PKU 1011 Sticks (11)

60天内阅读排行

  • 1. 湖大邀请赛(1464)
  • 2. 基本参数搜索(1323)
  • 3. 线性筛素数方法(1292)
  • 4. 杭州邀请赛 ACM在何方?(671)
  • 5. 【转】一页纸多一点的博士论文拿到诺奖(476)

Powered by: 博客园
模板提供:沪江博客
Copyright ©2008 oyjpart