OFFSET
1,2
COMMENTS
In this sequence there are no odd prime numbers, but there is even prime number 2.
EXAMPLE
2 is in this sequence because d(2)/1 = 2 is prime and sigma(2)/1 = 3 is prime, d(2)/2 = 1 is no prime and sigma(2)/2 = 3/2 is no prime, where 1, 2 are divisors of 2.
MATHEMATICA
okQ[n_] := Block[{d = Divisors[n], t0, t1=DivisorSigma[1, n]}, t0 = Length@ d; Length @Select[d, PrimeQ[ t0/#] &] == Length@ Select[d, PrimeQ[ t1/#] &]]; Select[Range@ 1000, okQ] (* Giovanni Resta, Mar 05 2017 *)
PROG
(PARI) is(n)=my(f=factor(n), d=numdiv(f), fd=factor(d)[, 1], s=sigma(f), fs=factor(s)[, 1]); sum(i=1, #fd, n%(d/fd[i])==0)==sum(i=1, #fs, n%(s/fs[i])==0) \\ Charles R Greathouse IV, Feb 27 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Feb 27 2017
EXTENSIONS
a(7)-a(44) from Charles R Greathouse IV, Feb 27 2017
a(45)-a(60) from Giovanni Resta, Mar 05 2017
STATUS
approved