今回は,キモすぎる公式,Ramanujan's master theorem(以下,ジャンマス,RMT)との出会いと,自力でガバ証明を与えるまでの葛藤を綴っていこうと思う.
RMTって略され方は英語圏の文章で見かける気がするけど,ジャンマスは勝手に俺がそう呼んでるだけだから注意されたし.
1年前ほど前(2022年5月),はてなで油そばを食いながら英wiki見漁っていたときに次の公式を発掘した.
$$ \int_0^\infty x^{s-1}\sum_{n=0}^\infty f(n)\frac{(-1)^n}{n!}x^n\mathrm dx =\Gamma(s)f(-s) $$
ラマヌジャンらしいキモい式だなと思った.一体何を食ったらこんな式思いつくんだ?
$f(x)=1$のとき,当然,和の部分は$e^{-x}$に等しく,$\Gamma$関数の定義式となるから,成り立ってそう.
俺は直感的に,この式の和の部分は$\Gamma$関数に持って行くために$e^{-x}$のマクローリン展開に重みづけをしているのだと感じた.
同時に,これを$\Gamma$関数ではなく,メリン変換ととらえると,$e^{-x}$のマクローリン展開のみならず,一般の関数のマクローリン展開の重みづけでも成り立ってそうと予想した.
$$ \int_0^\infty x^{s-1}\sum_{n=0}^\infty g(n)\frac{f^{(n)}(0)}{n!}x^n\mathrm dx =\mathcal M[f](s)g(-s) $$
右辺を思いついた経緯だが,普通のジャンマスが
$$ \int_0^\infty x^{s-1}\sum_{n=0}^\infty g(n)\frac{(-1)^n}{n!}x^n\mathrm dx =\cdots =\int_0^\infty x^{s-1}f(-s)\sum_{n=0}^\infty \frac{(-1)^n}{n!}x^n\mathrm dx =f(-s)\int_0^\infty x^{s-1}e^{-x}\mathrm dx =\Gamma(s)f(-s) $$
と示されるはずだから,
$$ \int_0^\infty x^{s-1}\sum_{n=0}^\infty g(n)\frac{f^{(n)}(0)}{n!}x^n\mathrm dx =\cdots =\int_0^\infty x^{s-1}g(-s)\sum_{n=0}^\infty \frac{f^{(n)}(0)}{n!}x^n\mathrm dx =g(-s)\int_0^\infty x^{s-1}f(x) dx =\mathcal M[f](s)g(-s) $$
となると予想した.実際,一般化の式で$f(x)=e^{-x}$としたら通常のジャンマスと一致している.
色々調べたら実際にこの直感は正しいことが分かったが,問題なのは,通常のジャンマスでさえ証明がどこにも載っていなかったことだ.
面倒なので,長い間放置していたが,自分の力ではこれ以上探せないため,自力で証明してみることにした.
予想式において,最初は$g(n)$をマクローリン展開して,和の順序を入れ替えるなどしたらうまくいきそうだと考えた.
実際はそうはいかなかったが,$n^k$という項が出てくるから,$x\longmapsto e^{-x}$の置換と部分積分で消すことを考えた.そして,部分積分の副産物で上手いこと$(-s)^k$がでてきた.総武線で中野駅を目指しているときだった.
ラマヌジャンはここからジャンマスを思いついたのだろうと感じた.
しかし,部分積分の計算過程で発散してうまく計算ができない箇所があった.
計算がうまくいかなかったが,すでに証明の一歩手前まで来ていたので,アプローチは変えなかった.
そこで,我々非数学屋の特権に手を出すことにした.に$\delta$関数を使ってしまおうという魂胆だ.
都合よく積分区間が$(-\infty,\infty)$(すなわち,原点を通っている)だったので微分に$\delta$関数を絡めてみたらうまく計算できた.ジムでバイクを漕ぎながらの出来事だった.思わずEureka!と叫んでしまいそうになった.
数学屋に怒られそうだったが,知見が欲しかったため,Twitterにこの証明を投稿した.悲しいことに,1つもいいねが付かなかった.
その証明が次である.
$$\begin{align}
n^k\int_{-\infty}^\infty e^{-sx}\cdot e^{-nx}\mathrm dx
&=n^k\int_{-\infty}^\infty e^{-sx}\int_{-\infty}^\infty\delta(y-x)e^{-ny}\mathrm dy\mathrm dx\\
&=\int_{-\infty}^\infty e^{-sx}\int_{-\infty}^\infty\delta(y-x)(-1)^k\frac{\mathrm d^ke^{-ny}}{\mathrm dy^k}\mathrm dy\mathrm dx\\
&=\int_{-\infty}^\infty e^{-sx}\int_{-\infty}^\infty\frac{\mathrm d^k\delta(y-x)}{\mathrm dy^k}e^{-ny}\mathrm dy\mathrm dx\\
&=\int_{-\infty}^\infty e^{-sx}\int_{-\infty}^\infty(-1)^k\frac{\mathrm d^k\delta(y-x)}{\mathrm dx^k}e^{-ny}\mathrm dy\mathrm dx\\
&=\int_{-\infty}^\infty e^{-ny}\int_{-\infty}^\infty\frac{\mathrm d^k\delta(y-x)}{\mathrm dx^k}e^{-sx}\mathrm dx\mathrm dy\\
&=\int_{-\infty}^\infty e^{-ny}\int_{-\infty}^\infty\delta(y-x)\frac{\mathrm d^ke^{-sx}}{\mathrm dx^k}\mathrm dx\mathrm dy\\
&=\int_{-\infty}^\infty e^{-ny}\int_{-\infty}^\infty\delta(y-x)(-s)^ke^{-sx}\mathrm dx\mathrm dy\\
&=(-s)^k\int_{-\infty}^\infty e^{-sy}\cdot e^{-ny}\mathrm dy
\end{align}$$
より,
$$\begin{align}
\therefore\int_0^\infty x^{s-1}\sum_{n=0}^\infty g(n)\frac{f^{(n)}(0)}{n!}x^n\mathrm dx
&=\int_0^\infty x^{s-1}\sum_{n,k=0}^\infty \frac{g^{(k)}(0)}{k!}n^k\frac{f^{(n)}(0)}{n!}x^n\mathrm dx\\
&=\int_{-\infty}^\infty e^{-sx}\sum_{n,k=0}^\infty \frac{g^{(k)}(0)}{k!}n^k\frac{f^{(n)}(0)}{n!}e^{-nx}\mathrm dx \ (x\longmapsto e^{-x})\\
&=\int_{-\infty}^\infty e^{-sx}\sum_{n,k=0}^\infty \frac{g^{(k)}(0)}{k!}(-s)^k\frac{f^{(n)}(0)}{n!}e^{-nx}\mathrm dx \\
&=\int_{-\infty}^\infty e^{-sx}\sum_{n=0}^\infty g(-s)\frac{f^{(n)}(0)}{n!}e^{-nx}\mathrm dx \\
&=\int_0^\infty x^{s-1}\sum_{n=0}^\infty g(-s)\frac{f^{(n)}(0)}{n!}x^n\mathrm dx\\
&=\int_0^\infty x^{s-1}f(x)\mathrm dx \ g(-s)\\
&=\mathcal M[f](s)g(-s).
\end{align}$$
ジャンマスそのものを用いた証明.
Riemann-Liouville 積分において,$f(x)$のメリン変換の収束条件から$a\to \infty$として,
$$\begin{align} \int_0^\infty x^{s-1}\sum_{n=0}^\infty g(n)\frac{f^{(n)}(0)}{n!}x^n\mathrm dx &=\Gamma(s)g(-s)f^{(-s)}(0)(-1)^{-s}\\ &=\Gamma(s)g(-s)\frac{1}{\Gamma(s)}\int_\infty^0 (-1)^{-s}(0-t)^{s-1}f(t)\mathrm dt\\ &=g(-s)\int_0^\infty t^{s-1}f(t)\mathrm dt\\ &=\mathcal M[f](s)g(-s). \end{align}$$
ジャンマスを既知のものとした証明.でないと循環論法になる.
先のツイートをする直前に反復積分に関する問題を解いていたので,そこから閃いた.
簡単のため$F(s)=\mathcal M[f(x)](s)$とおくと,メリン変換の倍公式より,$\mathcal M[f(a^kx)](s)=a^{-sk}F(s)$.
この両辺に$h^{(k)}(0)/k!$を掛けて$k=0,1,...$として辺々加えると,
$$\begin{align}
・(\mathrm{左辺})
&=\int_0^\infty x^{s-1}\sum_{n,k=0}^\infty \frac{h^{(k)}(0)}{k!}\frac{f^{(n)}(0)}{n!}a^{kn}x^n\mathrm dx\\
&=\int_0^\infty x^{s-1}\sum_{n=0}^\infty h(a^n)\frac{f^{(n)}(0)}{n!}x^n\mathrm dx\\
\end{align}$$
$$\begin{align}
・(\mathrm{右辺})
&=\sum_{k=0}^\infty \frac{h^{(k)}(0)}{k!}a^{-sk}F(s)\\
&=F(s)h(a^{-s})
\end{align}$$
$g(x)=h(a^x)$として,
$$ \therefore\int_0^\infty x^{s-1}\sum_{n=0}^\infty g(n)\frac{f^{(n)}(0)}{n!}x^n\mathrm dx =\mathcal M[f](s)g(-s).$$
任意の関数$g(x)$を$h(a^x)$と表せるかは微妙.