YadryaraНаписал. Даже AVX не пригодился, просто обычный асм х64. Скорость ещё почти в 800 раз выше! Это уже со счётом лишь половины (но числа в таблице сразу удвоил). Все числа 23#-43# выдаются правильно, 43# считается 2 секунды, 47# минуту, 53# полчаса, запустил 59# в три потока (в четвёртом считается 47# на PARI для образца), к обеду досчитается всё. Пока покажу 47# и 53# (c19-c50):
Код:
47#: 0, 0, 0, 0, 0, 0, 0, 0, 0, 24, 1406, 34278, 520308, 5590512, 44192112, 255924622, 1087592188, 3438714568, 8217471808, 14998221970, 20942301992, 22275622040, 17925228472, 10827430718, 4874885124, 1630911948, 406659102, 75530368, 9994728, 810792, 27480, 0, sum=107017666560
53#: 0, 0, 0, 0, 0, 0, 0, 0, 178, 16902, 460544, 7121330, 78869428, 657860428, 4116175026, 19280190180, 68067144514, 183171836012, 379227347730, 605935215088, 745332928752, 702113842442, 503399988380, 272878823040, 111161557082, 33924337564, 7758154760, 1318925414, 158040254, 11479392, 348600, 0, sum=3638600663040
Эх, замедлилось продвижение к c19, всего лишь c27. Похоже даже если ещё немного ускорить (раза в полтора-два, резервы ещё есть), всё равно c19>0 недостижимо ...