青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
隨筆:152 文章:0 評論:129 引用:0
Headacher
學習筆記,從一點一滴做起。
C++博客
首頁
發新隨筆
發新文章
聯系
聚合
管理
歐拉函數
E(x)表示比x小的且與x互質的正整數的個數。
*若p是素數,E(p)=p-1。
*E(p^k)=p^k-p^(k-1)=(p-1)*P^(k-1)
證:令n=p^k,小于n的正整數數共有n-1即(p^k-1)個,其中與p不質的數共[p^(k-1)-1]個(分別為1*p,2*p,3*p...p(p^(k-1)-1))。
所以E(p^k)=(p^k-1)-(p^(k-1)-1)=p^k-p^(k-1).得證。
*若ab互質,則E(a*b)=E(a)*E(b),歐拉函數是積性函數.
*對任意數n都可以唯一分解成n=p1^a1*p2^a2*p3^a3*...*pn^an(pi為素數).
則E(n)=E(p1^a1)*E(p2^a2)*E(p3^a3)*...*E(pn^an)
=(p1-1)*p1^(a1-1)*(p2-1)*p2^(a2-1)*...*(pn-1)*pn^(an-1)
=(p1^a1*p2^a2*p3^a3*...*pn^an)*[(p1-1)*(p2-1)*(p3-1)*...*(pn-1)]/(p1*p2*p3*...*pn)
=n*(1-1/p1)*(1-1/p2)*...*(1-1/pn)
* E(p^k) =(p-1)*p^(k-1)=(p-1)*p^(k-2)*p
E(p^(k-1))=(p-1)*p^(k-2)
->當k>1時,E(p^k)=E(p*p^(k-1))=E(p^(k-1))*p.
(當k=1時,E(p)=p-1.)
由上式: 設P是素數,
若p是x的約數,則E(x*p)=E(x)*p.
若p不是x的約數,則E(x*p)=E(x)*E(p)=E(x)*(p-1).
*快速求歐拉函數方法:
首先來回顧一下線性篩選素數方法:
code
for
(i
=
2
;i
<=
1000000
;i
++
)
{
if
(
!
c[i])prime[len
++
]
=
i;
for
(j
=
0
;j
<
len
&&
prime[j]
*
i
<=
1000000
;j
++
)
{
c[prime[j]
*
i]
=
1
;
//
不是質數
if
(i
%
prime[j]
==
0
)
break
;
//
}
}
}
然后求歐拉函數:
Phi1
phi[
1
]
=
1
;
for
(i
=
2
; i
<
10000
; i
++
) {
if
(
!
mark[i]) {
phi[i]
=
i
-
1
;
continue
;
}
for
(j
=
0
; j
<
size
&&
prime[j]
*
prime[j]
<=
i; j
++
) {
if
(i
%
prime[j]
==
0
) {
if
(i
/
prime[j]
%
prime[j]
==
0
)
phi[i]
=
prime[j]
*
phi[i
/
prime[j]];
else
phi[i]
=
(prime[j]
-
1
)
*
phi[i
/
prime[j]];
break
;
}
}
}
由以上思想,可以在篩選素數的過程中求出歐拉函數:
Phi
for
(i
=
2
;i
<=
limit;i
++
)
{
if
(mark[i]
==
0
)
{
prime[
++
prime[
0
]]
=
i;
E[i]
=
i
-
1
;
}
for
(j
=
1
;j
<=
prime[
0
]
&&
prime[j]
*
i
<=
limit;j
++
)
{
mark[prime[j]
*
i]
=
1
;
if
(i
%
prime[j]
==
0
)
{
E[i
*
prime[j]]
=
E[i]
*
prime[j];
break
;
}
else
E[i
*
prime[j]]
=
E[i]
*
(prime[j]
-
1
);
}
}
發表于 2009-07-19 12:39
Headacher
閱讀(3147)
評論(0)
編輯
收藏
引用
所屬分類:
數據結構和算法
只有注冊用戶
登錄
后才能發表評論。
相關文章:
POJ 2043 掃描 計算幾何
POJ 1113 凸包
POJ 3164 最小樹形圖 朱劉算法
POJ 2761 SBT 靜態數組實現
POJ 2778 自動機_矩陣乘法
HDU 2222 AC自動機
數位統計
無恥IO優化
哦哦
有上下界的可行流
網站導航:
博客園
IT新聞
BlogJava
博問
Chat2DB
管理
CALENDER
<
2008年11月
>
日
一
二
三
四
五
六
26
27
28
29
30
31
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
1
2
3
4
5
6
公告
留言簿
(8)
給我留言
查看公開留言
查看私人留言
隨筆分類
ACM-ICPC(7)
(rss)
操作系統
(rss)
計算機組成與體系結構(2)
(rss)
數據結構和算法(34)
(rss)
數據庫
(rss)
心情日記(20)
(rss)
隨筆檔案
2010年12月 (1)
2010年9月 (1)
2010年5月 (3)
2010年4月 (3)
2010年3月 (1)
2010年2月 (2)
2010年1月 (10)
2009年12月 (1)
2009年10月 (3)
2009年9月 (6)
2009年8月 (14)
2009年7月 (8)
2009年6月 (2)
2009年5月 (17)
2009年4月 (4)
2009年3月 (5)
2009年2月 (25)
2009年1月 (9)
2008年12月 (1)
2008年11月 (30)
2008年10月 (4)
2008年7月 (2)
ACM Teammates
Qinz
(rss)
SHFACM
(rss)
wudired
(rss)
The One
May
(rss)
搜索
積分與排名
積分 - 135369
排名 - 192
最新評論
1.?re: POJ 1379 run away 模擬退火算法[未登錄]
為何按你的代碼交會RE呢?
--zhang
2.?re: POJ 1947 樹狀dp[未登錄]
評論內容較長,點擊標題查看
--Sky
3.?re: 獨立集,覆蓋集,支配集,最大團,最大匹配
評論內容較長,點擊標題查看
--fly2best
4.?re: HDU HDOJ 1004 Let the Balloon Rise 字典樹[未登錄]
尼瑪 這就是個水題
--xxx
5.?re: nuaa 1017 最大0,1子矩陣[未登錄]
1 0 1 0 1
2 1 2 1 2
3 2 2 2 0
0 3 4 3 1
1 0 5 4 2 這個寫錯了吧
第三行第三列那個2應該為3才對
--hu
閱讀排行榜
1.?獨立集,覆蓋集,支配集,最大團,最大匹配(7952)
2.?原碼 補碼 反碼 移碼(6431)
3.?POJ 計算幾何入門題目推薦(轉)(5722)
4.?POJ 1379 run away 模擬退火算法(4431)
5.?數據的浮點數表示(3979)
評論排行榜
1.?POJ 1379 run away 模擬退火算法(12)
2.?我真是太笨了……(10)
3.?PKU POJ 2186 Popular Cows 強連通分量(5)
4.?PKU POJ 1679 The Unique MST 次小生成樹(4)
5.?HDU HDOJ 1005 Number Sequence(4)
Powered By:
博客園
模板提供
:
滬江博客
青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
亚洲电影av在线
|
99国产精品久久久久久久久久
|
性欧美暴力猛交69hd
|
韩国精品在线观看
|
国产精品国产亚洲精品看不卡15
|
久久久久免费观看
|
亚洲最黄网站
|
亚洲精品久久久久久下一站
|
亚洲中无吗在线
|
影视先锋久久
|
欧美日韩中国免费专区在线看
|
久久精品亚洲国产奇米99
|
一区二区三区精品视频在线观看
|
亚洲欧美韩国
|
亚洲欧美成人一区二区三区
|
在线观看国产成人av片
|
国产日本亚洲高清
|
国产一区二区电影在线观看
|
国产精品有限公司
|
欧美日韩在线不卡一区
|
国产精品v欧美精品v日韩精品
|
久久这里只有
|
久久婷婷色综合
|
欧美成人嫩草网站
|
国产精品社区
|
尤物yw午夜国产精品视频明星
|
亚洲国产成人tv
|
亚洲黄色成人网
|
国产一区视频在线观看免费
|
午夜精品久久
|
女仆av观看一区
|
久久不射电影网
|
久久一区中文字幕
|
欧美日韩高清在线
|
一区二区三区在线视频免费观看
|
国产精品国产精品
|
国语精品中文字幕
|
亚洲欧洲免费视频
|
亚洲少妇在线
|
欧美黑人多人双交
|
亚洲欧美亚洲
|
在线播放豆国产99亚洲
|
亚洲一区二区三区四区视频
|
久久精品99国产精品酒店日本
|
亚洲美女黄色片
|
欧美成黄导航
|
国产精品成人免费
|
亚洲三级观看
|
欧美一区中文字幕
|
久久综合伊人
|
亚洲精品国产精品乱码不99
|
久久riav二区三区
|
欧美国产激情
|
欧美激情亚洲
|
国产一区二区无遮挡
|
性欧美xxxx大乳国产app
|
亚洲一卡久久
|
在线电影国产精品
|
日韩亚洲成人av在线
|
国产原创一区二区
|
亚洲欧洲精品一区二区精品久久久
|
免费国产一区二区
|
欧美激情综合色
|
欧美一级艳片视频免费观看
|
久久久精品免费视频
|
99这里只有久久精品视频
|
亚洲一二区在线
|
亚洲精品黄色
|
久久蜜臀精品av
|
欧美一区二区三区在线视频
|
国产精品海角社区在线观看
|
欧美11—12娇小xxxx
|
国产精品久久一区主播
|
亚洲国产精品悠悠久久琪琪
|
国产精品vvv
|
亚洲精品一区二区三
|
亚洲人永久免费
|
久久综合久久88
|
亚洲人成网站在线播
|
国产精品国产
|
久久亚洲一区二区三区四区
|
亚洲性线免费观看视频成熟
|
亚洲黄一区二区三区
|
在线视频精品
|
亚洲区在线播放
|
亚洲精品久久视频
|
亚洲韩国精品一区
|
亚洲精品久久久久久一区二区
|
一区二区三区高清不卡
|
黑人极品videos精品欧美裸
|
亚洲一二三级电影
|
99av国产精品欲麻豆
|
中文国产一区
|
尤物yw午夜国产精品视频
|
欧美国产免费
|
亚洲欧美欧美一区二区三区
|
久热re这里精品视频在线6
|
黄色日韩在线
|
国产精品v欧美精品v日韩
|
久久国产精品久久精品国产
|
亚洲九九九在线观看
|
久久综合伊人77777蜜臀
|
日韩午夜在线观看视频
|
国产欧美精品日韩
|
欧美新色视频
|
欧美日产一区二区三区在线观看
|
欧美亚洲免费高清在线观看
|
亚洲午夜免费福利视频
|
亚洲人成网在线播放
|
欧美aⅴ99久久黑人专区
|
久久夜精品va视频免费观看
|
欧美与欧洲交xxxx免费观看
|
亚洲综合大片69999
|
亚洲性感美女99在线
|
在线亚洲伦理
|
亚洲欧美在线磁力
|
久久九九有精品国产23
|
久久精品国产亚洲aⅴ
|
欧美成人精品福利
|
国产一区视频网站
|
久久久久国内
|
午夜亚洲伦理
|
国产精品二区二区三区
|
久久久久网站
|
久久久亚洲国产美女国产盗摄
|
国产精品揄拍500视频
|
亚洲影院免费观看
|
欧美成人一区二区三区
|
国产视频在线观看一区
|
亚洲日本视频
|
另类专区欧美制服同性
|
久久久国产一区二区
|
精品不卡在线
|
国产精品任我爽爆在线播放
|
欧美美女喷水视频
|
国产精品高潮呻吟
|
亚洲国产成人av在线
|
亚洲视频在线观看免费
|
欧美一区二区三区在线看
|
久久精品卡一
|
欧美色综合网
|
亚洲国产婷婷
|
久久人人爽人人
|
亚洲图片欧美日产
|
欧美黄在线观看
|
黄网站免费久久
|
久久成人综合网
|
亚洲一区尤物
|
国产伦精品免费视频
|
欧美精品在线视频
|
国产伦精品一区二区三区照片91
|
91久久精品日日躁夜夜躁欧美
|
欧美一级久久久
|
亚洲少妇一区
|
国产精品老牛
|
亚洲欧美成人在线
|
亚洲午夜电影在线观看
|
欧美日韩极品在线观看一区
|
日韩午夜在线观看视频
|
亚洲精品三级
|
欧美日一区二区在线观看
|
一本色道88久久加勒比精品
|
欧美一区观看
|
久久精品国产69国产精品亚洲
|
国产一区二区看久久
|
久热精品视频在线免费观看
|
一本色道久久综合精品竹菊
|
欧美电影资源
|
欧美日韩mv
|
欧美综合第一页
|
久久精品夜色噜噜亚洲a∨
|
精品99一区二区
|
亚洲精品久久视频
|
国产美女精品在线
|
欧美激情第二页
|
国产精品视频99
|
欧美国产日韩一区二区
|
国产精品色午夜在线观看
|
欧美成人一二三
|
国产日韩精品一区观看
|
亚洲国产成人高清精品
|
国产精品一区毛片
|
亚洲成人在线网
|
在线观看三级视频欧美
|
亚洲综合视频1区
|
一本不卡影院
|
欧美成人午夜激情
|
蜜桃av一区二区在线观看
|
国产精品色午夜在线观看
|
最近中文字幕mv在线一区二区三区四区
|
国产精品一区二区在线观看网站
|
欧美成人国产
|
亚洲国产老妈
|
久久久久久久久久久久久9999
|
亚洲视频精选
|
国产精品www网站
|
亚洲视频1区2区
|
欧美一区二区三区免费观看视频
|
欧美日韩国产丝袜另类
|
在线视频日本亚洲性
|