モジュール:Sequence

提供: Xenharmonic Wiki
ナビゲーションに移動 検索に移動

このモジュールについての説明文ページを モジュール:Sequence/doc に作成できます

local getArgs = require('Module:Arguments').getArgs
local p = {}

function p.contains(seq, n)
	if n > seq[#seq] then
		return nil
	end
	for i = 1, #seq do
		if seq[i] == n then
			return true
		end
	end
	return false
end

function p.get_term(frame)
	local args = getArgs(frame)
	return p._get_term(args[1], tonumber (args[2]))
end

function p._get_term(seq, n)
	local sequences = {
		["odd_limit_diff"] = p.odd_limit_diff,
		["zeta_peak"] = p.zeta_peak,
		["zeta_integral"] = p.zeta_integral,
		["zeta_gap"] = p.zeta_gap,
		["square_superparticulars"] = p.square_superparticulars
	}
	return sequences[seq][n]
end

-- OEIS A072451
-- Number of odd terms in the reduced residue system of 2*n-1
-- Corresponds to the number of new interval pairs between the (2n-3)-odd-limit
-- and the (2n-1)-odd-limit, assuming the 1-odd-limit has 1 "pair" of intervals.
p.odd_limit_diff = {
	1, 1, 2, 3, 3, 5, 6, 4, 8, 9,
	6, 11, 10, 9, 14, 15, 10, 12, 18, 12,
	20, 21, 12, 23, 21, 16, 26, 20, 18, 29,
	30, 18, 24, 33, 22, 35, 36, 20, 30, 39,
	27, 41, 32, 28, 44, 36, 30, 36, 48, 30,
	50, 51, 24, 53, 54, 36, 56, 44, 36, 48,
	55, 40, 50, 63, 42, 65, 54, 36, 68, 69,
	46, 60, 56
}

-- OEIS A117536
p.zeta_peak = {
	0, 1, 2, 3, 4, 5, 7,
	10, 12, 19, 22, 27, 31, 41, 53, 72, 99,
	118, 130, 152, 171, 217, 224, 270, 342, 422, 441, 494, 742, 764, 935, 954,
	1012, 1106, 1178, 1236, 1395, 1448, 1578, 2460, 2684, 3395, 5585, 6079, 7033, 8269, 8539,
	11664, 14348, 16808, 28742, 34691,
	-- unconfirmed data from [[The Riemann zeta function and tuning #Zeta EDO lists]]
	36269, 57578, 58973, 95524,
	102557, 112985, 148418, 212147, 241200
}

p.zeta_peak_integer = {
	0, 1, 2, 3, 5, 7, 
	10, 12, 19, 22, 31, 41, 53, 87, 
	118, 130, 171, 224, 270, 311, 472, 494, 742, 
	1065, 1106, 1395, 1578, 2460, 2684, 3566, 4231, 4973, 5585, 8269, 8539, 
	14124, 14348, 16808, 28742, 30631, 34691, 36269, 57578, 58973
}

-- OEIS A117538
p.zeta_integral = {
	2, 5, 7,
	12, 19, 31, 41, 53, 72,
	130, 171, 224, 270, 764, 954,
	1178, 1395, 1578, 2684, 3395, 7033, 8269, 8539,
	14348, 16808, 36269, 58973
}

-- OEIS A117537
p.zeta_gap = {
	2, 3, 5, 7,
	12, 19, 31, 46, 53, 72,
	270, 311, 954,
	1178, 1308, 1395, 1578, 3395, 4190,
	-- unconfirmed data from [[The Riemann zeta function and tuning#Zeta EDO lists]]
	8539,
	14348, 58973, 95524
}

-- square superparticular indices by prime limit
-- counted by OEIS A117582
-- see https://github.com/lucasaugustus/oeis/blob/main/stormer.py
p.square_superparticulars = {
	[2] = {},
	[3] = {2, 3},
	[5] = {4, 5, 9},
	[7] = {6, 7, 8, 15, 49},
	[11] = {10, 11, 21, 55, 99},
	[13] = {
		12, 13, 14, 25, 26, 27, 64, 65,
		351
	},
	[17] = {
		16, 17, 33, 34, 35, 50, 51,
		120, 169, 441
	},
	[19] = {
		18, 19, 20, 39, 56, 76, 77,
		153, 170, 209, 324,
		2431
	},
	[23] = {
		22, 23, 24, 45, 69, 91,
		161, 208, 323, 391,
		2024
	},
	[29] = {
		28, 29, 57,
		115, 116, 144, 175, 231, 289, 494, 551, 783,
		1275, 2001, 3249, 9801,
		13311
	},
	[31] = {
		30, 31, 32, 63, 92,
		125, 154, 155, 341, 342, 495, 528, 714,
		1519, 1520,
		13455
	},
	[37] = {
		36, 37, 38, 75,
		111, 185, 186, 221, 260, 406, 407, 666, 703, 961,
		1000, 1331, 1444, 1665, 1701, 2737, 3626, 5291, 7105,
		17576
	},
	[41] = {
		40, 41, 81,
		124, 204, 246, 247, 286, 287, 288, 369, 493, 575,
		1024, 1025, 1189, 1518, 1681, 1682, 1886, 2295, 3773, 4960, 5083,
		29601, 32799,
		212381
	},
	[43] = {
		42, 43, 44, 85, 86, 87,
		129, 171, 216, 259, 300, 343, 344, 559, 560, 645, 989,
		1332, 1333, 1376, 1805, 1849, 2107, 2108, 2665, 2925, 4256, 4901, 5290,
		10879, 10880, 11439,
		134849
	},
	[47] = {
		46, 47, 48, 93, 94, 95,
		187, 188, 189, 329, 375, 376, 377, 517, 704, 799, 846, 987, 988,
		1035, 1364, 1457, 1599, 2209, 2255, 2584, 4465, 6579, 7568,
		16169, 17577, 19551, 21385,
		192511
	},
	[53] = {
		52, 53, 54,
		105, 160, 265, 319, 370, 371, 424, 476, 529, 637, 689, 741, 847, 900, 901, 902,
		1377, 2015, 2755, 2915, 5831, 11661, 1749, 4929, 6149,
		23001, 23374, 27455, 36518, 62699,
		247456, 248676,
		1205645
	},
	[59] = {
		58, 59,
		117, 118, 119, 176, 235, 295, 296, 530, 531, 532, 589, 649, 650, 944, 945,
		1121, 1887, 2184, 2300, 2419, 3009, 3363, 5075, 5781, 6137, 6727, 6784, 6785, 7315, 8555,
		10620, 13689, 17051, 19942, 27435, 31212, 37583, 46138, 55224, 63426,
		956449
	},
	[61] = {
		60, 61, 62,
		121, 122, 123, 183, 184, 243, 244, 245, 305, 550, 609, 610, 611,
		1036, 1160, 1219, 1220, 1221, 1463, 1768, 1769, 1770, 1891, 2014, 2256, 2379, 2500, 2623, 2624, 5247, 5734, 5796, 8601,
		11285, 14336, 16225, 18240, 22815, 74359, 88451,
		102601, 113399, 147499, 302499
	},
	[67] = {
		66, 67, 68,
		133, 134, 135,
		469, 671, 736, 737, 805,
		1072, 1272, 1273, 1274, 1407, 1475, 1540, 2145, 2210, 2278, 2279, 2479, 3484, 3552, 3751, 3886, 3887, 4557, 4624, 4959, 5829, 5830, 7503, 7905,
		17689, 18491, 19228, 20769, 25193, 26196, 27000, 31959, 41876, 43616, 44891, 45695, 48374, 50653, 98021,
		187265, 314432,
		2018979
	},
	-- using abc conjecture c < 100 * rad(abc)^2; may be incomplete
	[71] = {
		70, 71,
		141, 142, 143, 285, 425, 426, 496, 638, 639, 780, 781, 782, 851, 924,
		1065, 1420, 1633, 2058, 3267, 3478, 3479, 3480, 3550, 3551, 3905, 4047, 4757, 4900, 7525, 8449,
		10934, 11501, 16400, 26129, 26979, 30316, 39689, 52326,
		164151, 180481, 996557,
		1608575, 2456245, 3939649
	},
	-- using abc conjecture c < 100 * rad(abc)^2; may be incomplete
	[73] = {
		72, 73, 74,
		145, 146, 147, 220, 364, 365, 511, 512, 584, 657, 729, 730, 731, 804, 875,
		1023, 1240, 1241, 1680, 2116, 2262, 2773, 2774, 3212, 3431, 4088, 4599, 5184, 5329, 5475, 5547, 5985, 8177,
		10659, 12921, 17081, 19344, 21608, 22631, 31536, 31900, 35259, 35990, 37961, 39711, 43215, 46575, 48545, 70225, 70226,
		101761, 196882, 271999, 393471, 415151, 787815,
		1215449, 2307019
	},
	-- using abc conjecture c < 100 * rad(abc)^2; may be incomplete
	[79] = {
		78, 79, 80,
		159, 236, 237, 473, 474, 475, 552, 868, 869, 870, 949,
		1026, 1105, 1106, 1184, 1343, 1421, 1739, 1975, 2211, 2449, 2450, 2924, 3081, 3477, 3870, 3871, 4029, 4186, 4346, 4899, 5292, 5452, 6479, 7425, 7504, 8215, 8295, 9164,
		10585, 11375, 12167, 13431, 14535, 16353, 17458, 29784, 33496, 35075, 46137, 46530, 70151, 70152, 75999, 79236,
		118580, 122451, 146072, 163371, 209509, 525825,
		3607297, 4409069,
		15473809
	},
	-- using abc conjecture c < 100 * rad(abc)^2; may be incomplete
	[83] = {
		82, 83, 84,
		165, 248, 249, 414, 415, 497, 664, 665,
		1079, 1080, 1161, 1495, 1825, 1826, 2241, 2408, 2738, 2821, 2989, 3403, 3485, 4233, 4234, 4731, 4897, 4898, 5395, 5643, 6391, 6889, 7138, 7139, 8216,
		10044, 10126, 10374, 10625, 12615, 12616, 13114, 15105, 15687, 18425, 20501, 23903, 28304, 36519, 50049, 51128, 53534, 67229, 68475, 79764, 83249, 95201, 97525, 97526, 98686,
		114375, 144585, 154629, 171395, 218041, 275561, 318719, 336896, 830415,
		1572351, 3100215, 6571774, 7534241
	},
	-- using abc conjecture c < 100 * rad(abc)^2; may be incomplete
	[89] = {
		88, 89, 90,
		177, 266, 267, 355, 356, 533, 624, 711, 712, 713, 800,
		1156, 1246, 1247, 1334, 1425, 1691, 2047, 2491, 2848, 2849, 3025, 3381, 3649, 4095, 4806, 4896, 5074, 5428, 5695, 6231, 6319, 6320, 6497, 7030, 7743, 8722,
		10324, 10857, 12993, 15129, 15841, 18424, 19313, 23851, 24564, 28125, 28303, 29280, 29281, 29547, 35155, 46904, 56870, 68264, 71021,
		110449, 116144, 153791, 193843, 204611, 211641, 250001, 257299, 422304, 695267,
		1270565, 1313641, 3999216, 6167344,
		407498959
	},
	-- using abc conjecture c < 100 * rad(abc)^2; may be incomplete
	[97] = {
		96, 97, 98,
		195, 290, 291, 484, 485, 581, 582, 583, 775, 776, 874, 969,
		1066, 1067, 1456, 1551, 2133, 2134, 2135, 2232, 2813, 2911, 3008, 3783, 3976, 3977, 4558, 4559, 5625, 5916, 6498, 6499, 6888, 7371, 7565, 7566, 7567, 7955, 8051,
		10864, 13775, 13870, 13871, 14161, 14840, 16684, 17459, 25025, 28519, 37926, 40256, 46656, 57133, 62952, 88560,
		101269, 103500, 113296, 117078, 126294, 129108, 139195, 224847, 313600, 329121, 431649, 442224, 447849, 908600,
		1055361, 1231803, 1555008, 1584010, 7496644,
		43184401
	},
	-- using abc conjecture c < 100 * rad(abc)^2; may be incomplete
	[101] = {
		100, 101, 201, 202, 203, 304, 403, 404, 405, 505, 506,
		1312, 1313, 1616, 1716, 2323, 2324, 2625, 2626, 2627, 2727, 3233, 3535, 4544, 5251, 5453, 6060, 6161, 6364, 6665, 6969, 7372, 7474, 7475, 8281, 8585, 8788, 8990, 9797,
		10200, 12121, 12221, 12321, 12727, 13432, 16059, 17169, 17575, 20502, 21412, 21413, 21715, 44239, 61104, 61711, 75951, 80801, 83425, 84133, 84134, 89889, 92415, 94536,
		100595, 116654, 117363, 133825, 135642, 153217, 161601, 162811, 197456, 255529, 376831, 397536, 485001, 614384,
		2383095, 3444301, 5964959, 6530356, 9814169,
		14455826,
		116026274
	},
}

return p