Thursday, 30 May 2024

小説 君のクイズ

2022年 小川哲 朝日新聞出版

階差数列

生放送のクイズ番組の決勝戦に出たクイズプレーヤーの三島玲央が,まだ一言も問題が読まれていないのに正解を答えて優勝した対戦相手に対して不正を疑い,真相を解明しようとする話です.

僕は早押しクイズは数列と似ていると思っている。
1、2、4……と聞こえた時点で、僕は数列のルールがわかったと思ってボタンを押す。「この数列において10番目に来る数は何か」と問われているが、ボタンを押した時点で答えがわかっているわけではない。1、2、4の次は8だろう。この数列は、前の数を2倍にしていくも のなのだ ($a_n=2^{n-1}$)と考えて、10番目に何が来るのかを急いで計算する。 
実際に、それが正解ということもあるが、間違っている可能性もある。この数列はまだ確定していないからだ。
1、2、4の次に7が来る場合も考えられる。1から2は1増えている。2から4 は2増えている。4の次は3増えて、7になるかもしれない($a_n=\frac{n^2}{2}-\frac{n}{2}+1$)。この数列は階差数列かもしれない。

$1, 2, 4, 8, ......$ なら,ここまでは初項1,公比2の等比数列になると推測できます.そうすると一般項(第$n$項)は $2^{n-1}$ ですから,第10項は $2^{10-1}=2^9=512$ と計算できます.

$1, 2, 4, 7, ......$ なら一般項が$\frac{n^2}{2}-\frac{n}{2}+1$になることはすぐには分かりませんね.この数列を $\{a_n\}$ とし,階差数列を$\{d_n\}$とすると,$d_n=1, 2, 3, ......$ すなわち $d_k=k$ なので,$\{a_n\}$ の一般項は次のようになります.\begin{eqnarray}a_n&=&a_1+\displaystyle \sum_{k=1}^{n-1} d_k=1+\displaystyle \sum_{k=1}^{n-1} k&=&1+\frac{1}{2}n(n-1)=\frac{1}{2}\left(n^2-n+2\right)\tag{1}\end{eqnarray}ですから,その第10項は $\frac{1}{2}\left(10^2-10+2\right)=46$ になります.

この数列は「怠け仕出し屋の数列 (lazy caterer's sequence)」と呼ばれていて,$n=1, 2, 3, 4, ......$ のときに $1, 2, 4, 7, ......$ なら,円を$(n-1)$本の直線で切り分けたときにできる領域の最大数を表しています.(Wikipedia等のよくある説明では,円を$n$本の直線で切り分けるとしているので,$n=0, 1, 2, 3, ......$ のときに $1, 2, 4, 7, ......$ となり,一般項が $\frac{1}{2}\left(n^2+n+2\right)$となっています) (オンライン整数列大辞典A000124)

因みに,$1, 2, 4, 8, ......$ だからといって一般項が $2^{n-1}$ にならない例は他にもあります.

まず1つ目はこの数列です.$$\{b_n\} \quad 1, 2, 4, 8, 15, ...... $$第4項まで $2^{n-1}$ になっていますね.これの階差が $a_n=1, 2, 4, 7, ......$ すなわち「怠け仕出し屋の数列」=式$(1)$ になっていますから,$\{b_n\}$ の一般項は次のようになります.\begin{eqnarray}b_n&=&b_1+\displaystyle \sum_{k=1}^{n-1}a_k=1+\displaystyle \sum_{k=1}^{n-1} \frac{1}{2}\left(k^2-k+2\right)\\&=&\frac{1}{6}\left( n^3-3n^2+8n\right)\tag{2}\end{eqnarray}この数列は「ケーキ数 (cake number)」と呼ばれていて,$n=1, 2, 3, 4, 5, ......$ のときに $1, 2, 4, 8, 15, ......$ なら,円柱等の凸立体を$(n-1)$枚の平面で切り分けたときにできる領域の最大数を表しています.(Wikipedia等のよくある説明では,立体を$n$枚の平面で切り分けるとしているので,$n=0, 1, 2, 3, 4, ......$ のときに $1, 2, 4, 8, 15, ......$ となり,一般項が $\frac{1}{6}\left(n^3+5n+6\right)$となっています) (オンライン整数列大辞典A000125)

2つ目は次の数列です.$$\{c_n\} \quad 1, 2, 4, 8, 16, 31, ......$$第5項まで $2^{n-1}$ になっていますね.これの階差が $b_n=1, 2, 4, 8, 15, ......$ すなわち「ケーキ数」=式$(2)$ になっていますから,$\{c_n\}$ の一般項は次のようになります.\begin{eqnarray}c_n&=&c_1+\displaystyle \sum_{k=1}^{n-1} b_k=1+\displaystyle \sum_{k=1}^{n-1} \frac{1}{6}\left( k^3-3k^2+8k \right)\\&=&\frac{1}{24}\left( n^4-6n^3+23n^2-18n+24 \right) \tag{3}\end{eqnarray}これはモーザー数列 (Moser's sequence)」といって,円周上に $n$ 個の点があり, すべての2点を結ぶ弦で円を切り分けたときにできる領域の最大数を表しています.(オンライン整数列大辞典A000127

従って,$1, 2, 4, ......$ だからといって等比数列にならない例は,計算さえすれば次々とできます.

    $a_n= 1, 2, 4, 7, ......$            第3項まで$2^{n-1}$ 「怠け仕出し屋の数列」=式$(1)$
    $b_n= 1, 2, 4, 8, 15, ......$         第4項まで$2^{n-1}$ 「ケーキ数」=式$(2)$
    $c_n= 1, 2, 4, 8, 16, 31, ......$      第5項まで$2^{n-1}$ 「モーザー数列」=式$(3)$
    $A_n= 1, 2, 4, 8, 16, 32, 63, ......$  第6項まで$2^{n-1}$ 

3つ目,$\{A_n\}$の一般項も計算してみましょう!

これは階差数列が $1, 2, 4, 8, 16, 31, ...... $ すなわち「モーザー数列」=式(3) になっていますから,\begin{eqnarray}A_n&=&A_1+\displaystyle \sum_{k=1}^{n-1} c_k=1+\displaystyle \sum_{k=1}^{n-1} \frac{1}{24}\left( k^4-6k^3+23k^2-18k+24 \right)\\&=&\frac{1}{120}\left( n^5-10n^4+55n^3-110n^2+184n \right)  \end{eqnarray} この計算は大変でした.(オンライン整数列大辞典A006261

ついでにもう一つ,余裕のある人は計算してみてください.$1, 2, 4, 8, 16, 32, 64, 127, ...... $ の一般項は?
正解はこちらです.

まあ,実際のところ,早押しクイズに使えそうなのは公比2の等比数列だと思いますけどね.

余談ですが,「ストゥリクス・ウラレンスという学名を持ち,『森の番人』のイメージか......」という,主人公とそのライバルが答えられなかった問題で,正解がすぐに分かりました.問題の続きは「......ら,千葉駅前の交番のモチーフにもなっている生き物は何でしょう?」というものでした.(正解はこちら
 

Tuesday, 14 May 2024

小説 同志少女よ、敵を撃て

 2021年 逢坂冬馬 早川書房

ミル(mil=角度の単位)

第二次世界大戦中の1942年,旧ソ連の小さな村で,当時18歳の少女セラフィマの母親を含む村人全員をドイツ軍が惨殺し,さらに赤軍(ソ連地上軍)が村全体を燃やしたことから,生き残ったセラフィマが狙撃兵になって復讐しようとする話です.

狙撃兵訓練学校での座学の内容です.
基礎の基礎として、全員が「ミル」という単位を覚える。ミルとは射撃および砲撃の照準に用いられる角度の単位であり、周回360度を6000ミルとして定義する。すなわち正面に対して右に90度 は1500ミルであり、上に45度は750ミルである。

なぜこんな面倒な単位を用いるかといえば、「1000メートル先にある、幅1メートルのもの」がおおよそ1ミルであるからだ。この単位を用いることは照準に有利となる。

故にスコープを覗いて、「幅50センチと推定される物体が1ミルの幅に収まっている」という状態があるならば、彼我の距離は500メートルと計算できる。
もともと1ラジアン(rad)の1000分の1を1ミリラジアン(mrad)といい,1周360度=2$\pi$ラジアンは2$\pi$×1000≒6283(ミリラジアン)になります.

ラジアンの定義より,半径1の扇形の弧の長さが1になるときの中心角が1ラジアン(≒57.3度)ですから,半径1000mの扇形の弧の長さが1000mになるときの中心角も1ラジアンです.ということは半径1000mの扇形の弧の長さが1mになるときの中心角は$\frac{1}{1000}$ラジアン=1ミリラジアン(≒0.0573度)になります.

なので「1000メートル先にある、幅1メートルのもの」は1ミリラジアンになるわけですが,計算しやすいように旧ソ連では360度=6000ミル(日本では6400ミル)とし,「1000メートル先にある、幅1メートルのもの」を「おおよそ1ミル」としています.

射撃用の銃のスコープにはミル単位の目盛りがついていますが,ここでは角度というよりは標的の見かけの大きさをミルで測ります.

[例1] 上の例「幅50cm(0.5m)と推定される物体が1ミルの幅に収まっている」(下図の赤い長方形)なら,それは2倍大きく見えているので,距離は1000÷2で500mということになります.      
    
[例2] 高さ1.7mの標的が1ミルに見えたら,距離はその1000倍で1700mですが,4ミルの高さに見えたら(下図の赤い長方形),それは4倍大きく見えているので,距離は1700÷4で425mになります.
つまり,実際の標的の大きさを $a$ m,スコープ内での見かけの大きさを $m$ ミルとすると,次の式で距離 $d$ (m)を計算できます.$$d=\frac{1000a}{m}$$従って,標的の大きさが一定の場合,ミルと距離は反比例することになります.

[参考]
HB-PLAZA

Thursday, 9 May 2024

ドラマ 霊験お初 ~震える岩〜

2024年5月放送 テレビ朝日 宮部みゆき原作

算額

霊験(不思議な力)を持つ主人公のお初が連続殺人事件の解決に挑むというミステリー・ホラー時代劇です.算術の好きな与力見習いの右京之介が,お初に算額の問題を見せる場面(TVerで37分から38分の間)です.

お初「算額にはどのようなことが書かれているのですか?」
右京之介「お見せしましょうか.右の絵では外側の大きな円の長さ,左の絵では内側の2つの小円の径の長さを求めるのです」
お初「ひょうでもって測ったらどうかしら」
右京之介「それを算術で求めるところが学問なのです」
①の絵
①の絵では,外側の大きな円に甲,乙,丙,丁4種類の円が内接しています.

左右対称なので中心が直径上にある甲,丁の半径が分かればすぐに答えが出てしまいますから,甲,乙,丙または乙,丙,丁の半径が与えられていると推測できます.実際,3つの互いに外接する円があれば,その3円を内接させる大円がひとつ決まるので,乙,丙,丁の半径が与えられていると考えられます.

互いに外接する3円を内接させる大円はデカルトの円定理を使って求めることができます.3円の半径をa, b, cとし,大円の半径をrとすると,この定理で次の方程式が成り立ちます(証明難).$$2\left(\frac{1}{a^2}+\frac{1}{b^2}+\frac{1}{c^2}+\frac{1}{r^2}\right)=\left(\frac{1}{a}+\frac{1}{b}+\frac{1}{c}-\frac{1}{r}\right)^2$$例えば$a=1$,$b=2$,$c=3$のときは,$$2\left(\frac{1}{1^2}+\frac{1}{2^2}+\frac{1}{3^2}+\frac{1}{r^2}\right)=\left(\frac{1}{1}+\frac{1}{2}+\frac{1}{3}-\frac{1}{r}\right)^2$$これを解くと$r=6$となります.

ところで以前,これの類題を,各円の中心間を斜辺とする直角三角形に三平方の定理を適用して求めましたが,非常に大変な計算になりました.

②の絵
②の絵では,大円1つと小円2つが正三角形に内接し,大円と小円が外接しています.

問題文を読むことはできませんでしたが,正三角形の1辺と大円の中心の位置または半径が与えられていると思われます.計算をなるべく簡単にするため,正三角形の1辺を12として,底辺の中点が原点になるように置き,大円の中心を(0,  4)とします.小円の中心を(b,  c)とおき,半径cを求めます.

大円の半径DEは,点(0, 4)と直線 $y=-\sqrt{3}x+6\sqrt{3}$ との距離なので,$$DE=\frac{|0+4-6\sqrt{3}|}{\sqrt{1+3}}=3\sqrt{3}-2$$

直角三角形DGIに三平方の定理を当てはめると,$$(4-c)^2+b^2=(3\sqrt{3}-2+c)^2\tag{1}$$小円の半径 $c=IH$は,点(b, c)と直線 $y=-\sqrt{3}x+6\sqrt{3}$ の距離なので,$$c=\frac{|\sqrt{3}b+c-6\sqrt{3}|}{\sqrt{1+3}}\tag{2}$$連立方程式(1)(2)を解けば,cは求められます.複雑な計算を経て最後に次の方程式を解きます.$$3c^2-(18\sqrt{3}+4)c+21+12\sqrt{3}=0$$すると,小円の半径は次の値になります.$$c=\frac{2+9\sqrt{3}-2\sqrt{46}}{3}≒1.3413$$
 
(②の別解)  大円の中心より半径を与える方が計算が楽ではないかと思い,大円の半径を3,その中心を(0,  a),小円の中心を(b,  c)とおいて計算してみました.

大円の半径DEは,点(0, a)と直線 $y=-\sqrt{3}x+6\sqrt{3}$ との距離なので,$$DE=\frac{|0+a-6\sqrt{3}|}{\sqrt{1+3}}=3$$$$a=6\sqrt{3}-6$$直角三角形DGIに三平方の定理を適用すると,$$(6\sqrt{3}-6-c)^2+b^2=(3+c)^2\tag{1}$$小円の半径 $c=IH$は,点(b, c)と直線 $y=-\sqrt{3}x+6\sqrt{3}$ の距離なので,$$c=\frac{|\sqrt{3}b+c-6\sqrt{3}|}{\sqrt{1+3}}\tag{2}$$連立方程式(1)(2)を解けば,cは求められます.複雑な計算を経て最後に次の方程式を解きます.$$c^2+(2-8\sqrt{3})c+57-24\sqrt{3}=0$$すると,小円の半径は次の値になります.$$c=-1+4\sqrt{3}-2\sqrt{4\sqrt{3}-2}≒1.4883$$というわけで,初期条件を変えてもやはり大変な計算になりました.ああ,しんど(笑).

Friday, 3 May 2024

映画 Mean Girls (2024)

2024年 米国

関数の極限

直訳すると「意地悪な女の子たち」というタイトルの学園コメディで,2004年に公開された同名映画のリメイク版です.前作で主役だったリンジー・ローハンが学校対抗の数学コンテストの司会役でカメオ出演(特別出演)しています.このコンテストの最後,あと1問で勝敗が決まるという場面で,2004年の作品とまったく同じ次の問題が使われていました.$$\lim_{x\to 0} \frac{\ln(1-x)-\sin x}{1-\cos^2 x}$$

(ただし,$\ln x=\log_{e} x$)

相手校の生徒が「$-1$」と答えて不正解だった後,主人公のケイディが「極限は存在しない」と答えて勝利します.(前作と全く同じストーリーです)

なぜその答になるのか見てみましょう.まず上式を少し変形して,\begin{eqnarray}&=&\lim_{x\to 0} \frac{\ln(1-x)-\sin x}{\sin^2 x}\\&=&\lim_{x\to 0} \frac{\ln(1-x)-\sin x}{x^2} \cdot \frac{x^2}{\sin^2 x}\\ &=&\lim_{x\to 0} \frac{\ln(1-x)-\sin x}{x^2}\end{eqnarray}このままでは極限を求め難いので,分母も分子も微分可能で$\frac{0}{0}$(不定形)であることからロピタルの定理を使います.この分母子を微分すると,$$\lim_{x\to 0} \frac{\frac{-1}{1-x}-\cos x}{2x}$$定理より,これに極限があれば元の式の極限と一致するのですが,$$\lim_{x\to +0} \frac{\frac{-1}{1-x}-\cos x}{2x}=\frac{-2}{+0}=-\infty$$$$\lim_{x\to -0} \frac{\frac{-1}{1-x}-\cos x}{2x}=\frac{-2}{-0}=\infty$$となり,右極限と左極限が一致せず,極限が存在しないので,元の式も極限が存在しません.

因みに,グラフ描画アプリGeoGebraで確認するとこうなります.