Версия улучшенная, спасибо
Dmitriy40 за замечания
Код:
fpr(x)=my(f=factor(x),n=#f[,1]);print1("$",x,"_{10}=[");for(i=1,n,print1("P_{",primepi(f[i,1]),"}"));print("](",strjoin(f[,2],","),")$")
-- 10.08.2023, 13:38 --А использовать в публикациях эту прогу можно?
Да, пожалуйста.
и как на Вас ссылаться?
Там нет никакой интеллектуальной ценности, можно не ссылаться :)