%I#11 2022年9月8日08:46:13
%S 2,5,7,11,79,23,47769,5932831,831250931073329167733317912473,
%电话:227268435711263261133473593834848946675648947914706467503,
%电话:705492079724756358710211264577133669153313171049839863191449887
%形式为sigma(N)+N的除数乘积的N个素数。
%C或者:A259973中产生的底漆。
%H K.D.Bajpai,n的表,n的a(n)=1..10000</a>
%e a(5)=79;除数(8)={1,2,4,8};总和=1+2+4+8=15;乘积=1*2*4*8=64;15+64=79是质数。
%e a(8)=769;除数(27)={1,3,9,27};总和=1+3+9+27=40;乘积=1*3*9*27=729;40+729=769,这是质数。
%p with(numtheory):A260108:=n->(sigma(n)+convert(divisors(n),`*`)):选择(isprime,[seq((A260108(n)),n=1..800)));
%t选择[Table[DivisorSigma[1,n]+Times@@Divisors[n],{n,1,1000}],PrimeQ]
%o(PARI)(n=11000,d=除数(n));k=西格玛(n)+产品(i=1,#d,d[i]);如果(i质数(k),打印1(k,“,”));
%o(PARI)A007955(n)=如果
%o列表(lim)=v=列表([2]);对于素数(p=2,(lim-1)\2,if(isprime(2*p+1),listput(v,2*p/1)));对于素数(p=2,sqrtnint(lim\1,3),my(t=p^3+p^2+p+1);如果(t>lim,断裂);if(isprime(t),listput(v,t));对于复合材料(n=4,平方(lim\1),my(t=A007955(n)+σ(n));如果(t<=lim&&isprime(t),listput(v,t));集(v)\\_Charles R Greathouse IV_,2015年7月17日
%o(岩浆)[1..1000]中的k:n | IsPrime(k),其中k是(&*除数(n)+SumOfDivisors(n))]
%Y参考A0000203、A007955、A065512、A118369、A259973。
%K nonn公司
%O 1,1号机组
%A _K.D.Bajpai,2015年7月16日
|