%I#24 2022年1月15日21:40:25
%S 1,2,4,18,721302804958409009301216188026463183376217437,
%电话:21600294003065435207450245230266133105198012041641204320,
%电话:21934082494520254974532721783671570384337851667085398785870980589222708922480
%N数k,使非素数(k)除以k^2。
%第四个非素数是8,8除以4^2=16,所以8在序列中。
%从1do到n的p
%p如果n^2 mod A018252(n)=0,则
%p打印(n);
%p end if;
%截止日期:#R.J.Mathar_,2013年6月4日
%t序列={};np=n=0;p=2;当[n<10^5时,如果[++np==p,p=NextPrime@p,n++;如果[Mod[n^2,np]==0,打印@{n,np,n^2/np};附录[seq,n]]];seq(*乔瓦尼·雷斯塔,2013年5月27日*)
%t模[{nn=10^7,np,len},np=Complement[Range[nn],Prime[Range[PrimePi[nn]]];len=长度[np];转置[Select[Thread[{Range[len],np}],Divisible[First[#]^2,Last[#]]&]][1]](*哈维·P·戴尔,2014年9月10日*)
%Y参考A018252。
%K nonn,更少
%O 1,2号机组
%A _Gerasimov Sergey_,2013年5月27日
%E a(5)-a(38)摘自_Giovanni Resta_,2013年5月27日
|