😊出会いの泉😊

両極に次の2つの数字を配置しておけば、他の6組の和を=>での値（自然数)にする残りの数字がちょうど2度ずつ出現するような
6組の6個ずつの数字の組合せは山ほど構成可能となると思います。
1;1,2=>69
2;1,5=>68
3;1,8=>67
4;1,11=>66
5;1,14=>65
6;1,17=>64
7;1,20=>63
8;2,4=>68
9;2,7=>67
10;2,10=>66
11;2,13=>65
12;2,16=>64
13;2,19=>63
14;3,6=>67
15;3,9=>66
16;3,12=>65
17;3,15=>64
18;3,18=>63
19;4,5=>67
20;4,8=>66 (例の図のパターン)
21;4,11=>65
22;4,14=>64
23;4,17=>63
24;4,20=>62
25;5,7=>66
26;5,10=>65
27;5,13=>64
28;5,16=>63
29;5,19=>62
30;6,9=>65
31;6,12=>64
32;6,15=>63
33;6,18=>62
34;7,8=>65
35;7,11=>64
36;7,14=>63
37;7,17=>62
38;7,20=>61
39;8,10=>64
40;8,13=>63
41;8,16=>62
42;8,19=>61
43;9,12=>63
44;9,15=>62
45;9,18=>61
46;10,11=>63
47;10,14=>62
48;10,17=>61
49;10,20=>60
50;11,13=>62
51;11,16=>61
52;11,19=>60
53;12,15=>61
54;12,18=>60
55;13,14=>61
56;13,17=>60
57;13,20=>59
58;14,16=>60
59;14,19=>59
60;15,18=>59
61;16,17=>59
62;16,20=>58
63;17,19=>58
64;19,20=>57

1～９にこだわってみました。

(15768/3942)×(12345/9876)=13485/2697

(18534/9267)×(17469/5823)=34182/5697

(17469/5823)×(31689/4527)×(65934/1782)=748251/963=615384/792


なお昔ここに投稿されていた記事をメモしていたノートを見直していたら
たぶん2017年ごろ
DD++氏が
e≈(1+.2^(3^84))^(5^(9^(6*7)))
を投稿されていたと思います。
またこれを
(1+9^(-4^(6*7)))^(3^(2^85))
(1+2^(-76))^(4^38+.5)
等にもでき、円周率πも
π≈2^5^.4-.6-(.3^9/7)^.8^.1
(((2^7+8)/(90-1))^5.4+.6)*.3 (りらひい氏発見）
(8-(1+.6/(.2*.5+9))/(.3+7))*.4 (らすかる氏発見）
などの常連さんの驚くべき技が紹介されていました。

「小町算で無理数近似」の記事ですかね。

> "DD++"さんが書かれました:
> 「小町算で無理数近似」の記事ですかね。

この部分を読み返していたら
ネイピア数と小数点以下8368澗4289溝8906穣8425秭9438垓1759京0916兆4450億0188万7164桁まで一致する近似値
の精度がどうやって導けたのかの謎が
マクローリン展開が
1/e*(1+x)^(1/x)=1-1/2*x+11/24*x^2-7/16*x^3+2447/5760*x^4-959/2304*x^5+O(x^6)
これから
e-(1+x)^1/x≒1/2*e*x
これにx=.2^(3^84)=(1/5)^(3^84) なる微小な値を取ることで
両辺のlog[10]をとると右辺が
gp > log(exp(1)/2)/log(10)-3^84*log(5)/log(10)
%139 = -8368428989068425943817590916445001887164.5053429251
正にeとの誤差が
1/10^(8368428989068425943817590916445001887164)
つまり
8368澗4289溝8906穣8425秭9438垓1759京0916兆4450億0188万7164桁まで一致ということを示す。
という計算なのですね。

整理すると
e+o+r+z=0
e+n+o=1
o+t+w=2
2e+h+r+t=3
f+o+r+u=4
e+f+i+v=5
i+s+x=6
2e+n+s+v=7
e+g+h+i+t=8
e+i+2n=9
e+n+t=10
3e+l+n+v=11
2e+l+t+v+w=12
それぞれの文字が含まれる方程式の数を数えて個数の昇順にすると
1個: g,u,x,z
2個: f,h,l,s,w
3個: r
4個: i,o,v
5個: n,t
10個: e
g,u,x,zは一度しか登場しないので
e+o+r+z=0
f+o+r+u=4
i+s+x=6
e+g+h+i+t=8
の4つは後回しで残りは
e+n+o=1
o+t+w=2
2e+h+r+t=3
e+f+i+v=5
2e+n+s+v=7
e+i+2n=9
e+n+t=10
3e+l+n+v=11
2e+l+t+v+w=12
残った方程式で再度個数を数えると
1個: f,h,r,s
2個: i,l,o,w
4個: t,v
5個: n
8個: e
f,h,r,sは一度しか登場しないので
2e+h+r+t=3
e+f+i+v=5
2e+n+s+v=7
の3つは後回しで残りは
e+n+o=1
o+t+w=2
e+i+2n=9
e+n+t=10
3e+l+n+v=11
2e+l+t+v+w=12
残った方程式で再度個数を数えると
1個: i
2個: l,o,v,w
3個: t
4個: n
5個: e
iは一度しか登場しないので
e+i+2n=9
は後回しで残りは
e+n+o=1
o+t+w=2
e+n+t=10
3e+l+n+v=11
2e+l+t+v+w=12
残った方程式で再度個数を数えると
2個: l,o,v,w
3個: n,t
4個: e
2の式から1の式を引いてoを消去
t+w-e-n=1
12の式からこの式を引いてwを消去
3e+l+n+v=11
これは11の式と同じなので残った式は
e+n+t=10
3e+l+n+v=11
nとtを固定して
e=10-n-t
11の式に代入してeを消去
l+v=2n+3t-19
後回しにした式にvは登場しlは登場しないのでvを固定して
l=2n+3t-v-19
12の式のeにe=10-n-tとl=2n+3t-v-19を代入してwを算出すると
w=11-2t
2の式にw=11-2tを代入してoを算出すると
o=t-9
9の式から
i=t-n-1
7の式から
s=2n+2t-n-v-13
5の式から
f=2n-v-4
3の式に未知数h,rが同時に登場するので登場回数の多いrを固定して
h=2n-r+t-17
後は最初に後回しにした4式から
z=n-r-1
u=v-2n-r-t+17
x=v-3t+20
g=r-2t+16
以上から
n,r,t,vは固定
e=10-n-t
f=2n-v-4
g=r-2t+16
h=2n-r+t-17
i=t-n-1
l=2n+3t-v-19
o=t-9
s=n+2t-v-13
u=v-2n-r-t+17
w=11-2t
x=v-3t+20
z=n-r-1
元の式に代入すると0～12が出てすべて正しいので、後は
条件(-11～11)を満たすようにn,r,t,vを定めればよいのだが、
解は多数ありそうなので適当な解一つだけにする。
絶対値が最小になるように適当に値を決めると
(e,f,g,h,i,l,n,o,r,s,t,u,v,w,x,z)=
(1,4,5,-3,0,4,4,-4,-1,1,5,5,0,1,5,4)
※o+t=9なので絶対値を4以下にするのは不可能

もし異なるアルファベットには異なる整数(-11～11も含む)という条件が加わると
どれほどの組合せの可能性が発生するもんですか？

その場合は3通りですね。
(e,f,g,h,i,l,n,o,r,s,t,u,v,w,x,z)=
(-2,-6,0,-7,7,9,2,1,4,3,10,5,6,-9,-4,-3),
(-1,-4,5,-11,6,8,2,0,7,3,9,1,4,-7,-3,-6),
(3,9,6,1,-4,0,5,-7,-6,-1,2,8,-3,7,11,10)
-10～10ならば1通りでした。

16個の変数で13個の方程式から通常では変数の中の3つを固定して(定数とみる。)
13変数の連立方程式として解こうとしていきたいが、これを行列Mを利用して
いくとき、正にその係数を基とする行列がmatdet(M)=0 となってしまい、
またmatrank(M)を調べたとき12を返されたのはこの事だったのですね。
でもどの2つの式から、既存の式が産み出せるのかわからなかった。
全ての流れを詳しく示して頂き目的の組合わせが3つも知れたのはラッキーでした。

追伸
この技をいろいろ試していたら
[ E, F, G, H, I, L, N, O, R, S, T, U, V, W, X, Z]=
1;[3/4, 15/4, 3, -11/4, 5/4, 6, 7/2, -13/4, -3/2, 11/4, 23/4, 5, -3/4, -1/2, 2, 4]
2;[9/4, 17/4, 3, -17/4, 7/4, 5, 5/2, -15/4, -5/2, 13/4, 21/4, 6, -13/4, 1/2, 1, 4]
などの分数による対応でも可能な組合せも生まれてきました。

こちらのオハナシの続編ですか？

http://shochandas.xsrv.jp/mathbun/mathbun1164.html

200までの素数での全解(上にある解も含む)
5個組
[13,59,67,131,163] と [23,31,103,109,167]
[11,59,71,149,173] と [23,29,101,131,179]
[19,79,101,173,191] と [23,61,131,149,199]
[31,67,103,149,197] と [37,53,127,131,199]
6個組
[19,29,53,73,97,107] と [17,37,43,83,89,109]
[19,29,83,103,157,167] と [13,47,59,127,139,173]
[43,47,101,109,163,167] と [37,59,83,127,151,173]
[29,31,103,107,179,181] と [19,53,71,139,157,191]
[43,53,107,127,181,191] と [37,71,83,151,163,197]
ちなみに200を超えた次の解は
5個組
[61,79,151,197,227] と [67,71,157,191,229]
6個組
[19,53,89,157,193,227] と [17,67,73,173,179,229]

# 昔作ったプログラムがまだ残っていました。

過去この話題についての投稿があっていましたね。
一般にProuhet-Tarry-Escott problem と呼ばれることもあり、特に使う数を素数に限定するものが
何か特別に見えて面白いと思っていました。
と言うのも一般での整数では公式が存在できるので、幾つも組合せが発見できるが素数ではそうはいかなくなる。
ふと過去のノートを整理していたら、この素数に関する2組の解を見ていたら、他の解は有るんのか？
の疑問がわき実際にパソコンで探してみたら思ったものより多くの組合せが存在していることに驚いたのでした。
本に紹介されているのは、特にその中にある最小数での組合せとなっていることにしかないのかと認識でき
では探せるだけ探してみようと挑戦を始めてみたのが投稿の動機でした。
なにせ冪数が高まれば高まるだけ探索範囲が指数関数的に増大していき、ちょっとやそっとでは有限時間では
探し出せなくなって壁に突き当たってしまった状態になっております。

らすかるさんありがとうございます。
３乗までの発見できるプログラム（一分位での計算時間では終了した。)
の延長の意味での4乗での検索では3日計算し続けても一つも発見できずにいました。
200までの素数での結果を知るまでの時間はどれほどなんですか？

結果を点検していたら
5乗までの和を等しく6個組の最小組
S1=[19,29,53,73,97,107] と
S2=[17,37,43,83,89,109]
に対し
M1=[-22,-17,-5,5,17,22]
M2=[-23,-13,-10,10,13,23]
なる対照的配列と
q=2,r=63
の定数を選べば
S1[n]=q*M1[n]+r
S2[n]=q*M2[n]+r
(n=1,2,3,4,5,6)
の関係で結ばれるようです。

また
S1=[19,53,89,157,193,227]
S2=[17,67,73,173,179,229]
なら
M1=[-52,-35,-17,17,35,52]
M2=[-53,-28,-25,25,28,53]
q=2,r=123
となるようです。

時間測ってなかったので再度実行して確かめたところ、
5個組で20秒、6個組で3分半でした。

38123*45097=1719232931 ですか？

あと
56809*76529=4347535961
78623*91243=7173798389
78443*94079=7379838997

(1) 256/693
(2)1/2772
(3)11411/2520-2*log(2)-Pi
だと思います。

(1),(2)
は
一般に
∫[x=0,1](x*(1-x))^ndx=∫[x=0,1](1-x^2)^ndx/4^n
=Beta(n+1,n+1)
但し
Beta(s,t)=Γ(s)*Γ(t)/Γ(s+t)
(ベータ関数)
で繋がっていて
gp > bestappr(intnum(x=0,1,(1-x^2)^5))
%470 = 256/693
gp > bestappr(intnum(x=0,1,(1-x^2)^5)/4^5)
%471 = 1/2772
gp > bestappr(intnum(x=0,1,(x*(1-x))^5))
%472 = 1/2772
gp > Beta(s,t)=gamma(s)*gamma(t)/gamma(s+t);
gp > bestappr(Beta(6,6))
%475 = 1/2772


(3)は
∫[x=0,1]x^n/(1+x^2)dx=∫[t=0,Pi/4]tan(t)^ndt
の置換積分で
そこで
gp > (x*(1-x))^5
%476 = -x^10 + 5*x^9 - 10*x^8 + 10*x^7 - 5*x^6 + x^5
=(5*x^9+10*x^7+x^5)-(x^10+10*x^8+5*x^6)

ここで
I=∫[x=0,1](x*(1-x))^5/(1+x^2)dx
ここで
x=tan(t) と置くとdx=dt/cos(t)^2=dt*(1+tan(t)^2)=dt*(1+x^2)
x=0-->t=0 ; x=1-->t=Pi/4
よって
I=∫[t=0,Pi/4](5*tan(t)^9+10*tan(t)^7+tan(t)^5)dt
-∫[t=0,Pi/4](tan(x)^10+10*tan(t)^8+5*tan(t)^6)dt
=5*(1/2*log(2)-7/24)+10*(-1/2*log(2)+5/12)+(1/2*log(2)-1/4)
-(-1/4*Pi+263/315)-10*(1/4*Pi-76/105)-5*(-1/4*Pi+13/15)
=11411/2520-2*log(2)-Pi
(この計算はほとほとめんどくさい。)

gp > intnum(x=0,1,(x*(1-x))^5/(1+x^2))
%480 = 0.00028758846491931730606697697874715425500493507898685
gp > 11411/2520-2*log(2)-Pi
%481 = 0.00028758846491931730606697697874715425500493507898685

私もベータ関数で考えました。

(1)については、

∫[x=0→1](1-x^2)^5dx=(1/2)∫[x=-1→1]((1+x)*(1-x))^5dx
=∫[t=0→1](2t*2(1-t))^5dt=2^10*Β(6,6)=2^10*Γ(6)^2/Γ(12)
=2^10*(5!)^2/(11!)=1024/2772=256/693
となりますが、

(2)については、

∫[x=0→1](x*(1-x))^5dx=∫[x=0→1](x^5*(1-x)^5)dx
=Β(6,6)=Γ(6)^2/Γ(12)=(5!)^2/(11!)=(5*4*3*2*1)/(11*10*9*8*7*6)
=1/(11*2*3*2*7*3)=1/2772
となるのではないでしょうか。

(3)については、∫[x=0→1](x*(1-x))^n/(1+x^2)dx (n=1,2,3,4,5)を順次求めてみました。

(x(1-x))/(1+x^2)=(x-x^2)/(1+x^2)=(x+1)/(1+x^2)-1
より
∫[x=0→1](x*(1-x))/(1+x^2)dx=∫[x=0→1]((x+1)/(1+x^2)-1)dx
=ln(2)/2+π/4-1

(x(1-x))^2/(1+x^2)=((x+1)/(1+x^2)-1)*(x(1-x))
((x+1)x(1-x))/(1+x^2)=(x-x^3)/(1+x^2)=2x/(1+x^2)-x
より
∫[x=0→1](x*(1-x))^2/(1+x^2)dx=∫[x=0→1](2x/(1+x^2)-x-x(1-x))dx
=ln(2)-Β(2,1)-Β(2,2)=ln(2)-Γ(2)Γ(1)/Γ(3)-Γ(2)^2/Γ(4)
=ln(2)-(1!*0!)/2!-(1!)^2/(3!)=ln(2)-1/2-1/6=ln(2)-2/3

(x(1-x))^3/(1+x^2)=(2x/(1+x^2)-x-x(1-x))*(x(1-x))
(x^2(1-x))/(1+x^2)=(x^2-x^3)/(1+x^2)=(x-1)/(1+x^2)-(x-1)
より
∫[x=0→1](x*(1-x))^3/(1+x^2)dx
=∫[x=0→1](2(x-1)/(1+x^2)+2(1-x)-x^2(1-x)-x^2(1-x)^2)dx
=ln(2)-π/2+2Β(1,2)-Β(3,2)-Β(3,3)
=ln(2)-π/2+2Γ(1)Γ(2)/Γ(3)-Γ(3)Γ(2)/Γ(5)-Γ(3)^2/Γ(6)
=ln(2)-π/2+2*(0!1!)/2!-(2!1!)/4!-(2!)^2/5!
=ln(2)-π/2+1-1/12-1/30=ln(2)-π/2+53/60

(x(1-x))^4/(1+x^2)=(2(x-1)/(1+x^2)+2(1-x)-x^2(1-x)-x^2(1-x)^2)*(x(1-x))
-(x-1)^2*x/(1+x^2)=-x*(1-2x+x^2)/(1+x^2)=2x^2/(1+x^2)-x=-2/(1+x^2)+2-x
より
∫[x=0→1](x*(1-x))^4/(1+x^2)dx
=∫[x=0→1](-4/(1+x^2)+4-2x+2x(1-x)^2-x^3(1-x)^2-x^3(1-x)^3)dx
=-π+4Β(1,1)-2Β(2,1)+2Β(2,3)-Β(4,3)-Β(4,4)
=-π+4Γ(1)^2/Γ(2)-2Γ(2)Γ(1)/Γ(3)+2Γ(2)Γ(3)/Γ(5)-Γ(4)Γ(3)/Γ(7)-Γ(4)^2/Γ(8)
=-π+4(0!)^2/1!-2(1!0!)/2!+2(1!2!)/4!+(3!2!)/6!-(3!)^2/7!
=-π+4-1+1/6-1/60-1/140=-π+22/7

(x(1-x))^5/(1+x^2)=(-4/(1+x^2)+4-2x+2x(1-x)^2-x^3(1-x)^2-x^3(1-x)^3)*(x(1-x))
(x(1-x))/(1+x^2)=(x-x^2)/(1+x^2)=(x+1)/(1+x^2)-1
より
=∫[x=0→1](-4(x+1)/(1+x^2)+4+4x(1-x)-2x^2(1-x)+2x^2(1-x)^3-x^4(1-x)^3-x^4(1-x)^4)dx
=-2*ln(2)/2-π+4Β(1,1)+4Β(2,2)-2Β(3,2)+2Β(3,4)-Β(5,4)-Β(5,5)
=-2*ln(2)/2-π+4Γ(1)^2/Γ(2)+4Γ(2)^2/Γ(4)-2Γ(3)Γ(2)/Γ(5)+2Γ(3)Γ(4)/Γ(7)-Γ(5)Γ(4)/Γ(9)-Γ(5)^2/Γ(10)
=-2*ln(2)/2-π+4(0!)^2/1!+4(1!)^2/3!-2(2!1!)/4!+2(2!3!)/6!-(4!3!)/8!-(4!)^2/9!
=-2*ln(2)/2-π+4+2/3-1/6+1/30-1/280-1/630=-2*ln(2)-π+11411/2520

> "kuiperbelt"さんが書かれました:
> 私もベータ関数で考えました。

(3)もベータ関数に繋がるんですね。
参考にn=6,7,8,9,10
でやってみました。
I6=38429/13860-4*ln(2)
I7=2*π-4*ln(2)-421691/120120
I8=4*π-188684/15015
I9=8*ln(2)+4*π-17069771/942480
I10=16*ln(2)-1290876029/116396280

管理人さんの解答と数値が合わなくて悩んでましたが、やっぱり (2) は 1/2772 ですよね？

(1)と(2)はベータ関数なんて持ち出さなくても、5回部分積分すれば高校生でもすぐ答えられる問題ですね。

(1)
∫[x=0->1] (1-x^2)^5 dx
= (1/2) ∫[x=-1->1] (1+x)^5*(1-x)^5 dx
= (1/2)*(5/6) ∫[x=-1->1] (1+x)^6*(1-x)^4 dx
= (1/2)*(5/6)*(4/7) ∫[x=-1->1] (1+x)^7*(1-x)^3 dx
= (1/2)*(5/6)*(4/7)*(3/8) ∫[x=-1->1] (1+x)^8*(1-x)^2 dx
= (1/2)*(5/6)*(4/7)*(3/8)*(2/9) ∫[x=-1->1] (1+x)^9*(1-x) dx
= (1/2)*(5/6)*(4/7)*(3/8)*(2/9)*(1/10) ∫[x=-1->1] (1+x)^10 dx
= (1/2)*(5/6)*(4/7)*(3/8)*(2/9)*(1/10)*(1/11)*2^11
= 256/693

(2) (1) と同様にして
∫[x=0->1] (x*(1-x))^5 dx
= (5/6)*(4/7)*(3/8)*(2/9)*(1/10)*(1/11)*1^11
= 1/2772

(3)
まず、
(x*(1-x))^5 = (1+x^2)Q(x) + ax + b
とおき、x = ±i を代入すると、
-4 - 4i = b + ai
-4 + 4i = b - ai
となるので、a = b = -4

また、
∫[x=0->1] 1/(1+x^2) dx = π/4（x=tanθの置換積分による）
と
∫[x=0->1] 2x/(1+x^2) dx = log2
が成り立ちます。

よって、
∫[x=0->1] (x*(1-x))^5/(1+x^2) dx
= ∫[x=0->1] {x^5*(1-x)^5+4x+4}/(1+x^2) dx - 4 ∫[x=0->1] 1/(1+x^2) dx - 2 ∫[x=0->1] 2x/(1+x^2) dx
= Σ[n=0->∞] ∫[x=0->1] {x^(2n+5)*(1-x)^5+4*x^(2n+1)+4*x^(2n)}*(-1)^n dx - π - 2log2
= Σ[n=0->∞] {120/((2n+11)(2n+10)(2n+9)(2n+8)(2n+7)(2n+6))+4/(2n+2)+4/(2n+1)}*(-1)^n - π - 2log2
= Σ[n=0->∞] {-1/(2n+11)+5/(2n+10)-10/(2n+9)+10/(2n+8)-5/(2n+7)+1/(2n+6)+4/(2n+2)+4/(2n+1)}*(-1)^n - π - 2log2
= 4/1 + 4/2 - 4/3 - 4/4 + 4/5 + 5/6 - 9/7 + 5/8 - 1/9 - π - 2log2
= 11411/2520 - π - 2log2

最後の1行は手計算じゃ厳しい……。

11411/2520の正体が
4/1 + 4/2 - 4/3 - 4/4 + 4/5 + 5/6 - 9/7 + 5/8 - 1/9
であることには驚きました。
(x*(1-x))^6,(x*(1-x))^7
に現れる分数値
38429/13860、- 421691/120120
についてDD++さんの巧みな方法を参考に探すと
2 + 4/3 - 7/4 - 3/5 + 1/7 + 14/9 + 1/11 = 38429/13860
- 6 + 7/3 + 6/5 + 8/7 - 7/8 - 20/11 + 17/12 - 1/13 = - 421691/120120
によって構成されてくることになるんですね。
なおkuperbeltさんのベータ関数利用では
2 + 2/3 + 2/15 - 1/30 + 1/140 - 1/1260 - 1/2772 = 38429/13860
- 4 + 1/3 +2/15 + 1/35 - 1/140 + 1/630 - 1/5544 - 1/12012 = - 421691/120120
の構成になるようです。
しかしそもそも
∫[x=0,1]1/(x^2+1)dx=∑[n=0,oo]((-1)^n*∫[x=0,1]x^(2*n)dx)
∫[x=0,1]x/(x^2+1)dx=∑[n=0,oo](-1)^n*∫[x=0,1](x^(2*n+1)dx)
∫[x=0,1]x^2/(x^2+1)dx=∑[n=0,oo]((-1)^n*∫[x=0,1]x^(2*n+2)dx)
∫[x=0,1]x^s*(1-x)^sdx=∑[n=0,oo]((-1)^n*∫[x=0,1]x^(2*n+s)*(1-x)^sdx)
なんて式がどこから思いつけるんですか？

単純に、
1/(1+x^2) がいやだなあ……
→x か 1-x だけを因数に持つものの総和で書けたらなあ……
→初項 1 で公比 -x^2 の無限等比級数に展開すれば、それできるじゃん！
というだけです。

ところで、あとから気がついたんですが、普通に x^5*(1-x)^5 を 1+x^2 で割った商を普通に積分すれば同じ内容の分数の和に帰着するっぽいですね。
式変形頑張ったの、あんまり意味なかった疑惑。

19015*64926
でしょうか。

技術者に名前A、2、3、4、5、6、7、8、9、X、J、Q、Kを割り当てます。13枚の金貨を次のように名前付けします：

C_{A, 2, 6, Q}
C_{2, 3, 7, K}
C_{3, 4, 8, A}
C_{4, 5, 9, 2}
C_{5, 6, X, 3}
C_{6, 7, J, 4}
C_{7, 8, Q, 5}
C_{8, 9, K, 6}
C_{9, X, A, 7}
C_{X, J, 2, 8}
C_{J, Q, 3, 9}
C_{Q, K, 4, X}
C_{K, A, 5, J}

特性
1. 数学の概念である集合を用いて金貨の名前をインデックス付けしています。例えば、C_{K, A, 5, J}とC_{A, 5, K, J}は同じ金貨を指します。
2. 各技術者は、インデックスに自分の名前を含む4枚の金貨を独立して測定します。
3. 各金貨はちょうど4人の技術者によって測定されます。

測定プロセス
Tを、割り当てられた金貨について陽性結果を報告した技術者の集合とします。

Tの例
T = {A, 3, 4, 8, J, K}
T = {2, 4, 9, X, K}
T = {2, 8, X, Q}
T = {6, 8, 9}
T = {J, 9}

⬛⬜⬛⬜⬛⬜
上記の計測方法でうまくいくのだと
ハミング距離という概念の説明ぬきに
集合概念を知っている小学生に教えるにはどうしたらよいでしょうか？

計測の結果として得られた集合 T の要素となったおのおのの技術者が計測した４枚のコインに、技術者ごとに、偽物ポイントを 1 づつ与えます。

たとえば、T の要素に 2(技術者の名)があれば、
C_{A, 2, 6, Q}
C_{2, 3, 7, K}
C_{4, 5, 9, 2}
C_{X, J, 2, 8}
の４枚に偽物ポイントを 1 づつ付与します。

13 人の技術者のうち、放射線陽性の報告をあげた技術者の全てについて上記の操作をおこなった結果として、偽物ポイントが最大のものがユニークに定まります。それが偽金貨です。

⬛⬜⬛⬜⬛⬜

これがなぜうまくいくのか？
という疑問です。
小学生にもわかる説明はありますでしょうか？

申し遅れました。早見表は以下の通りです。

'A23456789XJQK'
'1100010000010', //C_{A, 2, 6, Q}
'0110001000001', //C_{2, 3, 7, K}
'1011000100000', //C_{3, 4, 8, A}
'0101100010000', //C_{4, 5, 9, 2}
'0010110001000', //C_{5, 6, X, 3}
'0001011000100', //C_{6, 7, J, 4}
'0000101100010', //C_{7, 8, Q, 5}
'0000010110001', //C_{8, 9, K, 6}
'1000001011000', //C_{9, X, A, 7}
'0100000101100', //C_{X, J, 2, 8}
'0010000010110', //C_{J, Q, 3, 9}
'0001000001011', //C_{Q, K, 4, X}
'1000100000101', //C_{K, A, 5, J}

最小ハミング距離は 6 となっています。

ハミング距離をよく知らない私が説明にトライ。
理論はわかってないので、こうすれば説明できるんじゃないかなと思うことを書いてみる。

全部のパターンをやってみて大丈夫だねというごり押し説明なので長いです。

------


金貨の名前が長いので、上から順に "あ","い",…,"す" としましょう。
もとの金貨の名前から、今度はそれぞれの技術者が測定する金貨の名前の集合を作ります。
ちょっと大変ですが、頑張れば次のようになります。

A:{あ,う,け,す}
2:{あ,い,え,こ}
3:{い,う,お,さ}
4:{う,え,か,し}
5:{え,お,き,す}
6:{あ,お,か,く}
7:{い,か,き,け}
8:{う,き,く,こ}
9:{え,く,け,さ}
X:{お,け,こ,し}
J:{か,こ,さ,す}
Q:{あ,き,さ,し}
K:{い,く,し,す}

ここでたとえば、金貨"い"を測定する技術者を集めてみます。
2:{あ,い,え,こ}
3:{い,う,お,さ}
7:{い,か,き,け}
K:{い,く,し,す}
すると、"い"以外の12枚の金貨はどれも、この4人の技術者(2,3,7,K)のなかのだれか1人だけが
測定していることがわかります。

同様に、ほかの金貨を測定する4人の技術者を見てみても、
残り12枚の金貨を4人のなかのだれか一人だけが測定するようになっています。
気になったらあとでそれぞれの金貨で確認してみてください。
このことがこの方法でうまくいく最大の理由です。



さて、まずは技術者全員が正確な報告をした場合を考えましょう。
この場合、4人が陽性(自分が測定した4枚中に偽物あり)と報告し、
残り9人が陰性(自分が測定した4枚中に偽物なし)と報告することになります。

たとえば金貨"い"が偽物だった場合、技術者{2,3,7,K}が陽性報告、他の技術者が陰性報告となります。
このとき、金貨"い"の偽物ポイントが4になり、残りの金貨の偽物ポイントは1になります。
なぜなら、先ほど確認したように、金貨"い"を測定する4人の技術者は、
"い"以外の12枚の金貨を重複することなくだれか1人だけが測定しているからです。

別の金貨が偽物だった場合でも同様に、偽物が4ポイント、その他の金貨が1ポイントになります。
よって、技術者全員が正確な報告をした場合は偽物を見分けられることがわかりました。



次に、不正確な報告をした技術者がいた場合を考えます。
これは次のように5パターンが考えられます。
・偽物を測定した4人のうちの1人が陽性(偽物あり)ではなく陰性(偽物なし)と報告した場合。
・偽物を測定した4人のうちの2人が陽性(偽物あり)ではなく陰性(偽物なし)と報告した場合。
・偽物を測定しなかった9人のうちの1人が陰性(偽物なし)ではなく陽性(偽物あり)と報告した場合。
・偽物を測定しなかった9人のうちの2人が陰性(偽物なし)ではなく陽性(偽物あり)と報告した場合。
・偽物を測定した4人のうちの1人が陽性(偽物あり)ではなく陰性(偽物なし)と報告し、
　　かつ、偽物を測定しなかった9人のうちの1人が陰性(偽物なし)ではなく陽性(偽物あり)と報告した場合。

順番に見ていきましょう。


偽物を測定した4人のうちの1人が陽性(偽物あり)ではなく陰性(偽物なし)と報告した場合は、
不正確な報告をした人が測定していた金貨の偽物ポイントが1ポイント減るので、
偽物金貨が3ポイント、その他の金貨は1ポイントまたは0ポイントになります。
よって偽物を見分けられます。

偽物を測定した4人のうちの2人が陽性(偽物あり)ではなく陰性(偽物なし)と報告した場合は、
不正確な報告をした人が共通して測定していた偽物のポイントが2ポイント減り、
不正確な報告をした人が測定していた他の金貨のポイントは1ポイント減るので、
偽物金貨が2ポイント、その他の金貨は1ポイントまたは0ポイントになり、偽物を見分けられます。


偽物を測定しなかった9人のうちの1人が陰性(偽物なし)ではなく陽性(偽物あり)と報告した場合は、
不正確な報告をした人が測定していた金貨の偽物ポイントが1ポイント足されることになります。
不正確な報告をした人が測定したのは偽物でない金貨なので、
偽物金貨は4ポイント、その他の金貨は1ポイントまたは2ポイントになります。
よって偽物を見分けられます。

偽物を測定しなかった9人のうちの2人が陰性(偽物なし)ではなく陽性(偽物あり)と報告した場合は、
不正確な報告をした人2人が共通して測定していた金貨は2ポイント追加され、
不正確な報告をした人のうちどちらか1人だけが測定していた金貨は1ポイントが追加されるので、
偽物金貨は4ポイント、その他の金貨は1ポイントまたは2ポイントまたは3ポイントになり、
偽物を見分けられます。


最後に、偽物を測定した4人のうちの1人が陽性(偽物あり)ではなく陰性(偽物なし)と報告し、
偽物を測定しなかった9人のうちの1人が陰性(偽物なし)ではなく陽性(偽物あり)と報告した場合を考えます。
陽性を陰性と報告した人が測定した金貨は1ポイント減少し、
陰性を陽性と報告した人が測定した金貨は1ポイント増加します。
これにより、偽物金貨は3ポイント、その他の金貨は1ポイントまたは0ポイントまたは2ポイントになるので
偽物を見分けられます。



以上により、不正確な報告を行った技術者が最大2人いたとしても、偽物を特定できることがわかります。


------

さて、小学生には難解すぎる説明になってしまったうえに、
これがなぜうまくいくのか？の「なぜ」の部分にちゃんと答えられていないことから、
Dengan kesaktian Indukmu さんの質問に答えたとはいえないでしょうね……。
すみません。

りらひいさんは「誰がどんなタイプの嘘をついたか事前にわかっている前提」で書いているような？
Dengan さんは、「誰がどんなタイプの嘘をついたかわかっていない前提」での話を尋ねているのではないかと思います。

説明はこんな感じでしょうか？


仮に全員が正直に報告した場合、本物は偽物ポイント1、偽物は偽物ポイント4になります。

ここで、研究者の誰かが嘘をついたとします。
もし、陽性を報告するはずの研究者が誤って陰性と報告した場合、該当する4枚の偽物ポイントが誤って1下がります。
逆に、もし、陰性を報告するはずの研究者が誤って陽性と報告した場合、該当する4枚の偽物ポイントが誤って1上がります。

すなわち、研究者が1人報告を誤るごとに、2枚のコインの偽物ポイントの差が誤って1変化する可能性があります。
……が、本物と偽物の偽物ポイント差は本来3ポイントあるので、誤りが2人までならその大小関係が逆転することはありません。
したがって、どんな場合でも偽物ポイントが最大であるものが偽物で確定です。

りらひいさん、DD++さん。
まことに有難うございました。おかげさまでスッキリいたしました。
容疑ポイントで逆転がおきえないことをスマートに説明することができるようになりました。

>残り12枚の金貨を4人のなかのだれか一人だけが測定するようになっています。

の視点も重要と存じます。
というのは、追加で14枚めのコインがあったとして、全技術者が必ずこのコインを計測する、その他のコインについてはここまでの従前通り、としても、
最小ハミング距離は 6 のままですので、やはり偽コインの特定は成功することとなります。
でも、この場合に今回の仕様の容疑者ポイントを使った方法では、なんだかとてもめんどくさいことになります。別の手が欲しいところですが、なかなかです。

考え直して、りらひいさんに誤った指摘をしてしまっていたことに気づきました。
証明はどんな嘘が混ざったかを前提にしてますが、最終的に判断する方法は「偽物ポイントが最大のもの」で共通しているので、嘘の混ざり方がわからないまま実行してOKな手順なんですねコレ。
失礼しました。

DD++さん

わたしも言葉が足りていなかったと反省しています。

技術者全員が正確な報告をした場合について書き始める前に、
「不正確な報告がどのような場合であっても偽物を特定できることを示すために、
　不正確な報告を分類してすべての場合で偽物を判別できるかを個別に考える。」
みたいな一文があればよかったのでしょう。
それから、「さて、まずは…」という書き始めも誤解を生んだ可能性がありますね。
言葉の選び方が悪かった面もあると思います。

DD++さんのお言葉から判断するに、
偽物の判断基準を私の投稿内で明示しなかったこともよくなかったみたいですね。
偽物の判断基準はDengan kesaktian Indukmuさんの投稿に示されていて
共通認識になっていると思い端折ってしまいました。
すみませんでした。


---------


さて、私が投稿した説明についてもう少し反省を……。

まず、後半の場合分け自体が過剰な回答でしたね。
DD++さんの説明が簡潔でわかりやすいと思います。

そして、前半の集合の再構築も特に不要でした。
「ある1つの金貨を測定する4人の技術者を集めたとき、
　残り12枚の金貨はそれぞれこの4人のなかのだれか一人だけが測定する」
ということを示すには各技術者が測定する金貨の名前の集合を作った方がわかりやすいと思ったのですが、
もとのDengan kesaktian Indukmuが書かれた金貨の名前のリストだけで十分に確認可能でした。
その説明は次のような感じです。
***
C_{2, 3, 7, K}を測定する4人の技術者2,3,7,Kが測定するほかの金貨について考える。
C_{2, 3, 7, K}以外の金貨の名前を見てみると、どれも2,3,7,Kのなかのどれか一つだけを含むことがわかる。
これは、C_{2, 3, 7, K}以外の12枚の金貨をこの4人のなかのだれか一人だけが測定することを意味する。
ほかの金貨の場合でも同様のことが成り立つ。
***
こんな感じの説明でいけました。
金貨に名前を付けなおして集合を作るくだりはいらなかったですね。


やっぱり思いついたことをそのまま書き連ねた後に大して推敲せずすぐに投稿すると、
後々直した方がいい箇所がいろいろ見つかりますね。

小学生向けの説明ではありませんが、13枚の金貨への13人の技術者の割り当ては、以下のようにGF(3)上の射影平面の9点に9人の技術者が対応し、残りの4人がGF(3)上の射影平面の無限遠点に対応しますね。

A:(0,0,1) 2:(1,0,1) Q:(2,0,1)
5:(0,1,1) 3:(1,1,1) X:(2,1,1)
J:(0,2,1) 7:(1,2,1) 4:(2,2,1)

6:(1,0,0),K:(0,1,0),8:(1,1,0),9:(2,1,0)は、直線A-2-Q,A-5-J,A-3-4,A-7-Xの延長上の無限遠点に対応します。

GF(3)上の射影平面の直線は
A-2-Q-6,5-3-X-6,J-7-4-6,
A-5-J-K,2-3-7-K,Q-X-4-K,
A-3-4-8,2-X-J-8,Q-5-7-8,
A-7-X-9,5-2-4-9,J-3-Q-9,
6-8-K-9

の13本ですが、13本の直線がそれぞれ13枚の金貨に相当します。
13人の技術者が全て真の報告をするのであれば、偽金貨が1枚のときは4人が陽性、9人が陰性の報告をすることになります。

13人の技術者のうち、最大2人が偽の報告をする場合、偽の報告を偽陽性、偽陰性とすると、13人からの報告は

①3人が陽性、1人が偽陰性、9人が陰性
②2人が陽性、2人が偽陰性、9人が陰性
③4人が陽性、1人が偽陽性、8人が陰性
④4人が陽性、2人が偽陽性、7人が陰性
⑤3人が陽性、1人が偽陰性、1人が偽陽性、8人が陰性

となります。

①②の場合
GF(3)上の射影平面上の直線は、2点が特定されると一意に決まるので、偽金貨を特定することができます。

③④の場合
GF(3)上の射影平面上の直線には4点が含まれ、陽性の報告をした技術者は同一直線上の4点に相当し、偽陽性の報告をした技術者は、当該直線外の点に相当します。偽陽性の報告をした技術者が2人までなので、直線外の点は2個しかなく、偽金貨に相当する直線以外に4点が並ぶことがないので偽金貨を特定することができます。

⑤の場合
陽性の報告をした技術者は同一直線上の3点に相当し、偽陽性の報告をした技術者は、当該直線外の1点に相当します。偽金貨に相当する直線上の3点のうち2点と当該直線外の1点を同時に通る直線というのは存在しないので、偽陽性を報告した技術者に相当する点を結ぶ直線上には、陽性あるいは偽陽性を報告した技術者に相当する点が2点しか存在しないことになります。こうして、陽性あるいは偽陽性の報告した技術者に相当する点のうち3点が並ぶ直線に相当する金貨が偽金貨と特定されます。

はい。
PG(2,3)および同じ構造ですが(2,3) Singerは以下の特徴を共有します：
点の数: 13点
直線の数: 13直線
各直線上の点の数: 4点
各点を通る直線の数: 4直線

もう、どこにぶら下げたらよいやら分からなくなりましたので

13個の金貨13人の技術者に投稿します。

技術者が13人いるのならば
金貨は88枚でいけるでしょ？
と教えられました。
まだ理解していません。

単に88×88行列で各行各列が1が13個ずつ配列できるものは構成可能ですが、これが限界という意味が私には理解できません。
(コンピュータによる出力のため2を0に読み替えて見てください。)

%01 = 2121221222122221222221222222122222212222222122222222122222222212222222222122222222222122
%02 = 2212122122212222122222122222212222221222222212222222212222222221222222222212222222222212
%03 = 2221212212221222212222212222221222222122222221222222221222222222122222222221222222222221
%04 = 1222121221222122221222221222222122222212222222122222222122222222212222222222122222222222
%05 = 2122212122122212222122222122222212222221222222212222222212222222221222222222212222222222
%06 = 2212221212212221222212222212222221222222122222221222222221222222222122222222221222222222
%07 = 2221222121221222122221222221222222122222212222222122222222122222222212222222222122222222
%08 = 2222122212122122212222122222122222212222221222222212222222212222222221222222222212222222
%09 = 2222212221212212221222212222212222221222222122222221222222221222222222122222222221222222
%10 = 2222221222121221222122221222221222222122222212222222122222222122222222212222222222122222
%11 = 2222222122212122122212222122222122222212222221222222212222222212222222221222222222212222
%12 = 2222222212221212212221222212222212222221222222122222221222222221222222222122222222221222
%13 = 2222222221222121221222122221222221222222122222212222222122222222122222222212222222222122
%14 = 2222222222122212122122212222122222122222212222221222222212222222212222222221222222222212
%15 = 2222222222212221212212221222212222212222221222222122222221222222221222222222122222222221
%16 = 1222222222221222121221222122221222221222222122222212222222122222222122222222212222222222
%17 = 2122222222222122212122122212222122222122222212222221222222212222222212222222221222222222
%18 = 2212222222222212221212212221222212222212222221222222122222221222222221222222222122222222
%19 = 2221222222222221222121221222122221222221222222122222212222222122222222122222222212222222
%20 = 2222122222222222122212122122212222122222122222212222221222222212222222212222222221222222
%21 = 2222212222222222212221212212221222212222212222221222222122222221222222221222222222122222
%22 = 2222221222222222221222121221222122221222221222222122222212222222122222222122222222212222
%23 = 2222222122222222222122212122122212222122222122222212222221222222212222222212222222221222
%24 = 2222222212222222222212221212212221222212222212222221222222122222221222222221222222222122
%25 = 2222222221222222222221222121221222122221222221222222122222212222222122222222122222222212
%26 = 2222222222122222222222122212122122212222122222122222212222221222222212222222212222222221
%27 = 1222222222212222222222212221212212221222212222212222221222222122222221222222221222222222
%28 = 2122222222221222222222221222121221222122221222221222222122222212222222122222222122222222
%29 = 2212222222222122222222222122212122122212222122222122222212222221222222212222222212222222
%30 = 2221222222222212222222222212221212212221222212222212222221222222122222221222222221222222
%31 = 2222122222222221222222222221222121221222122221222221222222122222212222222122222222122222
%32 = 2222212222222222122222222222122212122122212222122222122222212222221222222212222222212222
%33 = 2222221222222222212222222222212221212212221222212222212222221222222122222221222222221222
%34 = 2222222122222222221222222222221222121221222122221222221222222122222212222222122222222122
%35 = 2222222212222222222122222222222122212122122212222122222122222212222221222222212222222212
%36 = 2222222221222222222212222222222212221212212221222212222212222221222222122222221222222221
%37 = 1222222222122222222221222222222221222121221222122221222221222222122222212222222122222222
%38 = 2122222222212222222222122222222222122212122122212222122222122222212222221222222212222222
%39 = 2212222222221222222222212222222222212221212212221222212222212222221222222122222221222222
%40 = 2221222222222122222222221222222222221222121221222122221222221222222122222212222222122222
%41 = 2222122222222212222222222122222222222122212122122212222122222122222212222221222222212222
%42 = 2222212222222221222222222212222222222212221212212221222212222212222221222222122222221222
%43 = 2222221222222222122222222221222222222221222121221222122221222221222222122222212222222122
%44 = 2222222122222222212222222222122222222222122212122122212222122222122222212222221222222212
%45 = 2222222212222222221222222222212222222222212221212212221222212222212222221222222122222221
%46 = 1222222221222222222122222222221222222222221222121221222122221222221222222122222212222222
%47 = 2122222222122222222212222222222122222222222122212122122212222122222122222212222221222222
%48 = 2212222222212222222221222222222212222222222212221212212221222212222212222221222222122222
%49 = 2221222222221222222222122222222221222222222221222121221222122221222221222222122222212222
%50 = 2222122222222122222222212222222222122222222222122212122122212222122222122222212222221222
%51 = 2222212222222212222222221222222222212222222222212221212212221222212222212222221222222122
%52 = 2222221222222221222222222122222222221222222222221222121221222122221222221222222122222212
%53 = 2222222122222222122222222212222222222122222222222122212122122212222122222122222212222221
%54 = 1222222212222222212222222221222222222212222222222212221212212221222212222212222221222222
%55 = 2122222221222222221222222222122222222221222222222221222121221222122221222221222222122222
%56 = 2212222222122222222122222222212222222222122222222222122212122122212222122222122222212222
%57 = 2221222222212222222212222222221222222222212222222222212221212212221222212222212222221222
%58 = 2222122222221222222221222222222122222222221222222222221222121221222122221222221222222122
%59 = 2222212222222122222222122222222212222222222122222222222122212122122212222122222122222212
%60 = 2222221222222212222222212222222221222222222212222222222212221212212221222212222212222221
%61 = 1222222122222221222222221222222222122222222221222222222221222121221222122221222221222222
%62 = 2122222212222222122222222122222222212222222222122222222222122212122122212222122222122222
%63 = 2212222221222222212222222212222222221222222222212222222222212221212212221222212222212222
%64 = 2221222222122222221222222221222222222122222222221222222222221222121221222122221222221222
%65 = 2222122222212222222122222222122222222212222222222122222222222122212122122212222122222122
%66 = 2222212222221222222212222222212222222221222222222212222222222212221212212221222212222212
%67 = 2222221222222122222221222222221222222222122222222221222222222221222121221222122221222221
%68 = 1222222122222212222222122222222122222222212222222222122222222222122212122122212222122222
%69 = 2122222212222221222222212222222212222222221222222222212222222222212221212212221222212222
%70 = 2212222221222222122222221222222221222222222122222222221222222222221222121221222122221222
%71 = 2221222222122222212222222122222222122222222212222222222122222222222122212122122212222122
%72 = 2222122222212222221222222212222222212222222221222222222212222222222212221212212221222212
%73 = 2222212222221222222122222221222222221222222222122222222221222222222221222121221222122221
%74 = 1222221222222122222212222222122222222122222222212222222222122222222222122212122122212222
%75 = 2122222122222212222221222222212222222212222222221222222222212222222222212221212212221222
%76 = 2212222212222221222222122222221222222221222222222122222222221222222222221222121221222122
%77 = 2221222221222222122222212222222122222222122222222212222222222122222222222122212122122212
%78 = 2222122222122222212222221222222212222222212222222221222222222212222222222212221212212221
%79 = 1222212222212222221222222122222221222222221222222222122222222221222222222221222121221222
%80 = 2122221222221222222122222212222222122222222122222222212222222222122222222222122212122122
%81 = 2212222122222122222212222221222222212222222212222222221222222222212222222222212221212212
%82 = 2221222212222212222221222222122222221222222221222222222122222222221222222222221222121221
%83 = 1222122221222221222222122222212222222122222222122222222212222222222122222222222122212122
%84 = 2122212222122222122222212222221222222212222222212222222221222222222212222222222212221212
%85 = 2212221222212222212222221222222122222221222222221222222222122222222221222222222221222121
%86 = 1221222122221222221222222122222212222222122222222122222222212222222222122222222222122212
%87 = 2122122212222122222122222212222221222222212222222212222222221222222222212222222222212221
%88 = 1212212221222212222212222221222222122222221222222221222222222122222222221222222222221222

GAI さん。ありがとうございます。

小さなモデルでは次のようになります。
ガイガーカウンター管 で 11 回計測する。(技術者ひとりにつき１回、技術者の総数は 11 、うち０人~２人は嘘の報告をする)
このとき、22 枚の金貨のうち 1 枚の放射性物質を含む偽金貨を特定できる。
計測の準備には2通りの方法がある。
①少しづつ金貨を計測していき、その計測結果をもとにして次の計測の段取りをつける
②実際の計測の前に11人の技術者が計測する金貨を決定してしまう。

おとついみつけてたまげたのですが、下記のように計測すれば②で22枚中1枚を特定できます。

10111000100,
01011100010,
00101110001,
10010111000,
01001011100,
00100101110,
00010010111,
10001001011,
11000100101,
11100010010,
01110001001,
01000111011,
10100011101,
11010001110,
01101000111,
10110100011,
11011010001,
11101101000,
01110110100,
00111011010,
00011101101,
10001110110,

たまげた理由。
②のやりかたのシバリでは、kuiperbelt さんも、私も、最大 16 枚の金貨中から偽金貨1枚をみつける方法を過去に投稿して来ました。
符号長11で最小ハミング距離が 5 以上の系では、符号の最大総数は 16 だろうという雰囲気で。ところが上ではそれを6枚ぶん上回りました。
もっというと、all 0 と all 1 を勘定にいれれば 24 枚となるところに驚いております。

で、今回の 88 枚の件ですが、
技術者を13人とし、計測の方法は ① とする、ということとなります。嘘つきの数は変えません。少しづつ、計測しながら、その結果をみながら、次に測定する金貨を定めていく、しかも、いつ嘘をつかれるか不定である、そのような状況下で、88 枚も処理できるなんて！？？

この 88 は 先ほどの 22 の 4 倍になっています。なにか関係があるかどうか探っています。

wikipedia の プロトキン限界 の項目では
符号長 11 最小ハミング距離が 5 の二元コードの符号語数が 最大 24 となり、
やった！！！
みつけたぞ！！！
と大喜びしておりましたが、
これの限界の公式は、どうやら間違いのようで。ぐはっ

方法①ならこんな感じでどうしょう。
行き当たりばったりで組んだので汚いですが……。



金貨を4つのグループA,B,C,Dに振り分けながら測定を行っていく。

・グループA：偽金貨だった場合に現在までの全員が正しく報告をしている金貨の集合（偽金貨である可能性がある）
・グループB：偽金貨だった場合に現在までの一人が誤った報告をしている金貨の集合（偽金貨である可能性がある）
・グループC：偽金貨だった場合に現在までの二人が誤った報告をしている金貨の集合（偽金貨である可能性がある）
・グループD：偽金貨だった場合に現在までの三人が誤った報告をしている金貨の集合（偽金貨である可能性はない）

グループA,B,C,Dは、最上位:A→B→C→D:最下位となるように順位付けしておく。

各技術者の測定結果に応じて次のように金貨のグループを変更する。
陽性⇒測定した金貨は今のグループに残す。測定しなかった金貨を一つ下位のグループに移す(グループDの金貨を除く)。
陰性⇒測定した金貨を一つ下位のグループに移す。測定しなかった金貨は今のグループに残す。

最初は88枚すべてグループAである。
測定を繰り返してグループDが87枚になったら、残った1枚が偽金貨である。

測定方法は以下の通り(あくまでも一例である)。

丸数字は測定方法の場合分けであり、前の人までの測定結果によりどれかを選択していく。
丸数字のあとのA,B,C,Dの数字は測定前の各グループの金貨の枚数。
各グループからそれぞれ決められた枚数を選んでまとめて測定する。
グループ内でどの金貨を選ぶかは自由。

[1人目]
① A:88, B:0, C:0, D:0 ⇒ Aから44枚 測定
　・陽性/陰性 → 2人目①へ

[2人目]
① A:44, B:44, C:0, D:0 ⇒ Aから22枚、Bから22枚 測定
　・陽性/陰性 → 3人目①へ

[3人目]
① A:22, B:44, C:22, D:0 ⇒ Aから11枚、Bから22枚、Cから11枚 測定
　・陽性/陰性 → 4人目①へ

[4人目]
① A:11, B:33, C:33, D:11 ⇒ Aから6枚、Bから15枚、Cから12枚 測定
　・陽性 → 5人目①へ
　・陰性 → 5人目②へ

[5人目]
① A:6, B:20, C:30, D:32 ⇒ Aから3枚、Bから10枚、Cから15枚 測定
　・陽性/陰性 → 6人目①へ

② A:5, B:24, C:36, D:23 ⇒ Aから3枚、Bから11枚、Cから12枚 測定
　・陽性 → 6人目①へ
　・陰性 → 6人目②へ

[6人目]
① A:3, B:13, C:25, D:47 ⇒ Aから2枚、Bから5枚、Cから12枚 測定
　・陽性 → 7人目①へ
　・陰性 → 7人目②へ

② A:2, B:16, C:35, D:35 ⇒ Aから1枚、Bから9枚、Cから11枚 測定
　・陽性 → 7人目②へ
　・陰性 → 7人目③へ

[7人目]
① A:2, B:6, C:20, D:60 ⇒ Aから1枚、Bから3枚、Cから10枚 測定
　・陽性/陰性 → 8人目①へ

② A:1, B:10, C:18, D:59 ⇒ Aから1枚、Bから4枚、Cから7枚 測定
　・陽性 → 8人目①へ
　・陰性 → 8人目②へ

③ A:1, B:8, C:33, D:46 ⇒ Aから1枚、Bから4枚、Cから9枚 測定
　・陽性 → 8人目①へ
　・陰性 → 8人目③へ

[8人目]
① A:1, B:4, C:13, D:70 ⇒ Aから1枚、Bから1枚、Cから6枚 測定
　・陽性 → 9人目①へ
　・陰性 → 9人目②へ

② A:0, B:7, C:15, D:66 ⇒ Bから4枚、Cから5枚 測定
　・陽性 → 9人目②へ
　・陰性 → 9人目③へ

③ A:0, B:5, C:28, D:55 ⇒ Bから3枚、Cから12枚 測定
　・陽性 → 9人目③へ
　・陰性 → 9人目④へ

[9人目]
① A:1, B:1, C:9, D:77 ⇒ Aから1枚、Bから0枚、Cから3枚 測定
　・陽性 → 10人目①へ
　・陰性 → 10人目②へ

② A:0, B:4, C:8, D:76 ⇒ Bから2枚、Cから4枚 測定
　・陽性/陰性 → 10人目②へ

③ A:0, B:3, C:14, D:71 ⇒ Bから2枚、Cから5枚 測定
　・陽性 → 10人目②へ
　・陰性 → 10人目③へ

④ A:0, B:2, C:19, D:67 ⇒ Bから1枚、Cから10枚 測定
　・陽性 → 10人③目へ
　・陰性 → 10人④目へ

[10人目]
① A:1, B:0, C:4, D:83 ⇒ Aから1枚 測定
　・陽性 → 結論①へ
　・陰性 → 11人目①へ

② A:0, B:2, C:6, D:80 ⇒ Bから1枚、Cから3枚 測定
　・陽性/陰性 → 11人目①へ

③ A:0, B:1, C:11, D:76 ⇒ Bから1枚、Cから4枚 測定
　・陽性 → 11人目①へ
　・陰性 → 11人目③へ

④ A:0, B:1, C:10, D:77 ⇒ Bから1枚、Cから3枚 測定
　・陽性 → 11人目②へ
　・陰性 → 11人目③へ

[11人目]
① A:0, B:1, C:4, D:83 ⇒ Bから1枚、Cから1枚 測定
　・陽性 → 12人目①へ
　・陰性 → 12人目②へ

② A:0, B:1, C:3, D:84 ⇒ Bから1枚、Cから1枚 測定
　・陽性 → 12人目①へ
　・陰性 → 12人目③へ

③ A0:, B:0, C:8, D:80 ⇒ Cから4枚 測定
　・陽性/陰性 → 12人目②へ

[12人目]
① A:0, B:1, C:1, D:86 ⇒ Bから1枚 測定
　・陽性 → 結論②へ
　・陰性 → 13人目①へ

② A:0, B:0, C:4, D:84 ⇒ Cから2枚 測定
　・陽性/陰性 → 13人目①へ

③ A:0, B:0, C:3, D:85 ⇒ Cから2枚 測定
　・陽性 → 13人目①へ
　・陰性 → 結論③へ

[13人目]
① A:0, B:0, C:2, D:86 ⇒ Cから1枚 測定
　・陽性/陰性 → 結論③へ

[結論]
① A:1, B:0, C:0, D:87 ⇒ Aの1枚が偽金貨確定

② A:0, B:1, C:0, D:87 ⇒ Bの1枚が偽金貨確定

③ A:0, B:0, C:1, D:87 ⇒ Cの1枚が偽金貨確定

※12人目までに偽金貨が確定した場合、残りの技術者に偽金貨だけを測定させれば嘘つきを特定することもできる。


各測定の枚数は条件を満たすものをとりあえずで決めただけなので、ほかにも測定方法はいろいろありそうです。

りらひいさん！

これ、大変に興味深いです。！！！
ありがとうございます。
勉強させてください。
今日は病院待合が立て込んでおりまして
ボチボチとトレースしたいとぞんじます。


皆様への追伸:
皆様には当然のことですけれども一応。

■上界と上限の違い
上界 (Upper Bound): 集合 A の上界とは、集合のすべての要素がその数以下であるような数のことです。上界は複数存在する可能性があります。
上限 (Supremum): 上限は、上界の中で最小のものを指します。つまり、上限より小さい数は上界ではなくなります。上限は一意に定まります。

例として、集合 (0,1) の場合、1が上限ですが、2や3も上界です。

さて。
[2286]の拙投稿にて
バイナリコードにおける、符号長 11 、最小ハミング距離が 5 であって、符号語の数が 24 であるものの具体的な構成例をお示しいたしました。

昨夜にみつけたのですけれども、
バイナリコードにおける、符号長 11 、最小ハミング距離が 5 という条件のもとで、
符号語の数のひとつの上界として 24 が既に証明されているようです。
→　
Table of general binary codes
https://aeb.win.tue.nl/codes/binary-1.html

このページに記載の上界の一覧表では
符号長 n . 最小ハミング距離 d のもとで、符号語数 A₂(n,d) がまとめられています。ただし、d が偶数のときのみです。
d が奇数の時には以下の公式を使います。
A₂(n-1,2e-1) = A₂(n,2e)

そこで、n=12, e=3 とすると、
A₂(11,5) = A₂(12,6)
となります。
A₂(11,5)の値を知りたかったので一覧表からA₂(12,6)を引きます。
A₂(12,6) = 24
が得られました。すなわち
A₂(11,5) = 24
となります。
バイナリコードにおける、符号長 11 、最小ハミング距離が 5 という条件のもとで、
符号語の数のひとつの上界として 24 が判明したということとなります。

あらためまして先日みつけた符号を下記に再掲いたします。

00000000000,
10111000100,
01011100010,
00101110001,
10010111000,
01001011100,
00100101110,
00010010111,
10001001011,
11000100101,
11100010010,
01110001001,
01000111011,
10100011101,
11010001110,
01101000111,
10110100011,
11011010001,
11101101000,
01110110100,
00111011010,
00011101101,
10001110110,
11111111111,

24 という上界は上限でもあったと判明したこととなります(個人的に嬉しくてイキっております、ご寛恕を願います)

以下の問題に解があるとのお知らせをもらいました。このスレッドの冒頭に比べて金貨が51枚も増えています。

64 枚の金貨があり、そのうち1枚は放射性物質を含む偽物です。この偽物をガイガーカウンターを用いた一連の測定で特定する課題です。問題は以下のように構成されています：

金貨：64 枚の金貨があり、1 から 64 まで番号が振られています。ちょうど 1 枚が偽物です。

偽物の特性：偽の金貨は放射性物質を含んでおり、ガイガーカウンターで検出可能です。

技術者：13人の技術者が測定を行います。

測定プロセス：各技術者は 64 枚の金貨のうち一部を選んで測定します。技術者は選ばれた金貨をガイガーカウンターで同時に検査します。ガイガーカウンターは、選ばれた金貨の中に偽物が含まれている場合にのみ反応します。測定結果は、その技術者のみが知っています。

測定の制約：各技術者は正確に 1 回だけ測定を行い、合計 13 回の測定が行われます。

報告：各測定の後、技術者は「陽性」（放射線が検出された）または「陰性」（放射線が検出されなかった）と報告します。

信頼性の問題：技術者のうち最大 2 人までが、意図的な欺瞞や偶発的な誤りにより不正確な報告をする可能性があります。

目標：最大 2 件の不正確な報告があったとしても、偽の金貨を確実に特定することが目標です。

- - -

割とまじめに計算機を使いましたがせいぜいで 40 枚くらいまでで………
64 枚解を求めるためには本格的なバックトラッキングをしなくては解けない気がしてきまして慄いています。

[2299]の答です。
符号語数が64で符号長が13、最小ハミング距離が5となっています。

0000000000000
0000000011111
0000011100011
0000101101100
0000110110101
0000111011010
0001001110110
0001110001001
0010010101110
0010101010001
0011001001011
0011010010010
0011011111101
0011100100111
0011100111000
0011111000100
0100010111000
0100101000111
0101000101101
0101011001110
0101011010001
0101100010100
0101101111011
0101110100010
0110000110011
0110001011100
0110010000101
0110100001010
0110111101001
0110111110110
0111001100000
0111110011111
1000011010100
1000100101011
1001000110001
1001010000111
1001011101000
1001101000010
1001101011101
1001110111110
1010001100101
1010001111010
1010010011001
1010100010110
1010110100000
1010111001111
1011000001100
1011111110011
1100000100110
1100001001001
1100011111111
1100101110000
1100110001100
1100110010011
1101000011010
1101111100101
1110011000010
1110100111101
1111001010111
1111010101011
1111010110100
1111100000001
1111101101110
1111111011000

どうやっても私には独力では作れませんでした。ガックシ…… 符号語数が32までならスグに作れるのに。

ご参考。

OEIS の A005865
https://oeis.org/A005865
にてわかることとして
A005865(13, 5) = A005865(14, 6) = 64
となります。
嘘報告が最大2(すなわち5/2の整数部分)あっても13回測れば
64枚から1枚の偽金貨を特定できるということなので
その具体的な方法は？
というストーリーなのでした。