%I#56 2021年2月22日21:37:19

%S 2,3,5,7,11,13,23,29,31,43,47,53,59,61,67,71,79107131139157173，

%电话：21126326927728331173313473493673421431461463547，

%电话：557599643659661683691709733743787827853859863911941

%N素数p的形式为p=2^k*m+1，其中（i）m是无平方的奇数，（ii）所有除m的素数都在序列中，（iii）k是0、1或2。

%在定义中取m=1，我们得到素数2，3，5。

%C如果n在A226960中，则n是该序列项的乘积。

%如果k只允许为0或1，那么我们得到2、3、7、43以及更多_宋佳宁，2021年2月21日

%C也是A341858中术语的主要因素。据推测，这个序列是无限的_2021年2月22日，宋嘉宁

%H Ray Chandler，n的表，n=1..1000时的a（n）</a>

%t fa=系数整数；自由[n_]：=n==乘积[fa[n][[i，1]]，{i，长度[fa[n]}]；Os[b_，1]=真；Os[b_，b_]=真；Os[b_，n_]：=Os[b，n]=PrimeQ[n]&&free[（n-1）/b^IntegerExponent[n-1，b]]&&IntegerIndonent[n-1，b]<3&&Union@表[Os[b，fa[n-1][i，1]]]，｛i，Length[fa[n-1]]｝]==｛True｝；G[b_]：=选择[Prime[Range[1000]]，Os[b，#]&]；G[2]

%o（PARI）是（n）=如果（！isprime（n），return（0））；如果（n＜13，返回（1））；my（k=估值（n-1,2），m=n>>k，f）；如果（k>2，返回（0））；f=系数（m）；如果（lcm（f[，2]）>1，则返回（0））；对于（i=1，#f~，如果（！是（f[i，1]），返回（0））；2013年10月28日，查尔斯·R·格里特豪斯四世

%Y参见A227007、A226960、A007814、A229290、A229291、A289355、A341858。

%Y关于补码，请参见A289355。

%Y A066651的正确子序列。

%K nonn公司

%O 1,1号机组

%A _何塞·玛丽亚·格劳·里巴斯（JoséMaríA Grau Ribas），2013年10月5日

%E修订了2013年11月13日_Charles R Greathouse IV_的定义

%E术语由JoséMaría Grau Ribas更正，2013年11月14日