遊び tokidoki 仕事

数学と音楽と教育と遊び

| おしごと - きょういく - がくせい - ゼミ - イベント | すうがく - おんがく - 数理音楽 - DTM - かがく - scratch
| Art - photo - おきにー - Tips - ものもう - あれこれ | About - Top

バランスをとる数(1)―入り口は1+2=3


いまや当大学における絶滅危惧レッドリストの筆頭に挙げられる数学同人Sigmaでのある一幕.
https://aue-sigma.jimdo.com/
www.aichi-edu.ac.jp

久しぶりに教育実習から返ってきたメンバーが「実習中にこんなこと見つけた」と鼻息荒く話していた話題なんだが,それは例えば
\begin{align}
1+2&=3\\
4+5+6&=7+8\\
9+10+11+12&=13+14+15\\
16+17+18+19+20&=21+22+23+24
\end{align}
といった連続整数の和についての等式だった.つまり一般に

\displaystyle
\displaystyle\sum_{k=0}^a (a^2+k)=\sum_{k=a+1}^{(a+1)^2-1}(a^2+k)
ということだ.あるいは
\begin{equation}
\displaystyle2\sum_{k=0}^a (a^2+k)=\sum_{k=0}^{(a+1)^2-1}(a^2+k)
\end{equation}
と言い換えても良い.そしてこれは任意の自然数aについて確かに正しいことが計算すれば分かる.


さて,本当の遊びはここからだ.では一般に
\begin{equation}
\displaystyle2\sum_{k=a}^x k=\sum_{k=a}^b k
\end{equation}
となるような自然数a\le x\le bの組はどんなものだろうか?計算すれば
\begin{equation}
2(x-a+1)(x+a)=(b-a+1)(b+a)
\end{equation}
すなわち
\begin{equation}
2(2x+1)^2=(2a-1)^2+(2b+1)^2
\end{equation}
となる.X=2x+1,A=2a-1,B=2b+1と置けば奇数A,X,Bの組で
\begin{equation}
A^2+B^2=2X^2
\end{equation}
を満たすものを探せということになる.あるいはs=\frac{A}{X},t=\frac{B}{X}と置けば,
\begin{equation}
s^2+t^2=2
\end{equation}
を満たす有理点(s,t)で,分母分子共に奇数になっているものを探せ,ということになる.


そこでとりあえず有理点(s,t)を探すわけだが,例えば(s,t)=(\pm1,\pm1)などがすぐ見つかる.
図形的に見やすい(-1,-1)を足がかりに,有理数m>0を傾きとする直線 t=m(s+1)-1s^2+t^2=2 の交点を考えよう.
すでに(-1,-1)が交点,すなわちこの二式の連立の解なので,2次方程式
\begin{equation}
s^2+(m(s+1)-1)^2=2
\end{equation}
(s+1)で括られ,残りの解が
\begin{equation}
\displaystyle s=\frac{1+2m-m^2}{m^2+1},t=\frac{m^2+2m-1}{m^2+1}
\end{equation}
と得られる.あとはm=\frac{q}{p}と既約分数で表示したとき,s,t共に\frac{\text{奇数}}{\text{奇数}}となっているかを確認すればよい.計算すると
\begin{equation}
s=\displaystyle\frac{p^2+2pq-q^2}{p^2+q^2},t=\frac{q^2+2pq-p^2}{p^2+q^2}
\end{equation}
となり,p,q いずれか一方のみが奇数なら分母分子共に奇数,p,qいずれも奇数ならば
\begin{equation}
p^2+q^2\equiv p^2+2pq-q^2\equiv q^2+2pq-p^2\equiv 2\pmod{4}
\end{equation}
となるため,約分すればやはり分母分子共に奇数と分かる.


こうして任意の有理数mに対する組(s,t)が求めるもの全てとなり,更にp,qいずれか一方のみ奇数ならば
\begin{equation}
A=p^2+2pq-q^2,B=q^2+2pq-p^2,X=p^2+q^2,
\end{equation}
p,q共に奇数ならば
\begin{equation}
A=\displaystyle \frac{p^2-q^2}{2}+pq,B=\frac{q^2-p^2}{2}+pq,X=\frac{p^2+q^2}{2}
\end{equation}
と得られる.


もっとも,元の問題に立ち戻れば自然数和を考えていたので,
\begin{equation}
A=p^2+2pq-q^2=(p+q)^2-2q^2\ge 0
\end{equation}
\begin{equation}
B=q^2+2pq-p^2=(p+q)^2-2p^2\ge 0
\end{equation}
という条件,すなわち
\begin{equation}
(\sqrt{2}-1)p\le q\le (\sqrt{2}+1)p
\end{equation}
p,qに課されている.例えばこの条件下にあるq=p+1ならば,
\begin{equation}
A=2p^2-1,B=2(p+1)^2-1,X=2p^2+2p+1,
\end{equation}
すなわち
\begin{equation}
a=p^2,b=(p+1)^2-1,x=p^2+p
\end{equation}
となって,冒頭の学生が発見した等式に戻る.

また(p,q)=(3,5)とすると(A,B,X)=(7,23,17)すなわち(a,b,x)=(4,11,8)となり,新たな等式
\begin{equation}
4+5+6+7+8=9+10+11
\end{equation}
が見つかる.あるいは(p,q)=(4,7)とすると(A,B,X)=(23,89,65),すなわち(a,b,x)=(12,44,32)となって,別の新たな等式
\begin{equation}
12+13+\cdots+32=33+34+\cdots+44
\end{equation}
が得られる.q=2pも条件下にあり,このときは
\begin{equation}
A=p^2,B=7p^2,X=5p^2
\end{equation}
と綺麗な形になり,奇数pに対して新たな系列
\begin{equation}
a=\frac{p^2+1}{2},b=\frac{7p^2-1}{2},x=\frac{5p^2-1}{2}
\end{equation}
が得られる.



使っているのは高校数学程度なんだが,結構面白い.
おやおや,なんだかミイラ取りがミイラになりかかっているぞ...

今年もガッサガッサしてみた

この時期は毎年,担当講義である統計とコンピュータにてガッサガッサする実験をする.
tokidoki.hatenablog.jp

今年で3年目となるこの実験の元ネタは,啓林館のページにあった授業実践だった.
www.shinko-keirin.co.jp

ただし,気をつけないと隣のクラスからうるさいと言われかねない.
(実際,今年は隣で講義をされていた先生が何事かとそっと教室を見に来た.)
f:id:okiraku894:20150601112752j:plain

3年目ともなると結構な数の実験になる.一人50回×3年間実験参加者数272=13600回だ.
そうして得られた結果が以下.
f:id:okiraku894:20170626121052p:plain
ここまでくると コイン>ピン>キャップ のこの順が入れ替わるなんてことは
とんでもなく低い確率になろう.
実際,10×標準偏差分の幅を取ってやっとコインとピンが入れ替わる可能性が出てくる.
そしてその確率は7.7\times10^{-23}程度だ.
そんな評価ができるのも中心極限定理のおかげなんだが,
やっぱり学生にちっとも伝わらないのがこの現象.
もちろん今年も正方形ダーツボードの実験はしてもらうし,
なかなか軽い動作でシミュレーションできるこの回のEXCEL資料は
我ながら気に入っているのでまた掲載するよ↓
f:id:okiraku894:20150607135853g:plain
f:id:okiraku894:20150607145228p:plain

今年はこれに加えて「上か下か」のランダムウォークもEXCEL上で見せようと思った.
f:id:okiraku894:20170626124705p:plain
5本程度,1000ステップ分を描くと↑のようになるが,この程度ではふーんで終わる.
もっと本数を描いてみよう.例えば100本↓
f:id:okiraku894:20170626124706p:plain
こうなってくると形,いわゆる放物線がうっすら見えてくるのではなかろうか.
そしてこれらの軌跡は中心0に近いほど密度が高くなっている.
いるのだけどこれが平面では見せられない.
ということで今度は OpenProcessing にて立体的にこれを見せられるようにしてみた↓

上部は密度を立体的に表示し,その足元の影がランダムウォークの軌跡を表している.
因みに,パスの数が5000以下のときのみランダムウォークを表示するようにした.
それ以上の本数になるとヒストグラムと軌跡両方の表示はどうやら重すぎるので.


ところで全く関係ないが,とある事情によりイベントが急遽行われた.
この場を借りて,あざっす!
f:id:okiraku894:20170623122352j:plain

私の身の回りのべき乗則(2)―名字編

今クールNHKが放映している名字番組
「人名探究バラエティー 日本人のおなまえっ!」ってのがある.
www4.nhk.or.jp
帰宅がてら見るのにちょうど良い番組なのだけど,
レア名字編ではハンコタワー1つ5000種の印鑑で人口の80%がカバーできるとか,
10万種の印鑑の在庫を持つハンコ屋さんがあってそれでも人口の98%だとか,
って話を聞いてると,「おお!冪乗則!」ってなって,早速検索.
直ぐに名字ランキングサイトが見つかる.
以下のグラフはこのページから参照.
myoji-yurai.net
幸いにもその名字におよそ何人いるか,のデータも合わせて書かれてあったので,
冪乗則が調べられる.さてやってみたのが↓
f:id:okiraku894:20170609153627p:plain
EXCELが吐き出す近似曲線がy\fallingdotseq4\cdot 10^7x^{-1.14}で決定係数R^2が0.98を超える.
やはりなかなか良い冪乗則になっている様子.
\log-\log plot が↓
f:id:okiraku894:20170609153628p:plain

それでも何だか二次関数っぽくもあったので,\log_{10}-\log_{10} plotでの二次近似を行うと
X=\log_{10} x,Y=\log_{10} yとして
 Y= -0.189X^2 - 0.0647X + 6.1599,決定係数 R^2 = 0.9997 となった.
因みに,人数が一人,すなわち Y=0 となるのが X\fallingdotseq 5.54 で,
すなわちランキングとしては x=10^X\fallingdotseq 347014 となって,
名字が約29万種あると言われているところからすると,ちょっと多めの評価になる.
もっとも果たして「一人」の状態までカウントして良いのやら,というところでもある.
f:id:okiraku894:20170609173723p:plain


日本の名字は分家制で名字の派生が起こったり同じ読みのまま漢字を変えたり,と
名字の生成/消滅の過程は集団遺伝の力学系のようでもあって,
名字の力学モデルなんてのも考えても面白いのかもしれない.

私の身の回りのべき乗則

それは気のせいなだけなのかもしれないが,
何だかパソコンの動作がモッサリしてきたなぁ,
と思ったときにClean upした後デフラグをかけて,
モッサリがなくすという作業を定期的にしている.

で,そのとき分析結果として各ファイルの断片化数が報告されるわけだが,
それを分断数が多い順に並べ替えると
いわゆる「べき乗則」に従っているように以前から感じていた.
f:id:okiraku894:20170530130948p:plain

で,とある機会に実際にグラフにしてみたところ,なかなか見事な結果になった↓
f:id:okiraku894:20170526180652p:plain
ただし,デフラグのログが取り出せなかったので手打ちでデータを書き写したのだが,
分断数6以下のものは大量すぎて作業から外した.
EXCELによって与えられる近似曲線はy=5557.4x^{-0.95}で,
世にいうZipf則が成り立っているように見える.

それでも図に騙されている気がしたので両側対数軸 \log y/\log x で描いたのが↓
f:id:okiraku894:20170526180653p:plain
なかなか良い感じだ.

実はこの記事を書こうと思ったのは他でもない,
予算が削減されつつある大学の電子ジャーナルを契約更新するか否かの判断材料として
ジャーナルごとのアクセス数のデータをもらったことに始まる.
早速アクセス数の多いジャーナル順にソートして描いてみたら見事なべき乗則↓
f:id:okiraku894:20170526180655p:plain
f:id:okiraku894:20170526180654p:plain
おやおや.

この結果に気を良くして,では世間で言われているweb siteのアクセス数についても
べき乗則が見えるのかと,例えばこのブログの1日あたりアクセス数を
多い順にソートして並べてEXCELで近似曲線まで描かせると
どうもべき乗則ではなく指数分布のほうが決定係数が大きくなる.
(まぁ,決定係数が8種類もあってEXCELがどれを採用しているのか,
 とかいった厄介な話はこの際考えない.)
f:id:okiraku894:20170529121850p:plain

あれれ,何が起こっているのか?と思い両対数グラフを描いたら
\log x=2.5 あたりから急に落ち込んでいた↓
f:id:okiraku894:20170529121851p:plain

ならば\log x=2.5に該当する20アクセス未満の日を切り落とせば良いのか,と描いてみたのが↓
y=253.8x^{-0.408}と今度はべき乗則のほうが決定係数が大きくなった.
f:id:okiraku894:20170529121852p:plain
まぁ,不都合な部分を切り落としたのだから,そりゃべき乗則のほうが良くなるでしょ,
ともいえるのだけど.
そういえば最初のデフラグデータも分断数が小さいものは調べなかったから,
順位の大きい先のほうは切り落としたことになっているわけで,
世にいう「べき乗則」は順位の比較的小さい部分で起こっている話,と理解するものなのかなぁ.

Art Award IN THE CUBE 2017 @ 岐阜県美術館

f:id:okiraku894:20170503144428j:plain
ふらっと行ってきました,Art Award IN THE CUBE @ 岐阜県美術館.
art-award-gifu.jp
審査員には(私が勝手に)同士だと思っている高橋源一郎さんや
日々,折々の言葉で(勝手に)お世話になっている鷲田清一さん,
初年次演習ではかつて(勝手に)アレンジして遊んだことで(勝手に)お世話になった
逆シミュレーション音楽の三輪眞弘さんがいることも気になって行ってみた.
tokidoki.hatenablog.jp

f:id:okiraku894:20170503130918j:plain
応募者らはこの与えられたCUBEの中に作品を作るという寸法で,何でも今年が最初の試みとのこと.
古い世代はCUBEというとあの映画を思い出すのだけどね.

CUBEキューブ Blu-ray

CUBEキューブ Blu-ray

では,中を見ていこう.
f:id:okiraku894:20170503130755j:plain
三枝 愛「庭のほつれ」/O JUN賞.
美術館庭にぽつんと建てられたCUBE.何があるかと覗くところんと木が転がってる.
東日本大震災からシイタケ農家に必要な原木が不足したこと,
それによって生活の風景が一変したこと.
情報の激流によって薄皮でしか物を考えなくなった社会で,
敢えて考えるために立ち止まる人たち.
現代アートとはそういった人たちが居る場所でもあるのだよね.

f:id:okiraku894:20170503122832j:plain
f:id:okiraku894:20170503141404j:plain
安野 太郎「THE MAUSOLEUM ー大霊廟ー」/高橋源一郎賞.
館内に入ると何やら不協和音がブーブカブーブカ聞こえてくる.
その発生源がこれ.大小さまざまなリコーダーがポンプから送られる空気で鳴っている.
人類滅亡後,永遠に弔いを続けるロボットという想定らしい.
しかし,こんなデカいリコーダーもあるんだなぁ,と別の角度で感心.

f:id:okiraku894:20170503123047j:plain
柴山 豊尚「ニョッキ(如木)2017」/十一代大樋長左衛門賞.
板を重ねてそれを削ってできた,どこか別の惑星の景色のような空間.
この中で寝そべって想像すると楽しそうだ.

f:id:okiraku894:20170503123157j:plain
f:id:okiraku894:20170503124127j:plain
森 貞人「Mimesis Insect Cube」/中原 浩大賞.
ワサワサワサッとガラクタで作られたたくさんの昆虫.「ガラクタ虫」と呼んでるそうだ.
使い古された道具たちがこうして虫として生まれ変わった.
f:id:okiraku894:20170503141116j:plain:w400,left
f:id:okiraku894:20170503141127j:plain:w400,right

で,ひょいと目をやると何やら人の下半身が上の方で垂れ下がってる!
f:id:okiraku894:20170503123228j:plain
何々と中を覗くと人が倒れてる,っていうか手がカタツムリ!
f:id:okiraku894:20170503123312j:plain
f:id:okiraku894:20170503123356j:plain
耳のないマウス「移動する主体 (カタツムリ)」/三輪 眞弘賞.
小さな子が見たら夢に出てきそうな作品.これちょっとずつ動いてるんだ.
f:id:okiraku894:20170503123436j:plain
壁に貼りついてるヒトの指先だって,カタツムリの目のように動いてた.
こう,こっちに寄って来るとうわ~ん,てなるよ↓
f:id:okiraku894:20170503123506j:plain
あぁ,あれだ,ヒトがどんどんカタツムリになっていく
伊藤 潤二のマンガ「うずまき」を思い出したよ.

うずまき (2) (スピリッツ怪奇コミックス)

うずまき (2) (スピリッツ怪奇コミックス)

ちょっと落ち着こか.
f:id:okiraku894:20170503123702j:plain
宮原 嵩広「Missing matter」.
CUBE全体に敷き詰められた「アスファルト」,ただそれだけだ.
確かにアスファルトのニオイが充満して余りそこに立っていたくない感じだった.
やがて枯渇する化石燃料,大木や自然に感じる神性と同じ思いを
アスファルトに抱く時代がくるかもしれない,ってなことらしい.

f:id:okiraku894:20170503123752j:plain
f:id:okiraku894:20170503123808j:plain
三木 陽子「Conduit(導管)」.
そのとき全く分からなかったのだけどあとから調べたら陶芸作品とのこと.
建物の裏側なんかにありそうな導管なんだけど,これが蔓のように壁にまつわりつき,
したたかな生き物の様でもあり,
あぁ,うちの庭の草取りしなきゃなぁ,なんて思ったり.

f:id:okiraku894:20170503124554j:plain
f:id:okiraku894:20170503124616j:plain
松本 和子「透明の対話」.
どこか懐かしい感じはフレスコ画という表現方法にあるのだろうか.
朝の透明さがなんだか涼しかったりした.

f:id:okiraku894:20170503124808j:plain
ミルク倉庫+ココナッツ「cranky wordy things」/大賞.
寝言を言っている人の動画やら,突然音を立てるガンガンやら,
勝手にくるくる回り続ける机の上の小物.
その場は直ぐに立ち去ったが,家に帰ってから解説を読んだ.

物質という有限性を凌駕する何かしらに憑かれ、一時の『身体』がたち現れる

ああ,なるほど,だから何だかオカルティックだったのか.
そういえばケン・ウィルバーの「空像としての世界」を読んだとき,
僕らって要するに精神の受信機なんだ,と思ったのだっけ.
そう,僕らの身体はいつか土に還るまでの借り物なんだってことさ.

空像としての世界―ホログラフィをパラダイムとして (1984年)

空像としての世界―ホログラフィをパラダイムとして (1984年)

f:id:okiraku894:20170503124941j:plain
谷本 真理「この部屋とダンス」/田中 泯賞.
一見,現代アートにしばしばありがちなダダイズム的な何か,かと部屋を覗く.
田中 泯賞を取ったとのこと,そしてタイトルには「ダンス」.
あれ,確かに部屋の落書きは何か躍った後のような,全てのものが躍動しているような.
できればタイトル見ないでそれを感じ取りたかった.ちょいと悔しいような.

f:id:okiraku894:20170503125837j:plain
中村 潤「縫いの造形」.
CUBE全体が「縫う」という行為を表現していた.
ああ,できればその「縫う」行為の中に埋もれてみたかったかもしれない.

f:id:okiraku894:20170503125427j:plain
f:id:okiraku894:20170503125711j:plain
堀川 すなお「モノについて」.
実は入ったところに中が見えない箱が二つ置いてあり,
そこに手を入れて入っているものの様子を感じて,それから作品を見ると面白い.
子どもたちに同じことをしてもらって,感じたことを言葉にして,
その情報を元に再びデザインするといった行為が展開されていた.
とはいえ,かなりデザインが繊細で美しかったため,
だから子どもの言葉とどれくらい対話が行われて描かれたのだろうか,
と,かえって訝ってしまう自分がいた.

f:id:okiraku894:20170503125150j:plain
水無瀬 翔「DEMO DEPO イン・ザ・キューブ支店」/鷲田 清一賞.
後でじっくり解説を読んで,なるほど鷲田さんが選びそうな話だ,と思った.
ほら,例えば本来家庭が受け持っていた躾を始めとする様々な役割を
今どきの人たちは学校にアウトソーシングしているでしょ.
そのくせきちんと行われていないと学校を非難するでしょ?

この作品,本来どうしても社会に主張すべきことがあって
それを社会に知ってもらうために行うはずの「デモという行為」をアウトソーシングする,
っていうブラックジョークなわけだ.
しかもデモするのは人間ではなくロボット.
分業化が徹底的に進んで一生の営みがどんどんアウトソーシング化して,
果たして僕らのどこに僕らの人生があるんだろうね.
気付けばYahoo!知恵袋に課題の質問をする学生達ってのは,
学習自体もアウトソーシングしているわけで,
そんな風に上澄み液だけ啜って終わる人生って,楽しいのかね.

↓アハアハアハ,っていう機械的な笑いもアウトソーシングですか.
f:id:okiraku894:20170503125215j:plain

さてさて,しかしこの日をわざわざ選んで岐阜県美まで行った理由,それがこの作品↓
f:id:okiraku894:20170503123852j:plain
平野 真美「蘇生するユニコーン」.
この作品,ただのぬいぐるみではない.
架空の動物であるユニコーンを解剖学的見地を元に,
骨格から内臓,筋肉,表皮まで全てを実際に作ってあるのだ.
かつ,作家さん当人の話によれば口から肛門まできちんと一続きの仕組みにしてあり,
心臓も心房心室を作って外部からの「輸血」によって循環するように作り込み,
また口に伸びる管からは空気が送り込まれ,
それは確かに「肺」にとどけられているとのこと.
いや,実際胸のあたり,静かに上下しているのだ.
岐阜大獣医学科に足を運び獣医から学び,リアリティーを追及したとのこと.
本当は解剖にも立ち会う予定だったそうだが,
鳥インフルエンザ騒ぎで獣医学科がてんやわんやで話が流れてしまったのだそうだ.

しかし今回,その平野さんが会場に現れ「胃の交換手術」を行うのである.
f:id:okiraku894:20170503134905j:plain
↓赤いのが,交換前の「胃」.
f:id:okiraku894:20170503135859j:plain

f:id:okiraku894:20170503142544j:plain
ラテックス製の胃は随分と弾力性があった.
f:id:okiraku894:20170503142559j:plain
f:id:okiraku894:20170503143113j:plain
f:id:okiraku894:20170503143248j:plain
f:id:okiraku894:20170503143306j:plain
↓うわぁ,随分と大胆に胃を押し込んでいくよ.
f:id:okiraku894:20170503143331j:plain
f:id:okiraku894:20170503143335j:plain
f:id:okiraku894:20170503143340j:plain
↓右の黄色いのは取り出された方の胃.素材が違うのだそうだ.
f:id:okiraku894:20170503143924j:plain
手術の最後までは見届けなかったのだけど.
f:id:okiraku894:20170503143422j:plain
こうして思うのは,モノを一生懸命作る人というのは,やはり美しいということ.
ってお前,何を撮りに行ってるんだよ,というツッコミは無しの方向で.
いずれもILCE-6000+(SEL1670ZF4.0,SEL30F3.5Macro), Lightroomにて現像

おっと,そうそう,もちろん帰りの電車では高橋源一郎さんの本を読んでたよ.

丘の上のバカ ぼくらの民主主義なんだぜ2 (朝日新書)

丘の上のバカ ぼくらの民主主義なんだぜ2 (朝日新書)

実質日本会議である改憲勢力が2/3を超え維新を含めた大政翼賛会化してしまった国会と,
そしてその日暮らしのごく身の回りのことしか関心を持たない
羊よりも大人しい多数の同調主義者たち(これを「羊たちの沈没」という).
ひとたびこの国でテロが起これば,
あっという間に九条の精神は不安と怒りに呑みこまれてしまうだろうし,
それにかこつけて「個人」を「人間という動物」へと格下げする,
小学生の読書感想文にも劣る,しょーもない日本会議草案に流れてしまうのだろうな.
何てことあんまり書いていると,現代版治安維持法(テロ等準備罪)に引っかかったりして.
あたらしい憲法草案のはなし

あたらしい憲法草案のはなし

  • 作者: 自民党の憲法改正草案を爆発的にひろめる有志連合
  • 出版社/メーカー: 太郎次郎社エディタス
  • 発売日: 2016/07/02
  • メディア: 単行本(ソフトカバー)
  • この商品を含むブログ (5件) を見る