%我#19 2014年8月20日20:40:37
%S 1,2,6,30,36,54,90,96114120156174210216300330414510516546,
%电话:57666071472674480481083489493696610141044105611341170,
%电话:12961344135615001560158416261650160168168173417641770183618841926201020462064
%没有素数反除数的数字k。
%H Chai Wah Wu,n表,n(n)表示n=1..10000(n=1…167来自Michael De Vlieger)
%e a(3)=6,因为6有反除数4,所以它是复合的。
%e a(4)=30,因为30有反除数{4,12,20},并且没有素数。
%e所有整数6<k<30都至少有一个素数反除数,唯一没有素数反除数的整数k<6是k={1,2}。
%t primeAntiDivisors[n_]:=选择[Cases[Range[2,n-1],_?(Abs[Mod[n,#]-#/2]<1&)],PrimeQ];a241556[n_Integer]:=映射[Length[primeAntiDivisors[#]]&,Range[n]];压扁[位置[a241556[10^5],0]]
%o(Python)
%o从sympy导入isprime,divisors
%o A241557=[n代表范围(1,10**6)中的n,如果没有([i质数(x)代表x in
%o。。。。。。。。。。[2*d代表除数(n)中的d,如果n>2*d和n%(2*d)]+
%o。。。。。。。。。。[d代表除数(2*n-1)中的d,如果n>d>=2且n%d]+
%o。。。。。。。。。。[d代表除数(2*n+1)中的d,如果n>d>=2且n%d]])
%o#_Chai Wah Wu_,2014年8月19日
%Y参考A066272,A241556。
%K nonn公司
%O 1,2号机组
%2014年8月8日,A _Michael De Vlieger_
|