2021年8月1日日曜日

210801

Ruby


A282919

p(49*n+47) を出力してみた。

def s(k, n)
  s = 0
  (1..n).each{|i| s += i if n % i == 0 && i % k == 0}
  s
end

def A(ary, n)
  a_ary = [1]
  a = [0] + (1..n).map{|i| ary.inject(0){|s, j| s + j[1] * s(j[0], i)}}
  (1..n).each{|i| a_ary << (1..i).inject(0){|s, j| s - a[j] * a_ary[-j]} / i}
  a_ary
end

n = 15
a =
[7 ** 2 * 2546, 7 ** 4 * 48934, 7 ** 5 * 1418989, 7 ** 7 * 2488800, 7 ** 9 * 2394438,
 7 ** 11 * 1437047, 7 ** 12 * 4043313, 7 ** 15 * 161744, 7 ** 17 * 32136, 7 ** 18 * 31734,
 7 ** 20 * 3120, 7 ** 22 * 204, 7 ** 24 * 8, 7 ** 25]
ary = (0..13).map{|i| p [0] * i + A([[7, 4 * (i + 1)], [1, -4 * (i + 1) - 1]], n)}
(0..n).each{|i|
  print "#{(0..13).inject(0){|s, j| s + a[j] * ary[j][i]}} = "
  (0..13).each{|j|
    print " + " if j > 0
    print "#{a[j]} * #{ary[j][i]}"
  }
  puts
}

出力結果
[1, 5, 20, 65, 190, 506, 1265, 2986, 6745, 14645, 30767, 62745, 124706, 242110, 460337, 858673]
[0, 1, 9, 54, 255, 1035, 3753, 12483, 38701, 113193, 315013, 839802, 2155905, 5352252, 12894426, 30233558, 69160869]
[0, 0, 1, 13, 104, 637, 3276, 14820, 60697, 229360, 810498, 2705118, 8592857, 26134654, 76476816, 216174700, 592220696, 1576826355]
[0, 0, 0, 1, 17, 170, 1275, 7905, 42619, 206091, 912459, 3753328, 14500320, 53053498, 185046190, 618555931, 1990227519, 6186291009, 18633598578]
[0, 0, 0, 0, 1, 21, 252, 2233, 16170, 100926, 560945, 2837398, 13265679, 57989435, 239125579, 936702879, 3505361650, 12590400326, 43572202835, 145770820937]
[0, 0, 0, 0, 0, 1, 25, 350, 3575, 29575, 209405, 1312675, 7452201, 38939275, 189537775, 867436570, 3760131375, 15529994130, 61413915500, 233488417752, 856388420815]
[0, 0, 0, 0, 0, 0, 1, 29, 464, 5365, 49880, 394632, 2750969, 17296732, 99742368, 534126988, 2681856693, 12722233068, 57373155952, 247218913828, 1022189562610, 4070289420139]
[0, 0, 0, 0, 0, 0, 0, 1, 33, 594, 7667, 79101, 691119, 5299019, 36518791, 230122266, 1343028082, 7331536586, 37731144564, 184232285897, 857974579385, 3827695162667, 16420097827188]
[0, 0, 0, 0, 0, 0, 0, 0, 1, 37, 740, 10545, 119510, 1142338, 9548849, 71529474, 488650453, 3084466705, 18173253703, 100751920597, 529029597362, 2645187324766, 12651654794629, 58105915432081]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 41, 902, 14063, 173635, 1801745, 16300739, 131814181, 969824701, 6579564585, 41587633402, 246925024493, 1386436741480, 7402293438974, 37755020009290, 184685764132377]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 45, 1080, 18285, 244260, 2733804, 26606745, 230915656, 1819708110, 13198528010, 89041203249, 563420646090, 3366705675744, 19105222953420, 103448715353372, 536621238174195]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 49, 1274, 23275, 334425, 4015011, 41818315, 387605443, 3256150548, 25135003348, 180196297050, 1210028211210, 7663549175191, 46039891115155, 263630633610437, 1444741006154614]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 53, 1484, 29097, 447426, 5734918, 63638001, 627260142, 5594403499, 45779730871, 347453597091, 2466970932027, 16501339314082, 104588498225862, 631215364345159, 3642533720923593]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 57, 1710, 35815, 586815, 7997157, 94175267, 983458849, 9279004863, 80218101555, 642408637594, 4807304399931, 33855173217278, 225702273908048, 1431470152072364, 8673471170235715]
124754 = 124754 * 1 + 117490534 * 0 + 23848948123 * 0 + 2049633818400 * 0 + 96624210037866 * 0 + 2841511464047921 * 0 + 55964656476536913 * 0 + 767889588864220592 * 0 + 7475814197492884152 * 0 + 51676077116090188566 * 0 + 248951870848549443120 * 0 + 797603493910929561996 * 0 + 1532649851044531315208 * 0 + 1341068619663964900807 * 0
118114304 = 124754 * 5 + 117490534 * 1 + 23848948123 * 0 + 2049633818400 * 0 + 96624210037866 * 0 + 2841511464047921 * 0 + 55964656476536913 * 0 + 767889588864220592 * 0 + 7475814197492884152 * 0 + 51676077116090188566 * 0 + 248951870848549443120 * 0 + 797603493910929561996 * 0 + 1532649851044531315208 * 0 + 1341068619663964900807 * 0
24908858009 = 124754 * 20 + 117490534 * 9 + 23848948123 * 1 + 2049633818400 * 0 + 96624210037866 * 0 + 2841511464047921 * 0 + 55964656476536913 * 0 + 767889588864220592 * 0 + 7475814197492884152 * 0 + 51676077116090188566 * 0 + 248951870848549443120 * 0 + 797603493910929561996 * 0 + 1532649851044531315208 * 0 + 1341068619663964900807 * 0
2366022741845 = 124754 * 65 + 117490534 * 54 + 23848948123 * 13 + 2049633818400 * 1 + 96624210037866 * 0 + 2841511464047921 * 0 + 55964656476536913 * 0 + 767889588864220592 * 0 + 7475814197492884152 * 0 + 51676077116090188566 * 0 + 248951870848549443120 * 0 + 797603493910929561996 * 0 + 1532649851044531315208 * 0 + 1341068619663964900807 * 0
133978259344888 = 124754 * 190 + 117490534 * 255 + 23848948123 * 104 + 2049633818400 * 17 + 96624210037866 * 1 + 2841511464047921 * 0 + 55964656476536913 * 0 + 767889588864220592 * 0 + 7475814197492884152 * 0 + 51676077116090188566 * 0 + 248951870848549443120 * 0 + 797603493910929561996 * 0 + 1532649851044531315208 * 0 + 1341068619663964900807 * 0
5234371069753672 = 124754 * 506 + 117490534 * 1035 + 23848948123 * 637 + 2049633818400 * 170 + 96624210037866 * 21 + 2841511464047921 * 1 + 55964656476536913 * 0 + 767889588864220592 * 0 + 7475814197492884152 * 0 + 51676077116090188566 * 0 + 248951870848549443120 * 0 + 797603493910929561996 * 0 + 1532649851044531315208 * 0 + 1341068619663964900807 * 0
154043597379576030 = 124754 * 1265 + 117490534 * 3753 + 23848948123 * 3276 + 2049633818400 * 1275 + 96624210037866 * 252 + 2841511464047921 * 25 + 55964656476536913 * 1 + 767889588864220592 * 0 + 7475814197492884152 * 0 + 51676077116090188566 * 0 + 248951870848549443120 * 0 + 797603493910929561996 * 0 + 1532649851044531315208 * 0 + 1341068619663964900807 * 0
3617712763867604423 = 124754 * 2986 + 117490534 * 12483 + 23848948123 * 14820 + 2049633818400 * 7905 + 96624210037866 * 2233 + 2841511464047921 * 350 + 55964656476536913 * 29 + 767889588864220592 * 1 + 7475814197492884152 * 0 + 51676077116090188566 * 0 + 248951870848549443120 * 0 + 797603493910929561996 * 0 + 1532649851044531315208 * 0 + 1341068619663964900807 * 0
70593393646562135510 = 124754 * 6745 + 117490534 * 38701 + 23848948123 * 60697 + 2049633818400 * 42619 + 96624210037866 * 16170 + 2841511464047921 * 3575 + 55964656476536913 * 464 + 767889588864220592 * 33 + 7475814197492884152 * 1 + 51676077116090188566 * 0 + 248951870848549443120 * 0 + 797603493910929561996 * 0 + 1532649851044531315208 * 0 + 1341068619663964900807 * 0
1178875491155735802646 = 124754 * 14645 + 117490534 * 113193 + 23848948123 * 229360 + 2049633818400 * 206091 + 96624210037866 * 100926 + 2841511464047921 * 29575 + 55964656476536913 * 5365 + 767889588864220592 * 594 + 7475814197492884152 * 37 + 51676077116090188566 * 1 + 248951870848549443120 * 0 + 797603493910929561996 * 0 + 1532649851044531315208 * 0 + 1341068619663964900807 * 0
17229817230617210720599 = 124754 * 30767 + 117490534 * 315013 + 23848948123 * 810498 + 2049633818400 * 912459 + 96624210037866 * 560945 + 2841511464047921 * 209405 + 55964656476536913 * 49880 + 767889588864220592 * 7667 + 7475814197492884152 * 740 + 51676077116090188566 * 41 + 248951870848549443120 * 1 + 797603493910929561996 * 0 + 1532649851044531315208 * 0 + 1341068619663964900807 * 0
224282898599046831034631 = 124754 * 62745 + 117490534 * 839802 + 23848948123 * 2705118 + 2049633818400 * 3753328 + 96624210037866 * 2837398 + 2841511464047921 * 1312675 + 55964656476536913 * 394632 + 767889588864220592 * 79101 + 7475814197492884152 * 10545 + 51676077116090188566 * 902 + 248951870848549443120 * 45 + 797603493910929561996 * 1 + 1532649851044531315208 * 0 + 1341068619663964900807 * 0
2636785814481962651219075 = 124754 * 124706 + 117490534 * 2155905 + 23848948123 * 8592857 + 2049633818400 * 14500320 + 96624210037866 * 13265679 + 2841511464047921 * 7452201 + 55964656476536913 * 2750969 + 767889588864220592 * 691119 + 7475814197492884152 * 119510 + 51676077116090188566 * 14063 + 248951870848549443120 * 1080 + 797603493910929561996 * 49 + 1532649851044531315208 * 1 + 1341068619663964900807 * 0
28316911738879831363625420 = 124754 * 242110 + 117490534 * 5352252 + 23848948123 * 26134654 + 2049633818400 * 53053498 + 96624210037866 * 57989435 + 2841511464047921 * 38939275 + 55964656476536913 * 17296732 + 767889588864220592 * 5299019 + 7475814197492884152 * 1142338 + 51676077116090188566 * 173635 + 248951870848549443120 * 18285 + 797603493910929561996 * 1274 + 1532649851044531315208 * 53 + 1341068619663964900807 * 1
280403140023083872114273884 = 124754 * 460337 + 117490534 * 12894426 + 23848948123 * 76476816 + 2049633818400 * 185046190 + 96624210037866 * 239125579 + 2841511464047921 * 189537775 + 55964656476536913 * 99742368 + 767889588864220592 * 36518791 + 7475814197492884152 * 9548849 + 51676077116090188566 * 1801745 + 248951870848549443120 * 244260 + 797603493910929561996 * 23275 + 1532649851044531315208 * 1484 + 1341068619663964900807 * 57
2580469549453004933593920862 = 124754 * 858673 + 117490534 * 30233558 + 23848948123 * 216174700 + 2049633818400 * 618555931 + 96624210037866 * 936702879 + 2841511464047921 * 867436570 + 55964656476536913 * 534126988 + 767889588864220592 * 230122266 + 7475814197492884152 * 71529474 + 51676077116090188566 * 16300739 + 248951870848549443120 * 2733804 + 797603493910929561996 * 334425 + 1532649851044531315208 * 29097 + 1341068619663964900807 * 1710