No. 22≫ No.23 ≫No. 24
s_hskz
2016/01/21 11:44
>No.22 lengthさん
◆楽屋裏
307番目の素数が2017だと昨年末に気がついておりました。
今回の御出題を拝見して、すぐに、2016 を作るためには、2017から1を引けるといいなあと思いました。ですが、合成関数を作って自然数から1を引く真っ当な手段がみあたりませんでした。
先日、苦し紛れにMathematicaのドキュメントを読み、関数を片端から漁っていたところ、Digammaの無限遠点での級数展開の表記に目がとまり、ふと、こいつをExp[]に喰わせたらどうなのかと、何かが降りてきました。{※1}
さっそくWolframAlphaにて、Exp[Digamma[x]]を確認し、以下の通りとなりました。
(1) 正の範囲でグラフが一直線にかなり近い。
(2) x→∞ での級数展開が次のようになっていて感じがいい。
x - 1/2 + O(1/x)
…… うまくするとこれ、【ほぼ1/2を引くことになってないかなw】
……まあ無限遠点でのハナシだし…… Exp[Digamma[2]] は どうかな?→ 1.526... イイカンジ!!
===
《《《仮説登場》》》
n = Floor[Exp[Digamma[ n+1 ]]]
というわけでして、合成関数にて、自然数から1を引ける候補が生まれました。
今回の御出題における、ひとつの解としては、たとえば、2017から1を引けるかなど、ごく少数の例で減算がうまくできることを確認できればオッケーと腹をくくり、回答案をしたためました。任意の自然数でうまくいくか証明はしていませんが、ほぼ大丈夫??
●※1
Log[1+x]のx=無限での級数展開と、Digamma[x]のそれとで、主要項がなんとなく似ている……程度の野蛮な観察からです。
s_hskz 2016/01/21 11:44
>No.22 lengthさん
◆楽屋裏
307番目の素数が2017だと昨年末に気がついておりました。
今回の御出題を拝見して、すぐに、2016 を作るためには、2017から1を引けるといいなあと思いました。ですが、合成関数を作って自然数から1を引く真っ当な手段がみあたりませんでした。
先日、苦し紛れにMathematicaのドキュメントを読み、関数を片端から漁っていたところ、Digammaの無限遠点での級数展開の表記に目がとまり、ふと、こいつをExp[]に喰わせたらどうなのかと、何かが降りてきました。{※1}
さっそくWolframAlphaにて、Exp[Digamma[x]]を確認し、以下の通りとなりました。
(1) 正の範囲でグラフが一直線にかなり近い。
(2) x→∞ での級数展開が次のようになっていて感じがいい。
x - 1/2 + O(1/x)
…… うまくするとこれ、【ほぼ1/2を引くことになってないかなw】
……まあ無限遠点でのハナシだし…… Exp[Digamma[2]] は どうかな?→ 1.526... イイカンジ!!
===
《《《仮説登場》》》
n = Floor[Exp[Digamma[ n+1 ]]]
というわけでして、合成関数にて、自然数から1を引ける候補が生まれました。
今回の御出題における、ひとつの解としては、たとえば、2017から1を引けるかなど、ごく少数の例で減算がうまくできることを確認できればオッケーと腹をくくり、回答案をしたためました。任意の自然数でうまくいくか証明はしていませんが、ほぼ大丈夫??
●※1
Log[1+x]のx=無限での級数展開と、Digamma[x]のそれとで、主要項がなんとなく似ている……程度の野蛮な観察からです。