In [1]:
queen = [6, 14, 7, 121, 115, 48, 24, 26, 33, 45, 23, 37, 23, 12, 28, 33, 17, 12, 16, 12, 7, 19, 18, 22, 12, 39, 90, 103, 51, 70, 67, 64, 88, 55, 78, 104, 93, 97, 98, 119, 138, 79, 92, 91, 74, 64, 53, 73, 84, 69, 64, 92, 73, 61, 79, 77, 98, 106, 127, 107, 102, 111, 96, 75, 57, 42, 32, 19, 27, 41, 45, 60, 32, 58, 47, 44, 57, 67, 65, 70, 58, 58, 54, 66, 62, 74, 77, 63, 69, 36, 59, 69, 96, 59, 46, 61, 43, 62, 76, 77, 98, 62, 124, 277, 256, 325, 342, 279, 346, 353, 338, 374, 416, 480, 399, 471, 355, 392, 332, 376, 450, 358, 350, 317, 336, 421, 411, 466, 375, 414, 446, 468, 407, 472, 449, 418, 341, 358, 343, 385, 436, 457, 356, 414, 402, 411, 314, 512, 492, 352, 335, 413, 474, 416, 420, 452, 448, 498, 443, 402, 365, 363, 359, 262, 366, 283, 275, 297, 304, 370, 313, 318, 328, 371, 406, 512, 511, 406, 395, 422, 427, 473, 427, 515, 397, 404, 441, 496, 531, 642, 756, 641, 600, 627, 568, 730, 768, 856, 765, 698, 647, 650, 709, 714, 750, 608, 498, 467, 463, 562, 562, 580, 415, 457, 307, 290, 384, 625, 539, 436, 372, 561, 590, 765, 930, 915, 868, 709, 363, 376, 459, 445, 729, 597, 511, 602, 460, 602, 679, 564, 468, 473, 428, 467, 439, 553, 687, 623, 480, 451, 521, 489, 506, 479, 474, 469, 596, 1335, 954, 1587, 738, 1398, 2094, 1933, 1929, 2233, 2259, 2303, 1836, 1750, 2074, 3228, 4633, 6139, 5683, 4786, 4520, 4556, 4675, 4803, 4384, 4794, 4021, 4131, 4286, 4308, 3779, 3943, 4142, 4067, 4114, 4592, 6001, 6384, 7343, 6288, 5952, 5867, 5227, 4575, 4807, 4664, 4630, 4319, 4020, 3862, 4149, 4312, 5045, 5277, 5054, 4629, 5022, 4676, 5020, 5515, 5534, 4513, 4167, 5494, 4855, 4406, 4521, 6041, 5508, 4913, 3428, 3701, 3896, 4708, 3926, 3357, 3391, 3246, 3177, 3053, 3625, 3354, 2972, 3027, 3438, 4259, 3910, 3956, 3843, 3160, 3431, 3606, 3641, 4077, 4731, 4612, 4217, 4014, 4157, 4274, 4455, 4885, 4395, 3635, 3733, 3863, 3783, 3689, 3701, 4020, 3121, 3282, 3415, 3322, 3645, 3826, 3588, 3321, 3709, 3377, 2304, 2416, 2568, 2162, 1781, 1109, 1131, 1558, 1686, 2052, 2393, 2021, 1954, 2116, 2107, 2129, 2476, 2426, 1918, 1784, 1412, 1481, 1819, 2312, 2576, 2133, 1979, 1846, 1861, 2032, 2375, 2515, 2116, 1919, 1946, 2017, 2156, 2616, 2622, 2277, 2036, 2101, 2063, 2220, 2367, 2225, 1821, 1879, 1889, 1802, 1902, 2406, 2235, 2131, 1880, 1889, 2019, 2247, 2440, 2642, 2456, 2191, 2114, 2243, 2203, 2661, 2646, 2311, 2085, 2079, 2066, 2025, 2224, 2501, 2166, 2034, 1900, 2139, 2450, 2692, 2689, 1965, 2109, 1662, 2024, 2117, 2479, 2591, 2123, 1950, 1892, 1962, 1899, 2276, 2767, 2496, 2341, 2364, 2356, 2256, 2617, 2710, 2057, 1855, 2764, 3093, 2894, 2910, 3463, 2799, 1908, 2224, 2319, 2358, 2948, 3039, 2189, 1969, 2022, 1779, 1821, 2294, 2337, 1856, 1703, 1740, 1667, 1798, 1950, 2031, 1762, 1730, 1790, 2054, 2202, 1888, 1981, 1961, 1975, 2044, 2258, 2963, 1882, 2078, 1770, 1669, 1752, 1828, 1731, 2160, 2386, 1852, 1649, 1757, 1685, 1690, 2112, 2280, 1794, 1708, 1740, 1440, 1409, 1603, 1628, 1329, 1210, 1336, 1218, 1443, 1695, 1898, 1450, 1322, 1373, 1364, 1402, 2235, 2402, 1833, 1705, 1895, 1977, 1978, 2222, 2407, 1846, 1843, 1827, 1813, 1826, 2087, 2260, 1741, 1655, 1676, 1781, 1781, 2211, 2234, 1916, 1960, 2007, 2218, 8599, 4543, 4376, 3320, 2299, 2654, 3943, 3591, 3539, 3591, 3027, 2358, 2363, 2509, 1849, 2089, 2400, 2145, 1948, 1730, 1966, 2109, 2990, 2523, 2188, 2397, 2767, 2764, 4834, 4599, 4566, 4066, 3968, 3513, 4413, 8939, 10429, 12313, 11632, 14056, 10049, 9460, 9012, 10872, 11734, 9405, 8619, 9060, 14648, 18165, 21223, 19458, 16642, 15409, 15681, 16413, 19943, 23415, 26863, 23065, 21139, 21253, 20132, 19647, 20891, 23644, 20546, 18269, 18130, 18228, 21047, 22184, 22679, 19425, 18344, 17851, 18119, 18201, 20404, 20884, 18825, 17406, 17174, 20258, 33256, 42761, 47244, 39736, 34102, 33535, 34516, 37591, 44699, 48462, 40209, 34126, 38177, 39986, 35508, 35713, 40144, 35008, 29204, 29909, 30277, 32582, 38220, 44041, 36480, 31420, 31786, 32403, 33846, 38349, 41447, 33081, 29430, 29667, 31213, 32905, 34343, 36488, 32915, 31312, 31302, 31205, 33205, 37208, 38730, 32825, 31814, 32815, 29472, 34637, 39798, 42329, 35835, 32690, 33187, 27531, 35739, 39620, 42696, 40120, 38390, 35288, 31252, 37076, 42386, 43363, 31212, 36387, 33080, 35263, 36108, 37949, 42539, 36333, 31871, 28700, 28240, 21399, 24971, 24935, 19591, 18565, 20048, 24448, 25824, 27021, 27115, 24018, 21273, 23267, 22674, 26546, 29792, 30431, 26005, 26035, 25446, 26745, 29167, 30094, 20908, 18295, 16066, 15622, 19594, 30729, 34569, 35449, 30505, 28321, 30079, 32234, 34167, 40566, 44225, 36381, 32513, 33696, 34708, 35029, 41222, 45722, 36511, 31465, 36558, 40279, 47253, 51297, 38413, 38259, 32117, 34446, 36499, 39980, 46176, 51560, 42927, 37727, 36029, 33481, 41114, 46274, 51704, 40948, 35202, 38107, 39149, 39874, 47893, 50033, 41210, 35973, 37836, 37515, 35597, 38044, 47104, 44171, 45835, 38310, 38137, 40504, 48599, 55307, 46133, 39503, 40511, 42484, 49148, 53370, 55858, 46535, 42567, 42034, 37253, 42462, 51387, 56553, 47045, 43162, 43916, 62046, 49820, 52916, 57575, 45271, 40304, 43334, 45161, 44617, 51907, 56193, 46507, 41657, 42673, 44388, 44075, 52302, 56558, 46037, 40692, 43389, 44786, 47569, 50171, 51787, 44973, 42766, 45351, 46433, 50060, 50070, 45597, 43251, 43552, 45441, 45134, 50206, 68043, 48763, 47863, 44305, 46105, 48621, 50236, 58139, 61384, 49736, 43446, 39572, 42370, 46459, 49902, 55631, 45030, 39067, 41491, 42939, 49235, 56549, 58626, 48699, 43550, 43780, 46770, 55343, 59107, 61196, 50904, 43722, 34931, 33968, 47571, 40335, 42949, 33836, 31442, 33430, 33594, 35600, 39876, 42599, 35472, 32167, 32494, 32932, 35638, 42390, 44249, 36850, 32759, 32231, 32636, 34095, 39858, 42711, 34832, 31789, 31172, 31801, 33611, 39428, 42496, 34407, 32108, 32023, 33704, 34580, 37888, 43035, 34715, 32436, 31160, 32342, 33194, 38396, 41247, 35107, 31387, 32262, 35099, 35196, 37126, 41553, 35092, 30455, 30852, 31475, 33087, 37192, 38354, 33175, 30111, 30894, 31889, 32231, 36676, 38893, 31619, 29036, 30289, 31562, 32082, 36642, 38946, 33724, 33058, 31842, 32253, 34909, 36193, 33146, 29903, 27265, 26782, 26984, 29285, 32616, 37927, 32256, 28424, 29172, 30901, 31820, 36862, 39479, 33327, 29294, 30990, 31104, 36084, 43744, 62537, 46049, 33309, 34792, 34113, 37623, 41618, 48376, 42560, 42364, 76312, 64113, 59665, 59237, 58029, 51009, 46009, 45307, 46061, 48240, 52633, 56524, 47048, 42868, 45415, 46365, 46492, 51661, 53205, 46004, 41759, 43105, 45493, 43293, 50614, 51766, 43758, 42037, 42379, 42851, 46727, 49514, 48372, 39860, 37619, 37113, 42769, 47360, 48365, 46334, 39434, 37377, 37957, 39555, 40054, 38345, 44201, 38857, 38686, 38598, 43679, 40526, 47558, 48384, 42564, 41667, 42662, 43290, 44336, 48177, 48649, 42919, 41960, 42431, 42719, 43299, 46942, 46895, 40492, 39722, 40354, 40270, 41063, 45584, 45868, 39784, 37615, 36619, 37475, 38881, 43107, 44594, 40153, 37771, 38193, 37717, 39447, 46105, 46617, 40360, 38335, 42576, 41246, 40682, 47093, 49417, 41635, 39109, 41049, 40486, 43149, 50152, 58076, 94356, 310092, 153708, 83235, 70333, 73214, 73842, 60045, 53295, 52355, 56956, 62092, 71149, 74141, 58006, 52803, 52477, 53839, 62911, 73975, 77561, 60227, 56554, 57618, 55235, 55744, 73337, 75720, 62387, 52746, 53409, 53614, 54599, 64812, 66827, 55907, 48187, 47896, 46431, 46008, 54526, 58298, 49433, 43985, 41110, 46818, 49545, 58855, 63597, 49192, 45933, 45782, 48647, 49598, 56161, 61312, 51096, 48300, 48075, 47255, 49101, 55042, 61400, 53424, 49959, 46463, 48720, 49037, 55641, 61859, 50276, 44004, 44276, 46499, 49783, 58584, 62766, 52409, 48868, 50781, 63206, 164352, 97894, 90289, 67456, 55922, 55964, 55945, 55255, 59261, 61635, 50343, 44752, 46552, 46713, 49090, 53523, 58829, 48974, 42934, 49831, 48780, 48157, 54500, 57689, 47330, 44267, 45905, 57616, 66929, 70746, 75997, 71772, 74749, 67095, 67052, 68738, 77326, 115148, 79808, 68916, 61989, 63817, 68560, 78570, 79940, 67596, 62640, 64680, 67462, 69071, 79662, 85542, 69993, 68812, 82977, 77421, 74516, 84513, 93809, 74590, 67714, 67866, 68738, 72827, 89207, 93555, 69082, 61494, 60490, 60951, 63668, 72094, 80181, 64782, 58037, 63075, 74326, 89055, 96900, 92786, 78974, 71777, 79511, 78742, 79597, 91417, 93892, 77329, 73059, 71527, 71533, 73786, 82908, 89398, 73484, 82706, 77404, 75231, 78340, 86637, 89418, 73777, 69880, 70043, 71372, 71427, 85804, 90965, 73516, 65388, 68666, 69225, 69230, 77914, 85123, 76978, 71083, 69138, 72496, 76741, 84963, 93588, 76185, 67810, 69099, 71616, 109804, 188856, 167012, 114210, 93405, 95505, 105155, 105733, 104540, 99818, 89426, 85923, 80597, 78361, 82854, 87287, 94991, 78793, 71288, 70039, 72593, 76856, 90229, 101922, 87270, 75758, 76399, 91224, 101718, 118736, 121830, 98491, 93438, 87381, 80691, 84821, 93223, 99647, 85512, 76838, 78454, 79159, 83028, 93808, 110245, 92664, 77599, 79286, 82068, 84885, 88907, 89379, 78419, 70508, 78119, 103175, 83837, 91723, 95874, 85715, 78090, 76076, 80315, 83066, 91572, 98116, 85287, 74677, 76313, 78541, 81760, 88000, 91460, 80954, 70421, 79619, 79588, 81823, 90819, 94462, 79194, 71748, 70962, 73019, 72875, 81250, 81112, 68942, 67653, 70678, 66113, 69524, 77446, 78227, 65923, 64469, 67137, 66840, 67855, 75318, 79308, 68472, 70060, 132772, 110192, 96196, 107548, 97042, 96109, 91155, 84522, 75839, 67047, 63144, 80759, 61505, 67793, 69217, 68506, 71518, 101851, 107598, 81595, 82447, 81350, 79153, 77531, 90453, 93054, 73589, 71272, 71413, 69959, 75009, 82058, 82778, 108010, 113402, 105150, 95877, 92169, 98317, 100312, 84710, 80412, 78882, 81261, 85958, 98206, 99996, 86313, 82099, 80405, 79639, 82226, 88072, 92175, 80125, 75845, 77888, 143308, 92955, 95825, 96418, 85738, 77076, 76512, 78971, 87692, 95967, 100848, 80681, 72108, 70088, 72820, 79514, 90493, 97923, 76087, 69188, 72808, 71754, 74010, 85876, 95423, 82799, 75640, 80019, 82789, 107055, 104386, 111869, 87817, 79940, 80198, 82190, 83842, 96271, 105085, 84171, 78414, 75002, 79429, 84385, 96951, 106194, 85602, 78151, 81185, 102391, 108426, 130495, 128693, 106336, 88307, 87279, 87935, 98271, 110380, 119668, 104619, 89184, 83221, 90235, 90585, 102877, 125725, 96461, 89336, 91737, 92896, 94875, 104702, 110547, 93378, 86757, 91602, 95336, 94109, 117706, 150592, 107447, 85301, 89776, 78443, 83011, 97130, 105852, 86704, 77974, 78459, 85340, 94564, 109451, 102964, 83485, 75491, 81160, 82612, 85561, 93173, 105024, 81461, 71861, 82350, 84835, 81600, 91698, 101460, 90262, 94608, 92524, 76627, 80209, 87913, 93060, 89697, 134432, 80855, 75369, 72858, 86698, 94353, 78328, 71949, 76227, 82088, 88726, 90971, 98283, 85251, 83249, 80001, 81987, 87127, 103294, 125874, 97489, 91151, 83665, 101882, 98384, 101889, 118176, 94885, 90583, 89581, 86162, 88016, 103717, 114702, 90509, 85950, 87466, 86555, 87873, 101139, 110065, 92498, 80723, 77654, 80403, 86682, 98837, 114383, 84111, 74372, 55427, 400, 115, 146, 178, 123, 101, 111, 119, 107, 128, 132, 148, 93, 86, 92, 113, 107, 144, 135, 32313, 55494, 45201, 43969, 62411, 79947, 59779, 50892, 51562, 63349, 66713, 72569, 83049, 54493, 44948, 50849, 76248, 83853, 90368, 97443, 94380, 83684, 80276, 80023, 86785, 92427, 102389, 81286, 79652, 83938, 87610, 92545, 103125, 113060, 89897, 77764, 79557, 87354, 85885, 101394, 107225, 92624, 88258, 88196, 89619, 82947, 95889, 105657, 87181, 74855, 82576, 86709, 85687, 99779, 113129, 98900, 81395, 88328, 89847, 91368, 99008, 107261, 88029, 73457, 73571, 78715, 81131, 96143, 97880, 79740, 75902, 73944, 79450, 82547, 90968, 96150, 80236, 72938, 77525, 82815, 83968, 90172, 102266, 82264, 75890, 73887, 78808, 80817, 85991, 97041, 83933, 93880, 87113, 80791, 81744, 90653, 98654, 81699, 73660, 87397, 77840, 83336, 96701, 93341, 74868, 72151, 78227, 83382, 89438, 93948, 90380, 75046, 72098, 75238, 72818, 76239, 77885, 78838, 67334, 65455, 69391, 69240, 75511, 81344, 91360, 74685, 72834, 70710, 78975, 79943, 88896, 87267, 76240, 73783, 76471, 74643, 76783, 86895, 85654, 72621, 71799, 75722, 64416, 59621, 64504, 64043, 57308, 53909, 56921, 57419, 58628, 67370, 68733, 58313, 55054, 58189, 59098, 60490, 65569, 69452, 57210, 56106, 58927, 58827, 60822, 70033, 73885, 61716, 58700, 60666, 59420, 62253, 68544, 72388, 60454, 57377, 61731, 63209, 109891, 81890, 78541, 63788, 54614, 52614, 54986, 57144, 64572, 69461, 60103, 52166, 51836, 52125, 54008, 62734, 67868, 58451, 55118, 54985, 54578, 55870, 64880, 69950, 57548, 48903, 53283, 54820, 57128, 64671, 67928, 57350, 54207, 55522, 56844, 60799, 67110, 74191, 58856, 54497, 52039, 54126, 56965, 67229, 72334, 57288, 53182, 54886, 59484, 66230, 71403, 73881, 62304, 57806, 60512, 60520, 69977, 76193, 78112, 71572, 56718, 58486, 59145, 61717, 71834, 77440, 65931, 59786, 58627, 58974, 62393, 71816, 81122, 63642, 55674, 53871, 62712, 64741, 73533, 84409, 96442, 70022, 63617, 60310, 68319, 77714, 83611, 65657, 66286, 84953, 77168, 78526, 78975, 86815, 85016, 83582, 68029, 66980, 68027, 78963, 79250, 61531, 55945, 62707, 62685, 66903, 72866, 70904, 59857, 61237, 71871, 74891, 62862, 66642, 66468, 59510, 66669, 108167, 61116, 63441, 65540, 71114, 60142, 66775, 65241, 70522, 77819, 79707, 81364, 67111, 67058, 72040, 74715, 74975, 75151, 79900, 65672, 63833, 67505, 65954, 70109, 78877, 79741, 65513, 59650]
In [2]:
queen_d0 = 1217548800000
In [3]:
beeges = [4, 17, 135, 349, 190, 327, 265, 257, 386, 305, 332, 338, 315, 824, 733, 1312, 1257, 1249, 1450, 1625, 1345, 817, 1177, 1212, 883, 1122, 1359, 1235, 1108, 1349, 1280, 880, 1050, 1310, 621, 400, 688, 432, 396, 451, 482, 475, 404, 415, 642, 436, 293, 307, 215, 205, 197, 264, 698, 597, 723, 639, 501, 701, 521, 691, 737, 1078, 783, 916, 631, 648, 953, 1157, 1237, 1085, 1016, 939, 1283, 1472, 1639, 2547, 1649, 1442, 1687, 1486, 1945, 2548, 3032, 2141, 1710, 1485, 1537, 1592, 2055, 2025, 1393, 1155, 1480, 1497, 1675, 1929, 2222, 1910, 1805, 2010, 1838, 1953, 2530, 2759, 1669, 1609, 1717, 1731, 1935, 2716, 2365, 1825, 1415, 1581, 2154, 3614, 3230, 3964, 2986, 2575, 2790, 2721, 2445, 3235, 6115, 2284, 1991, 1848, 1793, 2109, 2490, 2364, 1757, 1445, 1571, 1854, 2058, 2574, 2753, 2242, 2086, 2096, 2312, 2044, 2531, 2925, 3169, 2650, 2998, 2193, 1321, 1889, 3321, 3116, 3628, 2404, 2427, 2614, 3696, 4125, 2780, 2231, 2657, 2340, 2043, 2832, 3408, 2340, 2516, 2564, 2770, 2933, 3311, 3769, 3223, 2658, 2776, 3328, 2728, 3254, 3714, 2966, 2315, 2520, 3150, 3181, 4587, 6095, 4932, 4483, 4273, 4418, 4316, 4992, 5577, 4507, 3480, 3535, 4516, 5825, 6749, 7072, 5833, 5178, 5072, 4933, 5904, 6312, 7105, 6100, 5137, 5499, 5364, 4555, 4818, 5781, 5209, 4140, 3699, 3841, 3831, 4758, 5439, 4606, 3514, 3779, 3987, 4143, 4723, 5120, 4203, 3894, 3859, 4276, 4030, 4455, 4518, 4215, 3805, 3800, 3747, 3994, 4586, 4820, 3750, 3575, 3757, 3695, 3884, 4556, 5333, 4385, 3939, 3828, 3038, 3318, 3672, 3917, 3187, 3050, 3462, 2894, 2898, 2965, 3784, 2816, 2712, 3149, 3016, 2912, 3591, 3565, 2861, 2739, 3167, 2764, 2453, 2751, 2944, 2628, 2537, 2842, 2787, 2574, 2350, 3305, 2354, 2123, 2463, 2856, 2817, 3354, 3548, 2881, 2582, 2719, 3005, 3049, 3449, 2543, 2054, 2035, 1923, 2267, 4883, 6136, 6699, 5613, 5968, 6730, 6887, 7225, 7954, 8653, 7517, 7164, 6952, 7358, 7587, 8266, 10397, 8055, 6243, 7607, 7892, 8981, 9587, 11085, 13890, 11679, 12282, 12770, 14500, 17402, 19384, 15740, 14536, 13753, 11339, 13993, 16645, 18784, 16048, 15163, 13912, 13616, 14749, 16730, 17579, 13504, 12570, 13110, 12950, 12939, 15483, 17853, 14916, 13331, 14176, 13383, 14516, 17670, 20251, 17619, 15527, 15489, 16240, 18352, 20348, 21402, 17788, 16116, 15687, 11030, 8936, 7537, 9259, 7757, 6703, 6962, 7413, 7948, 8079, 8628, 6847, 6203, 6253, 6061, 6668, 7165, 8740, 6949, 5871, 6017, 6475, 5953, 6712, 7613, 6127, 5534, 5340, 5428, 5918, 6812, 7486, 6130, 5431, 5748, 6000, 6384, 6991, 6704, 6138, 6217, 6563, 6667, 7786, 10214, 7538, 7140, 6043, 6435, 6742, 6831, 7682, 8648, 7436, 5938, 5631, 5650, 14331, 17285, 20030, 16233, 14371, 14267, 14319, 14992, 17828, 20565, 16175, 14565, 15182, 15698, 16678, 21207, 27731, 22370, 18536, 16750, 16760, 21263, 21030, 23641, 17929, 16265, 16575, 16969, 19431, 21205, 23423, 18608, 16677, 16578, 16470, 20197, 23966, 27082, 22239, 20610, 19811, 20671, 21857, 25681, 28784, 24869, 25022, 23150, 21620, 23840, 27518, 30453, 24144, 21268, 21008, 20680, 22028, 24503, 28198, 23451, 20493, 20646, 20929, 22289, 25864, 28885, 23332, 20478, 20810, 22435, 22323, 27060, 30278, 25424, 19887, 18367, 19430, 22116, 25204, 26922, 22607, 20220, 21579, 21048, 21488, 26117, 28907, 22930, 19884, 22219, 21698, 22542, 25608, 28821, 23574, 20931, 22236, 21749, 23486, 24001, 23906, 23018, 21413, 20993, 21023, 21722, 22863, 26728, 22603, 18556, 19427, 19654, 20575, 23579, 26096, 22320, 19896, 20678, 20621, 22646, 25587, 26713, 22258, 20016, 19247, 19715, 19824, 23573, 30724, 24035, 20940, 19816, 20436, 20642, 24296, 25507, 21803, 20763, 20171, 20730, 22042, 23384, 24615, 20786, 18717, 18786, 19760, 20520, 24735, 26227, 21719, 19740, 20354, 21456, 21888, 26667, 27782, 23217, 21153, 21062, 21199, 22811, 24299, 26563, 21618, 19125, 19352, 19862, 21092, 22678, 23607, 20301, 18886, 19161, 19317, 19973, 23434, 24237, 20677, 18805, 19777, 20806, 20975, 22948, 24477, 21585, 20605, 21692, 21961, 22325, 23928, 25761, 22048, 20313, 20863, 21281, 22053, 24465, 25827, 22339, 19383, 20443, 20970, 21240, 23737, 23806, 21663, 20715, 20414, 20636, 24699, 26671, 26070, 22181, 21122, 20863, 20836, 21968, 24866, 26523, 23085, 20735, 22169, 22542, 21538, 24978, 27963, 23184, 20123, 21113, 20741, 22321, 23802, 26180, 23548, 21172, 19990, 21653, 21359, 25327, 27512, 22539, 20384, 21772, 21912, 22455, 26134, 29928, 24926, 23105, 21738, 21449, 23238, 25428, 28338, 21675, 19464, 21256, 22035, 21731, 25570, 28147, 23714, 20106, 21596, 22116, 22290, 26963, 30183, 24353, 21461, 20225, 21029, 22138, 27189, 28859, 23973, 20958, 19548, 21754, 23251, 25472, 28148, 22199, 21328, 24519, 25033, 24660, 27871, 31731, 26145, 23901, 23910, 23327, 23887, 27413, 31118, 25969, 21808, 22050, 23166, 24278, 27935, 29892, 24675, 21475, 21473, 23088, 24136, 27889, 30780, 28044, 26163, 23621, 24630, 26601, 28009, 30002, 24946, 21256, 22862, 22903, 23126, 25585, 28617, 24237, 20936, 20851, 21437, 22490, 26528, 28148, 22432, 20017, 20290, 20898, 20810, 24022, 26258, 21392, 20435, 21198, 20838, 23336, 23181, 26513, 22388, 20226, 19596, 20648, 21892, 24480, 42682, 23058, 20903, 19398, 23487, 27250, 28969, 29206, 25624, 21270, 22457, 22690, 24352, 28026, 31437, 25936, 23896, 24773, 24124, 25353, 30045, 31620, 25384, 22434, 23567, 24136, 24577, 27450, 29445, 23733, 21384, 21113, 21942, 22802, 25865, 29595, 23446, 20947, 21242, 21933, 21998, 25888, 28680, 23735, 20820, 21438, 21931, 23378, 26996, 28622, 25214, 22511, 20884, 21454, 21333, 23984, 27603, 22448, 20129, 20973, 20647, 22175, 25997, 27105, 21868, 19434, 19013, 21447, 22175, 24949, 28078, 22377, 19698, 19893, 21640, 22946, 25336, 26467, 22430, 20019, 20631, 20651, 23957, 25563, 28712, 23007, 20887, 22015, 21961, 22398, 25408, 27249, 23130, 20142, 20604, 21029, 23296, 23909, 26732, 25083, 22761, 22260, 22436, 22812, 26188, 69389, 45959, 28786, 27010, 30980, 30284, 33849, 37496, 38451, 31085, 27754, 26445, 27214, 28876, 30880, 27410, 25371, 25664, 22667, 23707, 26103, 29205, 24572, 22133, 22599, 23310, 23712, 25654, 28668, 23939, 21316, 22209, 22383, 27303, 29081, 28093, 167462, 414565, 157669, 89010, 70009, 60390, 58159, 47435, 39120, 38103, 36793, 36626, 40736, 44222, 37296, 31750, 30393, 29792, 31143, 36666, 36736, 30804, 27708, 27244, 26972, 27848, 30664, 33228, 28216, 24242, 24728, 25681, 25636, 27509, 29549, 24630, 22845, 23308, 21800, 22270, 26426, 27089, 22304, 20870, 21633, 21137, 21245, 23870, 25436, 21969, 20622, 21119, 21012, 23168, 30394, 34181, 28911, 26150, 26011, 25712, 23324, 23489, 31334, 22656, 24753, 26382, 25577, 26527, 35014, 34745, 28097, 25302, 25822, 24601, 25038, 27398, 29778, 25747, 24258, 24043, 24631, 25420, 27986, 29976, 28079, 27171, 29742, 28804, 28863, 32310, 34141, 29869, 27820, 27640, 28506, 30346, 35273, 36591, 30356, 26723, 27689, 27612, 30927, 36576, 42613, 36564, 32919, 32261, 30936, 31006, 34470, 35068, 30491, 26565, 26699, 27715, 25277, 27432, 27402, 21227, 18719, 18120, 19136, 23009, 24183, 27479, 21779, 19023, 19520, 19617, 20212, 23580, 25944, 22079, 20173, 21333, 22286, 23515, 27713, 29000, 23514, 21054, 20412, 21321, 22459, 25587, 28467, 23234, 20408, 20440, 21510, 21967, 25457, 28279, 23822, 20659, 21454, 21613, 21929, 25869, 28529, 22800, 20539, 20929, 21242, 23789, 27859, 29594, 25464, 22176, 21335, 21602, 21713, 24948, 29440, 23614, 21113, 22099, 22687, 23391, 25475, 26923, 22213, 20894, 21675, 21988, 22425, 24547, 29784, 23799, 21044, 21325, 21012, 20178, 23796, 25435, 20829, 18853, 19695, 21729, 22780, 24944, 25404, 21102, 19420, 20472, 21796, 21385, 24732, 25848, 19786, 19613, 20052, 20638, 25706, 31155, 25954, 23065, 25811, 21820, 19661, 20269, 22337, 25196, 25724, 46201, 22209, 19287, 18864, 22130, 24519, 20417, 18340, 19890, 19484, 20538, 23893, 28726, 22743, 20172, 20361, 20408, 21457, 25785, 27620, 22519, 20979, 22542, 23604, 24147, 27691, 30166, 23794, 20520, 23736, 22459, 23661, 26730, 29629, 24160, 22647, 22647, 22332, 22697, 26452, 29613, 24259, 21375, 20822, 20893, 22732, 25164, 28177, 22757, 20769, 20366, 20068, 22000, 25979, 29937, 23693, 19883, 23676, 23590, 24882, 28125, 32403, 23546, 22964, 22794, 23880, 25208, 30719, 33163, 27676, 24764, 23512, 23388, 24023, 24724, 32513, 25340, 22399, 20454, 21638, 20467, 18285, 19228, 20337, 16585, 17554, 16205, 17265, 19484, 19506, 17663, 14563, 15608, 15407, 16885, 19212, 23179, 26334, 17327, 18989, 21258, 21724, 39729, 28080, 20546, 17516, 17560, 18789, 18203, 21729, 25053, 19200, 17779, 17580, 19180, 19357, 21272, 22910, 18685, 16066, 18195, 18720, 17810, 21143, 23125, 18095, 17301, 18715, 19044, 20807, 22889, 24512, 19670, 17455, 18333, 19431, 19232, 22894, 24423, 18858, 20230, 18943, 19536, 20785, 23031, 24237, 20150, 18818, 19366, 19860, 21618, 22493, 24618, 20590, 18761, 17878, 17112, 18974, 20116, 22417, 19319, 17579, 17155, 16990, 17982, 19762, 22219, 19269, 21964, 17037, 16081, 17058, 18646, 20074, 17172, 16160, 17378, 18119, 21393, 25541, 27368, 22115, 21372, 21858, 21766, 28865, 25588, 24794, 21057, 28412, 32949, 32557, 33976, 36809, 41140, 35393, 31787, 33921, 33193, 33768, 37116, 37976, 32733, 29849, 32605, 31773, 32890, 37122, 39468, 32851, 30898, 31617, 33455, 34732, 39240, 38894, 35237, 30952, 32774, 34052, 34351, 40703, 41918, 35296, 32802, 33462, 33827, 34588, 37516, 40237, 33642, 30925, 32042, 32553, 34581, 38927, 44679, 36847, 32479, 33647, 33002, 35035, 40153, 43119, 38421, 32532, 33362, 34126, 36111, 37376, 43965, 37881, 35324, 34482, 35785, 37868, 45546, 46223, 38517, 32832, 32927, 32421, 35028, 41370, 45267, 39007, 33213, 33678, 33018, 35106, 41827, 45959, 40119, 34845, 35841, 37098, 40899, 44373, 48000, 39456, 35013, 37713, 40061, 41302, 46192, 52749, 43302, 39138, 37954, 41136, 41994, 47469, 52192, 41451, 36577, 36243, 38546, 41285, 48605, 54202, 44600, 38391, 38972, 41110, 43647, 47737, 54029, 47243, 39383, 40472, 42037, 44235, 46902, 52036, 42790, 36688, 35288, 37977, 38282, 43293, 46884, 39777, 33707, 34442, 33895, 35720, 44848, 48322, 38010, 34219, 35443, 34783, 38281, 40387, 45566, 35461, 31320, 33764, 36071, 37418, 42062, 45352, 37059, 33427, 33140, 31564, 35347, 41657, 46845, 37733, 32452, 34207, 34868, 38747, 41578, 46375, 49053, 38710, 50559, 42814, 39313, 41050, 46962, 43142, 44254, 86015, 43099, 38260, 40542, 44234, 39830, 37307, 36498, 39074, 40398, 47298, 60641, 49537, 44478, 44406, 44019, 43561, 48604, 56084, 45787, 39411, 40645, 39893, 42082, 50117, 54212, 43857, 39518]
In [4]:
beeges_d0 = 1257724800000
In [5]:
ny = [10, 11, 9, 9, 15, 6, 2, 12, 5, 18, 20, 11, 23, 23, 12, 5, 13, 15, 17, 15, 21, 6, 9, 13, 21, 15, 21, 19, 14, 10, 18, 10, 12, 13, 14, 14, 16, 8, 20, 36, 21, 16, 6, 17, 18, 29, 20, 20, 16, 6, 7, 23, 18, 17, 17, 23, 22, 11, 16, 33, 32, 22, 9, 12, 24, 22, 40, 22, 22, 29, 43, 41, 26, 59, 179, 80, 54, 47, 52, 43, 27, 46, 53, 38, 28, 37, 71, 64, 45, 20, 32, 35, 33, 44, 49, 36, 49, 48, 31, 75, 115, 40, 55, 83, 34, 48, 33, 33, 39, 27, 98, 42, 34, 34, 53, 85, 61, 62, 52, 57, 52, 77, 55, 67, 53, 52, 58, 64, 106, 70, 120, 134, 88, 113, 107, 102, 92, 76, 101, 58, 50, 146, 248, 140, 157, 166, 142, 101, 147, 134, 109, 95, 100, 101, 87, 96, 105, 84, 107, 94, 111, 68, 89, 93, 126, 164, 120, 157, 166, 171, 185, 126, 137, 115, 93, 92, 127, 130, 81, 70, 69, 70, 51, 40, 46, 85, 111, 91, 272, 206, 157, 74, 98, 98, 300, 209, 217, 151, 124, 104, 168, 126, 242, 277, 262, 100, 79, 144, 116, 88, 77, 82, 76, 189, 131, 126, 71, 107, 211, 177, 143, 206, 240, 235, 211, 300, 265, 304, 251, 208, 184, 171, 201, 200, 202, 215, 288, 272, 141, 213, 239, 200, 302, 267, 194, 179, 173, 150, 175, 1645, 482, 269, 232, 251, 207, 195, 197, 179, 211, 210, 243, 321, 332, 277, 257, 273, 254, 213, 229, 187, 197, 172, 126, 153, 175, 223, 239, 287, 298, 265, 148, 122, 103, 184, 113, 101, 103, 86, 139, 207, 173, 119, 150, 129, 69, 104, 193, 152, 113, 103, 119, 196, 263, 173, 189, 94, 142, 253, 252, 215, 118, 115, 138, 138, 212, 250, 309, 230, 212, 258, 170, 191, 143, 219, 195, 208, 228, 302, 262, 278, 295, 290, 287, 325, 234, 348, 266, 262, 488, 683, 587, 493, 362, 377, 352, 298, 334, 292, 246, 270, 295, 544, 332, 368, 336, 335, 301, 351, 414, 468, 408, 389, 453, 359, 339, 546, 401, 763, 478, 323, 449, 400, 331, 353, 365, 383, 415, 401, 503, 620, 430, 471, 525, 578, 431, 465, 470, 472, 444, 439, 410, 381, 448, 508, 465, 422, 433, 410, 404, 389, 483, 467, 485, 388, 553, 693, 796, 874, 809, 793, 783, 825, 779, 931, 1004, 1179, 1302, 1217, 871, 882, 649, 781, 868, 769, 638, 600, 566, 548, 592, 529, 662, 582, 607, 573, 572, 640, 558, 633, 561, 514, 467, 447, 493, 516, 593, 497, 613, 438, 529, 541, 588, 465, 458, 447, 513, 526, 802, 595, 595, 498, 464, 434, 368, 419, 527, 661, 583, 629, 535, 517, 475, 563, 529, 516, 520, 436, 424, 544, 561, 519, 552, 521, 464, 439, 508, 481, 624, 559, 572, 524, 449, 528, 519, 638, 565, 530, 530, 461, 486, 522, 507, 662, 506, 424, 428, 519, 570, 633, 660, 663, 561, 519, 696, 938, 694, 686, 680, 603, 556, 654, 761, 666, 717, 804, 672, 671, 740, 767, 678, 633, 624, 656, 719, 969, 714, 635, 817, 671, 635, 610, 713, 642, 897, 889, 676, 764, 726, 787, 760, 711, 619, 599, 648, 709, 746, 705, 790, 767, 668, 593, 658, 853, 782, 729, 622, 760, 767, 663, 825, 871, 852, 798, 745, 877, 785, 811, 813, 750, 660, 791, 697, 759, 717, 700, 766, 721, 785, 1596, 1329, 1203, 1979, 1512, 1099, 1007, 792, 928, 1112, 929, 836, 849, 1016, 857, 938, 992, 1165, 805, 734, 722, 664, 751, 911, 961, 907, 962, 1126, 1025, 993, 1117, 1041, 1037, 982, 998, 977, 893, 1154, 1000, 1039, 1031, 877, 861, 1062, 1106, 1153, 1305, 1013, 1272, 778, 614, 800, 926, 910, 1029, 922, 825, 765, 975, 955, 903, 918, 946, 1076, 1043, 1029, 1092, 1159, 1021, 949, 940, 948, 1038, 1047, 1151, 883, 885, 1174, 1314, 1200, 1123, 1160, 1076, 1070, 1052, 1137, 1205, 1086, 978, 1073, 953, 890, 1149, 1170, 1390, 1380, 1352, 1327, 1406, 1494, 1735, 1438, 1416, 1186, 1934, 1292, 1059, 1158, 1131, 1040, 1106, 8959, 6408, 3697, 3305, 3835, 2310, 1715, 1419, 1287, 1310, 1271, 1069, 1123, 1125, 1110, 820, 1291, 1072, 1009, 880, 811, 727, 926, 992, 1023, 922, 989, 969, 911, 918, 960, 1097, 1136, 1315, 1323, 1314, 1282, 2353, 3904, 4085, 2725, 2012, 1910, 1644, 1369, 1555, 1344, 1297, 1174, 1114, 1174, 1127, 1048, 1203, 1025, 1021, 1038, 1047, 1005, 1008, 1078, 852, 924, 965, 873, 933, 939, 929, 901, 900, 879, 876, 864, 898, 894, 931, 942, 930, 1014, 989, 1049, 1062, 1026, 1055, 1063, 1094, 1024, 1045, 1164, 1263, 1240, 1204, 1237, 1213, 1210, 1196, 1136, 1117, 1292, 1109, 1108, 1364, 1231, 1024, 1142, 1054, 944, 954, 1282, 1188, 985, 971, 1017, 975, 975, 1014, 953, 953, 943, 990, 1080, 983, 946, 1056, 867, 945, 909, 904, 811, 771, 851, 902, 882, 792, 862, 722, 795, 862, 942, 791, 894, 845, 820, 812, 924, 811, 852, 905, 980, 828, 824, 920, 940, 1273, 1273, 970, 792, 800, 903, 1108, 1125, 969, 956, 801, 845, 875, 775, 797, 776, 773, 713, 949, 821, 800, 829, 845, 864, 853, 780, 860, 982, 1342, 1212, 916, 1004, 895, 941, 954, 841, 802, 820, 808, 931, 895, 907, 890, 837, 1022, 1010, 1616, 1067, 914, 982, 1100, 1007, 1089, 1098, 997, 988, 1014, 970, 897, 1034, 966, 812, 881, 805, 844, 847, 865, 909, 875, 830, 856, 823, 863, 879, 885, 835, 802, 847, 832, 881, 822, 1013, 779, 731, 794, 812, 823, 793, 810, 854, 870, 814, 840, 810, 1040, 920, 809, 810, 839, 803, 879, 1033, 1009, 915, 882, 878, 797, 796, 861, 867, 821, 917, 710, 689, 861, 986, 935, 862, 824, 840, 852, 816, 857, 971, 918, 957, 827, 1186, 889, 974, 1003, 871, 797, 762, 918, 898, 822, 903, 811, 882, 819, 872, 823, 894, 963, 868, 1052, 1011, 965, 856, 738, 704, 799, 945, 989, 921, 988, 858, 786, 968, 1009, 977, 947, 918, 850, 810, 703, 775, 1108, 1313, 1120, 1018, 1074, 1011, 1094, 1131, 1341, 1385, 1376, 1513, 1414, 1250, 1209, 1523, 1391, 1498, 1406, 1132, 1246, 1030, 1053, 1001, 909, 986, 978, 940, 974, 952, 999, 1027, 1252, 1014, 949, 926, 1013, 937, 1065, 1054, 989, 861, 847, 918, 1009, 940, 967, 973, 1023, 1090, 926, 839, 884, 959, 918, 919, 1080, 1312, 1206, 1174, 1177, 1088, 1191, 1045, 990, 945, 915, 983, 945, 966, 924, 886, 890, 823, 930, 803, 827, 2474, 913, 870, 884, 886, 854, 874, 833, 838, 920, 870, 869, 817, 821, 827, 751, 912, 905, 843, 751, 1230, 960, 858, 850, 885, 798, 771, 856, 815, 741, 757, 841, 693, 768, 720, 779, 844, 706, 702, 714, 679, 698, 733, 660, 668, 680, 763, 695, 687, 708, 660, 648, 624, 726, 806, 694, 697, 717, 754, 605, 612, 681, 655, 658, 691, 729, 1048, 836, 714, 723, 646, 681, 784, 698, 709, 620, 708, 746, 778, 725, 674, 651, 746, 725, 733, 672, 692, 674, 710, 724, 706, 653, 653, 688, 620, 673, 961, 855, 791, 901, 981, 805, 749, 780, 624, 619, 681, 607, 574, 595, 549, 635, 878, 871, 742, 465, 561, 519, 591, 548, 485, 513, 447, 538, 510, 479, 557, 554, 487, 461, 422, 533, 519, 521, 572, 516, 521, 530, 453, 484, 467, 499, 464, 384, 415, 467, 432, 435, 445, 888, 477, 503, 532, 491, 508, 549, 486, 424, 461, 469, 507, 554, 465, 474, 493, 475, 537, 687, 497, 585, 840, 704, 782, 814, 702, 777, 634, 739, 625, 775, 751, 652, 741, 684, 730, 775, 797, 700, 698, 672, 807, 745, 730, 866, 836, 742, 781, 798, 714, 765, 753, 706, 704, 786, 679, 775, 804, 838, 730, 725, 844, 724, 693, 708, 760, 1063, 1315, 806, 828, 822, 1258, 1082, 829, 666, 801, 814, 738, 744, 772, 1198, 5003, 5211, 3755, 5442, 4952, 5134, 4512, 3682, 2673, 2284, 2055, 1950, 2161, 1790, 1806, 1615, 1546, 1411, 1435, 1497, 1375, 1300, 1238, 1421, 1331, 1217, 1241, 1411, 1269, 1230, 1258, 1150, 1216, 1144, 1064, 1299, 1185, 1125, 1092, 1016, 1084, 918, 743, 857, 1012, 1618, 2562, 1168, 1082, 976, 992, 1055, 1000, 1053, 1008, 1115, 979, 1031, 1007, 977, 1112, 6407, 1799, 1097, 1169, 919, 1020, 1065, 1028, 1013, 1029, 1030, 1042, 1048, 1132, 1151, 967, 953, 977, 989, 1145, 1174, 1168, 964, 958, 996, 949, 1002, 906, 1043, 925, 866, 891, 809, 817, 933, 917, 807, 804, 887, 804, 954, 1109, 795, 685, 632, 590, 558, 574, 543, 855, 859, 887, 838, 843, 879, 957, 895, 871, 799, 857, 842, 869, 889, 906, 824, 892, 753, 797, 768, 789, 842, 851, 769, 867, 789, 737, 777, 706, 708, 703, 768, 772, 759, 890, 993, 764, 768, 854, 752, 844, 692, 763, 739, 813, 733, 763, 770, 683, 804, 771, 705, 688, 718, 836, 779, 822, 835, 831, 761, 698, 768, 716, 779, 742, 717, 829, 713, 749, 705, 821, 718, 742, 691, 830, 705, 813, 1089, 918, 784, 760, 905, 823, 831, 641, 683, 758, 689, 880, 741, 724, 861, 889, 714, 701, 654, 685, 638, 653, 689, 732, 602, 736, 611, 574, 615, 588, 598, 555, 509, 664, 590, 577, 678, 668, 639, 716, 632, 648, 590, 694, 659, 645, 550, 590, 524, 560, 617, 627, 546, 630, 590, 722, 825, 753, 738, 644, 628, 1194, 760, 650, 620, 567, 664, 745, 636, 550, 573, 661, 581, 649, 582, 607, 524, 602, 646, 616, 596, 593, 566, 577, 561, 641, 536, 525, 557, 545, 563, 572, 626, 599, 619, 657, 1794, 852, 707, 939, 1097, 788, 828, 3762, 2034, 1264, 1987, 2295, 1515, 958, 1925, 966, 847, 1063, 2028, 1728, 1207, 1077, 1035, 1130, 976, 795, 862, 741, 738, 756, 771, 712, 1427, 3530, 2021, 1118, 887, 1135, 858, 826, 758, 757, 788, 725, 825, 782, 873, 695, 962, 844, 1014, 859, 737, 758, 799, 849, 859, 749, 789, 827, 1030, 781, 834, 837, 777, 875, 773, 772, 709, 703, 4037, 1469, 1084, 1165, 784, 800, 760, 739, 801, 799, 767, 779, 1057, 952, 775, 832, 1013, 801, 786, 869, 731, 762, 741, 873, 738, 688, 856, 692, 639, 626, 772, 667, 649, 749, 706, 973, 829, 764, 707, 658, 609, 723, 701, 736, 693, 630, 682, 611, 645, 658, 683, 693, 650, 600, 579, 559, 590, 579, 589, 661, 684, 540, 784, 651, 675, 667, 613, 716, 738, 802, 770, 706, 686, 700, 647, 737, 670, 747, 834, 673, 713, 664, 721, 664, 703, 714]
In [6]:
ny_d0 = 1238198400000
In [7]:
gangnam_style = [520578, 429666, 397952, 505910, 405054, 256381, 649459, 535209, 524792, 462833, 565211, 687343, 636699, 603342, 813458, 832641, 924225, 970433, 1058284, 1789240, 1467022, 1361674, 1580239, 1585645, 1649960, 1581240, 1763891, 1505520, 1341269, 1544174, 1874544, 2031506, 2077288, 2214882, 2111487, 2087280, 2300856, 2672590, 3057718, 3212088, 3594824, 3607928, 3295138, 3589528, 3929102, 4219268, 4768364, 5018013, 5357443, 5019866, 5173064, 5444334, 5581247, 6202701, 7006687, 7032648, 6919060, 6719197, 6842544, 7212990, 7910952, 9664447, 9108651, 8820094, 9146978, 8418257, 8473949, 9550893, 11018902, 11736860, 10416584, 9738119, 9483200, 9321228, 8909267, 11481444, 12784802, 10849805, 9454005, 8803530, 8835570, 9069084, 11005339, 11615236, 9848207, 8449565, 8017559, 8126286, 8133914, 9467998, 10859777, 9832122, 7273529, 7614319, 7739443, 7633689, 9214745, 10957737, 9175074, 7587595, 7764514, 8484355, 8414463, 9836827, 11189149, 9377162, 7906863, 8065552, 8051642, 8601761, 9574078, 10454058, 8760752, 6828075, 6940728, 6758396, 7475981, 9105561, 11206549, 9367574, 7667185, 7120974, 7332373, 7178476, 8518937, 10511076, 8755494, 6932099, 7023241, 7389004, 7887420, 8565469, 11487589, 10209032, 7512689, 6682767, 6494478, 6549704, 7957156, 9311429, 7476852, 5808087, 5671298, 5951764, 6036078, 7156204, 8888429, 6897801, 5419944, 5476154, 5673707, 6040946, 7060284, 8138955, 6665434, 5369141, 8624024, 6962725, 8956125, 14924298, 11541865, 9492325, 9921559, 9019929, 7682540, 7052114, 7316434, 7819302, 7296451, 11085546, 7287532, 6389752, 5914194, 6215161, 7246611, 5928877, 5080581, 4849609, 5043994, 5336680, 5874136, 6878945, 5797948, 4385026, 4584479, 4541108, 4655647, 5267185, 6306349, 5376396, 4308772, 3681413, 4489439, 4271413, 5032795, 5953710, 5049273, 3962326, 4012856, 4048600, 4162806, 4786452, 5786704, 4960040, 3794802, 3764029, 3957399, 4014468, 4707340, 5617656, 4848504, 3889920, 3797241, 3733985, 3762549, 4454666, 5103411, 4363097, 3582828, 3470084, 3544961, 3638580, 4401248, 5426250, 4224641, 3164109, 3468998, 3306348, 3157505, 3746716, 5144933, 4329939, 3242996, 3108167, 3206673, 3408942, 3726897, 4575528, 3766935, 3022835, 3021680, 3015885, 3118800, 3735088, 4648755, 3934828, 2959366, 2901021, 2873277, 2889656, 3390074, 4060047, 3522396, 2772557, 2779615, 2880541, 2942589, 3232238, 3370352, 3114496, 2736286, 2486267, 2574656, 2719661, 2746360, 3346764, 2809162, 2240907, 2218082, 2329051, 2445492, 2926574, 3990257, 4220378, 3631999, 3322043, 2839298, 2688346, 3550368, 3893333, 2973921, 2137509, 2307979, 2349896, 2240535, 2501047, 2853962, 2509205, 2009597, 2143262, 2423911, 2316451, 2310104, 2656082, 2238667, 1823209, 1759979, 1847124, 1946184, 2144452, 2439730, 2003199, 1606718, 1598226, 1661170, 1667831, 1899119, 2238454, 2116519, 1752631, 1645825, 1642035, 1742668, 1989617, 2290706, 2013891, 1692240, 1625839, 1669896, 1740956, 1834506, 2139404, 1856432, 1614242, 1539186, 1503506, 1565075, 1729917, 1887569, 1660447, 1394235, 1444947, 1416540, 1390298, 1505303, 1669974, 1479979, 1287720, 1299237, 1276124, 1300069, 1381207, 1522580, 1305555, 1215816, 1235000, 1272230, 1374422, 1394487, 1529056, 1346645, 1260346, 1252824, 1255996, 1295045, 1403945, 1477193, 2191564, 1395833, 1372207, 2463147, 1341864, 1927184, 1650021, 2623358, 2947237, 1684224, 1189390, 1213568, 1326633, 1329957, 1171079, 1127822, 1124372, 1113470, 1100733, 1205090, 1270241, 1127974, 1094596, 1087548, 1053101, 1054306, 1214791, 1217006, 1090662, 1034406, 1039773, 1047610, 1068348, 1121260, 1174848, 1085700, 1033750, 1024127, 1062159, 1037714, 1082550, 1125624, 1042973, 1009021, 1076690, 1062314, 1069452, 1131626, 1201404, 1308454, 939452, 951912, 940037, 959168, 1080650, 1154721, 1055020, 919882, 880427, 905276, 896710, 998658, 1163679, 984427, 840402, 833984, 848721, 837183, 975047, 1205779, 1008897, 767616, 771686, 756693, 780916, 929299, 1073402, 912294, 709954, 705218, 705834, 718649, 850786, 1014891, 892357, 720660, 725652, 748325, 741509, 868118, 1066284, 905775, 719367, 735545, 803910, 803072, 943124, 1105894, 975401, 818138, 791833, 812249, 835844, 961647, 1109471, 935630, 768193, 793144, 818822, 829571, 962540, 1124023, 1005348, 844987, 839410, 785031, 873256, 1025542, 1168894, 1053445, 812220, 796252, 821785, 811054, 947204, 1111561, 1034968, 853981, 810984, 823659, 819858, 938461, 1152921, 969425, 763053, 772884, 799418, 820273, 971502, 1152222, 966040, 786763, 810022, 863079, 915460, 999834, 1155809, 970919, 792316, 801611, 854061, 860599, 969760, 1147332, 970812, 785217, 776486, 860520, 795485, 962033, 1146396, 959024, 767085, 766319, 800528, 854124, 935753, 1070980, 974719, 932175, 1135712, 1034452, 983544, 1062462, 1099100, 1042753, 1056354, 1687223, 1085504, 1021200, 1059459, 1114790, 1001356, 883433, 891419, 876013, 873431, 974177, 1136040, 970624, 792575, 806275, 845833, 821561, 940386, 1131729, 973537, 816732, 815515, 827487, 821041, 940088, 1095140, 959687, 779763, 804907, 833333, 841005, 912807, 1072332, 922285, 742791, 740753, 808256, 794750, 895617]
In [8]:
gangnam_style_d0 = 1339718400000
In [9]:
flashdance_d0 = 1265155200000
In [10]:
flashdance = [9, 27, 19, 31, 60, 204, 317, 354, 318, 338, 3318, 5105, 5235, 5483, 5426, 5275, 6459, 7322, 6347, 5018, 6322, 6133, 6690, 10798, 15256, 14565, 7026, 6671, 7309, 7253, 7409, 8210, 6862, 7238, 7229, 7332, 6998, 8590, 5174, 4660, 3974, 4212, 3857, 3584, 4168, 5403, 4931, 4443, 4757, 4542, 4471, 4991, 5032, 4303, 3804, 6213, 6859, 5191, 5435, 5526, 8448, 8280, 6787, 6366, 6086, 6848, 8112, 7644, 7215, 6839, 7018, 6461, 7461, 7431, 6885, 6020, 6326, 6513, 6953, 8284, 8906, 9247, 7864, 7776, 8309, 8347, 9118, 10026, 9438, 8520, 8464, 8546, 8610, 9293, 10601, 9506, 8956, 9760, 9784, 10418, 10354, 13088, 10628, 9553, 9764, 9822, 10509, 10380, 11209, 12451, 11258, 9973, 10212, 9605, 10415, 10886, 9895, 8985, 9224, 8900, 8893, 8797, 9823, 8957, 8206, 8481, 8579, 8664, 8505, 9236, 9045, 8411, 8653, 8792, 8704, 9481, 10548, 9930, 9038, 8200, 8098, 8170, 8790, 9548, 8576, 8553, 8286, 8712, 8901, 9114, 8876, 8818, 8675, 8346, 7837, 8538, 9174, 8676, 7617, 8036, 8590, 7511, 8475, 9193, 9949, 8677, 7874, 8353, 8632, 11880, 13177, 10659, 9213, 8567, 8420, 8352, 8831, 9358, 9432, 8456, 8419, 8343, 8213, 8857, 8739, 8600, 7967, 7311, 7631, 8222, 7924, 8829, 9163, 8210, 7945, 7783, 8134, 8108, 9222, 9245, 8991, 7968, 8199, 8279, 8314, 10021, 10070, 9778, 8282, 7976, 8850, 8754, 9649, 11149, 9524, 8583, 8032, 8046, 11667, 10065, 10409, 8762, 7668, 8450, 8342, 8518, 9288, 10394, 9055, 8396, 8649, 8984, 8717, 9763, 11121, 9471, 8512, 8494, 8595, 8711, 10402, 11627, 9392, 8188, 8840, 8980, 9249, 10453, 11274, 10021, 9110, 9310, 8752, 10021, 13443, 12927, 11454, 9607, 9369, 10142, 10473, 11642, 12451, 11118, 9668, 9709, 10022, 9922, 11139, 11913, 10066, 10005, 9664, 9851, 9943, 10880, 12457, 11089, 9420, 9924, 10682, 11377, 12197, 13005, 11029, 9831, 9582, 9847, 9631, 10430, 12226, 11201, 10393, 9591, 9957, 10098, 11470, 13055, 11681, 9453, 9422, 9667, 9790, 11020, 11618, 10129, 9182, 9278, 9252, 9663, 10899, 13396, 11244, 9598, 9861, 10032, 9350, 10548, 11056, 9941, 8983, 9471, 9529, 9235, 8485, 8368, 9035, 8932, 9180, 9730, 10066, 12593, 9245, 9113, 8328, 8445, 9452, 9341, 10371, 11666, 10748, 8904, 9905, 10193, 10256, 10838, 12810, 10644, 9088, 8982, 6800, 9426, 11064, 12340, 10888, 9314, 9529, 9733, 9882, 11184, 12244, 11124, 9375, 9504, 9301, 10708, 12071, 13399, 12734, 11086, 11158, 11206, 12475, 12532, 14610, 12615, 10202, 10980, 10890, 11431, 12472, 15210, 15018, 11958, 11564, 11832, 11933, 13564, 17231, 15118, 13115, 12492, 12228, 12575, 14482, 15815, 13766, 12283, 11943, 12032, 12086, 12809, 15309, 14505, 13285, 12065, 12010, 11895, 14610, 17113, 13830, 11684, 11088, 11553, 12213, 12449, 16061, 12407, 11441, 11915, 11781, 11996, 12459, 13475, 11830, 10529, 11087, 12081, 12168, 12488, 13832, 11996, 11569, 12028, 11834, 12417, 13006, 16373, 13953, 11772, 11534, 11260, 12407, 12445, 12297, 11388, 10404, 10336, 10420, 10418, 10613, 11785, 10921, 10023, 10479, 10549, 10868, 11578, 12210, 10665, 10113, 10598, 11225, 11416, 12333, 14456, 13101, 11202, 11971, 11353, 12620, 13243, 14332, 12051, 11162, 11648, 11643, 11639, 13653, 14451, 12967, 13159, 13521, 12567, 13358, 13143, 16074, 13527, 11996, 11838, 11940, 11720, 12961, 14413, 12000, 11320, 12654, 10966, 10864, 12866, 13070, 12505, 12410, 13010, 13211, 12766, 13521, 13849, 11311, 10283, 9992, 10509, 11425, 11747, 12042, 10069, 9634, 9524, 9855, 9623, 10791, 10623, 9763, 9575, 10224, 11361, 11576, 12383, 12014, 11966, 11210, 11287, 10493, 11026, 11247, 11396, 10195, 9857, 10506, 10657, 10540, 10839, 10802, 9572, 9192, 9015, 10169, 10019, 10690, 10674, 9557, 10021, 10063, 9644, 9721, 10604, 10975, 11052, 9995, 9378, 9066, 9405, 10631, 10971, 9659, 9282, 9635, 9938, 10471, 11830, 13298, 12124, 10046, 10003, 9589, 9464, 10353, 11354, 10312, 9262, 9002, 9346, 9256, 10286, 11529, 10716, 9516, 9752, 10305, 10395, 12252, 13647, 11999, 10455, 10381, 10072, 9804, 11657, 12182, 10488, 8897, 8892, 9511, 9569, 10961, 12122, 10987, 9199, 9437, 9671, 9630, 11031, 11735, 10670, 9291, 9215, 9239, 9329, 11455, 12233, 11286, 9944, 8959, 10014, 9970, 11563, 13540, 11723, 10728, 15486, 15022, 13828, 16029, 16133, 12983, 11294, 11581, 11890, 12467, 13288, 14137, 12940, 10712, 10879, 11080, 11697, 12885, 13494, 11498, 10056, 10120, 11089, 11484, 11839, 12913, 11726, 9503, 9828, 10056, 10187, 10912, 12057, 10133, 9110, 9330, 9825, 9544, 10088, 11437, 9665, 8513, 8768, 9142, 9389, 9931, 10844, 9964, 9027, 9396, 9098, 9682, 10422, 11814, 10689, 11602, 12013, 10048, 9655, 9360, 9107, 9169, 9275, 9481, 10296, 10573, 10908, 13678, 9158, 9110, 9136, 12910, 18714, 16053, 14957, 12638, 12253, 12839, 13046, 12509, 13105, 13631, 14636, 11514, 11117, 11173, 11997, 12738, 16605, 14953, 12743, 11000, 10862, 11471, 12575, 13540, 11545, 9959, 9744, 10157, 10693, 11300, 12343, 10648, 9526, 9377, 10121, 10120, 11388, 12713, 11817, 10086, 10285, 15521, 14347, 13305, 13065, 11244, 10635, 10897, 11088, 11387, 13548, 15435, 12435, 11133, 11182, 11413, 11557, 13003, 13866, 12322, 10660, 10479, 10682, 11323, 12727, 14617, 13066, 11016, 11080, 11082, 10920, 11687, 15592, 12761, 10132, 9763, 9817, 9951, 11374, 12090, 10170, 9754, 9601, 9484, 10041, 11281, 11809, 10420, 9336, 9720, 15282, 12675, 12381, 12227, 10986, 10497, 10103, 9992, 9990, 10986, 12008, 10541, 9440, 10719, 10681, 10956, 11846, 13015, 12912, 10675, 10946, 10793, 10924, 11868, 13114, 11789, 10958, 11263, 10025, 10260, 11281, 12787, 11268, 10374, 11912, 11922, 12358, 13572, 13652, 11484, 10602, 12687, 13744, 16818, 17023, 16073, 13373, 11750, 11733, 11795, 11636, 12362, 12573, 12131, 11705, 10914, 11194, 11786, 13846, 15280, 13244, 11595, 11893, 12583, 12459, 13508, 14389, 12891, 10842, 11759, 11833, 11697, 12635, 13000, 11253, 10055, 10284, 11272, 11058, 11567, 12558, 12197, 10985, 10703, 10114, 10282, 11674, 12139, 10867, 10570, 10463, 10104, 10124, 11011, 11532, 10510, 9762, 9796, 9732, 10179, 10582, 11096, 9905, 9146, 8978, 8886, 8863, 9811, 7901, 9347, 9095, 9756, 9959, 9493, 10125, 8740, 8261, 7916, 7755, 7670, 7666, 8164, 8543, 8280, 8244, 7854, 7926, 8436, 8430, 8740, 8467, 8068, 8991, 9214, 8692, 8500, 10895, 9283, 8398, 8075, 8294, 8859, 9336, 9689, 9371, 8087, 8241, 8190, 8424, 9132, 10034, 8971, 8748, 7969, 8100, 8266, 8870, 8919, 8991, 7319, 7211, 7319, 7184, 7180, 7400, 7367, 5983, 6025, 5866, 6344, 7537, 7977, 6558, 5755, 5851, 6225, 6358, 6856, 7507, 7301, 7672, 7840, 7751, 7952, 9092, 10246, 8681, 9947, 8523, 8005, 8161, 9660, 9989, 8413, 7937, 7397, 7837, 8134, 9501, 10211, 8532, 7617, 7782, 8525, 8412, 9522, 10405, 9087, 7981, 8170, 8034, 9204, 11475, 11612, 10639, 8917, 8666, 8941, 9018, 9850, 12025, 10476, 9178, 8629, 9168, 8922, 9246, 10061, 8570, 7989, 8272, 8424, 9230, 10874, 15204, 11355, 8575, 8328, 8395, 8166, 9058, 9370, 8751, 7742, 7883, 8326, 8312, 9979, 10830, 11167, 9314, 9874, 9568, 8934, 7996, 9560, 7119, 6067, 6625, 6830, 6600, 7132, 9043, 8321, 7476, 6995, 7163, 6992, 7619, 8898, 9545, 12634, 8008, 7044, 7315, 7865, 8454, 7952, 7192, 7185, 7453, 7564, 8141, 9159, 7998, 6910, 6913, 7062, 7409, 8401, 9369, 8276, 7592, 7845, 8011, 8258, 9550, 10492, 9275, 8735, 8330, 7989, 8134, 9243, 9932, 8787, 7860, 7823, 8087, 8294, 9049, 9888, 8772, 7080, 7151, 7178, 7533, 7648, 9448, 7818, 7395, 7117, 7127, 7597, 8598, 9669, 8634, 7855, 8558, 8267, 8300, 9256, 9795, 7882, 6622, 6705, 7325, 7295, 8583, 13763, 10994, 8121, 8207, 8255, 8143, 10040, 11077, 9634, 8952, 7995, 7948, 8249, 8431, 9253, 9483, 7279, 7951, 7944, 8078, 8220, 8995, 8565, 8166, 7313, 7323, 7712, 7612, 8803, 7647, 6950, 7597, 8463, 8736, 10214, 12139, 11029, 9606, 8878, 9353, 9043, 9185, 9875, 9170, 8044, 8095, 8331, 8922, 9478, 10663, 9334, 8064, 8465, 9112, 7528, 9477, 10001, 8233, 7345, 8185, 8553, 8798, 9471, 10087, 8576, 7238, 8008, 8234, 7613, 8629, 9458, 7872, 7927, 7687, 7817, 8153, 9229, 10005, 8536, 8082, 8895, 8740, 8537, 9263, 10672, 9192, 8510, 8330, 8145, 8070, 8632, 9321, 8544, 7515, 7630, 7429, 8104, 8228, 8804, 7619, 6887, 7142, 7070, 7236, 7683, 8373, 7676, 6914, 7099, 7707, 8771, 9175, 9414, 8185, 7731, 7925, 9329, 9224, 9768, 9567, 8424, 7788, 7821, 8062, 8255, 9011, 9295, 8218, 7540, 7428, 7627, 7906, 8309, 8686, 7843, 7533, 7645, 8044, 7923, 8738, 9508, 8285, 8247, 8305, 8936, 8452, 9221, 9355, 9046, 8331, 8640, 8597, 9024, 9700, 10513, 9141, 8663, 8801, 8920, 8985, 9439, 10134, 9139, 9475, 9334, 8982, 9459, 10234, 11630, 10868, 8955, 9477, 9190, 10670, 11262, 11146, 10018, 8936, 9916, 9299, 9204, 9742,11132, 10317, 9496, 10963, 10040, 10072, 11188, 11660, 10581, 9040, 9432, 9017, 9472, 9893, 12266, 11310, 9354, 9591, 9423, 9703, 10866, 12101, 9997, 8283, 8263, 8828, 9095, 10098, 11296, 9963, 8855, 9167, 9799, 9806, 10494, 12413, 11298, 10043, 9204, 10645, 10259, 11690, 12432, 11637, 9943, 10124, 10302, 10488, 11328, 12501, 10969, 9676, 9861, 9354, 9710, 11278, 11982, 10319, 9529, 9133, 9237, 9659, 11338, 12138, 10867, 9423, 9506, 9000, 9729, 11036, 11877, 10371, 9242, 9311, 9698, 9764, 10196, 10756, 10248, 9414, 8830, 8979, 9282, 9756, 10943, 12915, 8427, 8631, 8380, 8449, 9286, 10202, 9084, 8941, 8109, 7756, 8800, 9443, 11025, 9977, 8131, 7403, 7762, 9032, 9410, 9285, 9408, 7951, 8321, 7874, 8530, 8770, 10097, 9692, 10045, 13777, 8388, 8504, 9473, 10109, 8773, 7469, 8006, 8979, 8516, 11349, 12440, 9644, 8648, 8406, 11043, 11522, 9797, 11849, 9776, 9245, 9502, 9573, 10340]
In [11]:
mj = [1, 1, 0, 3, 3, 4, 2, 3, 2, 1, 0, 8, 3, 0, 2, 1, 0, 0, 3, 2, 1, 1, 0, 1, 1, 1, 3, 2, 0, 6, 4, 12, 3485, 13884, 17711, 13614, 16923, 21907, 16639, 10998, 11137, 12429, 15204, 21099, 25443, 20448, 19432, 19192, 20851, 21148, 18898, 17503, 16963, 16732, 18130, 21634, 29555, 21493, 20257, 15601, 11529, 11203, 12122, 15030, 18276, 22234, 16300, 10176, 10522, 11276, 12958, 17454, 21340, 17606, 9104, 6340, 6305, 8742, 17208, 23277, 16640, 15045, 13436, 13108, 14749, 18994, 26205, 20598, 16159, 16018,
16313, 17430, 24328, 32169, 24847, 15779, 19034, 17414, 17804, 23301, 30526, 24199, 21375, 17585, 17061, 1980, 2222, 3152, 2489, 1734, 1602, 1654, 1967, 2637, 3489, 2647, 1879, 1933, 10148, 14778, 19442, 27692, 19676, 13208, 14246, 13470, 14764, 19972, 26723, 18407, 13773, 13975, 14039, 12990, 16372, 20710, 17270, 12135, 11521, 12068, 16789, 20399, 21144, 17659, 13957, 14472, 13643, 15135, 18429, 18207, 14417, 12641, 11517, 11844, 12781, 15533, 21318, 16640, 13831, 14263, 13065, 12150, 14876, 22869, 17628, 12196, 12525, 13742, 14183, 17944, 22918, 17796, 12083, 13775, 13599, 14427, 17509, 23090, 18058, 13328, 14004, 15150, 16909, 20999, 27405, 20125, 14083, 14202, 15045, 14822, 14681, 18691, 16149, 12184, 17079, 17440, 18131, 21445, 23368, 18406, 13809, 13966, 15686, 18950, 23556, 27284, 23988, 21406, 18221, 18490, 18693, 20123, 26953, 23462, 18810, 19938, 19583, 20698, 20848, 24598, 22192, 20131, 17962, 17921, 20659, 23282, 25247, 23241, 22674, 24553, 23718, 28584, 59251, 44241, 33749, 29994, 24683, 23628,
25605, 25141, 24540, 21784, 21743, 21727, 21540, 21885, 23122, 23793, 21432, 21755, 20767, 19087, 20218, 21275, 21986, 18784, 20193, 20529, 21499, 40308, 29451, 40029, 35258, 30579, 33174, 34227, 29228, 25461, 35877, 29900, 27971, 18269, 31345, 27753, 30571, 34381, 31012, 25844, 29338, 27550, 17588, 19645, 20625, 17954, 16002, 16852, 17245, 18363, 18354, 20075, 17959, 16114, 15210, 15366, 15544, 17551, 21810, 20718, 15220, 13895, 13483, 14977, 18067, 19709, 17850, 15847, 12285, 12672, 14894, 18555, 24626, 20239, 13874, 14247, 16000, 16771, 20818, 26004, 20894, 14001, 14639, 15029, 17076, 21731, 23307, 20671, 15478, 15300, 17599, 17394, 19304, 25621, 23034, 15088, 15226, 17639, 18858, 24248, 28612, 23963, 19767, 17225, 14643, 21768, 27006, 32884, 25658, 21344, 23364, 27336, 29079, 40459, 47549, 38252, 33685, 33025, 39208, 54066, 103990, 126484, 149206, 36993, 26563, 22641, 22981, 27266, 33648, 27108, 18370, 18563, 19382, 23333, 25199, 30705, 23435, 16025, 18616, 19423, 17846, 21793, 28716, 24047, 18476, 21442,
26695, 31199, 26916, 27275, 20671, 15368, 15826, 17710, 31443, 22950, 28741, 22070, 16614, 15844, 18587, 21073, 25901, 35056, 28427, 21760, 21481, 21662, 23081, 25540, 31514, 25281, 22704, 22875, 24785, 26485, 27061, 23373, 21990, 23386, 23939, 24698, 26814, 35349, 32391, 25467, 20291, 19864, 20688, 17648, 26547, 36601, 27058, 19679, 20368, 19950, 21838, 27307, 34926, 31142, 25831, 17693, 19916, 20089, 28168, 34778, 27000, 19087, 22670, 22538, 24786, 28417, 37446, 28795, 20454, 24213, 25022, 26653, 33195, 42072, 38545, 29716, 27904, 27286, 28175, 35505, 45351, 33939, 22774, 22350, 24921, 23889, 30431, 38190, 32655, 26934, 21334, 21035, 22605, 28461, 34375, 24970, 19009, 18813, 19646, 21068, 26343, 40516, 34240, 21972, 20426, 22334, 25900, 28976, 42155, 30475, 23888, 22603, 21148, 19796, 23231, 28629, 23063, 17164, 15939, 18599, 21774, 28512, 36914, 29007, 20430, 22080, 23569, 22674, 26016, 31664, 24518, 18992, 18540, 19111, 20749, 26384, 31586, 30018, 57401, 67782, 62210, 63195, 75029, 88190, 67054, 58051,
58808, 64424, 77571, 87612, 87913, 75854, 65428, 62402, 62877, 65986, 73257, 85215, 70609, 56749, 56873, 60037, 62218, 72932, 85299, 67529, 54348, 56093, 58338, 63007, 77801, 98163, 80096, 61639, 64725, 64471, 70518, 82824, 96338, 78233, 63570, 65042, 67213, 71397, 85946, 98406, 82708, 68295, 71349, 71537, 75622, 80182, 89693, 76922, 71468, 74979, 76252, 77149, 87369, 99764, 83182, 71840, 74112, 77958, 80391, 98874, 109236, 95292, 83439, 78121, 84589, 87547, 112798, 182020, 120286, 118117, 95728, 92760, 91512, 95161, 93366, 2286, 9383, 68284, 75233, 75077, 62907, 288, 5, 12349, 64895, 74623, 76884, 84790, 94288, 79696, 75480, 79374, 77401, 78524, 85481, 89687, 79149, 75216, 74211, 71878, 86347, 80631, 83752, 71971, 73175, 76648, 78678, 77380, 77505, 82619, 73378, 71855, 70332, 70595, 71699, 69379, 85110, 73706, 72229, 67868, 70178, 63301, 76632, 70447, 71596, 65372, 66064, 66878, 68264, 76609, 90461, 77653, 81532, 70807, 65735, 67427, 77601, 90478, 85257, 72175, 66187, 68965, 63865, 80077, 98087, 76906,
58557, 60339, 70226, 76062, 91235, 108087, 88491, 66488, 67381, 70568, 73266, 88467, 100656, 77766, 62323, 67501, 76202, 75723, 88898, 106207, 90198, 69086, 72204, 82057, 82377, 102427, 122711, 100510, 77297, 75512, 81456, 80456, 100173, 116726, 94136, 77460, 75243, 88183, 92346, 111171, 129617, 101712, 98567, 107389, 122701, 136391, 196627, 252805, 203350, 407464, 149768, 107500, 96787, 106543, 124088, 102453, 92182, 79905, 78129, 80873, 96813, 111954, 87476, 65377, 64321, 60432, 64524, 83136, 102114, 81152, 63323, 65945, 72927, 79941, 84031, 96088, 75161, 56864, 62041, 69453, 67387, 75365, 89974, 70662, 56080, 56046, 57393, 61238, 65962, 76758, 60365, 49703, 52964, 54890, 57226, 68155, 77381, 64046, 56366, 59039, 60598, 64330, 67291, 76478, 62120, 60958, 62190, 64756, 62332, 67223, 104084, 61898, 56526, 50715, 51878, 54289, 62776, 71015, 57386, 46645, 50047, 57956, 63866, 69184, 77819, 78259, 67349, 52907, 52710, 53568, 62780, 72207, 58788, 47222, 47322, 49167, 79520, 75759, 74653, 52916, 42118, 42175,
44444, 46251, 53279, 64275, 51946, 41558, 41824, 46790, 45610, 62605, 75369, 67169, 54003, 54457, 57903, 66598, 76020, 81652, 67658, 59421, 53483, 51204, 52476, 52781, 49190, 43876, 37151, 36234, 38048, 39868, 46731, 52114, 44689, 33415, 37717, 40247, 41864, 49411, 54434, 44825, 36339, 36917, 38985, 41791, 49743, 54286, 51653, 42885, 41513, 43746, 52912, 61061, 71274, 55962, 46655, 45230, 46742, 50803, 63273, 69756, 57020, 48989, 49064, 51910, 57782, 60240, 61983, 56080, 49474, 46586, 44566, 43078, 46478, 55089, 47702, 41237, 41350, 39750, 42712, 50704, 57560, 49884, 40198, 40294, 39103, 40689, 48627, 55426, 48402, 43932, 45091, 39201, 42473, 50915, 55532, 46841, 38811, 40357, 39119, 45045, 51752, 57892, 47954, 39697, 39789, 41532, 43647, 48611, 52336, 48052, 39804, 41875, 43389, 43725, 51306, 59378, 52207, 45928, 43602, 45375, 48171, 56256, 64341, 55234, 46886, 48371, 49509, 47218, 53587, 53249, 44158, 35839, 39948, 42145, 41500, 46020, 53017, 47240, 41856, 43210, 44176, 44051, 48405, 60263, 60712, 81114,
58069, 50789, 47838, 51063, 53878, 47514, 43924, 45538, 44455, 43852, 54871, 49704, 44788, 45602, 48889, 47720, 49435, 53287, 55174, 49743, 46298, 46529, 46634, 48583, 43925, 53721, 42972, 44593, 48261, 49693, 50737, 53261, 53570, 48106, 43552, 44495, 44627, 44860, 49619, 50798, 44935, 41576, 41819, 43412, 44300, 47658, 50009, 44032, 41503, 45609, 46181, 48074, 50462, 53907, 49471, 45510, 46007, 46157, 47252, 53272, 61430, 54585, 47224, 48668, 66468, 56474, 59026, 66456, 57807, 49971, 45153, 43563, 46019, 52086, 61123, 53985, 43991, 43938, 54577, 63194, 72676, 84680, 69920, 61781, 56443, 59544, 62161, 74652, 89149, 71615, 54595, 56152, 62412, 63624, 77699, 93380, 77699, 62883, 68293, 70265, 72806, 88498, 107165, 85814, 71237, 69358, 72789, 79249, 92069, 109620, 88606, 73652, 72941, 80388, 86215, 104627, 123399, 99706, 83579, 89301, 99084, 110804, 141850, 180545, 129063, 117025, 158023, 468998, 154947, 135004, 134321, 101834, 73915, 69014, 68925, 70547, 85047, 102542, 81098, 65094, 62820, 67491, 65892, 73530,
92532, 76338, 63487, 63685, 70439, 76775, 79248, 89178, 70226, 55720, 57836, 55550, 59209, 107004, 92700, 70944, 59416, 54751, 55003, 48866, 54420, 76249, 62494, 49863, 55870, 57161, 57726, 65791, 75773, 63564, 50208, 53177, 57937, 59294, 67648, 75676, 72224, 85789, 77104, 66424, 64753, 68889, 76332, 79081, 125846, 69927, 60152, 61191, 73821, 81495, 66688, 53270, 52514, 56130, 59816, 71843, 88717, 73827, 57173, 54853, 57706, 63668, 68517, 71324, 61370, 50712, 50411, 52530, 59823, 82862, 102346, 80630, 62386, 63725, 65084, 66176, 81850, 99832, 82834, 64098, 63651, 66367, 68767, 79947, 96176, 77001, 63292, 59576, 59388, 63810, 74183, 89874, 74772, 64949, 61101, 62193, 65873, 62524, 75256, 60568, 47822, 53797, 52634, 56058, 64596, 85005, 63565, 46705, 50949, 59486, 75571, 85084, 106951, 87628, 68213, 69395, 74285, 73645, 78993, 113004, 94012, 69426, 67953, 73304, 68162, 70832, 94042, 81921, 62787, 63933, 65838, 73524, 87734, 89997, 79034, 63559, 59771, 68185, 73024, 80171, 92280, 73081, 56323, 60141, 65481,
70746, 79507, 98581, 79904, 58267, 61085, 65311, 65532, 80431, 92767, 74550, 55942, 57585, 60887, 68033, 79404, 95721, 80880, 63129, 68237, 77773, 68654, 79931, 99849, 79959, 63684, 64125, 71540, 74570, 84001, 107553, 88439, 66271, 71204, 73893, 78320, 89626, 107447, 96667, 72852, 73214, 73094, 76689, 89969, 110041, 91026, 75946, 74040, 75616, 75881, 89980, 109071, 92093, 71909, 73941, 79363, 84303, 88662, 105551, 88998, 68002, 72266, 72391, 74821, 87036, 95610, 82854, 67875, 69292, 68453, 79952, 91844, 98194, 95336, 80593, 138750, 101592, 96330, 109544, 114056, 92107, 83928, 87457, 89180, 94768, 102388, 113371, 93825, 84501, 84385, 84618, 85125, 92976, 101860, 89694, 80875, 75821, 78167, 82174, 89826, 97342, 83539, 77392, 77578, 79677, 85934, 92681, 98093, 82973, 77607, 83530, 80352, 82322, 93741, 96755, 89445, 79035, 80962, 81677, 86795, 93639, 99616, 90115, 80559, 81352, 84422, 85581, 94084, 101635, 86379, 76320, 79058, 80043, 80854, 90508, 104129, 89681, 71704, 72322, 76427, 110910, 103073, 109505, 97758,
78419, 74392, 76983, 80199, 87834, 109041, 91921, 69288, 72775, 75823, 78164, 99917, 114888, 97605, 70754, 71978, 71400, 74032, 90019, 112417, 91486, 69556, 71051, 74574, 78055, 94001, 119125, 101282, 75231, 78228, 85510, 89069, 111705, 133673, 104720, 81453, 90981, 96396, 103398, 125892, 147559, 123677, 105059, 99548, 105119, 115716, 137562, 168867, 134765, 115472, 124386, 138602, 156848, 195127, 243984, 192341, 166765, 184674, 253809, 750403, 301988, 258883, 210095, 119303, 112451, 112541, 118896, 140470, 172767, 132778, 100829, 92738, 97590, 100068, 120299, 153104, 122347, 92352, 90115, 98221, 102805, 121894, 150466, 117080, 90284, 92906, 102111, 110684, 107602, 132547, 106744, 103330, 91935, 78685, 79913, 93038, 118363, 103357, 81748, 80523, 88073, 91921, 116627, 137915, 109397, 85771, 86714, 97928, 106425, 116741, 131268, 118496, 113253, 145706, 129998, 112025, 121826, 128285, 118141, 122200, 210764, 122733, 112554, 119241, 127501, 108910, 95798, 92548, 94580, 99801, 124608, 154362, 125078, 92656,
94464, 99522, 101454, 123085, 153796, 132949, 105149, 98081, 100137, 106850, 124602, 145628, 115731, 93033, 97626, 101547, 102305, 123840, 152432, 112648, 82245, 87790, 103838, 106761, 124169, 159409, 127303, 97966, 102912, 107975, 109081, 130608, 162148, 133751, 107647, 105736, 106733, 105497, 129935, 156676, 125495, 96858, 101507, 104736, 113212, 135700, 160555, 131530, 107590, 104519, 98901, 105411, 132055, 154629, 120278, 96742, 96113, 95245, 105169, 124174, 160073, 133896, 100198, 93924, 105215, 105400, 127264, 169279, 136750, 103826, 104864, 108976, 106478, 124892, 149093, 119067, 90531, 100052, 105902, 116322, 136127, 159161, 131242, 95534, 96898, 102070, 94665, 127499, 138827, 112800, 100746]
In [12]:
from collections import Callable
import matplotlib.pylab as pylab
#pylab.rcParams['figure.figsize'] = 8, 4
pylab.rcParams['figure.figsize'] = 16, 8
from matplotlib import pyplot as plt
from phoenix import models
import lmfit
import pandas as pd
import numpy as np
np.seterr(all='warn')
%matplotlib inline
In [13]:
def fit_and_plot(data, day_zero):
import datetime as dt
data = np.asanyarray(data)
model = models.WavePhoenixR(threshold=np.inf, score_func='mdl')
fit_time = int(data.shape[0])
params = model.fit(data)
for key in model.parameters:
print(key, model.parameters[key])
deltas = model(data.shape[0])
day = 60 * 60 * 24
first_date = (day_zero / 1000)
timestamps = [first_date + i * day for i in xrange(len(data))]
dates = [dt.datetime.fromtimestamp(ts) for ts in timestamps]
data_series = pd.Series(data, index=dates)
model_series = pd.Series(deltas, index=dates)
data_series.plot(style='wo', label='data')
model_series.plot(style='k-', label='PhoenixR')
plt.tight_layout(pad=0)
plt.legend(loc='upper left')
plt.show()
In [14]:
fit_and_plot(mj, 1)
In [15]:
fit_and_plot(queen, queen_d0)
In [16]:
fit_and_plot(flashdance, flashdance_d0)
In [17]:
fit_and_plot(ny, ny_d0)
In [18]:
fit_and_plot(gangnam_style, gangnam_style_d0)
In [19]:
fit_and_plot(beeges, beeges_d0)