У меня с функцией
valuation такой казус приключился. Я рискнул добавить последовательность производную от простых экспонет Вагстафа
A000978.
Описание было такое:
"a(n) is the positive integer x such that 3^((W-1)/(2*p)) == +-2^x (mod W), where p > 3 is prime, W=(2^p+1)/3 is the n-th Wagstaff prime and x < p."
Её одобрили и опубликовали. Но через несколько дней я нашёл в ней ошибку, расстроился и попросил админов удалить последовательность, они удалили.
Последовательность была такая:
Код:
3, 4, 1, 5, 15, 11, 11, 4, 36, 57, 44, 66, 5, 103, 150, 145, 287, 267, 329, 476, 2363, 2915, 156, 3641, 8305, 2689, 10266, 7662
Ошибка была в третьем члене, должно быть 0 вместо 1. На этом члене степень двойки вырождается в единицу.
Сейчас ничего лучше такого исправленного кода придумать не могу
Код:
wgs()=
{
A000978=[3, 5, 7, 11, 13, 17, 19, 23, 31, 43, 61, 79, 101, 127, 167, 191, 199, 313, 347, 701, 1709, 2617, 3539, 5807, 10501, 10691, 11279, 12391, 14479, 42737, 83339, 95369, 117239, 127031, 138937, 141079, 267017, 269987, 374321, 986191, 4031399];
for(n=2, #A000978,
p= A000978[n];
W=(2^p+1)/3;
s= (W-1)/2/p;
T= Mod(3, W)^s;
if(T==1,
print1(0, ", ")
,
x= valuation(lift(T), 2);
if(x, print1(x, ", "), x= valuation(lift(-T), 2); print1(x, ", "))
)
)
};
Возможно ли проще сделать код? И добавляют ли такие последовательности с одним диссонирующим членом?