%I#13 2020年7月2日13:05:11

%S 1,1,12,4240，4813441923840128033792307255910430080245760，

%电话：49152167116809830405229772802752512201211105280173015040，

%电话：144703488062914560229008998480179968140928614446976204811676942330407689065201664

%N a（N）=分母（b_N（x）），其中b_N是A335947中定义的多项式。

%C通过Kellner和Sondow（2019）的方法，也可以在不参考伯努利多项式的情况下计算序列（最终得益于von Staudt-Clausen定理）。比较SageMath程序。

%H Bernd C.Kellner和Jonathan Sondow，<a href=“https://arxiv.org/abs/1902.10672“>关于Carmichael和多边形数、Bernoulli多项式和以p为底的数字之和，arXiv:1902.10672[math.NT]，2019。

%F a（n）=min{m|m*（[x^k]b（n，x））是所有k=0..n}的整数。

%F（n）的奇数部分是无平方的（A000265）。

%F a（n）和A144845（n）具有相同的奇素因子。

%F a（n）/A144845（n）=4^楼层（n/2）/2，对于n>=1。

%Fα（n）/拉德（a（n））=A158302（n+1），（拉德=A007947）。

%o（SageMath）

%o定义A335949（n）：

%o a=集合（素数除数（n+1））-集合（[2]）

%o b=(

%o p公司

%p在prime_range中为o（3，（n+2）//（2+n%2））

%如果不是p，则为o。除以（n+1）和和（n+1。数字（基数=p）>=p

%o）

%o p=列表（a.union（集合（b）））

%o返回4^（n//2）*mul（p）

%o打印（[A335949（n）代表范围（31）内的n]）

%Y参考A335947/A335948、A144845、A158302。

%K nonn公司

%0、3

%A _彼得·卢什尼，2020年7月1日