%I#25 2022年9月8日08:46:13
%S 1,5,11,17,23,25,29,41,47,53,55,59,71,83,85,89101107113115121,
%电话:125131137145149167173179187191197205227233235239251,
%电话:25325726326526927528293929531317317335359383
%N只能被6*m-1形式的素数整除的数。
%H Charles R Greathouse IV,n表,n=1..10000的a(n)</a>
%e 275是按顺序排列的,因为275=5*5*11,其中5=6-1和11=6*2-1。
%tokQ[n_]:=n==1||与@@(Mod[#,6]==5&/@FactorInteger[n][[All,1]]);选择[Range[400],okQ](*_Jean-François Alcover_,2015年7月2日*)
%o(岩浆)[1..400]|forall{d:d in PrimeDivisors(n)|d mod 6 eq 5}];
%o(PARI)list(lim)=my(v=list([1]),mn,mx,t);对于素数(p=5,lim=1,if(p%6==5,listput(v,p)));如果(lim<25,返回(Vec(v)));对于素数(p=5,平方(lim),如果(p%6<5,next));mx=1;而(v[mx+1]*p<=lim,对于(i=mn=mx+1,mx=v,t=p*v[i];如果(t>lim,中断);listput(v,t)));集(v)\\_Charles R Greathouse IV_,2018年1月11日
%Y A007310的子序列。
%Y参见A004611(数字只能被6*m+1形式的素数整除),A007528。
%K nonn公司
%O 1,2号机组
%A _布鲁诺·贝塞利_,2015年6月30日
|