「利用者:Nayuta Ito/素微分友愛数の研究の進展/独自研究」の版間の差分

提供: 数学を愛する会Wiki
ナビゲーションに移動 検索に移動
147行目: 147行目:
  
 
よって\( 30 \)番目の素数階乗まで調べればよい。わんだほーい!!!
 
よって\( 30 \)番目の素数階乗まで調べればよい。わんだほーい!!!
 +
 +
===実験===
 +
 +
====コード====
 +
\( (p_n#)' \)を計算するPython3のコード
 +
 +
<pre>
 +
# coding: utf-8
 +
# Your code here!
 +
 +
import sympy
 +
 +
n = 10
 +
 +
primes = []
 +
 +
for i in range(n*n):
 +
    if sympy.isprime(i):
 +
        primes.append(i)
 +
    if len(primes) == n:
 +
        break
 +
 +
total = 0
 +
 +
for p1 in primes:
 +
    term = 1
 +
    for p2 in primes:
 +
        if p1 != p2:
 +
            term *= p2
 +
    total += term
 +
 +
print(total)
 +
print(sympy.factorint(total))
 +
</pre>
 +
 +
====判定法====
 +
 +
\( 30 \)番目までの素数の逆数和が\( 2 \)に満たないこと、\( \frac{1}{7}+\frac{1}{11}+\cdots+\frac{1}{23} < 0.5 \)であることから、\( 5# \)以降について、素因数が6個以下であれば素微分友愛数ではないことがわかる。
 +
 +
====結果====
 +
 +
# (2#)'=1=1
 +
# (3#)'=5=5
 +
# (5#)'=31=31
 +
# (7#)'=247=13×19
 +
# (11#)'=2927=2927
 +
# (13#)'=40361=40361
 +
# (17#)'=716167=19×37693
 +
# (19#)'=14117683=37×381559
 +
# (23#)'=334406399=43×163×47711
 +
# (29#)'=9920878441=269×36880589
 +
# (31#)'=314016924901=61×4253×1210397
 +
# (37#)'=11819186711467=79×82913×1804421
 +
# (41#)'=492007393304957=2683561×183341237
 +
# (43#)'=21460568175640361=47×456607833524263
 +
# (47#)'=1021729465586766997=1559×655374897746483
 +
# (53#)'=54766551458687142251=1701437×32188409831623
 +
# (59#)'=3263815694539731437539=1087×3002590335363138397
 +
# (61#)'=201015517717077830328949=201015517717077830328949
 +
# (67#)'=13585328068403621603022853=13585328068403621603022853
 +
# (71#)'=972416614407737400870501653=653×24310877×61254579539162813
 +
# (73#)'=71544353681891529224514036059=79×905624730150525686386253621
 +
# (79#)'=5692733621468679832887230172131=5692733621468679832887230172131
 +
# (83#)'=475714535349241099037539188841003=6991×17117×31727×125299823719985846887
 +
# (89#)'=42605658161771733665696611824842057=128427649157×331748330218886478419701
 +
# (97#)'=4156517583588203716343221884611037839=11887×349669183443106226663011851990497
 +
 +
26以降は後で書く

2022年5月14日 (土) 18:54時点における版

\( p_n \)で\( n \)番目の素数を表す。また、\( \log \)は全て自然対数である。

予想: 素数の2倍が素微分友愛数となることはない

pを素数とし、2pが素微分友愛数であるとする

(2p)'=p+2である

p+2が約数がすごく多くて微分すると2pに戻ってくる

pを6で割った余りで場合分けする

p=6a+1のとき

p+2=6a+3=3(2a+1)なのでこれは3の倍数になる

2a+1=bとおくと2p=(2p)''=(2a+1)+3b'となる

これが偶数ということはb'は奇数

b'=2c+1とすると

a+3c+2=pである

p=6a+1と仮定したのでa%3=2

a=3d+2とする(このときp=18d+13である)

3d+3c+4=pである

このとき明らかにc+dは奇数

p=6a-1のとき

p+2=6a+1は奇数

予想: 素数の5倍が素微分友愛数となることはない

pを素数とし、5pが素微分友愛数であるとする

(5p)'=p+5である

p+5が約数がすごく多くて微分すると5pに戻ってくる

かなりえぐい

pを5で割った余りで場合分けする

p=6a+1のとき

p+5=6a+6=2・3(a+1)

(p+2)'=3(a+1)+2(a+1)+6(a+1)'=5(a+1)+6(a+1)'=5a+5+6(a+1)'

すなわち6(a+1)'=a+1である

a+1は6の倍数なのでa=6b-1とする

このとき(a+1)'=(6b)'=5b+6b'であるから

6(a+1)'=30b+36b'>6b=a

よって矛盾が生じた

p=6a-1のとき

p+5=2(3a+2)は偶数

3a+2=bとおくと(5p)''=(3a+2)+2b'となる

左辺は奇数なのでaは奇数

左辺は3で割って1余るのでb'は3で割って1余る

a=2c+1とおく(このときp=12c+5である)

5p=(5p)''=(6c+5)+2b'

mod5で考えるとc+2b'は5の倍数

予想: 素数階乗は素微分友愛数にならない

素数定理

Rosser's Theoremによると、\( n \)番目の素数は\( n \log{n} \)より大きい。

また、\( n \geq 6 \)のとき、\( n \)番目の素数は\( n (\log{n} + \log{\log{n}}) \)より小さいことが知られている。

\( n (\log{n} + \log{\log{n}}) < 2n \log{n} \)であることは容易にわかる。また、数値計算と合わせて\( p_n < 2n \log{n} \)が\( n \geq 3 \)で成り立つことが分かる。

証明

10番目(29#)までの素数階乗は実験により素微分友愛数でないことが確かめられる

\( n \geq 11 \)とする

N=(n番目の素数階乗)について考えその大きさを評価する

\( N < (2n \log{n})^n \)である

\( \begin{align*} & N' \\ <& \frac{1}{2} \cdot (2n \log{n})^n \cdot n \\ =& 2^{n-1} \cdot (\log{n})^n \cdot n^{n+1} \\ =& n^{n+1 + (n-1)\log_n{2} + n\log_n{\log{n}}} \\ <& n^{n+1 + n\log_n{(2\log{n})}} \\ <& n^{n+1 + n\log_n{(2n^{0.37})}} (\because \log{n} \leq n^{\frac{1}{e}}) \\ <& n^{n+1 + 0.37n + n\log_n{2}} \\ <& n^{n+0.1n + 0.37n + 0.29n} \\ <& n^{1.76n} \end{align*} \)

\( N' \)の素因数は\( n \)より大きいので、素因数は最大でも\( 1.76n \)個しか持たない

よって\( N' \)の素因数の逆数の総和を\( T \)とすると

\( 0 < T < \frac{1.76n}{n\log{n}} = \frac{1.76}{\log{n}} \)

が成り立つ

また\( N \)の素因数の逆数の総和を\( S \)とすると

\( \begin{align*} & S \\ <& \frac{1}{2} + \frac{1}{3} + \cdots + \frac{1}{p_n} \\ \leq & \sum_{k=1}^{10}\frac{1}{p_k} + \sum_{k=11}^n \frac{1}{n\log{n}} \\ <& \sum_{k=1}^{10}\frac{1}{p_k} + \int_{10}^n \frac{1}{x\log{x}} dx (\because \frac{1}{x\log{x}}は単調減少) \\ =& \sum_{k=1}^{10}\frac{1}{p_k} + \log{\log{n}} - \log{\log{10}} \\ <& \log{\log{n}} + 0.7 \\ \end{align*} \)

であるから

\( S \cdot T < \frac{1.76}{\log{n}}(\log{\log{n}} + 0.7) \)

である。

\( f(x) = \frac{1.76}{\log{x}}(\log{\log{x}} + 0.7) \)

とおくと、

\( \frac{d}{dx}f(x) = - \frac{2(\log{\log{x}})}{x(\log{x})^2} \)

であるから、\( x > e \)で\( f(x) \)は単調減少である。

数値計算によると\( x \geq 30 \)で\( f(x) < 1 \)となるから矛盾。

よって\( 30 \)番目の素数階乗まで調べればよい。わんだほーい!!!

実験

コード

\( (p_n#)' \)を計算するPython3のコード

# coding: utf-8
# Your code here!

import sympy

n = 10

primes = []

for i in range(n*n):
    if sympy.isprime(i):
        primes.append(i)
    if len(primes) == n:
        break

total = 0

for p1 in primes:
    term = 1
    for p2 in primes:
        if p1 != p2:
            term *= p2
    total += term

print(total)
print(sympy.factorint(total))

判定法

\( 30 \)番目までの素数の逆数和が\( 2 \)に満たないこと、\( \frac{1}{7}+\frac{1}{11}+\cdots+\frac{1}{23} < 0.5 \)であることから、\( 5# \)以降について、素因数が6個以下であれば素微分友愛数ではないことがわかる。

結果

  1. (2#)'=1=1
  2. (3#)'=5=5
  3. (5#)'=31=31
  4. (7#)'=247=13×19
  5. (11#)'=2927=2927
  6. (13#)'=40361=40361
  7. (17#)'=716167=19×37693
  8. (19#)'=14117683=37×381559
  9. (23#)'=334406399=43×163×47711
  10. (29#)'=9920878441=269×36880589
  11. (31#)'=314016924901=61×4253×1210397
  12. (37#)'=11819186711467=79×82913×1804421
  13. (41#)'=492007393304957=2683561×183341237
  14. (43#)'=21460568175640361=47×456607833524263
  15. (47#)'=1021729465586766997=1559×655374897746483
  16. (53#)'=54766551458687142251=1701437×32188409831623
  17. (59#)'=3263815694539731437539=1087×3002590335363138397
  18. (61#)'=201015517717077830328949=201015517717077830328949
  19. (67#)'=13585328068403621603022853=13585328068403621603022853
  20. (71#)'=972416614407737400870501653=653×24310877×61254579539162813
  21. (73#)'=71544353681891529224514036059=79×905624730150525686386253621
  22. (79#)'=5692733621468679832887230172131=5692733621468679832887230172131
  23. (83#)'=475714535349241099037539188841003=6991×17117×31727×125299823719985846887
  24. (89#)'=42605658161771733665696611824842057=128427649157×331748330218886478419701
  25. (97#)'=4156517583588203716343221884611037839=11887×349669183443106226663011851990497

26以降は後で書く