%I#32 2015年10月16日03:38:08
%S 7,8,13,19,20,24,25,28,33,36,37,40,43,49,50,52,55,56,63,64,66,67,68,
%电话74,75,79,85,88,98100103108109112113116117123124126131,
%U 132133134136140143145150153156159160163164167168
%N N-d(N)从不取这些值,其中d(N)=A000005(N)。
%A236562的C补体_雅罗斯拉夫·克里泽克,2014年2月9日
%C A060990中零的位置,边关系A049820(child)=父级生成的树中的叶节点_Antti Karttunen,2015年10月6日
%C由于A000005(x)<=1+x/2,如果没有x<=2*(n+1)且n=x-d(x),则n在序列中_罗伯特·伊斯雷尔(Robert Israel),2015年10月12日
%C这可以改进为:如果没有x<=n+A002183(2+A261100(n)),n=x-d(x),则n在序列中。另请参阅A070319、A262686_Antti Karttunen,2015年10月12日
%H Antti Karttunen,n的表,n=1..10000的a(n)</a>
%p N:=1000:#获取所有项<=N
%p排序(转换({$1..N}减去{seq(x-numtheory:-tau(x),x=1..2*(1+N))},list));#_罗伯特·伊斯雷尔(Robert Israel),2015年10月12日
%t lim=10000;取[Complement[Range@lim,Sort@DeleteDuplicates@Table[n-DivisorSigma[0,n],{n,lim}]],57](*_Michael De Vlieger_,2015年10月13日*)
%o(PARI)
%o分配((2^31)+(2^30));
%o上传=36756720+640;\\=A002182(53)+A002183(52)。
%o v060990=矢量(uplim);
%o表示(n=3,uplim,v060990[n-numdiv(n)]++);
%o A060990=n->如果(!n,2,v060990[n]);
%o uplim2=36756720;
%o n=0;k=1;while(n<=uplim2,如果(0==A060990(n),则写入(“b045765_big.txt”,k,“”,n);k++);n++;);
%o\\_Antti Karttunen_,2015年10月9日
%o(方案)
%o(定义A045765(零位置1 1 A060990))
%o;;同时使用安蒂·卡图宁的IntSeq-library,2015年10月6日
%Y A262898的顶行。
%Y参见A000005、A002183、A049820、A060990、A070319、A236562、A236565、A261100、A262511、A262686、A262901、A26902、A26290、A2629、A263081。
%Y参考A263091(此序列中的素数),A263095(正方形)。
%Y参见A259934(给出了同一棵树的无限树干,推测是唯一的)。
%K nonn公司
%O 1,1号机组
%A·热心的W·威尔逊_
|