A111267号:=proc(nmin)局部a,k,qn,dvs;a:=[1];qn:=5;而nops(a)<nmin-do-dvs:=排序([op(numtheory[除数](qn))]);对于k从1到nops(dvs),如果不是a中的dvs[k],则a:=[op(a),dvs[k]];qn:=qn+ithprime(nops(a)+1);断裂;fi;od;od;返回(a);结束时间:A111272号:=proc(nmin)局部a,a111267,n,i;a:=[];a111267:=A111267号(nmin);n:=1;成员(ithprime(n),a111267,'i')执行a:=[操作(a),i];n:=n+1;od;返回(a);结束时间:A111272号(1000) ; - Richard J.Mathar(Mathar(AT)strw.leidenuniv.nl),2007年8月20日