[b]このワークシートは[url=https://www.geogebra.org/m/twxxx3yq]Math by Code[/url]の一部です。[br][/b][br]集合というと、外延的定義と内包的定義の2種類の説明だったり、[br]ベン図をかいてその交わりや合併を見える化してみたり、[br]何か、止まったイメージがあるね。[br][br]でも、[color=#0000ff][b]関数を使うことで、集合と集合がダイナミックにつながったり、生成したりするイメージ[/b][/color]が[br]わいてくる。[br][br][b][size=150]<部分集合を作る関数>[br][/size][/b][br]たとえば、部分集合というのがあるね。[br]4つの要素の集合があるとしたら、[br]その部分集合は、各要素を入れる/入れないを4回選べるので、[br]2の4乗通りの部分ができることは基本的なことだね。[br]極端でつまらない部分集合、自明な部分集合が2つあったね。[br][br]4つとも入らない∅[br]4つとも入る全集合[br][br]つまらないといっても、極端なことは数学では重要になることがある。[br]この2つが入ることで、きれいに2の要素数乗の部分集合があるといえるしね。[br]この部分集合たちを1つの集合と見ると、その個数は2のベキ乗になる。[br]そんなところから、この集合(部分集合の集合)を[b]ベキ集合[/b]と呼ぶこともあるね。[br][br][color=#9900ff][b][u][size=150]質問:集合から部分集合を作り出す関数はどうやって作ったらよいですか?[br][/size][/u][/b][/color][br]juliaで作るとしたらどうしよう。[br]たとえば、4要素の集合xsの部分集合を作りたかったら、[br]ターゲット(target)は、はじめは∅だけのリストにしよう。[br]xsの先頭要素itemを取り出す。[br][br][color=#0000ff][b]・そのターゲットのコピー(参照ではなくて、別メモリ上)res_sameを作る[br]・そのターゲットの各要素に、itemをつけたして上書きする。[br]・更新されたtargetに対して、もとのコピー(res_same)を連結して、倍要素数のリストにする。[br][/b][/color][br]xsの次の要素itemを取り出す。[br] 上記3項目をやる。[br][br]これをxsの要素数length(xs)だけ繰り返そう。[br][br][color=#0000ff][IN]julia[br]#============================================================[br]function makeSubList(xs)[br] target=[[]][br] for id in 1:length(xs)[br] res_same = deepcopy(target) #copyや[=]では参照のみコピー[br] for item in target[br] push!(item, xs[id])[br] end[br] target=append!(target,res_same)[br] println(target)[br] end[br] return target[br]end[br]start=[1,2,3,4][br]makeSubList(start)[br]#============================================================[br][/color][out][br]Vector{Any}[[1], []][br]Vector{Any}[[1, 2], [2], [1], []][br]Vector{Any}[[1, 2, 3], [2, 3], [1, 3], [3], [1, 2], [2], [1], []][br]Vector{Any}[[1, 2, 3, 4], [2, 3, 4], [1, 3, 4], [3, 4], [1, 2, 4], [2, 4], [1, 4], [4], [1, 2, 3], [2, 3], [1, 3], [3], [1, 2], [2], [1], []][br][br][18]:[br]16-element Vector{Vector{Any}}:[br] [1, 2, 3, 4][br] [2, 3, 4][br] [1, 3, 4][br] [3, 4][br] [1, 2, 4][br] [2, 4][br] [1, 4][br] [4][br] [1, 2, 3][br] [2, 3][br] [1, 3][br] [3][br] [1, 2][br] [2][br] [1][br] [][br][br]ただし、集合といいながら、順番のあるリストというデータ構造を利用したコードなので、[br]本物の集合では、残念ながら動かない。[br][br]start=Set([1,2,3,4])[br]makeSubList(start)[br]MethodError: no method matching getindex(::Set{Int64}, ::Int64)[br][br]Stacktrace:[br] [1] makeSubList(xs::Set{Int64})[br] @ Main ./In[18]:6[br] [2] top-level scope[br] @ In[21]:2[br][br][color=#0000ff]start=[1,2,3,4][br]res=makeSubList(start)[br]Set(map(x->Set(x),res))[br][/color]集合として扱いたいならば、[br]こうして、結果を集合の集合に直しとけばいいね。[br][br][b][br][/b]
集合というデータ構造は[color=#0000ff][b]順番がない[/b][/color]ので、順次的な処理には向かない。[br]しかし、いいところがある。[br]それは、[color=#0000ff][b]要素がユニーク[/b][/color]だということ。[br]そして、処理によってリストが変化していくときに、それを集合に直すことで、[br]使った[color=#0000ff][b]要素を引き算[/b][/color]して、残りのまだ使える要素を取り出したりできることだ。[br][br][color=#9900ff][b][u][size=150]質問:リストを集合にするにはどうしたらいい?[br][/size][/u][/b][/color][br]juliaではSet()という関数でできます。[br][color=#0000ff][size=150]#[IN]julia[br]#================================[br]lst=[1,1,2,2,3,3,4,5][br]set=[b]Set[/b](lst)[br]#================================[br][/size][/color][OUT][br]Set{Int64} with 5 elements:[br] 5[br] 4[br] 2[br] 3[br] 1[br][br][color=#9900ff][b][u][size=150]質問:使った要素を引き算して、残りの要素を取り出したいときはどうしますか?[br][/size][/u][/b][/color][br]たとえば、1から100までの整数の集合lstを作ります。[br]次に、1から100までの整数の目がでるサイコロを300回ふって出た目の集合usedを作ります。[br]2つの集合lst, usedの差setdiff()は、まだサイコロを300回ふっても出ていない目です。[br][br][size=150][color=#0000ff]#[IN]julia[br]#================================[/color][color=#0000ff][br]lst=[b]Set[/b]([x for x in 1:100])[br]used = [b]Set[/b](rand(1:100, 1,300))[br][b]setdiff[/b](lst,used)[br]#================================[br][/color][/size][OUT][br]Set{Int64} with 6 elements:[br] 32[br] 91[br] 84[br] 76[br] 49[br] 38[br][br]毎回、結果が変わりますね。[br][br][b][size=150]<集合演算の確認>[br][/size][/b][br]2つのAとBの集合演算の関数を確認しておこう。[br]juliaではリストに対しても集合関数が使える。[br]ということは、リスト操作関数のときにも、思い出して使えるといいね。[br][br][color=#0000ff]A∩B:積集合(重なり)の関数[br][/color]julia> [b]intersect[/b]([1, 2, 3], [3, 4, 5])[br]1-element Vector{Int64}:[br] 3[br]julia> (0, 0.0) ∩ (-0.0, 0)[br]1-element Vector{Real}:[br] 0[br][color=#0000ff]A∪B:和集合(合併)の関数[br][/color][b]julia> union[/b]([1, 2], [3])[br]3-element Vector{Int64}:[br] 1[br] 2[br] 3[br]julia> (0, 0.0) ∪ (-0.0, NaN)[br]3-element Vector{Real}:[br] 0[br] -0.0[br] NaN[br][br][color=#0000ff]AーB:差集合の関数[br][/color]julia> [b]setdiff[/b]([1,2,3], [3,4,5])[br]2-element Vector{Int64}:[br] 1[br] 2[br][color=#0000ff]A△B = (A−B)∪(B−A):対称差(重なりを引いた残り)の関数[br][/color]julia> [b]symdiff[/b]([1,2,3], [3,4,5], [4,5,6])[br]3-element Vector{Int64}:[br] 1[br] 2[br] 6[br][br][color=#0000ff]x∈ A:xが集合(リストでもよい)Aの要素である[/color]関数[br]julia> a = 1:3:20[br]1:3:19[br]julia> 4 in a[br]true[br]julia> x∈(4, a)[br]true[br][color=#0000ff]x∉A:xが集合(リストでもよい)Aの要素でない関数[br][/color]∉(x, A) [br]juliaはユニコードも演算子や関数に使えるので、数学記号で動く。[br][br][color=#0000ff]A ⊆ B:AがBの部分集合である関数[br][/color]julia> issubset([1, 2], [1, 2, 3])[br]true[br][br]julia> [1, 2, 3] ⊆ [1, 2][br]false[br][br][br]
[b][size=150]<覆面算を集合関数で解こう>[br][/size][/b]覆面算[br][br][color=#0000ff][b][size=200]read + write + talk = skill[br][/size][/b][/color][br]のすべての解を求めたいとしましょう。[br][color=#9900ff][b][u][size=150]質問:どんなところに目をつけてコードを作りますか。[br][/size][/u][/b][/color][br]文字に数字を入れる処理と入れた結果で等式が成り立つかチェックする処理にわけると[br]かきやすいです。[br][br]・覆面算では、同じ文字は同じ数字、別の文字は別の数字が入ります。[br]・覆面算でも、最上位r,w,t,sは0がこないけど、他の位は0がこれます。[br]・a,e,t,r,kは2回出てくるので、先に入れてしまいたい。[br]・一の位に目をつけると、d+e+k≡l(mod 10)だから、lとの関係でdも先にいれる。[br]・さらに考えると数や和の範囲はもっと絞られるけど、[color=#0000ff]覆面算を正面から解くわけではない[/color]ので、[br]あとはコードで調べるのもよいでしょう。[br][br]pythonやjuliaでは、順列や組み合わせのループ用のパッケージがあるので、コードがすっきりしますよ。[br]aetrkdliwsの10文字に0から9の10数字を割り当てるから順列をすべて割り当ててもよい。[br]それだと、ループを10重にするよりはループが1つになりきれいにかけるけど、時間がかかりすぎるでしょう。無駄な調査を減らすために、ループを2段階にしてみるのもありだ。[br][br]・第1ループ[br]使える数字usable=[0,1,2,3,4,5,6,7,8,9]の集合setU=Set(usable)から6個を取り出した[color=#0000ff][b]順列vの6要素[/b][/color]をタプルにいれる。[b](a,e,t,r,k,d) = v[br][/b]そして、tとrが0でないなら、lをd+e+kの一の位として、vの数字集合Set(v)とset([l])が被らない。[br]・第2ループ[br]残った数字、[b]diffs=setdiff(setU,Set([a,e,t,r,k,d,l]))から残り3数字i,w,sをsetからArrayに直してから、順列pで埋めていく。[/b][br]埋めたら等式チェックをして、成立したらカウントして表示することにしよう。[br][color=#0000ff]#[IN]julia[br]#=================================================[br][/color]# REPL画面でCombinatoricsをインポートする必要がある。[br]# Pkg> add Combinatorics[br]# 順列も、このパッケージに入っている。[br][color=#0000ff]using Combinatorics[br][/color][br][color=#0000ff][b]usable=[0,1,2,3,4,5,6,7,8,9][br]setU=Set(usable)[br]cnt=0[br][/b]function [b]check[/b](a,e,t,r,k,d,l,i,w,s)[br] global cnt[br][/color] A=r*1000+e*100+a*10+d[br] B=w*10000+r*1000+i*100+t*10+e[br] C=t*1000+a*100+l*10+k[br] D=s*10000+k*1000+i*100+l*11[br] if A+B+C==D[br] cnt +=1[br] println("$A + $B + $C = $D")[br] end[br]end[br] [br][b][color=#9900ff]for v in Combinatorics.permutations(usable,6)[br] (a,e,t,r,k,d) = v[br][/color][color=#0000ff] if t!=0 && r!=0[br] l=(d+e+k) % 10[br] if ∉(l, Set([a,e,t,r,k,d]))[br] diffs=setdiff(setU,Set([a,e,t,r,k,d,l]))[br] iws =[ x for x in diffs][br][/color][color=#9900ff] for p in Combinatorics.permutations(iws)[br] (i,w,s) = p[br][/color][color=#0000ff] if w!=0 && s!=0[br] check(a,e,t,r,k,d,l,i,w,s)[br] end[br] end[br] end[br] end[br]end[br][/color][/b][br]print("$cnt 通り")[br]#=================================================[br][OUT][br]4905 + 24689 + 8017 = 37611[br]9728 + 19467 + 6205 = 35400[br]1632 + 41976 + 7380 = 50988[br]2543 + 72065 + 6491 = 81099[br]5270 + 85132 + 3764 = 94166[br]5180 + 65921 + 2843 = 73944[br]5094 + 75310 + 1962 = 82366[br]5096 + 35710 + 1982 = 42788[br]7092 + 37510 + 1986 = 46588[br]7092 + 47310 + 1986 = 56388[br]10 通り[br][br]