完全数(perfect number)とは、自分自身を除く正の約数の和が自分自身に等しい正の整数、すなわち約数の和 $\sigma(N)$ が $\sigma(N)=2N$ を満たす $N$ のことである。$6=1+2+3$、$28$、$496$、$8128$ がその最初の四つである。偶数の完全数は Euclid–Euler の定理により $2^p-1$ が素数(Mersenne 素数)となる $p$ に対する $2^{p-1}(2^p-1)$ に限られ、Mersenne 素数と 1 対 1 に対応する。奇数の完全数は一つも知られておらず存在しないことも証明されていないが、存在すれば $N=p^{\alpha}M^2$($p\equiv\alpha\equiv1\pmod 4$)の形をもち、$10^{1500}$ より大きいなど多くの必要条件が知られている。
真の約数の和と $N$ の大小で正の整数は三つに分かれる。$\sigma(N)>2N$ のとき $N$ を過剰数、$\sigma(N)<2N$ のとき不足数といい、完全数はその境目 $\sigma(N)=2N$ にある数である。$1$ は真の約数をもたず $\sigma(1)=1<2$ なので不足数であり、素数 $p$ は $\sigma(p)=p+1<2p$ なので不足数である。完全数と過剰数・不足数の区別は $\sigma(N)/N$ の値が $2$ に等しいか、$2$ より大きいか、小さいかの区別であり、prop-perfect-number-product-criterion のように素因数分解から読み取れる。
完全数は「自分自身を除く約数を全部足すと自分に戻る」数であり、$6=1+2+3$、$28=1+2+4+7+14$ がその最初の二つである。約数の和 $\sigma(N)$ は乗法的関数なので、$\sigma(N)=2N$ という条件は素因数分解を通して各素因数の寄与の積 $\prod\sigma(p^{e})/p^{e}=2$ に書き直せる。偶数の完全数についてはこの条件が完全に解けており、$2^p-1$ が素数となる $p$ ごとに完全数 $2^{p-1}(2^p-1)$ がちょうど一つ対応する(Euclid–Euler の定理)。一方、奇数の完全数は一つも見つかっておらず、存在しないことも証明されていない。奇数の完全数が満たすべき必要条件(形・法 $4$ や $12$ での剰余・素因数の個数・大きさ)は多数知られており、それらが本記事の後半の主題である。
$6$, $28$, $496$, $8128$ は完全数である。実際、
正の整数 $m,n$ が互いに素ならば $\sigma(mn)=\sigma(m)\sigma(n)$ である。
$d$ が $mn$ の正の約数であるとき、$d_1:=\gcd(d,m)$、$d_2:=\gcd(d,n)$ とおく。$m$ と $n$ が互いに素なので $d_1$ と $d_2$ も互いに素であり、$d_1d_2$ は $d$ を割る。また $d$ の任意の素因数 $q$ は $m$ または $n$ の一方だけを割り、たとえば $q\mid m$ なら $q\nmid n$ と $d\mid mn$ から $d$ における $q$ の指数は $m$ における指数以下なので、$d_1=\gcd(d,m)$ における $q$ の指数は $d$ における指数に等しい($q\mid n$ のときは $d_2$ について同様)。よって $d=d_1d_2$ である。逆に $m$ の約数 $d_1$ と $n$ の約数 $d_2$ に対し $d_1d_2$ は $mn$ の約数であり、$\gcd(d_1d_2,m)=d_1$、$\gcd(d_1d_2,n)=d_2$ である。よって $d\mapsto(d_1,d_2)$ は $mn$ の正の約数の集合から $m$ の正の約数と $n$ の正の約数の組の集合への全単射であり、
$$\sigma(mn)=\sum_{d\mid mn}d=\sum_{d_1\mid m}\sum_{d_2\mid n}d_1d_2=\sigma(m)\sigma(n)$$
となる。$\square$
素数 $p$ と正の整数 $e$ に対し
$$\sigma(p^e)=1+p+\cdots+p^e=\frac{p^{e+1}-1}{p-1}$$
であり、$N=p_1^{e_1}\cdots p_r^{e_r}$($p_i$ は相異なる素数、$e_i\ge1$)と素因数分解するとき
$$\frac{\sigma(N)}{N}=\prod_{i=1}^r\frac{1-p_i^{-(e_i+1)}}{1-p_i^{-1}}$$
が成り立つ。したがって $N$ が完全数であることは
$$\prod_{i=1}^r\frac{1-p_i^{-(e_i+1)}}{1-p_i^{-1}}=2$$
と同値である。また各因子は $1<\dfrac{1-p_i^{-(e_i+1)}}{1-p_i^{-1}}<\dfrac{p_i}{p_i-1}$ を満たす。
$p^e$ の正の約数は $1,p,\dots,p^e$ に限るので、$\sigma(p^e)$ は等比数列の和であり、$(p-1)(1+p+\cdots+p^e)=p^{e+1}-1$ から第 1 の等式を得る。lem-perfect-number-sigma-multiplicative を繰り返し用いると $\sigma(N)=\prod_i\sigma(p_i^{e_i})$ であり、
$$\frac{\sigma(p^e)}{p^e}=\frac{p^{e+1}-1}{p^e(p-1)}=\frac{1-p^{-(e+1)}}{1-p^{-1}}$$
から第 2 の等式を得る。完全数の条件 $\sigma(N)=2N$ は $\sigma(N)/N=2$ と同値である。最後に、$1-p^{-(e+1)}$ は $1-p^{-1}$ より大きく $1$ より小さいので、因子は $1$ より大きく $1/(1-p^{-1})=p/(p-1)$ より小さい。$\square$
正の整数 $n\ge2$ について、$2^n-1$ が素数ならば $n$ は素数である。
$n=ab$($a,b\ge2$)と分解できたとすると、$x^b-1=(x-1)(x^{b-1}+\cdots+x+1)$ に $x=2^a$ を代入して $2^n-1=(2^a-1)(2^{a(b-1)}+\cdots+2^a+1)$ となる。$a\ge2$ より $2^a-1\ge3$ であり、第 2 因子は $2^a+1\ge5$ 以上なので、$2^n-1$ は $1$ でも自身でもない約数 $2^a-1$ をもち、素数でない。$\square$
$p$ が素数で $2^p-1$ も素数であるとき、$2^p-1$ を Mersenne素数 という。lem-perfect-number-mersenne-exponent の逆は成り立たない($2^{11}-1=23\cdot89$)。
正の整数 $N$ について、次は同値である。
(2 ⇒ 1) $q:=2^p-1$ が素数であるとする。$q$ は奇数なので $2^{p-1}$ と $q$ は互いに素であり、lem-perfect-number-sigma-multiplicative と prop-perfect-number-product-criterion により
$$\sigma(N)=\sigma(2^{p-1})\sigma(q)=(2^p-1)(q+1)=(2^p-1)2^p=2\cdot2^{p-1}(2^p-1)=2N$$
となる。$p\ge2$ なので $N$ は偶数である。
(1 ⇒ 2) $N$ を偶数の完全数とし、$N=2^{k}M$($k\ge1$、$M$ は奇数)と書く。$2^k$ と $M$ は互いに素なので
$$2^{k+1}M=2N=\sigma(N)=\sigma(2^k)\sigma(M)=(2^{k+1}-1)\sigma(M)$$
である。$2^{k+1}-1$ は奇数なので $2^{k+1}$ と互いに素であり、$2^{k+1}-1$ は $M$ を割る。$M=(2^{k+1}-1)M'$ とおくと、両辺を $2^{k+1}-1$ で割って $\sigma(M)=2^{k+1}M'$ を得る。ここで $M'$ と $M$ は $M$ の相異なる正の約数であり($2^{k+1}-1\ge3$ なので $M'< M$)、
$$M+M'=(2^{k+1}-1)M'+M'=2^{k+1}M'=\sigma(M)$$
である。$\sigma(M)$ は $M$ の正の約数すべての和であり、その中の二つ $M$ と $M'$ だけで和全体に達しているので、$M$ の正の約数は $M$ と $M'$ の二つしかない。正の約数が二つしかない数は素数であり、その約数は $1$ と自身なので $M'=1$、$M=2^{k+1}-1$ は素数である。$p:=k+1$ とおくと lem-perfect-number-mersenne-exponent により $p$ は素数であり、$N=2^{p-1}(2^p-1)$ となる。
一意性:$N=2^{p-1}(2^p-1)$ において $2^p-1$ は奇数なので、$p-1$ は $N$ の素因数分解における $2$ の指数であり、$N$ から定まる。$\square$
$p\mapsto2^{p-1}(2^p-1)$ は、Mersenne 素数 $2^p-1$ を与える素数 $p$ の集合から偶数の完全数の集合への全単射である。特に、偶数の完全数が無限に存在することと Mersenne 素数が無限に存在することとは同値である。
thm-perfect-number-euclid-euler により写像は偶数の完全数の集合への全射であり、一意性の主張により単射である。$\square$
Mersenne 素数が無限に存在するかどうかは未解決である。$2^p-1$ が素数となる $p$ は小さい順に $2,3,5,7,13,17,19,31,\dots$ であり、Great Internet Mersenne Prime Search(GIMPS)の報告によれば、2024 年 10 月の $2^{136279841}-1$ の発見時点で 52 個の Mersenne 素数が知られている(GIMPS24)。したがって同じ個数の偶数の完全数が知られている。偶数の完全数が $2^{p-1}(2^p-1)$ の形をもつことは Euclid『原論』第 9 巻命題 36 にあり、逆にすべての偶数の完全数がこの形であることは Euler が示した(HW08 §16.8、Dic05 Chapter I)。
奇数の完全数は本記事執筆時点(2026 年)で一つも知られておらず、存在しないことも証明されていない(HW08 §16.8、OchemRao2012)。以下は、奇数の完全数が存在するとすれば満たさなければならない条件である。最初のものは Euler による(Dic05 Chapter I)。
$N$ が奇数の完全数ならば、素数 $p$、正の整数 $\alpha$、$p$ で割り切れない正の整数 $M$ により
$$N=p^{\alpha}M^2,\qquad p\equiv\alpha\equiv1\pmod 4$$
と書ける。すなわち、$N$ の素因数分解 $N=p^{\alpha}q_1^{2b_1}\cdots q_s^{2b_s}$ において、指数が奇数の素因数 $p$ はちょうど一つであり、その指数 $\alpha$ と $p$ はともに $4$ を法として $1$ に合同(合同式)である。
$N=p_1^{a_1}\cdots p_r^{a_r}$ と素因数分解する。$N$ は奇数なので各 $p_i$ は奇素数であり、$2N\equiv2\pmod 4$ である。prop-perfect-number-product-criterion により $2N=\sigma(N)=\prod_i\sigma(p_i^{a_i})$ である。
$\sigma(p^a)=1+p+\cdots+p^a$ は $a+1$ 個の奇数の和なので、$a$ が偶数のとき奇数、$a$ が奇数のとき偶数である。積 $\prod_i\sigma(p_i^{a_i})$ は偶数なので、$a_i$ が奇数となる $i$ が少なくとも一つある。もし二つ以上あれば積は $4$ で割り切れ、$2N\equiv2\pmod 4$ に反する。よって指数が奇数の素因数はちょうど一つであり、それを $p=p_1$、$\alpha=a_1$ とおくと、残りの指数はすべて偶数なので $N=p^{\alpha}M^2$($M:=\prod_{i\ge2}p_i^{a_i/2}$、$p\nmid M$)と書ける。さらに $\sigma(p^{\alpha})\cdot(\text{奇数})=2N\equiv2\pmod 4$ であるから、$\sigma(p^{\alpha})$ は $4$ で割り切れず、$\sigma(p^{\alpha})\equiv2\pmod 4$ である。
$p\equiv3\pmod 4$ と仮定すると $p^k\equiv(-1)^k\pmod 4$ なので、$\sigma(p^{\alpha})\equiv\sum_{k=0}^{\alpha}(-1)^k\pmod 4$ となる。$\alpha+1$ は偶数なので右辺は $0$ であり、$\sigma(p^{\alpha})\equiv0\pmod 4$ となって矛盾する。よって $p\equiv1\pmod 4$ である。このとき $p^k\equiv1\pmod 4$ なので $\sigma(p^{\alpha})\equiv\alpha+1\pmod 4$ であり、$\alpha+1\equiv2\pmod 4$、すなわち $\alpha\equiv1\pmod 4$ を得る。$\square$
$N=p^{\alpha}q_1^{2b_1}\cdots q_s^{2b_s}$ を thm-perfect-number-euler-odd-form の形の奇数の完全数とする。このとき $\sigma(p^{\alpha})/2$ は整数であり、$\sigma(p^{\alpha})/2$ および各 $\sigma(q_i^{2b_i})$ はいずれも $N$ を割り切る。さらに $N\equiv1\pmod 4$ である。
prf-perfect-number-euler-odd-form で見たように $\sigma(p^{\alpha})\equiv2\pmod 4$ なので $\sigma(p^{\alpha})/2$ は整数であり、$\sigma(q_i^{2b_i})$ は奇数である。乗法性により
$$N=\frac{\sigma(N)}{2}=\frac{\sigma(p^{\alpha})}{2}\prod_{i=1}^s\sigma(q_i^{2b_i})$$
であり、右辺の各因子は $N$ を割り切る。最後に、$p\equiv1\pmod 4$ より $p^{\alpha}\equiv1\pmod 4$ であり、$M=q_1^{b_1}\cdots q_s^{b_s}$ は奇数なので $M^2\equiv1\pmod 4$ である。よって $N=p^{\alpha}M^2\equiv1\pmod 4$ である。$\square$
$N\equiv1\pmod 4$ は Stern による(Stern1886)。より強い合同条件は Touchard が Jacobi のtheta関数の理論を用いて証明し(Touchard1953)、その後 Satyanarayana(Satyanarayana1959)、Raghavachari(Raghavachari1966)、Holdener(Holdener2002)がそれぞれ独立に初等的な証明を与えた。
$N$ が奇数の完全数ならば $N\equiv1\pmod{12}$ または $N\equiv9\pmod{36}$ である。やや強く、$N\equiv1\pmod{12}$、$N\equiv81\pmod{324}$、$N\equiv117\pmod{468}$ のいずれかが成り立つ。
前半の証明は Holdener2002 に短い初等的証明がある(原証明は Touchard1953)。後半の精密化は Roberts による(Roberts2008)。
thm-perfect-number-euler-odd-form の素因数分解 $N=p^{\alpha}q_1^{2b_1}\cdots q_s^{2b_s}$ に現れる指数 $b_1,\dots,b_s$($M=q_1^{b_1}\cdots q_s^{b_s}$ の指数)にも、多くの制限が知られている。それらの議論の出発点は、$2b_i+1$ の約数が $\sigma(q_i^{2b_i})$ の因数を生み、その因数が prop-perfect-number-odd-divisors により $N$ を割り切るという次の観察である。
素数 $q$ と正の整数 $b$ について、正の整数 $d$ が $2b+1$ を割り切るならば、$\sigma(q^{d-1})$ は $\sigma(q^{2b})$ を割り切る。したがって $N=p^{\alpha}q_1^{2b_1}\cdots q_s^{2b_s}$ を thm-perfect-number-euler-odd-form の形の奇数の完全数とするとき、$d$ が $2b_i+1$ を割り切るならば $\sigma(q_i^{d-1})$ は $N$ を割り切る。特に $b_i\equiv1\pmod 3$ ならば $q_i^2+q_i+1$ は $N$ を割り切る。
要点:$2b+1=de$ とおくと、$q^0,\dots,q^{2b}$ を $d$ 項ずつまとめて $\sigma(q^{2b})=\sigma(q^{d-1})\sum_{l=0}^{e-1}q^{dl}$ と因数分解でき、prop-perfect-number-odd-divisors と合わせて主張が従う。
$2b+1=de$ とおく。$0\le k\le de-1$ を満たす整数 $k$ は、$d$ で割った商と余りにより $k=j+dl$($0\le j\le d-1$、$0\le l\le e-1$)と一意に書けるので
$$\sigma(q^{2b})=\sum_{k=0}^{de-1}q^k=\Bigl(\sum_{j=0}^{d-1}q^j\Bigr)\Bigl(\sum_{l=0}^{e-1}q^{dl}\Bigr)=\sigma(q^{d-1})\sum_{l=0}^{e-1}q^{dl}$$
となり、右辺の和は整数である。よって $\sigma(q^{d-1})\mid\sigma(q^{2b})$ である。prop-perfect-number-odd-divisors により $\sigma(q_i^{2b_i})\mid N$ なので、$\sigma(q_i^{d-1})\mid N$ を得る。最後に、$2$ は $3$ を法として可逆なので、$b_i\equiv1\pmod 3$ と $2b_i+1\equiv0\pmod 3$ は同値であり、$d=3$ とすると $\sigma(q_i^2)=q_i^2+q_i+1$ が $N$ を割り切る。$\square$
たとえば $q=5$、$b=4$ のとき $2b+1=9$ は $d=3$ で割り切れ、$\sigma(5^8)=488281=31\cdot15751$ は $\sigma(5^2)=31$ で割り切れる。第 2 因子は $1+5^3+5^6=15751$ である。
$N=p^{\alpha}q_1^{2b_1}\cdots q_s^{2b_s}$ を thm-perfect-number-euler-odd-form の形の奇数の完全数とする。
thm-perfect-number-odd-exponents の各項と、その他の主な結果を「条件・結論・文献」の形で並べる。定理の 2 以外の証明は各論文に譲る。
| 条件 | 結論 | 文献 |
|---|---|---|
| すべての $b_i=1$ | 起こらない | Steuerwald1937 |
| すべての $b_i=2$ | 起こらない | Kanold1941 |
| $t$ がすべての $2b_i+1$ を割り切る | $t^4\mid N$(定理の 1) | Kanold1941 |
| すべての $b_i\equiv1\pmod 3$ | 起こらない(定理の 2) | McDaniel1970 |
| すべての $b_i=3$ | 起こらない | HagisMcDaniel1972 |
| すべての $b_i=\beta$、$\beta\in\{5,12,17,24,62\}$ | 起こらない | McDanielHagis1975 |
| すべての $b_i\equiv17\pmod{35}$($35\mid2b_i+1$) | 起こらない | McDanielHagis1975 |
| すべての $b_i=\beta$、$\beta\in\{6,8,11,14,18\}$ | 起こらない | CohenWilliams1985 |
| 各 $b_i\equiv1\pmod 3$ または $b_i\equiv2\pmod 5$ | 最小の素因数 $r$ は $10^8 < r < 10^{1000}$(定理の 3) | FletcherNielsenOchem2012 |
| すべての $b_i=\beta$ | $\omega(N)\le2\beta^2+6\beta+3$(定理の 4) | Yamada2021 |
$N$ がどの素数の $6$ 乗でも割り切れないときは各 $b_i$ が $1$ か $2$ なので定理の 3 の仮定が成り立ち、この場合には $5\nmid N$ であることも示されている(FletcherNielsenOchem2012)。すべての $b_i$ が等しい値 $\beta$ をとる場合については、以上を含む結果から $\beta\ge9$、$\beta\not\equiv1\pmod 3$、$\beta\not\equiv2\pmod 5$ が知られている(Yamada2021 の序論にまとめがある)。この形の奇数の完全数は存在しないと予想されている(McDanielHagis1975)。これらの証明の多くは、定理の 1 を出発点として、lem-perfect-number-odd-exponent-divisor のように指数から生じる $\sigma(q_i^{d-1})$ の素因数を次々に追い、ある素数の冪が $N$ を割り切る回数の上限に矛盾させるという形をとる(McDanielHagis1975)。定理の 2 は定理の 1 と 3 から導け、同じ議論ですべての $b_i\equiv2\pmod 5$ とはならないことも従う。
すべての $b_i\equiv1\pmod 3$ と仮定する。$2$ は $3$ を法として可逆なので $b_i\equiv1\pmod 3$ と $3\mid2b_i+1$ は同値であり、$3$ はすべての $2b_i+1$ を割り切る。定理の 1 により $3^4\mid N$ となり、$N$ の最小の素因数 $r$ は $r\le3$ である。一方、各 $b_i$ は $b_i\equiv1\pmod 3$ を満たすので定理の 3 の仮定が成り立ち、$r>10^8$ となって矛盾する。同様に、$2$ は $5$ を法として可逆なので $b_i\equiv2\pmod 5$ と $5\mid2b_i+1$ は同値であり、すべての $b_i\equiv2\pmod 5$ とすると定理の 1 により $5^4\mid N$、定理の 3 により $r>10^8$ となって矛盾するから、すべての $b_i\equiv2\pmod 5$ とはならない。この導き方は定理の 3 に依存しており、McDaniel1970 の証明とは別のものである。
定理の 1 だけから、次のような制限がすぐに従う(この系は定理の 1 に依存し、定理の 1 の証明は Kanold1941 に譲る)。
$N=p^{\alpha}q_1^{2b_1}\cdots q_s^{2b_s}$ を thm-perfect-number-euler-odd-form の形の奇数の完全数とする。すべての $2b_i+1$ を割り切る整数 $t>1$ について、$\sigma(t) < 2t$ である。
要点:定理の 1 により $t$ は $N$ の真の約数であり、$\sigma(n)/n=\sum_{d\mid n}1/d$ は約数が増えると大きくなるので $\sigma(t)/t < \sigma(N)/N=2$ となる。
正の整数 $n$ について、$d\mapsto n/d$ は $n$ の正の約数の集合の全単射なので
$$\frac{\sigma(n)}{n}=\sum_{d\mid n}\frac{d}{n}=\sum_{d\mid n}\frac1d$$
である。$t>1$ がすべての $2b_i+1$ を割り切るとすると、thm-perfect-number-odd-exponents の 1 により $t^4\mid N$、特に $t\mid N$ であり、$t< t^4\le N$ である。$t$ の正の約数はすべて $N$ の正の約数であり、$N$ 自身は $t$ の約数でないので
$$2=\frac{\sigma(N)}{N}=\sum_{d\mid N}\frac1d>\sum_{d\mid t}\frac1d=\frac{\sigma(t)}{t}$$
となる。よって $\sigma(t) < 2t$ である。$\square$
たとえば $945=3^3\cdot5\cdot7$ は $\sigma(945)=40\cdot6\cdot8=1920>1890$ を満たす奇数なので、すべての $2b_i+1$ が $945$ で割り切れる(すべての $b_i\equiv472\pmod{945}$ となる)奇数の完全数は存在しない。
$N$ を奇数の完全数とし、$\omega(N)$ を $N$ の相異なる素因数の個数、$\Omega(N)$ を重複を込めた素因数の個数(素因数の個数)とする。
1 の $\omega(N)\ge10$ は Nielsen2015、$3\nmid N$ のときの $\omega(N)\ge12$ は Nielsen2007 による。2 は OchemRao2012、3 は OchemRao2014 による。4 は Clayton–Hansen(ClaytonHansen2023)によるもので、Zelinsky による $3\nmid N$ のとき $\Omega(N)\ge(8\omega(N)-7)/3$、$3\mid N$ のとき $\Omega(N)\ge(21\omega(N)-39)/8$(Zelinsky2018)、および $3\nmid N$ のとき $\Omega(N)\ge(302\omega(N)-286)/113$、$3\mid N$ のとき $\Omega(N)\ge(66/25)\omega(N)-5$(Zelinsky2021)を改良したものである。いずれも証明は各論文に譲る。
1 は Ochem–Rao による計算機を用いた結果である(OchemRao2012)。2 は Dickson の定理であり、証明は Dickson1913 に譲る。奇数の完全数が有限個か無限個かは(そもそも存在するかを含めて)未解決である。
完全数の研究は古代ギリシアに遡り、$6,28,496,8128$ は古くから知られていた(Dic05 Chapter I)。完全数の条件 $\sigma(N)=2N$ を $\sigma(N)=kN$ に一般化した数はk倍完全数と呼ばれ、また二つの数が互いに他方の真の約数の和になっているものは友愛数と呼ばれる。
Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する