妖精現実



最新記事 妖精の森 ♌ ペル方程式の夏(2020-12-27)

[ 新着記事 | 数学・プログラミング | 天文・暦
| シリア語・Unicode・詩 | ジョーク | 漫画・アニメ
| 字幕 | 哲学・ファンタジー | 全記事 | チラ裏 ]

チラ裏

2021-03-05 ∫sec x dx (その3)

【9】 −π/2 < θ < π/2 の範囲において t = sec θ の原始関数が
  t = arsinh (tan θ) 【A】
であることを確かめたが、それは
  sinh t = tan θ 【ハ】
  cosh t = sec θ 【ヒ】
の【ハ】を使ったもので、表現の一例にすぎない。【ヒ】の両辺の arcosh を取れば
  t = arcosh (sec θ) 【B】
とも書けるし、【ハ】を【ヒ】で辺々割れば
  tanh t = sin θ 【フ】
  t = artanh (sin θ) 【C】
とも書ける(【フ】の右辺は、tan θ / sec θ = (sin θ/cos θ) / (1/cos θ) = sin θ / 1 による)。

さらに、例えば【ハ】の両辺の逆数を考えれば
  csch t = cot θ
  t = arcsch (cot θ) 【D】
のような表現も可能だし、同様に【ヒ】【フ】の逆数から、次の式が生じる:
  t = arsech (cos θ) 【E】
  t = arcoth (csc θ) 【F】

これら6種類のどれも、本質的に同じ意味の逆グーデルマン関数であり、従って、多かれ少なかれ sec の原始関数を表している。そのうち【A】については、前回一通り見た。最初の作図の点 S のy座標を利用したものだった。

一方、【B】は点 S のx座標を利用したものだが、このx座標は常に正で(作図から一目瞭然…双曲線の右の枝を見ているのだから、当然そうなる)、θ の正負の情報が失われてしまう(0 ≤ θ < π/2 の範囲に限り、そのまま有効)。言い換えると、cosh は1対1の関数でないため、arcosh には主値の問題が付きまとう。面倒なので後回しにして、今回は【C】を考えてみたい。

【10】 artanh の対数表現。y = tanh t = sinh t / cosh t = (et − e−t) / (et + e−t) について、A = et と置くと:
  y = (A − 1/A) / (A + 1/A) = (A2 − 1) / (A2 + 1)
  分母を払って y(A2 + 1) = A2 − 1, A2 − 1 = A2y + y
  A2(1 − y) = 1 + y
  A2 = (1 + y) / (1 − y)
y = tanh t の値域は (−1, 1) なので上記右辺は正、A = et も常に正、従って:
  A = et = [(1 + y) / (1 − y)]1/2
  t = log [(1 + y) / (1 − y)]1/2 = (1/2) log [(1 + y) / (1 − y)]

上の式が t = artanh y の対数表現。【C】は y = sin θ に当たる:
  t = artanh (sin θ) = (1/2) log [(1 + sin θ) / (1 − sin θ)] 【ヘ】

【A】経由で得られる
  log | sec θ + tan θ | 【ホ】
と比べると全然形は異なるが、【ホ】【ヘ】は、それぞれ、グーデルマン関数の2通りの表記【A】【C】の対数表現に当たる(従って、どちらも sec θ の原始関数)。三角関数の範囲で【ヘ】⇔【ホ】を直接変換することもできるが、例によって、ゴチャゴチャする。【ヘ】自体、シンプルな【C】と比べると、複雑で分かりにくい。

【ヘ】の log の引数に絶対値記号を付けている資料もあるが、(1 + sin θ) / (1 − sin θ) は負にならないので、絶対値記号は(無害だが)必要ない。

実は【C】バージョン t = artanh (sin θ) は、そのままの形で −π/2 < θ < π/2 の範囲外でも有効…という素晴らしい性質を持つ。sec の原始関数の決定版といえるかもしれない。

【11】 念のため【C】の導関数が sec θ であることを、直接的に確かめておこう。y = artanh x の導関数を得るため、まず x = tanh y を y について微分:
  dx/dy = (sinh y / cosh y)′ = (cosh y cosh y − sinh y sinh y) / cosh2 y
  = 1 / cosh2 y = 1 − tanh2 y
  (最後の等号は、恒等式 1 = cosh2 y − sinh2 y の両辺を cosh2 y で割ったもの)
  従って dy/dx = 1 / (1 − tanh2 y) = 1 / (1 − x2) 【マ】

y = artanh x の導関数は【マ】になることが分かった。従って【C】の導関数は:
  d[artanh (sin θ)] / dθ = [1 / (1 − sin2 θ)] cos θ
  = (1 / cos2 θ) cos θ = 1 / cos θ = sec θ

逆に言えば、sec θ を積分すると【C】になり、それを書き換えると【ヘ】になる(積分定数省略)。

2021-03-04 ∫sec x dx (その2)

【8】 グーデルマン関数 θ = arctan (sinh t) の逆関数
  t = arsinh (tan θ) 【ナ】
は、sec θ の原始関数。それを確かめるため、前回は素直に、逆関数の微分法を使った。今回は【ナ】を直接 θ で微分してみる。

【ナ】を微分するには arsinh を微分しなければならない。y = arsinh x の導関数は何だろうか。x = sinh y, dx/dy = cosh y なので:
  dy/dx = 1 / cosh y 【ニ】
【ニ】の右辺を x の式にすれば、その答えが得られる。そのためには、分母の cosh y を sinh y の式に変換すればいい(sinh y = x だから)。
  恒等式 cosh2 y − sinh2 y = 1 から
  cosh2 y = 1 + sinh2 y
  cosh y = ±1 + sinh2 y
実変数の cosh は常に正なので、複号のプラスを選んで【ニ】に代入すると:
  dy/dx = 1 / 1 + sinh2 y = 1 / 1 + x2 【ヌ】

y = arsinh x の導関数は【ヌ】になることが分かった。

さっそく【ヌ】を使って本題の【ナ】を θ で微分する。(tan θ)′ = sec2 θ なので*2:
  (arsinh (tan θ))′ = (1 / 1 + tan2 θ) sec2 θ = (1 / sec2 θ) sec2 θ 【ネ】
考えている範囲 −π/2 < θ < π/2 では sec は正なので、符号の心配なく【ネ】の根号を外せる:
  = (1 / sec θ) sec2 θ = sec θ

これで【ナ】を微分すると sec θ になること、すなわち sec θ の原始関数が(積分定数を無視して)【ナ】であることが直接、示された!

このやり方の方が、前回の方法より少し簡単かもしれない。

*2 証明は、(tan θ)′ = (sin θ / cos θ)′ に商の微分法を適用するだけ。
  [cos θ cos θ − sin θ (−sin θ)] / cos2 θ = 1 / cos2 θ = sec2 θ

付記 グーデルマン関数の逆関数を使って sec の積分を考える方法は、比較的マイナーとはいえ、ネット上で検索すれば、有益な資料が見つかる。けれど【ナ】のバリエーションとして arcosh (sec θ) を採用しているテキストもある(例=ウィキペディアのロシア語版、Wolf Holzmann のPDFファイル)。arcosh の主値は負の値を取れないので、arcosh を使うなら θ が負の場合、arcosh を −arcosh に置き換えるべきだろう。

<補足> 通常のグーデルマン関数では、値域が −π/2 < θ < π/2 に制限されている。その逆関数【ナ】の定義域も、当然その範囲。従って、今のところ、その範囲では sec θ の原始関数は【ナ】…としか言えない。その範囲をはみ出して、例えば π/2 をまたぐような区間で sec を積分することはできないので(θ = π/2 では sec θ は定義されない)、その意味では、この制限には何の実害もない。一方、不連続点をまたがない形で、上記範囲外の積分区間を選ぶことは可能であり、それについては別途考える必要がある。arsinh 表現の場合、θ が任意だと、本文の「符号の心配なく根号を外せる」が成り立たなくなり、符号の調整が必要。これは、上の「付記」とは別の問題。

2021-03-03 ∫sec x dx グーデルマニアンの威力

sec x つまり 1 / cos x の積分は、いろいろな形で書き表されるものの、伝統的・教科書的な方法では、どれもこれも場当たり的・天下り的だった。グーデルマン関数の立場から眺めると、全てが明快に整理される。必要なのは、微積の初歩と好奇心だけ。

【4】 前回見たように、双曲線上の任意の点 S は、(sec θ, tan θ) と (cosh t, sinh t) の2通りに媒介変数表示される。もちろん、どちらも同一の点の別表記にすぎないので:
  sec θ = cosh t  【カ】
  tan θ = sinh t  【キ】

特に【キ】では、t と θ が1対1に対応する(双射)。そこで【キ】の両辺の arctan を取って、θ を t の関数 θ = f(t) と考えよう:
  θ = arctan (sinh t)  【ク】
定義域は −∞ < t < ∞。arctan の値域は、通常の主値 −π/2 < θ < π/2 とする。

【キ】の両辺の arsinh を取れば、【ク】の逆関数 t = f−1(θ) を考えることもできる:
  t = arsinh (tan θ)  【ケ】

【ケ】の意味は、前回の作図から一目瞭然。θ が上記の範囲(±90°)を動くとき、黄緑の点 Q が縦線上を自在に動く。それに対応して、青い点 S が双曲線上を動き、水色の面積の半分に当たる t が定まる。イメージの分かりやすさがいい。

以下で見るように、【ケ】の導関数は sec θ、従って【ケ】は sec θ の原始関数。三角関数の範囲で処理しようとすると不透明なトリックが必要になるが、グーデルマン関数を使うと、分かりやすい一本道になる。

【5】 双曲線関数の解析的な定義は cosh t = (et + e−t)/2, sinh t = (et − e−t)/2、従って:
  (sinh t)′ = cosh t, (cosh t)′ = sinh t  【サ】
  cosh2 t − sinh2 t = 1, 1 + sinh2 t = cosh2 t  【シ】

【サ】は (et)′ = et, (e−t)′ = −e−t による。【シ】が成り立つのは当然といえば当然だが(cosh t と sinh t は 双曲線 x2 − y2 = 1 の x と y なので)、定義からも容易に確認可能

y = sinh t = (et − e−t)/2 の逆関数を t = arsinh y と書く。A = et と置くと:
  y = (A − A−1)/2
  両辺を 2A 倍して 2yA = A2 − 1, A2 − 2yA − 1 = 0
  2次方程式を解いて A = y ± y2 + 1
  つまり et = y + y2 + 1 (下記の理由で複号のマイナスを無視)
  両辺の対数を取って arsinh y = t = log (y + y2 + 1)  【ス】

−∞ < t < ∞ は実数なので A = et は正の実数であり、複号のマイナスは題意に適さない(y も実数)。

【6】 【ケ】の導関数 dt / dθ が sec θ であることを示したい。直接計算してもいいのだが、逆関数の微分なので、まずは素直に逆関数の微分法で。次の基本公式*1については既知とする:
  d(arctan x) / dx = 1 / (1 + x2)  【タ】

とりあえず【ク】を t について微分。【タ】【サ】を使うと:
  dθ / dt = [1 / (1 + sinh2 t)] cosh t (合成関数の微分法)
  = (1 / cosh2 t) cosh t (【シ】より)
  つまり dθ / dt = 1 / cosh t   ……… (*)
  従って dt / dθ = cosh t = sec θ (なぜなら【カ】)

レガシーな方法と比べると、あっけないほど簡単に、自然に結論が得られた。何よりいいこととして、「双曲線は (cosh t, sinh t) または (sec θ, tan θ) で表される」という認識さえあれば(前回の作図から、鮮明にイメージできる)、【カ】【キ】【ク】【ケ】は間違いようのない一本道。双曲線関数についての簡単な知識【5】や、逆関数の微分法【6】は、応用範囲の広い事柄であり、この方法なら「1/cos x の積分のやり方のトリック」を個別的に暗記する必要がない。

<参考> (*)によると、グーデルマン関数【ク】は、1 / cosh t = sech t の原始関数。グーデルマン関数のフォーマルな定義
  gd u = 0t sech u du
は、この事実に対応している(t = 0 のとき【ク】の値が 0 であることから、積分定数は 0)。

【7】 【ケ】などを使って、sec θ の積分の伝統的表現(log | sec θ + tan θ | + C など)を導くことは、難しくない。【ケ】は arsinh の対数表現【ス】において y = tan θ としたものなので:
  arsinh (tan θ) = log (tan θ + tan2 θ + 1) = log (tan θ + sec θ)

最後の等号については、前回の(☆☆☆)参照。…意味不明にも思える教科書の式 log | sec θ + tan θ | + C の正体は、【ケ】の arsinh (tan θ) を対数表現したものだった!

ただし、グーデルマン関数では θ の範囲が (−π/2, π/2) に制限されている。その範囲外のことは別途考える必要がある。

*1 検索すれば、説明ページがたくさん見つかる。
Derivative of arctan(x) - MIT OpenCourseWare (PDF)
Derivative of the Inverse of a Function - MIT OpenCourseWare (PDF)

2021-03-01 グーデルマン関数の可視化 円と双曲線 虹の懸け橋

【1】 cos θ の逆数 1 / cos θ を sec θ と書く(sec は secant の略)。cos θ が正なら sec θ も正。

PNG画像

図のように、単位円 x2 + y2 = 1 上の点 P を考える。三角関数の基礎によれば、P の座標は (cos θ, sin θ)、AQ の長さは tan θ。これは点 Q のy座標に当たる:
  y = tan θ  (☆)

恒等式 cos2 θ + sin2 θ = 1 の両辺を cos2 θ で割ると:
  1 + tan2 θ = sec2 θ  (☆☆)

さて、θ が −π/2 < θ < π/2 の範囲を動くとき、Q は黄緑の縦線上を動くが、「Q とy座標が等しい」点 S が、双曲線 x2 − y2 = 1 の右の枝(青い曲線)の上にも存在する。この点 S のx座標は 1 以上の正の数で、
  x2 = 1 + y2
  x = 1 + y2
を満たす。(☆)(☆☆)を使うと:
  x = 1 + tan2 θ = sec2 θ = sec θ  (☆☆☆)

考えている範囲では sec θ は正の数なので、上記のように、単純に根号を外していい。

従って、双曲線 x2 − y2 = 1 の右の枝の上にある任意の点 S は、何らかの角度 θ を使って (sec θ, tan θ) と表される。ちなみに OQ の長さは sec θ なので、点 B は OQ を半径とする円上にある(図解)。

θ が −π/2 < θ < π/2 の範囲を動くとき、Q は黄緑の縦線上の全範囲を動き、それに対応して S は青い双曲線上の全範囲を動く。

一方、双曲角 t(それは水色の領域 SOA の面積の2倍に等しい)を使うと、点 S は (cosh t, sinh t) とも表される(水色の領域がx軸より下になる場合、その面積を負、つまり t を負とする)。

双曲角 t は、普通の意味での ∠SOA とは意味が異なる。普通の意味での ∠SOA 自体、θ = ∠QOA とは異なる角度。では θ と t は、どういう関係にあるのか? 点 S の座標を (sec θ, tan θ) と書くことも (cosh t, sinh t) と書くこともでき、それらは同じ値なのだから:
  sec θ = cosh t   【ア】
  tan θ = sinh t   【イ】

【イ】を使うと、θ の値を定めれば、自動的に t の値が定まる。逆に、t を定めれば θ が定まる。例えば【イ】の両辺の arctan を考えると:
  θ = arctan (sinh t)   【ウ】

これによって、θ を t の関数 f(t) = arctan (sinh t) と考えることができる(グーデルマン関数と呼ばれる)。

【ア】を使ってもほぼ同様だが、【ア】は S のx座標なので常に正になり、それだけでは θ ないし t の正負を区別できない。【イ】は S のy座標なので θ と t が1対1対応(一方が正なら他方も正、一方が負なら他方も負)。だから、とりあえず【イ】を使うのが便利。

【2】 上記の関数 f は、sec の原始関数と関係している。自分の目で確かめてみよう。上の図では t = 1、つまり水色の面積は 1/2。このとき AQ の長さ(言い換えれば S のy座標)が 1 より少し大きいのは作図から明白。すなわち θ は π/4 = 45° より少し大きい。もし仮に θ = π/3 = 60° なら AQ = tan θ = √3 = 1.732… だが、図の AQ は、1.732… よりはだいぶ小さい。45° より少し大きいが、60° よりだいぶ小さいのだから、ざっと 50° くらいだろう。

さて、sec の原始関数のグラフ(下図の赤い曲線)を見ると、確かに縦座標が 1 のとき、横座標は π/4 より少し大きい。

PNG画像

グラフの縦線の間隔は π/12。一方、上記の目分量による 50° = 50π/180 = 5π/18 は π/12 を単位とすると 10/3 = 3.333… に当たる(5π/18 ÷ π/12 より)。縦座標が 1 の場所の赤曲線をじっくり見ると、確かに横座標 3.3 × π/12 くらい。50° というのは目分量なので誤差はあるが、「t = 1 のとき θ ≈ 3.3π/12」という対応は、sec の原始関数のグラフとよく一致している。…グーデルマン関数について、少し具体的な感触が得られた。

正確な計算によると【ウ】から
  θ = arctan (sinh 1) = arctan 1.1752011… = 0.8657694… = 3.3069958… × π/12
この θ に対して
  sec θ + tan θ = 2.7182818… = e
なので、教科書的な sec の原始関数を使っても:
  log | sec θ + tan θ | = log e = 1

つまり上記の角度 θ と双曲角 t = 1 は対応していて、その対応関係は sec の原始関数で表される。このこと(グーデルマン関数の逆関数が sec の原始関数)については、次回、証明する。

【3】 これを考えるメリットは大きい:

  1. 多くの人を悩ませる ∫√x2 − 1 dx のようなもの(因習的には、それは三角関数を使った置換積分でゴチャゴチャ計算される)。双曲線関数を使うと、それがはるかに簡単になる。y = x2 − 1 は双曲線の第1象限の部分なので、双曲線関数を使って自然に表されることは当然だが、【1】で見たように、sec と tan でも表現可能。…2種類の置換の関係を見通すことができれば、視野が広がり理解も深まるだろう。
  2. 上記のような問題を普通の方法で考えると、sec x や sec3 x の積分が現れる。sec x の積分は、極端に面倒なわけではないが、多少トリッキーな操作を要求され、結果も意味不明な形に思える。【2】で見たように、この意味不明に見える関数の正体は【ウ】の逆関数であり、その観点から考えると、全てがすっきりする。この関数は、ネット上の地図(メルカトル図法)の縦座標として、日常的にも役立っている。
  3. 6種類の双曲線関数について、単純明快な視覚化が可能になる。

いくらメリットが大きいからといって、積分を使って教科書的にグーデルマンを定義すると、その定義自体が天下り的に思える。上記の作図はこの欠点を補い、グーデルマン関数の意味を可視化したもの。アイデアの一部は The Weird World of the Gudermannian に基づく。内容を簡単化し、半径 OQ の円を考えることで OB = sec θ をイメージしやすくした。

チラ裏より

チラ裏」は、きちんとまとまった記事ではなく、断片的なメモです。


新着記事

妖精の森 ♌ ペル方程式の夏(2020-12-27) e
x2 − 79y2 = 5 を満たす整数 (xy) は存在しません。その証明は意外と難しく、しかも隠された深い意味を持っています。この種の問題を扱います。ハイライトは、2020年夏に発見されたばかりの「改良版コンラッドの不等式」。中学生向け。 〔最終更新=v3: 2021年1月10日〕
ばびっと数え歌 でかい数編 (2019-09-01) l
31桁の 1,000,000,000,000,000,000,000,000,000,000(=100穣)までの数え歌。日本語・英語・SI接頭辞・2進数付き。
(za)b = zab の成立条件(2019-06-09) m
(za)b = zab は一般には不成立。ではどういう条件で、この等式が成り立つか。(za)bzab は、どういう関係にあるのか。「巻き戻しの数」(unwinding number)は、この種のモヤモヤをすっきりさせるための便利なコンセプト。 〔最終更新=v5: 2021年1月24日〕
ツンデレ素数 (チラ裏3題)(2019-05-13)
べ、別に、あんたを分解したくて楕円曲線法を実装したんじゃないんだからね。任意精度演算ライブラリを作り過ぎちゃって、もったいないから、ついでに分解しただけなんだから…。勘違いしないでよね!
−1 の 3/2 乗? オイラーの公式(その2)(2019-03-03)
(−1)3/2 って ((−1)3)1/2 = (−1)1/2 = i なのか、((−1)1/2)3 = i3 = −i なのか、それとも…? exp zez が同じという根拠は? 〔最終更新=v7: 2021年1月24日〕
60%他の生物【人体の細胞】100%星くず(2019-02-24)
ヒトの体は約25兆の細胞から成るが、体には65兆の細胞が…。本人以外の40兆は何なんでしょ? 〔v8: 2019年4月18日〕
曇りなきオイラーの公式 微分を使わない直接証明(2019-02-17)
exp ix = cos xi sin x のこんな証明。目からうろこが落ちまくる! 〔v11: 2020年12月23日〕

数学・プログラミング

まあるい緑の単位円 (三角関数覚え歌)(2017-12-24)
まあるい緑の単位円/半径 斜辺の三角形/「高さ」の「さ」の字はサインの「サ」/サインは 対辺 高さ
アルファとベータが角引いた (加法定理・図解の歌)(2017-12-24)
「ごんべさんの赤ちゃん」のメロディーで。「アルファさんとベータさんが麦畑」でもOK。 〔最終更新: 2018年1月28日〕
cos 36° 魔法のにおい(2018-01-14)
五角形を使った解法も優雅だが、代数的に… 〔最終更新: 2019年9月30日〕
cos π/7 正七角形の七不思議(2018-01-28)
日頃めったに見掛けない正七角形。その作図不可能性は、有名な「角の3等分問題」に帰着する。コンパス・定規・「角度3等分」器があれば、360° を7等分できる! 〔最終更新: 2018年12月30日〕
覚えやすさを重視した3次方程式の解法(2018-02-11)
分数なくして、すっきり。語呂合わせ付き。 〔v8: 2019年3月17日〕
3次方程式の奥(2018-03-04)
3次方程式は奥が深い。「判別式の図形的解釈」は1990年代の新発見だという。 〔v14: 2021年1月12日〕
3次方程式の判別式(2018-03-18)
いろいろな判別式。Qiaochu Yuan による恐ろしくエレガントな解法。 〔最終更新: 2018年12月30日〕
3次方程式と双曲線関数 ☆ 複素関数いじっちゃお(2019-02-17)
定義から始めてのんびり進むので、双曲線関数の予備知識は不要。3次方程式も別記事で初歩から解説。三角・指数関数なら知ってるという探検気分のあなたへ。複素関数プチ体験。 〔v7: 2021年2月19日〕
cos i = ?
曇りなきオイラーの公式 微分を使わない直接証明(2019-02-17)
exp ix = cos xi sin x のこんな証明。目からうろこが落ちまくる! 〔v11: 2020年12月23日〕
−1 の 3/2 乗? オイラーの公式(その2)(2019-03-03)
(−1)3/2 って ((−1)3)1/2 = (−1)1/2 = i なのか、((−1)1/2)3 = i3 = −i なのか、それとも…? exp zez が同じという根拠は? 〔最終更新=v7: 2021年1月24日〕
(za)b = zab の成立条件(2019-06-09)
(za)b = zab は一般には不成立。ではどういう条件で、この等式が成り立つか。(za)bzab は、どういう関係にあるのか。「巻き戻しの数」(unwinding number)は、この種のモヤモヤをすっきりさせるための便利なコンセプト。 〔最終更新=v5: 2021年1月24日〕
フェルマーのクリスマス定理で遊ばせて!(2018-12-23)
1640年のクリスマスの日、フェルマーはメルセンヌに宛てた手紙の中で、こう言った。「4の倍数より1大きい全ての素数は、ただ一通りの方法で、2個の平方数の和となります」 〔v5: 2020年12月27日〕
「西暦・平成パズル」を解くアルゴリズム(2016-03-27)
整数28と四則演算で2016を作るには、最小でも9個の28が必要。
2016 = (28+28+28)×[28−(28+28+28+28)/28]
一見全数検索は大変そうだが、50行程度の平易なスクリプトで高速に解決される。ES6 の Map の長所、splice より速い要素挿入法も紹介。 〔最終更新: 2016年4月10日〕
[JS] 100行のプチ任意精度ライブラリ(2016-05-08)
JavaScript 用に最小構成的な「任意精度整数演算」ライブラリを作ってみた。 〔最終更新: 2019年6月23日〕
[JS] メルセンヌ数の分類と分解(2016-06-05)
数千万桁のメルセンヌ素数が脚光を浴びるが、その裏では、たった数百桁のメルセンヌ合成数が分解できない。 〔v6: 2019年5月5日〕
楕円曲線で因数分解(2016-08-14)
楕円曲線を使って、巨大整数に含まれる数十桁の因数を検出できる。計算は、曲線上の勝手な点を選んで整数倍するだけ。ステージ1、モンゴメリー形式、標準版ステージ2、素数ペアリングについて整理した。 〔最終更新: 2019年6月4日〕
楕円曲線の位数: 点の擬位数に基づく計算法(2016-10-02)
元の位数を考えると群の位数計算が高速化されるが、それには高速な素因数分解が必要。「擬位数」はどの教科書にも載ってないような概念だが、ハンガリー人数学者 Babai László によって研究された。 〔最終更新: 2016年10月23日〕
アルカンの異性体の数の公式・第1回 小さなパズルと不思議な解(2015-09-20)
異性体の数は難しいが、炭素数12くらいまでなら素朴な計算ができる。中学数学くらいの予備知識で気軽に取り組めて、めちゃくちゃ奥が深い。(全9回予定だが第6回の途中で止まっている。そのうち気が向いたら完結させたい)
「マイナス×マイナス=プラス」は証明できるか?(2014-08-03)
数学的に正しい質問は、「なぜマイナス×マイナス=プラスか?」ではなく「いつマイナス×マイナス=プラスか?」 〔最終更新: 2019年9月29日〕
平方剰余の相互法則(2003-03-26)
「バニラ素数とチョコレート素数」という例えを用いた「お菓子な」説明。
楕円曲線暗号(2003-11-28)
最初歩から具体例で。書き手も手探りというライブ感あふれる記事6本。手探りだからエレガントではないが、JavaScriptでは世界初の実装? 実装はダサいが、内容(ロジック)は正しい。
触って分かる公開鍵暗号RSA(2004-02-04)
理論的説明でなく、実地に体験。JavaScriptで実現したので結構注目され、大学の授業などの参考資料としても使われたらしい。ダサい実装だが、ちゃんと動作する。
デスノートをさがして: 論理パズル(2006-04-10)
真神・偽神・乱神。間違いだらけの乱神探し。

天文・暦

13日は金曜になりやすく31日は水曜になりにくい(2017-09-03)
曜日は「日月火…」の繰り返しだから各曜日は均等のようだが、「毎月1日の曜日」「13日の曜日」のように「特定の日にちが何曜になるか」を考えると、曜日分布に偏りが… 〔v6: 2019年4月21日〕
「春夏秋冬」は「夏秋冬春」より長い(2017-11-26)
「春分→夏→秋→冬→春分」と「夏至→秋→冬→春→夏至」は、どっちも春・夏・秋・冬1回ずつなのに、前者の方が長い。素朴な図解(公転最速理論?)、簡易計算、そして精密な解析解。春分間隔から春分年へ… 〔最終更新: 2018年12月30日〕
公式不要の明快な曜日計算(2016-10-23)
公式や表を使わず、何も覚えていない状態で、手軽に任意の年月日の曜日を暗算。
ぼくの名前は冥王星(2013-09-30)
いいもん、いいもん! これからは小惑星になって、ジュノーちゃんやベスタちゃんと遊ぶから! …と思っていたら、「おまえは小惑星でもないんだよ」と言われてしまった。そんなー。ぼくのアイデンティティーは粉々さ。 〔v6: 2019年3月24日〕
さよなら第9惑星・冥王星 カイパーベルト終着駅(2019-03-24)
海王星~海王星~。目蒲めかま線はお乗り換えです。
第9惑星・追悼演説(2019-03-24)
我々は一つの惑星を失った。しかし、これは「終わり」を意味するのか? 否、始まりなのだ!
ケプラー方程式(微積・三角公式を使わないアプローチ)(2018-01-14)
微積分を使わず、算数的にケプラー方程式を導く。倍角・半角などの公式を使わずに、離角の関係を導く。特別な予備知識は不要。 〔最終更新: 2018年2月4日〕
ケプラー方程式・2 エロい感じの言葉(2018-01-28)
「ケプラー方程式(微積・三角公式を使わないアプローチ)」の別解・発展。 〔最終更新: 2018年2月4日〕

シリア語・Unicode・詩

シリア語: カラバシ注解(2013-12-01)
カラバシ『読み方のレッスン』はシリア語文語・西方言の教科書。ウェブ上で公開されている。その魅力を紹介し、第1巻全21課に注釈を付けた。 〔最終更新: 2016年5月8日〕
ばびっと数え歌 シリア語編(2014-02-09)
「シリア語の数詞の1~10」を覚えるための数え歌。「ごんべさんの赤ちゃん」のメロディーでも歌えます。 〔最終更新: 2017年12月24日〕
ペシタ福音書における「女性聖霊・男性聖霊」の混在について(2014-12-14)
キリスト教の「聖霊」はイエス自身の言語では女性だったが、後に男性イメージに変化した。この変化は興味深いが、そこに注目し過ぎると中間期の状況が正しく理解できない。3種類のシリア語聖書とギリシャ語聖書を比較し「叙述トリック」を検証。 〔最終更新: 2018年11月4日〕
少年と雲 (シリア語の詩)(2017-12-24)
雲さん、どこから来たんだい?/背中に何をしょってるの?/そんなに顔を曇らせて/空から何を見ているの?
黙示録の奇妙な誤訳: 楽しいシリア語の世界(2018-04-15)
「南の子午線を飛ぶハゲタカ」が、なぜか「尾が血まみれのハゲタカ」に…。誤訳の裏にドラマあり。 〔最終更新: 2018年5月6日〕
ターナ文字入門: 表記と発音(2013-01-16)
以前公開していた記事を全面改訂。ターナ文字は、インドの南、南北1000キロにわたって散らばる島々で使われる文字。 〔最終更新: 2014年5月4日〕
HTML5 の bdi 要素と Unicode 6.3 の新しい双方向アルゴリズム(2012-12-04)
ブログのコメント欄で起きる身近な例を出発点に、双方向性が絡む問題と解決法を探る。HTML の dir 属性は落とし穴が多い。HTML5 の <bdi> は役立つ。近い将来、「ユーザー入力欄などの語句は、このタグで隔離」が常識になるかも。 〔最終更新: 2014年4月27日〕
空は青くて真白くて(2014-11-23)
「わたしの心は躍り上がる」(ワーズワース)/「空は青くて白くて」(フィンランド民謡)

ジョーク

未来の水 フリーズドライ ☆ 粉末乾燥水(2012-04-01)
宇宙旅行のお供に/非常時の備えに… 場所を取らない超軽量・携帯用のインスタントお水です。
イヤ~な「金縛り」を強制解除 ☆ 全自動かなほど機(2019-04-01)
睡眠中の金縛り。嫌なものですね…。そこでご紹介するのが、この「かなほど機」。金縛りになったとき、ワサビの匂いで身体を自動リセットする未来の製品です。
さよなら第9惑星・冥王星 カイパーベルト終着駅(2019-03-24)
海王星~海王星~。目蒲めかま線はお乗り換えです。
漢詩と唐代キリスト教 「日本の影響」説も(2019-04-01)
客舍かくしゃ青青せいせい 柳色りゅうしょく新たなり」仏教徒でもあった唐の大詩人・王維(おうい)。彼がキリスト教とも関わっていたことは、ほとんど知られていない。(エイプリルフールのジョーク記事)
円周率は12個の2 スパコンで判明/ほか 3題(2016-04-01)
三原則ロボットおちょくられて仕返し?/円周率は12個の2 スパコンで判明/人間を模倣する学習AI 学習し過ぎ?
ISOとJISによる「ハッカー」の正式な定義(2005-02-19)
JIS規格では「ハッカー」という言葉が定義されてる。
ヒマワリをふてくされさせる実験(2005-02-20)
お花はとってもデリケート。
「確信犯」たちの「開発動機」(2005-09-23)
ストラビンスキー「ファゴット奏者を苦しめてやろうとしてやった。苦しそうな音なら何でも良かった」
「水からの伝言」の世界(2006-08-21)
水さん、ちょっと漏れ過ぎです。
脳内ディベート大会(2009-07-31)
応援団を応援することは正しいか。タンポポの綿毛を吹いて飛ばしていいか。

漫画・アニメ

大島弓子の漫画 (チラ裏3題)(2019-04-28)
バナブレは「漫画で何ができるのか?」という世界の枠組みそのものを変えた。綿国(わたくに)は、漫画・アニメ史上「猫耳の発明」という意味も持つ。もともとは「自分は半分人間だと思っている子猫」の主観的世界を表す絶妙な表現。
ラピュタ滅びの呪文は波動砲かフェーザー砲か?(2006-01-28)
ムスカは、ジブリ作品では珍しい悪役と評されるが、ラピュタ文字の解読は、現実世界ならノーベル賞もの。
勇者よ、侵略者から東京を守れ(2006-01-22)
「ブジュンブラにキメラアニマが現れたわ!」 お気に入りのネタだが、アニオタ以外の一般人には意味不明かも。
チラ裏
アニメ関係の小ネタも多い。イタリアのアニメ事情もあるよ。

字幕

MKV埋め込み字幕用フォントのMIME問題 (2019-10-20)
字幕用フォントが、ロードされない事例が起きている。問題の背景・対策・対応状況。
SSA入門 中級編(2004-08-27)
二つの入門編(音声タイミング・基本スタイリング)に続くフレーム・タイミング関連の内容。古い記事で使用ツールは時代遅れだが、考え方は依然参考になるかも。
[SSA/ASS] 高品質のフェイドイン・フェイドアウト(2005-12-21)
単純な fad() は濁りやすい。各種の代替手段を紹介。
ASS: 縁ワイプと縦カラオケ(2006–2009)
字幕と音声のずらし方/縁ワイプ/字幕のリップシンク/縦カラオケ/他。古い記事だが参考までに。

哲学・ファンタジー

60%他の生物【人体の細胞】100%星くず(2019-02-24)
ヒトの体は約25兆の細胞から成るが、体には65兆の細胞が…。本人以外の40兆は何なんでしょ? 〔v8: 2019年4月18日〕
至るところ青山 (チラ裏3題)(2019-04-14)
3丁目が見えない理由(先行きの不安)は、1丁目にいるからで、2丁目まで行けば自然と選択肢は狭まる。
不死でないから星は輝く (チラ裏3題)(2019-04-14)
「核融合には燃料が必要。燃料を使い果たせば反応は止まる」という当たり前のことを言い換えると「いつかは終わるから今輝いている」。
猫のしっぽを思い切り引っ張ることは十戒のどれに違反するか?(2014-11-23)
南泉は言った。「この猫の命が惜しければ、禅を一言で語れ。さもないと猫を斬り殺す」 〔最終更新: 2019年4月24日〕
神から見た「主の祈り」(2004-10-04)
「天にましますわれらの父よ」 神「はい?」 — へリング牧師は、ジョークのような設定で深い問題を提示した。 〔最終更新: 2013年10月2日〕
「無断コピー以外」を禁止するライセンス(2004-10-04)
人間の心理的困難があまりに大きいようなので、 それに対抗するために、次のような新しいライセンス形態を思いつくほどだ。いわく…
妖精物語 3題(2005-07-02)
王様の赤いばらと白いばら。
「反辞書」の著者フレッド・レスラー(2009-02-03)
Urban Dictionary というサイトをご存じでしょうか。 ウィキペディアみたいな、でもそれよりずっと砕けた新語辞典…

お知らせ

メールアドレス<PNG画像>