約数関数

提供: testwiki
2024年3月15日 (金) 04:54時点におけるimported>阿賀 雲斗による版 (注釈)
(差分) ← 古い版 | 最新版 (差分) | 新しい版 → (差分)
ナビゲーションに移動 検索に移動
nの約数の個数を表す
σ0(n)≡d(n) のグラフ(n≦250)
nの約数の総和を表す
σ1(n)≡σ(n) のグラフ(n≦250)

約数関数(やくすうかんすう、テンプレート:Lang-en-short)は、自然数 n変数とする関数で、n の全ての約数を整数乗した数の総和を値にとるものである。

定義

自然数 テンプレート:Mvar に対して、約数関数 テンプレート:Math とは、テンプレート:Mvar の約数 テンプレート:Mvarテンプレート:Mvar 乗和を値に取る関数である:

σx(n)=d|ndx

特に、テンプレート:Math2 のとき テンプレート:Mathテンプレート:Mvar の約数の個数を表し、テンプレート:Mathテンプレート:Math と表されることもある。テンプレート:Math2 のとき テンプレート:Mathテンプレート:Mvar の約数の総和であり、単に省略して テンプレート:Math と表す場合もある。

また、約数関数 テンプレート:Math の[[写像の反復|テンプレート:Mvar階反復]]を

σxk(n):=σx((σxk(n))

と書く。例えば σx2(n)=σx(σx(n)),σx3(n)=σx(σx(σx(n))) である。

k = 1 、x = 1 のときはどちらもそれぞれ省略して、σ(n) = σテンプレート:Sub(n)(k=x=1の場合)、σテンプレート:Sup(n)(k=2,x=1の場合)などと表記する場合もある。

概要

テンプレート:Math の値は、小さい順に次のようになる:

1, 2, 2, 3, 2, 4, 2, 4, 3, 4, 2, 6, 2, 4, 4 …(テンプレート:OEIS

テンプレート:Math の値は、小さい順に次のようになる:

1, 3, 4, 7, 6, 12, 8, 15, 13, 18, 12, 28, 14, 24, 24 …(テンプレート:OEIS

テンプレート:Math の値は、小さい順に次のようになる:

1, 5, 10, 21, 26, 50, 50, 85, 91, 130, 122, 210, 170, 250, 260 …(テンプレート:OEIS

計算例

例えば n = 15 では、

d(15) = σテンプレート:Sub(15) = 1テンプレート:Sup + 3テンプレート:Sup + 5テンプレート:Sup + 15テンプレート:Sup = 4,
σ(15) = σテンプレート:Sub(15) = 1テンプレート:Sup + 3テンプレート:Sup + 5テンプレート:Sup + 15テンプレート:Sup = 24,
σテンプレート:Sub(15) = 1テンプレート:Sup + 3テンプレート:Sup + 5テンプレート:Sup + 15テンプレート:Sup = 260

特徴

テンプレート:Mvar素数とすると、テンプレート:Mvar の約数は テンプレート:Mathテンプレート:Mvar の 2個のみであるから テンプレート:Math2 となる。また、テンプレート:Mvar を自然数とすると、テンプレート:Mvar の約数は テンプレート:Math2テンプレート:Math2個なので テンプレート:Math2 となる。

テンプレート:Math および テンプレート:Mathテンプレート:Math2 のとき最小値 1 をとる。テンプレート:Math2 の解は テンプレート:Math2 の 2 個のみであり、テンプレート:Math2 の解や テンプレート:Math2 の解は テンプレート:Math2 のみである。テンプレート:Math2 では テンプレート:Math2 が成り立つ。

約数関数 テンプレート:Math乗法的関数テンプレート:Lang-en-short)であるが、テンプレート:仮リンクではない。

gcd(a,b)=1σx(ab)=σx(a)σx(b).

テンプレート:Mvar素因数分解して以下の式の形で表す。

n=i=1rpiai

ここで rn素因子の個数、pテンプレート:Sub はその中で i 番目に小さい素因子、aテンプレート:Sub は素因数分解で現れる各素因子の指数部である。ここから

σx(n)=i=1rpi(ai+1)x1pix1(x0)

が導かれる。これは

σx(n)=i=1rj=0aipijx=i=1r(1+pix+pi2x++piaix)

同値である。x = 0 のときは

σ0(n)d(n)=i=1r(ai+1)

となる。例えば n = pqp, q は素数)とすると、σ(n) = (1 + p)(1 + q) = n + p + q + 1, d(n)=(1 + 1)(1 + 1) = 4 となる。

  • 約数関数から導き出される数列 an=σ(an1) はその初期値によって異なる発散の仕方をする。( aテンプレート:Sub = 1 を除く)
例. aテンプレート:Sub = 2 のとき 2, 3, 4, 7, 8, 15, 24, 60, 168, 480, … (テンプレート:OEIS)
aテンプレート:Sub = 5 のとき 5, 6, 12, 28, 56, 120, 360, 1170, 3276, … (テンプレート:OEIS)
aテンプレート:Sub = 16 のとき 16, 31, 32, 63, 104, 210, 576, 1651, 1792, … (テンプレート:OEIS)
この初期値は 2, 5, 16, 19, 27, 29, 33, 49, 50, 52, 66, 81, 85, 105,… (テンプレート:OEIS)

その他の公式

オイラーは約数関数が以下のように表されることを示した。[1]

  σ1(n)=σ1(n1)+σ1(n2)σ1(n5)σ1(n7)+σ1(n12)+σ1(n15)...=i=1(1)i+1(σ1(n12(3i2i))+σ1(n12(3i2+i)))

なおこの数式で、n<0 のとき σ1(n)=0 とし、σ1(0)=n とする。

約数関数の母関数はテンプレート:仮リンクである。

n=1naxn1xn=n=1m=1naxmn=n=1xnσa(n)

約数関数は以下の三角関数を用いた式で表すこともできる。 テンプレート:Indent

またゼータ関数 ζ(s) とは テンプレート:Indent という関係式をもつ。

σ(n)の増加の割合は以下の式で表される。 テンプレート:Indent γ はオイラー定数である。

また、d(n)の増加の割合は以下の式で表される。 テンプレート:Indent 実際、左辺の上極限記号内の分数の値が最大となるのは n=6983776800 のときで、その値は 1.0660186 であることが知られている[2]。 特に、任意の ε > 0 に対して、d(n) = o(nε) が成り立つ。

σ(n)<eγnloglogn (n > 5040)

が真であるならリーマン予想も真であることが証明されている。つまりこの不等式を満たさない最大の数が 5040 であり[3]、5041 以上の全ての自然数がこの不等式を満たすならばリーマン予想は真である。もしリーマン予想が偽ならこの不等式を満たさない n は無数に存在する。


約数関数の値

x=0~21についてのσx(n)の値はオンライン整数列大辞典に数列として掲載されている。

オンライン整数列大辞典に掲載されている約数関数
x 約数関数 σx(n) 値のリスト
0 σ0(n) (テンプレート:OEIS) Table of n, a(n) for n = 1..10000
1 σ1(n) (テンプレート:OEIS) Table of n, a(n) for n = 1..10000
2 σ2(n) (テンプレート:OEIS) Table of n, a(n) for n = 1..10000
3 σ3(n) (テンプレート:OEIS) Table of n, a(n) for n = 1..10000
4 σ4(n) (テンプレート:OEIS) Table of n, a(n) for n = 1..10000
5 σ5(n) (テンプレート:OEIS) Table of n, a(n) for n = 1..10000
6 σ6(n) (テンプレート:OEIS) Table of n, a(n) for n = 1..1000
7 σ7(n) (テンプレート:OEIS) Table of n, a(n) for n = 1..10000
8 σ8(n) (テンプレート:OEIS) Table of n, a(n) for n = 1..1000
9 σ9(n) (テンプレート:OEIS) Table of n, a(n) for n = 1..1000
10 σ10(n) (テンプレート:OEIS) Table of n, a(n) for n = 1..10000
11 σ11(n) (テンプレート:OEIS) Table of n, a(n) for n = 1..10000
12 σ12(n) (テンプレート:OEIS) Table of n, a(n) for n = 1..10000
13 σ13(n) (テンプレート:OEIS) Table of n, a(n) for n = 1..10000
14 σ14(n) (テンプレート:OEIS) Table of n, a(n) for n = 1..10000
15 σ15(n) (テンプレート:OEIS) Table of n, a(n) for n = 1..10000
16 σ16(n) (テンプレート:OEIS) Table of n, a(n) for n = 1..10000
17 σ17(n) (テンプレート:OEIS) Table of n, a(n) for n = 1..10000
18 σ18(n) (テンプレート:OEIS) Table of n, a(n) for n = 1..10000
19 σ19(n) (テンプレート:OEIS) Table of n, a(n) for n = 1..10000
20 σ20(n) (テンプレート:OEIS) Table of n, a(n) for n = 1..10000
21 σ21(n) (テンプレート:OEIS) Table of n, a(n) for n = 1..10000

σ(n) < 2n を満たす n不足数、σ(n) = 2n を満たす n完全数、σ(n) > 2n を満たす n過剰数という。

6, 28, 496 などが完全数として知られている。偶数の完全数全体はメルセンヌ素数 2p − 1 に対して 2p−1(2p − 1) と表されるもの全体と一致することが知られている。奇数の完全数が存在するかどうかは古くからの数論の未解決問題として有名である。

このほかにも、約数関数、特に約数の和の関数 σ(n) の値に関しては多くの概念が考察され、多くの未解決問題が提示されている。いくつかの例を挙げる。

  • σ(n) = 2n − 1 を満たす n概完全数といい、σ(n) = 2n + 1 を満たす n準完全数という。概完全数は 2 の累乗(1 も含む)が知られているが、それ以外に存在するかどうか知られていない。準完全数は存在するかどうか未だに分かっていない。準完全数が存在するならば、それは奇数の平方数でなければならないことが知られている。
  • σ(n) = kn (k:整数) を満たす nk-倍完全数という。例えば 120 は3倍完全数である。現在知られている倍積完全数n = 1(このとき、k = 1)を除いて全て偶数である。1 以外に奇数の倍積完全数が存在するか否かは知られていない。
  • σ(σ(n)) = 2n を満たす n超完全数という。偶数の超完全数メルセンヌ素数 2p − 1 に対して、2p−1 と表されるもの全体と一致することが知られている。奇数の超完全数が存在するか否かは知られていない。奇数の超完全数が存在するならば、それは平方数で少なくとも2つの相異なる素因数を持たなければならないことが知られている。
  • σ(n) = σ(m) = n + m を満たす相異なる数 n, m の組を友愛数という。(n, m) = (220, 284)などがそれである。友愛数が無限に存在するか否かは知られていない。
    • それと対照的に、 σ(n) = σ(m) = n + m + 1 を満たす相異なる数 n, m の組を婚約数という。(n, m) = (140, 195)などがそれである。婚約数は無限に存在するか否かは証明されていない。
  • d(n) = d(n + 1) を満たす n は無数に存在することが証明されている。 (例 : n = 2, 14, 21, 26, 33, 34, ・・・(テンプレート:OEIS))

関連項目

注釈

  1. テンプレート:Cite arXiv
  2. J. L. Nicolas et G. Robin, Majorations explicites pour le nombre de diviseurs de $N$, Canad. Math. Bull. 26 (1983), 485--492.
  3. σ(5040) = 19344, eγ ・ 5040 log log 5040 = 19237.84...

de:Teilersumme hu:Osztóösszeg-függvény pl:Funkcja σ テンプレート:Divisor classes