約数関数

ファイル:Divisor.svg ファイル:Sigma function.svg 約数関数（やくすうかんすう、英: divisor function）は、自然数 n を変数とする関数で、n の全ての約数を整数乗した数の総和を値にとるものである。

定義

約数関数 σ_x(n) は自然数 n の約数 d の x 乗の総和を値として持つ関数であり、式では

$\sigma_x (n)=\sum_{d|n} d^x\,\!$

のように表される。x = 0 のとき σ₀(n) は n の約数の個数を表し、d(n)と表記されることもある。x = 1 のとき σ₁(n) は n の約数の総和であり、単に σ(n) と表記する場合もある。

計算例と特徴

例えば n = 15 では d(15) = σ₀(15) = 1⁰ + 3⁰ + 5⁰ + 15⁰ = 4, σ(15) = σ₁(15) = 1¹ + 3¹ + 5¹ + 15¹ = 24, σ₂(15) = 1² + 3² + 5² + 15² = 260

素数 p を変数とすると、p の約数は 1 と p の 2 個のみであるから d(p) = 2, σ(p) = p + 1 となる。また pⁿ の約数は 1, p, p², …, pⁿ の n + 1 個なので d(pⁿ) = n + 1, σ(pⁿ) = (p^{n + 1} - 1)/(p - 1) となる。

d(n) および σ(n) は n = 1 のとき最小値 1 をとる。d(n) = n の解は n = 1, 2 の 2 個のみであり、σ(n) = n の解や d(n) = σ(n) の解は n = 1 のみである。n ≧ 3 では 2 ≦ d(n) < n < σ(n) が成り立つ。

nを素因数分解して以下の式の形で表す。

$n= \prod_{i=1}^r {p_i}^{a_i}$

ここで r は n の素因子の個数、p_i はその中で i 番目に小さい素因子、a_i は素因数分解で現れる各素因子の指数部である。ここから

$\sigma_x (n)= \prod_{i=1}^r \frac{{p_i}^{(a_i +1)x} -1}{{p_i}^x -1} \quad (x \ne 0)$

が導かれる。これは

$\sigma_x (n) = \prod_{i=1}^r \sum_{j=0}^{a_i} {p_i}^{jx} = \prod_{i=1}^r (1+ {p_i}^x + {p_i}^{2x} + \cdots + {p_i}^{a_i x} )$

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

$\sigma_0(n)\equiv d(n)=\prod_{i=1}^{r} (a_i+1)$

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

• 約数関数から導き出される数列 $a_n=\sigma(a_{n-1})$ はその初期値によって異なる発散の仕方をする。( a₁ = 1 を除く)

例. a₁ = 2 のとき 2, 3, 4, 7, 8, 15, 24, 60, 168, 480, … (オンライン整数列大辞典の数列 A007497)

a₁ = 5 のとき 5, 6, 12, 28, 56, 120, 360, 1170, 3276, … (オンライン整数列大辞典の数列 A051572)

a₁ = 16 のとき 16, 31, 32, 63, 104, 210, 576, 1651, 1792, … (オンライン整数列大辞典の数列 A257349)

この初期値は 2, 5, 16, 19, 27, 29, 33, 49, 50, 52, 66, 81, 85, 105,… (オンライン整数列大辞典の数列 A257348)

その他の公式

約数関数は以下の三角関数を用いた式で表すこともできる。

$\sigma_x (n)= \sum_{\mu=1}^n \mu^{x-1} \sum_{\nu=1}^{\mu} \cos \frac{2\pi\nu n}{\mu}$

またゼータ関数 ζ(s) とは

$\sum_{n=1}^{\infty} \frac{\sigma_a (n)}{n^s} =\zeta(s) \zeta(s-a)$

という関係式をもつ。

σ(n)の増加の割合は以下の式で表される。

$\limsup_{n\rightarrow\infty}\frac{\sigma(n)}{n\ \log \log n}=e^\gamma$

γ はオイラー定数である。

また、d(n)の増加の割合は以下の式で表される。

$\limsup_{n\rightarrow\infty}\frac{\log d(n) \log\log n}{\log n}=\log 2$

実際、左辺の上極限記号内の分数の値が最大となるのは $n=6983776800$ のときで、その値は $1.0660186\ldots$ であることが知られている。 特に、任意の ε > 0 に対して、d(n) = o(n^ε) が成り立つ。

$\sigma (n)<e^{\gamma} n\log \log n \,$　(n > 5040)

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

また、実用性のある式ではないが、

$\lim_{n\to\infty}\sum_{k=2}^n\left(\frac{\frac{\sigma\left(\prod_{m=1}^np(m)^n\right)}{\prod_{m=1}^np(m)^n}}{\frac{\sigma\left(\prod_{m=1}^np(m)^k\right)}{\prod_{m=1}^np(m)^k}}-1\right)=1$

である。ただし、p(m)はm番目の素数。

約数関数の値

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

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

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

• σ(n) = 2n - 1 を満たす n を概完全数といい、σ(n) = 2n + 1 を満たす n を準完全数という。概完全数は 2 の累乗（1 も含む）が知られているが、それ以外に存在するかどうか知られていない。準完全数は存在するかどうか未だに分かっていない。準完全数が存在するならば、それは奇数の平方数でなければならないことが知られている。
• σ(n) = kn (k:整数) を満たす n を k-倍完全数という。例えば 120 は3倍完全数である。現在知られている倍積完全数は n = 1(このとき、k = 1)を除いて全て偶数である。1 以外に奇数の倍積完全数が存在するか否かは知られていない。
• σ(σ(n)) = 2n を満たす n を超完全数という。偶数の超完全数はメルセンヌ素数 2^p - 1 に対して、2^{p - 1} と表されるもの全体と一致することが知られている。奇数の超完全数が存在するか否かは知られていない。奇数の超完全数が存在するならば、それは平方数で少なくとも2つの相異なる素因数を持たなければならないことが知られている。
  • 超完全数は 2, 4, 16, 64, 4096, 65536, 262144,…である。 (オンライン整数列大辞典の数列 A019279)
• σ(n) = σ(m) = n + m を満たす相異なる数 n, m の組を友愛数という。(n, m) = (220, 284)などがそれである。友愛数が無限に存在するか否かは知られていない。
• d(n) = d(n + 1) を満たす n は無数に存在することが証明されている。 (例 : n = 2, 14, 21, 26, 33, 34, ･･･(オンライン整数列大辞典の数列 A005237))

関連項目

• ゼータ関数
• リーマン予想
• 不足数-完全数-過剰数
• 友愛数