%I#18 2024年2月27日07:36:16

%S 1,5,7,1,11,13,17,19,1,23,1,1,29,31,1,1,37,1,41,43,1,47,1,1,53,1,1，

%电话59,61,1,1,67,1,71,73,1,1,79,1,83,1,1,89,1,1,1,97,1101103,1107109，

%U 1113,1,1,1,1127,1131,1,1137139,1,1,11,1149151,1,1157,11163

%N a（N）=A145609（N）的平方部分。

%C关于无方形零件，请参见A145738。A128059是一个非常相似的序列。

%如果2n+1是素数，F a（n）=2n+1，否则，对于n>1。

%F来自_Gary Detlefs_，2011年10月18日：（开始）

%F a（n）=分母（n！*（和{k=1..n}k^3）/（和{k=1..n{k^2））

%F=分母（n！*3*n*（n+1）/（2*（2*n+1）））。（结束）

%p序列（分母（n！*3*n*（n+1）/（2*（2*n+1））），n=1..81）；#_Gary Detlefs，2011年10月18日

%t m=1；aa={}；做[k=0；做[k=k+m^（2r+1-d）/d，{d，1，2r}]；b=Sqrt[分子[k]]/。Sqrt[_]->1；附加到[aa，b]，{r，1，137}]；aa（*_阿图尔·贾辛斯基_*）

%o（Python）

%o来自sympy import isprime

%o def A145737（n）：返回a if isprime（a:=（n<<1）+1）和n>1其他1#_Chai Wah Wu_，2024年2月26日

%Y参考A008833、A128059、A145609、A145738。

%K nonn公司

%O 1,2号机组

%A _阿特尔·贾辛斯基_，2008年10月17日