令和7年3月30日

[流れ星]

  第452数学的な連続応募解答

   <解答募集期間:3月2日〜3月30日>

単位分数]

 分子が1の分数を単位分数といい、分数を異なる単位分数の和で表すことを単位分数分解という。

実は、単位分数分解の方法は無限通りあり、証明することも可能です。

後述にある参考を参照ください。

 

単位分数分解の方法は、2通りあります。

1 分解したい分数から、最大の単位分数をひいていく方法

2 不等式によって、分母にあてはまる自然数をしぼっていく方法

 

単位分数分解の方法1

@ もとの分数から、それより小さい最大の単位分数をひく。

A その結果が単位分数なら終了。単位分数でないなら、その分数より小さい最大の単位分数をひく。

B Aを繰り返す。

* 有限回の計算で必ず単位分数が登場する。

例    の単位分数分解

@   より小さい最大の単位分数は、  なので、

        ∴ 

A  は単位分数のため、分解終了。

 

単位分数分解の方法2

@ もとの分数     (x<yとおく。

 

A x にあてはまる可能性のある自然数を求める。

 

B  x をそれぞれ代入し、そのときのy の値が自然数(が単位分数)なら終了し、そうでなければその分数について(1)から同じ作業を行っていく。

 

 x にあてはまる可能性のある自然数は、n より大きく、2n より小さい分数となります。

 

ここで、問題です

問題  のように3個の単位分数に分解できます。

 

1を4個の単位分数に分解してください。

 

余力問題 1を5個の異なる単位分数に分解すると。何組になりますか。

紙の上で、解を見つけるには、相当な時間がかかり、途中で諦めました。 

誰か、PCで求めて頂けませんか。

      

参考文献 BLUEBACKS 大学受験で語る数論の世界 清水健一 講談社

 

参考に

単位分数の和

単位分数分解の帰納法

追加問題(出題者は「ジョーカー」) 新作シリーズ11,12問目

 今回で終了です。

問題1

 

問題2

 

NO1「三角定規」 03/02        2142分  受信 更新 3/30

[問題452] (余力問題) 11/a+1/b+1/c+1/d+1/e (a<b<c<d<e)につきまして

参考資料もキチンと読んでおらず, (問題1 も 余力問題も) まだキチンとは考えていないのですが,

> 誰か、 PC で求めて頂けませんか。

とあったので,エクセルでマクロを組み,解の例のみ求めてみました。

添付 エクセルファイル「1を5個の単位分数の和で.xlsmをご覧下さい。

e100 42個の,e200 57個の 解の例があるようです。

「三角定規」 03/15        0937分  受信 更新 3/30

「三角定規」 03/16        0918分  受信 更新 3/30

 

前回から修正された解答

 

NO2「ジョーカー」     03/04     1613       受信  更新 3/30

「ジョーカー」     03/07     2131       受信  更新 3/30

前回から予力問題の解法が追加され解答

 

NO3「日曜数学者」     03/05       1232分      受信  更新 3/30

[解答]

 

     (1/x)+(1/y)+(1/z)+(1/w)=1 --------- (1)

 

を満たす4個の自然数x,y,z,w (0<x<y<z<w)を求める。

 

(1)および0<x<y<z<wより、

     (1/x)<(1/x)+(1/y)+(1/z)+(1/w)=1

よって、

     x>1 -------- (2)

である。

 

また、仮に、x>=3ならば、

     (1/x)+(1/y)+(1/z)+(1/w)<=1/3+1/4+1/5+1/6=19/20<1

より、(1)は成立しない。

よって、x<3 ------- (3)

自然数x(2),(3)を満たすので、x=2である。

 

(1)x=2を代入して、

     (1/y)+(1/z)+(1/w)=1/2 ---------- (4)

2<y<z<wより、仮に、y>=6ならば

     (1/y)+(1/z)+(1/w)<=1/6+1/7+1/8=73/168<1/2

となるので、(4)は成立しない。

よって、2<y<6である。

 

(i)y=3のとき

(4)y=3を代入して、

     (1/z)+(1/w)=1/6 -------- (5)

     6*(z+w)=z*w

     (z-6)*(w-6)=36

(5)より

     z>6

(z-6),(w-6)は、36の約数で、0<x-6<w-6なので、

     (z-6,w-6)=(1,36),(2,18),(3,12),(4,9)

     (z,w)=(7,42),(8,24),(9,18),(10,15)

よって、

     (x,y,z,w)=(2,3,7,42),(2,3,8,24),(2,3,9,18),(2,3,10,15)

 

(ii)y=4のとき

 

(4)y=4を代入して、

     (1/z)+(1/w)=1/4 -------- (6)

     4*(z+w)=z*w

     (z-4)*(w-4)=16

(6)より

     z>4

(z-4),(w-4)は、16の約数で、0<z-4<w-4なので、

     (z-4,w-4)=(1,16),(2,8)

     (z,w)=(5,20),(6,12)

よって、

     (x,y,z,w)=(2,4,5,20),(2,4,6,12)

 

(iii)y=5のとき

 

(4)y=5を代入して、

     (1/z)+(1/w)=3/10 -------- (7)

     10*(z+w)=3*z*w

     (3*z-10)*(3*w-10)=100

5<z<wより、整数(3*z-10),(3*w-10)は、100の約数で、0<3*z-10<3*w-10なので、

    (3*z-10,3*w-10)=(1,100),(2,50),(4,25),(5,20)

    (z,w)=(11/3,110/3),(4,20),(14/3,35/3),(5,10)

この中で、z,wが自然数で、z>5となるものは存在しない。

 

(i),(ii),(iii)より、

 

     (x,y,z,w)=(2,3,7,42),(2,3,8,24),(2,3,9,18),(2,3,10,15),(2,4,5,20),(2,4,6,12)

 

 

(1)の自然数解(x,y,z,w)0<x<y<z<wを満たすものは、6個である。

 

     1/2+1/3+1/7+1/42=1

     1/2+1/3+1/8+1/24=1

     1/2+1/3+1/9+1/18=1

     1/2+1/3+1/10+1/15=1

     1/2+1/4+1/5+1/20=1

     1/2+1/4+1/6+1/12=1

 

 

[参考]

 

1n個の分母の異なる単位分数の和で表す方法の数をPari/GPのプログラムを作成して、

計算させてみました。

結果は、以下の通りです。

 

自然数n>=2に対して、

 

S(n)={[x_1,x_2,x_3,...,x_n] :(1/x_1)+(1/x_2)+(1/x_3)+...+(1/x_n)=1,

      x_1,x_2,x_3,...,x_nは自然数, 0<x_1<x_2<x_3<...<x_n }

とする。

 

 

#S(2)=0

#S(3)=1

#S(4)=6

#S(5)=72

#S(6)=2320     (添付ファイル参照)

#S(7)=249264      (添付ファイル参照)

#S(8)=???  (157万以上)

 

 

[Pari/GPによる計算結果]

 

                         GP/PARI CALCULATOR Version 2.15.0 (released)

                 amd64 running netbsd (x86-64/GMP-6.1.2 kernel) 64-bit version

                    compiled: Oct 15 2022, gcc version 5.5.0 (nb2 20180327)

                                   threading engine: single

                        (readline v7.0 enabled, extended help enabled)

 

                            Copyright (C) 2000-2022 The PARI Group

 

PARI/GP is free software, covered by the GNU General Public License, and comes WITHOUT ANY WARRANTY WHATSOEVER.

 

Type ? for help, \q to quit.

Type ?18 for how to get moral (and possibly technical) support.

 

parisize = 6144000000, primelimit = 100000000

(08:03) gp > sss2a(1)

%1 = 0

(08:03) gp > sss3a(1)

1:[2, 3, 6]

%2 = 1

(08:04) gp > sss4a(1)

1:[2, 3, 7, 42]

2:[2, 3, 8, 24]

3:[2, 3, 9, 18]

4:[2, 3, 10, 15]

5:[2, 4, 5, 20]

6:[2, 4, 6, 12]

%3 = 6

(08:04) gp > sss5a(1)

1:[2, 3, 7, 43, 1806]

2:[2, 3, 7, 44, 924]

3:[2, 3, 7, 45, 630]

4:[2, 3, 7, 46, 483]

5:[2, 3, 7, 48, 336]

6:[2, 3, 7, 49, 294]

7:[2, 3, 7, 51, 238]

8:[2, 3, 7, 54, 189]

9:[2, 3, 7, 56, 168]

10:[2, 3, 7, 60, 140]

11:[2, 3, 7, 63, 126]

12:[2, 3, 7, 70, 105]

13:[2, 3, 7, 78, 91]

14:[2, 3, 8, 25, 600]

15:[2, 3, 8, 26, 312]

16:[2, 3, 8, 27, 216]

17:[2, 3, 8, 28, 168]

18:[2, 3, 8, 30, 120]

19:[2, 3, 8, 32, 96]

20:[2, 3, 8, 33, 88]

21:[2, 3, 8, 36, 72]

22:[2, 3, 8, 40, 60]

23:[2, 3, 8, 42, 56]

24:[2, 3, 9, 19, 342]

25:[2, 3, 9, 20, 180]

26:[2, 3, 9, 21, 126]

27:[2, 3, 9, 22, 99]

28:[2, 3, 9, 24, 72]

29:[2, 3, 9, 27, 54]

30:[2, 3, 9, 30, 45]

31:[2, 3, 10, 16, 240]

32:[2, 3, 10, 18, 90]

33:[2, 3, 10, 20, 60]

34:[2, 3, 10, 24, 40]

35:[2, 3, 11, 14, 231]

36:[2, 3, 11, 15, 110]

37:[2, 3, 11, 22, 33]

38:[2, 3, 12, 13, 156]

39:[2, 3, 12, 14, 84]

40:[2, 3, 12, 15, 60]

41:[2, 3, 12, 16, 48]

42:[2, 3, 12, 18, 36]

43:[2, 3, 12, 20, 30]

44:[2, 3, 12, 21, 28]

45:[2, 3, 14, 15, 35]

46:[2, 4, 5, 21, 420]

47:[2, 4, 5, 22, 220]

48:[2, 4, 5, 24, 120]

49:[2, 4, 5, 25, 100]

50:[2, 4, 5, 28, 70]

51:[2, 4, 5, 30, 60]

52:[2, 4, 5, 36, 45]

53:[2, 4, 6, 13, 156]

54:[2, 4, 6, 14, 84]

55:[2, 4, 6, 15, 60]

56:[2, 4, 6, 16, 48]

57:[2, 4, 6, 18, 36]

58:[2, 4, 6, 20, 30]

59:[2, 4, 6, 21, 28]

60:[2, 4, 7, 10, 140]

61:[2, 4, 7, 12, 42]

62:[2, 4, 7, 14, 28]

63:[2, 4, 8, 9, 72]

64:[2, 4, 8, 10, 40]

65:[2, 4, 8, 12, 24]

66:[2, 4, 9, 12, 18]

67:[2, 4, 10, 12, 15]

68:[2, 5, 6, 8, 120]

69:[2, 5, 6, 9, 45]

70:[2, 5, 6, 10, 30]

71:[2, 5, 6, 12, 20]

72:[3, 4, 5, 6, 20]

time = 1 ms.

%4 = 72

(08:04) gp > sss6c(1)

time = 23 ms.

%5 = 2320

(08:04) gp > sss7c(1)

time = 8,337 ms.

%6 = 249264

 

 

<水の流れ:[Pari/GPのプログラム(n=4,5の場合)]は割愛させていただきます。>

<水の流れ:この解答は異なる単位分数の和を考えています。>

「日曜数学者」     03/06       1423分      受信  更新 3/30

[参考2]

 

単位分数(分母の重複を許す)の和で表す場合についても、計算しました。

 

14個の単位分数(分母の重複を許す)の和で表す方法は、以下の14通りである。

 

     1/2+1/3+1/7+1/42=1

     1/2+1/3+1/8+1/24=1

     1/2+1/3+1/9+1/18=1

     1/2+1/3+1/10+1/15=1

     1/2+1/3+1/12+1/12=1

     1/2+1/4+1/5+1/20=1

     1/2+1/4+1/6+1/12=1

     1/2+1/4+1/8+1/8=1

     1/2+1/5+1/10+1/10=1

     1/2+1/6+1/6+1/6=1

     1/3+1/3+1/4+1/12=1

     1/3+1/3+1/6+1/6=1

     1/3+1/4+1/4+1/6=1

     1/4+1/4+1/4+1/4=1

 

自然数n>=2に対して、

 

T(n)={[x_1,x_2,x_3,...,x_n] :(1/x_1)+(1/x_2)+(1/x_3)+...+(1/x_n)=1,

      x_1,x_2,x_3,...,x_nは自然数, 0<x_1<=x_2<=x_3<=...<=x_n }

とする。

 

#T(2)=1

#T(3)=3

#T(4)=14

#T(5)=147

#T(6)=3483      (添付ファイル参照)

#T(7)=294314      (添付ファイル参照)

#T(8)=???

 

<水の流れ:6個、7個でできる単位分数の和が1となる場合が4桁、6桁となるとは驚きです・。

 

[Pari/GPによる計算結果]

(09:24) gp > sss2(1)

1:[2, 2]

%2 = 1

 

(09:24) gp > sss3(1)

1:[2, 3, 6]

2:[2, 4, 4]

3:[3, 3, 3]

%3 = 3

 

(09:31) gp > sss4(1,1)

1:[2, 3, 7, 42]

2:[2, 3, 8, 24]

3:[2, 3, 9, 18]

4:[2, 3, 10, 15]

5:[2, 3, 12, 12]

6:[2, 4, 5, 20]

7:[2, 4, 6, 12]

8:[2, 4, 8, 8]

9:[2, 5, 5, 10]

10:[2, 6, 6, 6]

11:[3, 3, 4, 12]

12:[3, 3, 6, 6]

13:[3, 4, 4, 6]

14:[4, 4, 4, 4]

%2 = 14

 

(09:31) gp > sss5(1,1)

1:[2, 3, 7, 43, 1806]

2:[2, 3, 7, 44, 924]

3:[2, 3, 7, 45, 630]

4:[2, 3, 7, 46, 483]

5:[2, 3, 7, 48, 336]

6:[2, 3, 7, 49, 294]

7:[2, 3, 7, 51, 238]

8:[2, 3, 7, 54, 189]

9:[2, 3, 7, 56, 168]

10:[2, 3, 7, 60, 140]

11:[2, 3, 7, 63, 126]

12:[2, 3, 7, 70, 105]

13:[2, 3, 7, 78, 91]

14:[2, 3, 7, 84, 84]

15:[2, 3, 8, 25, 600]

16:[2, 3, 8, 26, 312]

17:[2, 3, 8, 27, 216]

18:[2, 3, 8, 28, 168]

19:[2, 3, 8, 30, 120]

20:[2, 3, 8, 32, 96]

21:[2, 3, 8, 33, 88]

22:[2, 3, 8, 36, 72]

23:[2, 3, 8, 40, 60]

24:[2, 3, 8, 42, 56]

25:[2, 3, 8, 48, 48]

26:[2, 3, 9, 19, 342]

27:[2, 3, 9, 20, 180]

28:[2, 3, 9, 21, 126]

29:[2, 3, 9, 22, 99]

30:[2, 3, 9, 24, 72]

31:[2, 3, 9, 27, 54]

32:[2, 3, 9, 30, 45]

33:[2, 3, 9, 36, 36]

34:[2, 3, 10, 16, 240]

35:[2, 3, 10, 18, 90]

36:[2, 3, 10, 20, 60]

37:[2, 3, 10, 24, 40]

38:[2, 3, 10, 30, 30]

39:[2, 3, 11, 14, 231]

40:[2, 3, 11, 15, 110]

41:[2, 3, 11, 22, 33]

42:[2, 3, 12, 13, 156]

43:[2, 3, 12, 14, 84]

44:[2, 3, 12, 15, 60]

45:[2, 3, 12, 16, 48]

46:[2, 3, 12, 18, 36]

47:[2, 3, 12, 20, 30]

48:[2, 3, 12, 21, 28]

49:[2, 3, 12, 24, 24]

50:[2, 3, 13, 13, 78]

51:[2, 3, 14, 14, 42]

52:[2, 3, 14, 15, 35]

53:[2, 3, 14, 21, 21]

54:[2, 3, 15, 15, 30]

55:[2, 3, 15, 20, 20]

56:[2, 3, 16, 16, 24]

57:[2, 3, 18, 18, 18]

58:[2, 4, 5, 21, 420]

59:[2, 4, 5, 22, 220]

60:[2, 4, 5, 24, 120]

61:[2, 4, 5, 25, 100]

62:[2, 4, 5, 28, 70]

63:[2, 4, 5, 30, 60]

64:[2, 4, 5, 36, 45]

65:[2, 4, 5, 40, 40]

66:[2, 4, 6, 13, 156]

67:[2, 4, 6, 14, 84]

68:[2, 4, 6, 15, 60]

69:[2, 4, 6, 16, 48]

70:[2, 4, 6, 18, 36]

71:[2, 4, 6, 20, 30]

72:[2, 4, 6, 21, 28]

73:[2, 4, 6, 24, 24]

74:[2, 4, 7, 10, 140]

75:[2, 4, 7, 12, 42]

76:[2, 4, 7, 14, 28]

77:[2, 4, 8, 9, 72]

78:[2, 4, 8, 10, 40]

79:[2, 4, 8, 12, 24]

80:[2, 4, 8, 16, 16]

81:[2, 4, 9, 9, 36]

82:[2, 4, 9, 12, 18]

83:[2, 4, 10, 10, 20]

84:[2, 4, 10, 12, 15]

85:[2, 4, 12, 12, 12]

86:[2, 5, 5, 11, 110]

87:[2, 5, 5, 12, 60]

88:[2, 5, 5, 14, 35]

89:[2, 5, 5, 15, 30]

90:[2, 5, 5, 20, 20]

91:[2, 5, 6, 8, 120]

92:[2, 5, 6, 9, 45]

93:[2, 5, 6, 10, 30]

94:[2, 5, 6, 12, 20]

95:[2, 5, 6, 15, 15]

96:[2, 5, 7, 7, 70]

97:[2, 5, 8, 8, 20]

98:[2, 5, 10, 10, 10]

99:[2, 6, 6, 7, 42]

100:[2, 6, 6, 8, 24]

101:[2, 6, 6, 9, 18]

102:[2, 6, 6, 10, 15]

103:[2, 6, 6, 12, 12]

104:[2, 6, 7, 7, 21]

105:[2, 6, 8, 8, 12]

106:[2, 6, 9, 9, 9]

107:[2, 7, 7, 7, 14]

108:[2, 8, 8, 8, 8]

109:[3, 3, 4, 13, 156]

110:[3, 3, 4, 14, 84]

111:[3, 3, 4, 15, 60]

112:[3, 3, 4, 16, 48]

113:[3, 3, 4, 18, 36]

114:[3, 3, 4, 20, 30]

115:[3, 3, 4, 21, 28]

116:[3, 3, 4, 24, 24]

117:[3, 3, 5, 8, 120]

118:[3, 3, 5, 9, 45]

119:[3, 3, 5, 10, 30]

120:[3, 3, 5, 12, 20]

121:[3, 3, 5, 15, 15]

122:[3, 3, 6, 7, 42]

123:[3, 3, 6, 8, 24]

124:[3, 3, 6, 9, 18]

125:[3, 3, 6, 10, 15]

126:[3, 3, 6, 12, 12]

127:[3, 3, 7, 7, 21]

128:[3, 3, 8, 8, 12]

129:[3, 3, 9, 9, 9]

130:[3, 4, 4, 7, 42]

131:[3, 4, 4, 8, 24]

132:[3, 4, 4, 9, 18]

133:[3, 4, 4, 10, 15]

134:[3, 4, 4, 12, 12]

135:[3, 4, 5, 5, 60]

136:[3, 4, 5, 6, 20]

137:[3, 4, 6, 6, 12]

138:[3, 4, 6, 8, 8]

139:[3, 5, 5, 5, 15]

140:[3, 5, 5, 6, 10]

141:[3, 6, 6, 6, 6]

142:[4, 4, 4, 5, 20]

143:[4, 4, 4, 6, 12]

144:[4, 4, 4, 8, 8]

145:[4, 4, 5, 5, 10]

146:[4, 4, 6, 6, 6]

147:[5, 5, 5, 5, 5]

time = 1 ms.

%5 = 147

 

(09:31) gp > sss5(1,0)

time = 1 ms.

%6 = 147

 

(09:31) gp > sss6(1,0)

time = 18 ms.

%7 = 3462

 

(09:32) gp > sss7(1,0)

time = 9,275 ms.

%8 = 294314

 

<水の流れ:[Pari/GPのプログラム(n=4,5の場合)]は割愛させていただきます。>

NO4「スモークマン」 03/06 2003分  受信 更新 3/30

問題 単位分数3 or 4個の和が1

1/a+1/b+1/c+1/d=1

a<b<c<d

4/a>1, 4/d<

Soa3以下、d5以上

(1)     a=3

2/3=1/b+1/c+1/d

3/b>2/39/2>b・・・b4以下、

b=4

2/3-1/4=5/12=1/c+1/d

5/12<2/c24/5>cc4以下・・・b=c=4で適さない

(2)     a=2

1/2=1/b+1/c+1/d

3/b>1/26>b・・・b5, 4, 3

(21) b=5

1/2-1/5=3/10<2/cc<20/3・・・c6以下

3/10-1/6=(18-10)/60・・・X

(22) b=4

1/2-1/4=1/4<2/c・・・c7以下

1/4-1/7=(7-4)/28・・・X

1/4-1/6=(6-4)/24=1/12

1/4-1/5=1/20

(23) b=3

1/2-1/3=1/6<2/c・・・c12未満

1/6-1/11=(11-6)/66・・・X

1/6-1/10=(10-6)/60=1/15

1/6-1/9=(9-6)/54=1/18

1/6-1/8=(8-6)/48=1/24

1/6-1/7=1/42

結局…

1/2+1/3+1/10+1/15

1/2+1/3+1/9+1/18

1/2+1/3+1/8+1/24

1/2+1/3+1/7+1/42

1/2+1/4+1/6+1/12

1/2+1/4+1/5+1/20

 

同様にして

1/a+1/b+1/c+1/d+1/e=1

5/a>1・・・a4以下

(1)     a=4

3/4<4/bb<16/3・・・b5以下

3/4-1/5=(15-4)/20=11/20<3/cc/60/11・・・c5以下でX

(2)     a=3

2/3<4/bb<6・・・b5以下

(21)b=5

2/3-1/5=(10-3)/15=7/15<3/cc<45/7・・・c6以下

7/15-1/6=(42-15)/90=27/90=3/10<2/dd<20/3・・・d6以下でX

(22)b=4

2/3-1/4=(8-3)/12=5/12<3/cc<36/5・・・cは7以下

(31)c=5

5/12-1/5=(25-12)/60=13/60<2/dd<120/13・・・dは9以下

13/60-1/9=(117-60)/540=57/540・・・X

13/60-1/8=(104-60)/480=44/480・・・X

13/60-1/7=(91-60)/420=31/420・・・X

13/60-1/6=(78-60)/360=18/360=1/20

1/3+1/4+1/5+1/6+1/20

(32)c=6

5/12-1/6=(30-12)/72=18/72=1/4<2/dd<8・・・d7以下

1/4-1/7=(7-3)/28・・・X

(33)c=7

5/12-1/7=(35-12)/84=23/94<2/dd<188/23・・・d8以下

23/94-1/8=(184-94)/(8*94)=90/(8*94)・・・X

(3)     a=2

1/2<4/b・・・b7以下

(31)b=3

1/2-1/3=1/6<3/c・・・c17以下

1/6-1/17=11/102<2/d・・・d<204/11・・・d18以下

11/102・・・なし

1/6-1/16=10/(6*16)=5/48<2/d・・・d<96/5・・・d19以下・・・適するものはなし

1/6-1/15=1/10<2/d・・・d19以下

ここもなし

1/6-1/14=8/84=2/21<2/d・・・d20以下

2/21-1/15=1/35

1/2+1/3+1/14+1/15+1/35

1/6-1/13=7/78<2/d・・・d<156/7・・・d22以下

ここはなし

1/6-1/12=1/12<2/d・・・d23以下

1/12-1/21=1/28

1/2+1/3+1/12+1/21+1/28

1/12-1/20=1/30

1/2+1/3+1/12+1/20+1/30

1/12-1/18=1/36

1/2+1/3+1/12+1/18+1/36

1/12-1/16=1/48

1/2+1/3+1/12+1/16+1/48

1/12-1/15=1/60

1/2+1/3+1/12+1/15+1/60

1/12-1/14=1/84

1/2+1/3+1/12+1/14+1/84

1/12-1/13=1/156

1/2+1/3+1/12+1/13+1/156

1/6-1/11=5/66<2/dd<132/5・・・d26以下

5/66-1/22=1/33

1/2+1/3+1/11+1/22+1/33

5/66-1/15=1/110

1/2+1/3+1/11+1/15+1/110

5/66-1/14=1/231

1/2+1/3;1/11+1/14+1/231

1/6-1/10=1/15<2/d・・・d29以下

1/15-1/24=1/40

1/2+1/3+1/10+1/24+1/40

1/15-1/20=1/60

1/2+1/3+1/10+1/20+1/60

1/15-1/18=1/90

1/2+1/3+1/10+1/18+1/90

1/15-1/16=1/240

1/2+1/3+1/10+1/16+1/240

1/6-1/9=1/18<2/d・・・d35以下

1/18=1/30=1/45

1/2+1/3+1/9+1/30+1/45

1/18-1/27=1/54

1/2+1/3+1/9+1/27+1/54

1/18-1/24=1/72

1/2+1/3+1/9+1/24+1/72

1/18-1/22=1/99

1/2+1/3+1/9+1/22+1/99

1/18-1/21=1/126

1/2+1/3+1/9+1/21+1/126

1/18-1/20=1/180

1/2+1/3+1/9+1/20+1/180

1/18-1/19=1/342

1/2+1/3+1/9+1/19+1/342

1/6-1/8=1/24<2/d・・・d47以下

1/24-1/42=1/56

1/2+1/3+1/8+1/42+1/56

1/24-1/40=1/60

1/2+1/3+1/8+1/40+1/60

1/24-1/36=1/72

1/2+1/3+1/8+1/36+1/72

1/24-1/33=1/88

1/2+1/3+1/8+1/33+1/88

1/24-1/32=1/96

1/2+1/3+1/8+1/32+1/96

1/24-1/30=1/120

1/2+1/3+1/8+1/30+1/120

1/24-1/28=1/168

1/2+1/3+1/8+1/28+1/168

1/24-1/27=1/216

1/2+1/3+1/8+1/27+1/216

1/24-1/26=1/312

1/2+1/3+1/8+1/26+1/312

1/24-1/25=1/600

1/2+1/3+1/8+1/25+1/600

1/6-1/7=1/42<2/d・・・d83以下

1/42-1/78=1/91

1/2+1/3+1/7+1/78+1/91

1/42-1/70=1/105

1/2+1/3+1/7+1/70+1/105

1/42-1/63=1/126

1/2+1/3+1/7+1/63+1/126

1/42-1/60=1/140

1/2+1/3+1/7+1/60+1/140

1/42-1/56=1/168

1/2+1/3+1/7+1/56+1/168

1/42-1/54=1/189

1/2+1/3+1/7+1/54+1/189

1/42-1/51=1/238

1/2+1/3+1/7+1/51+1/238

1/42-1/49=1/294

1/2+1/3+1/7+1/49+1/294

1/42-1/48=1/336

1/2+1/3+1/7+1/48+1/336

1/42-1/46=1/483

1/2+1/3+1/7+1/46+1/483

1/42-1/45=1/630

1/2+1/3+1/7+1/45+1/630

1/42-1/44=1/924

1/2+1/3+1/7+1/44+1/924

1/42-1/43=1/1806

1/2+1/3+1/7+1/43+1/1806

 

(32)b=4

1/2-1/4=1/4<3/c・・・c11以下

1/4-1/11=7/44<2/d・・・d<88/7・・・d12以下

7/44-1/12=5/66・・・X

1/4-1/10=6/40=3/20<2/d・・・d<40/3・・・d13以下

3/20-1/12=1/15

1/2+1/4+1/10+1/12+1/15

1/4-1/9=5/36<2/d・・・d14以下

5/36-1/12=1/18

1/2+1/4+1/9+1/12+1/18

1/4-1/8=1/8<2/d・・・d15以下

1/8-1/12=1/24

1/2+1/4+1/8+1/12+1/24

1/8-1/10=1/40

1/2+1/4+1/8+1/10+1/40

1/8-1/9=1/72

1/2+1/4+1/8+1/9+1/72

1/4-1/7=3/28<2/d・・・d18以下

3/28-1/14=1/28

1/2+1/4+1/7+1/14+1/28

3/28-1/12=1/42

1/2+1/4+1/7+1/12+1/42

3/28-1/10=1/140

1/2+1/4+1/7+1/10+1/140

1/4-1/6=1/12<2/d・・・d23以下

1/12-1/21=1/28

1/2+1/4+1/6+1/21+1/28

1/12-1/20=1/30

1/2+1/4+1/6+1/20+1/30

1/12-1/18=1/36

1/2+1/4+1/6+1/18+1/36

1/12-1/16=1/48

1/2+1/4+1/6+1/16+1/48

1/12-1/15=1/60

1/2+1/4+1/6+1/15+1/60

1/12-1/14=1/84

1/2+1/4+1/6+1/14+1/84

1/12-1/13=1/156

1/2+1/4+1/6+1/13+1/156

1/4-1/5=1/20<2/d・・・d39以下

1/20-1/36=1/45

1/2+1/4+1/5+1/36+1/45

1/20-1/30=1/60

1/2+1/4+1/5+1/30+1/60

1/20-1/28=1/70

1/2+1/4+1/5+1/28+1/70

1/20-1/25=1/100

1/2+1/4+1/5+1/25+1/100

1/20-1/24=1/120

1/2+1/4+1/5+1/24+1/120

1/20-1/22=1/220

1/2+1/4+1/5+1/22+1/220

1/20-1/21=1/420

1/2+1/4+1/5+1/21+1/420

(33)b=5

1/2-1/5=3/10<3/c・・・cは9以下

3/10-1/9=17/90<2/dd<180/17・・・d10以下

3/10-1/8=7/40<2/dd<80/7・・・d11以下

3/10-1/7=11/70<2/dd<140/11・・・d12以下

3/10-1/6=4/30=2/15<2/d・・・d14以下

2/15-1/12=1/20

1/2+1/5+1/6+1/12+1/20

2/15-1/10=1/30

1/2+1/5+1/6+1/10+1/30

2/15-1/9=1/45

1/2+1/5+1/6+1/9+1/45

2/15-1/8=1/120

1/2+1/5+1/6+1/8+1/120

(34)b=6

1/2-1/6=1/3<3/c・・・c8以下

1/3-1/8=5/24<2/d・・・dは9以下

適するものなし

1/3-1/7=4/21<2/d・・・d10以下

適するものなし

(35)

1/2-1/7=5/14<3/cc<42/5・・・c8以下

5/14-1/8=(40-14)/(8*14)=13/56<2/dd<112/13・・・d8以下・・・なし

 

PCでも

1/2+1/3+タイプ・・・45

1/2+1/4+タイプ・・・22

1/2+1/5+タイプ4

1/3+1/4+1/5+1/6+1/20・・・1

全部で72

So…上のものですべてのはず

 

「スモークマン」 03/07   2121分  受信 更新 3/30

 

 

452sumoku2

NO5「二度漬け白菜」    03/16       1959分      受信  更新 3/30

452回の応募問題の解答:

[単位分数]

14個の単位分数に分解する方法は次の14通り.()

1=(1/2)+(1/3)+(1/7)+(1/42),

=(1/2)+(1/3)+(1/8)+(1/24),

=(1/2)+(1/3)+(1/9)+(1/18),

=(1/2)+(1/3)+(1/10)+(1/15),

=(1/2)+(1/3)+(1/12)+(1/12),

=(1/2)+(1/4)+(1/5)+(1/20),

=(1/2)+(1/4)+(1/6)+(1/12),

=(1/2)+(1/4)+(1/8)+(1/8),

=(1/2)+(1/5)+(1/5)+(1/10),

=(1/2)+(1/6)+(1/6)+(1/6),

=(1/3)+(1/3)+(1/4)+(1/12),

=(1/3)+(1/3)+(1/6)+(1/6),

=(1/3)+(1/4)+(1/4)+(1/6),

=(1/4)+(1/4)+(1/4)+(1/4).

また,14個の単位分数に,単位分数分解する表す方法は次の6通り.

1=(1/2)+(1/3)+(1/7)+(1/42),

=(1/2)+(1/3)+(1/8)+(1/24),

=(1/2)+(1/3)+(1/9)+(1/18),

=(1/2)+(1/3)+(1/10)+(1/15),

=(1/2)+(1/4)+(1/5)+(1/20),

=(1/2)+(1/4)+(1/6)+(1/12).

 

1=(1/a)+(1/b)+(1/c)+(1/d), abcd

を満たすような正整数a,b,c,dの組(a,b,c,d)

すべて求めればよい.

1-(1/a)=(1/b)+(1/c)+(1/d)>0 より,2a.

(4/a)(1/a)+(1/b)+(1/c)+(1/d)=1 より,a4.

よって,2a4.

1-(1/a)-(1/b)=(1/c)+(1/d)>0 より,floor(1/(1-(1/a)))+1b.

(3/b)(1/b)+(1/c)+(1/d)=1-(1/a) より, bfloor(3/(1-(1/a))).

よって,max(a,floor(1/(1-(1/a)))+1)bfloor(3/(1-(1/a))) ---(1)

1-(1/a)-(1/b)-(1/c)=(1/d)>0 より,floor(1/(1-(1/a)-(1/b)))+1c.

(2/c)(1/c)+(1/d)=1-(1/a)-(1/b) より, cfloor(2/(1-(1/a)-(1/b))).

よって,max(b,floor(1/(1-(1/a)-(1/b)))+1)cfloor(2/(1-(1/a)-(1/b))) ---(2)

 

a=2の場合:

(1)より,3b6.

(a,b)=(2,3)のとき,(2)より,7c12.

このうち,(d)=1/(1-(1/a)-(1/b)-(1/c)) 7 以上の整数となるような(a,b,c)は,

(a,b,c)=(2,3,7),(2,3,8),(2,3,9),(2,3,10),(2,3,12).

(a,b)=(2,4)のとき,(2)より,5c8.

このうち,(d)=1/(1-(1/a)-(1/b)-(1/c)) 5 以上の整数となるような(a,b,c)は,

(a,b,c)=(2,4,5),(2,4,6),(2,4,8).

(a,b)=(2,5)のとき,(2)より,5c6.

このうち,(d)=1/(1-(1/a)-(1/b)-(1/c)) 5 以上の整数となるような(a,b,c)は,

(a,b,c)=(2,5,5)のみ.

(a,b)=(2,6)のとき,(2)より,c=6.

このとき,(d)=1/(1-(1/a)+(1/b)+(1/c))=6 となる.

 

a=3の場合:

(1)より,3b4.

(a,b)=(3,3)のとき,(2)より,4c6.

このうち,(d)=1/(1-(1/a)-(1/b)-(1/c)) 4 以上の整数となるような(a,b,c)は,

(a,b,c)=(3,3,4),(3,3,6).

(a,b)=(3,4)のとき,(2)より,c=4.

このとき,(d)=1/(1-(1/a)-(1/b)-(1/c))=6 となる.

 

a=4の場合:

(1)より,b=4.

(a,b)=(4,4)のとき,(2)より,c=4.

このとき,(d)=1/(1-(1/a)-(1/b)-(1/c))=4 となる.

 

以上より,次の14組が解となる.

(a,b,c,d)

=(2,3,7,42),(2,3,8,24),(2,3,9,18),(2,3,10,15),(2,3,12,12),

(2,4,5,20),(2,4,6,12),(2,4,8,8),(2,5,5,10),(2,6,6,6),

(3,3,4,12),(3,3,6,6),(3,4,4,6),(4,4,4,4).

 

[余力問題]

1=(1/a)+(1/b)+(1/c)+(1/d)+(1/e), a<b<c<d<e

を満たすような正整数a,b,c,d,eの組(a,b,c,d,e)の総数を

数えればよい.この組の総数は 72 ()

 

1-(1/a)=(1/b)+(1/c)+(1/d)+(1/e)>0 および

(5/a)>(1/a)+(1/b)+(1/c)+(1/d)+(1/e)=1 より,

2a4.

1-(1/a)-(1/b)=(1/c)+(1/d)+(1/e)>0 および    

(4/b)>(1/b)+(1/c)+(1/d)+(1/e)=1-(1/a) より,

max(a+1,floor(1/(1-1/a))+1)bceil(4/(1-1/a))-1.

1-(1/a)-(1/b)-(1/c)=(1/d)+(1/e)>0 および    

(3/c)>(1/c)+(1/d)+(1/e)=1-(1/a)-(1/b) より,

max(b+1,floor(1/(1-(1/a)-(1/b)))+1)cceil(3/(1-(1/a)-(1/b)))-1.

1-(1/a)-(1/b)-(1/c)-(1/d)=(1/e)>0 および    

(2/d)>(1/d)+(1/e)=1-(1/a)-(1/b)-(1/c) より,

max(c+1,floor(1/(1-(1/a)-(1/b)-(1/c)))+1)dceil(2/(1-(1/a)-(1/b)-(1/c)))-1.

ここで,

e=1/(1-(1/a)-(1/b)-(1/c)-(1/d))

としたとき,eが整数か否かをチェックすればよい.

以上をmaximaで計算させるプログラムが,次です.

 

count:0$

for a:2 thru 4 do (

for b:max(a+1,floor(1/(1-1/a))+1) thru ceiling(4/(1-1/a))-1 do ( for c:max(b+1,floor(1/(1-1/a-1/b))+1) thru ceiling(3/(1-1/a-1/b))-1 do ( for d:max(c+1,floor(1/(1-1/a-1/b-1/c))+1) thru ceiling(2/(1-1/a-1/b-1/c))-1 do ( s:[a,b,c,d,1/(1-1/a-1/b-1/c-1/d)],

if floor(s[5])+floor(-s[5])=0 then (count:count+1,print(count,":",s))

))));

 

このプログラムの実行結果 → https://fpseries.exblog.jp/33549829/

 

OEISには,1n個の単位分数に単位分数分解する方法の総数 a(n)

の値が掲載されていますね.

https://oeis.org/A006585

この数え上げ問題は,かなりの難問のようです.

a(9)以降の値は未登録です.

 

単位分数分解の方法が無限にあることは,次の等式からもわかります.

2以上の任意の整数 n, m に対して,

(1/n) = 1/(n+1) + 1/(n+1)^2 + 1/(n+1)^3 + + 1/(n+1)^m + 1/(n*(n+1)^m)

この等式から,(1/n)は,(m+1)個の異なる単位分数の和で表されることがわかります.

 

また,次の等式からも単位分数分解の方法が無限にあることがわかります.

(1/n) = 1/(n*(m+1)!) + Σ[k=1m]1/(n*(k+1)!/k).

 

[追加問題 1]

(1) (AE/CE)=(1-t)/2 ()

(2) CE=(2*a)/(5-4*t+t^2), AE=(1-t)*a*(2/(5-4*t+t^2))  ()

BC=xとおく.

線分CEの中点を M とし,∠CBM=β とおく.

MBCEMFCEであるから,3 B,M,F は同一直線上にあることがわかる.

BFC=(α/2) である.

ΔBCFに正弦定理を使って,

FC=sin(β)*(x/sin(α/2)) ---(1)

DCF=BCF-(π/2)=(π-(β+α/2))-(π/2)=(π/2)-(β+α/2).

CDF=π-(DCF+CFD)=(π/2+β-α/2).

ΔCDFに正弦定理を使って,

FC=sin(π/2+β-α/2)*(x/sin(α)), すなわち,

FC=cos(β-α/2)*(x/sin(α)) ---(2)

(1),(2)より,

sin(β)*(x/sin(α/2))=cos(β-α/2)*(x/sin(α)).

よって,

sin(α)*sin(β)=sin(α/2)*cos(β-α/2).

両辺をsin(α/2)で割って,

2*cos(α/2)*sin(β)=cos(β-α/2).

よって,

2*cos(α/2)*sin(β)=cos(β)*cos(α/2)+sin(β)*sin(α/2).

両辺を cos(α/2)*sin(β) で割って,

2=cot(β)+tan(α/2) ---(3)

線分AEの中点を N とする.

AE=2*AN=2*BA*sin(ABN)=2*x*sin((π/2-2*β)/2)

=2*x*sin(π/4-β)=(2^(1/2))*x*(cos(β)-sin(β)).

CE=2*x*sin(β).

よって,

(AE/CE)

=(2^(-1/2))*(cot(β)-1)

=(2-tan(α/2)-1)/2 ((3)より,cot(β)=2-tan(α/2))

=(1-tan(α/2))/2

=(1-t)/2.

(sin(β))^2

=1-(cos(β))^2

=1-1/(1+(tan(β))^2)

=(tan(β))^2/(1+(tan(β))^2)

=1/((cot(β))^2+1)

=1/((2-t)^2+1)

=1/(5-4*t+t^2).

よって,

sin(β)=1/(5-4*t+t^2).

CE=2*x*sin(β)=(2*a)/(5-4*t+t^2).

AE=((1-t)/2)*CE=(1-t)*a*(2/(5-4*t+t^2)).

 

[追加問題 2]

BC=(1/2)*(5*3^(1/2)-11^(1/2)) ()

ΔOAC ∽ ΔCAD である.

(まず明らかに,∠OAC=CAD.

 OA=OCであることから,∠OAC=OCA.

 CA=CDであることから,∠CAD=CDA.

よって,∠OCA=CDA.  )

OA=xとおく.

ΔOAC ∽ ΔCAD であることから,

2:x=1:2. よって,x=4. OD=OA-AD=4-1=3.

AOC=2*α,∠AOB=2*β とおく.

sin(α)=(1/2)*(AC)/(OA)=1/x=(1/4).

cos(α)=(1-(sin(α))^2)^(1/2)=(1/4)*15^(1/2).

cos(2*β)=(1/2)*(OD)/(OB)=(1/2)*(3)/x=(3/8).

cos(β)=((1+cos(2*β))/2)^(1/2)=(1/4)*11^(1/2).

sin(β)=((1-cos(2*β))/2)^(1/2)=(1/4)*5^(1/2).

よって,

BC=2*(OB)*sin(β-α)

=2*4*(sin(β)*cos(α)-cos(β)*sin(α))

=(1/2)*(5*3^(1/2)-11^(1/2)).

 

(以上)

<水の流れ:以下の情報ありがとうございます。数列として何かの規則性を見出そうとして考えましたが、断念した次第です。助かった。>

OEISには,1n個の単位分数に単位分数分解する方法の総数 a(n)

の値が掲載されていますね.

https://oeis.org/A006585

この数え上げ問題は,かなりの難問のようです.

a(9)以降の値は未登録です.

 

NO6kasama          03/21       0126分     受信  更新 3/30

寄せされた解答

NO7「よふかしのつらいおじさん」   03/23  1605   受信  更新 3/30  

  寄せされた解答

 

「水の流れ」                                          更新 3/30

参考:実は過去の第291回の問題と同じでした。