%I#12 2019年12月7日12:18:25
%S 3,4,5,12,24,41,60,8411314418026426531248076110121104,
%电话:1105120114051624174017411861211232133281415304,
%电话:5512572459406385661273207564756585808845938596609619941
%N个反除数和为素数的数。
%C关于反除数的定义,请参见A066272。
%H Chai Wah Wu,n的表格,n=1..1000的a(n)</a>
%H Jon Perry<a href=“http://www.users.globalnet.co.uk/~perry/methemats/antidivisor.htm“>反除数。
%H乔恩·佩里<a href=“/A066272/a066272a.html”>反除数</a>[缓存副本]
%H乔恩·佩里<a href=“/A066272/A066272.html”>反除数:更多反除数</a>[缓存副本]
%761的反除数是{2,3,9,13,39,117,169,507},这些和是859,一个素数,所以761在序列中。
%t反除数[n_]:=并集[Drop[Divisors[2*n-1],1],-1],映射[2*n/#&,Drop[Select[Divisor[2*n],OddQ],1]],Drop[Cdivisors[2],1];Do[p=Plus@@AntiDivisors[n];如果[PrimeQ[p],打印[n]],{n,2,10^4}]
%o(Python)
%o从sympy导入isprime,divisors
%o A109350=[n代表范围(1,10**6)中的n,如果是素数(如果n>2*d和n%(2*d),则求和([2*d代表除数(n)中的d,如果n>d>=2和n%d]+[d代表除法中的d(2*n+1),如果n>d>=2和n%d)]
%o#_Chai Wah Wu_,2014年8月12日
%K nonn公司
%O 1,1号机组
%2005年8月21日,A_Ryan Propper_
|