cp's OEIS Frontend

A286095 Composite numbers n such that tau(n) (number of divisors of n) is prime and sigma(n) (sum of divisors of n) is not prime.

%I A286095 #30 Jun 05 2017 03:48:43
%S A286095 49,81,121,169,361,529,625,841,961,1024,1369,1849,2209,2809,3721,4489,
%T A286095 5329,6241,6889,9409,10609,11449,11881,12769,14641,16129,18769,19321,
%U A286095 22201,22801,24649,26569,32041,32761,36481,37249,38809,39601,44521,49729,51529,52441,54289
%N A286095 Composite numbers n such that tau(n) (number of divisors of n) is prime and sigma(n) (sum of divisors of n) is not prime.
%C A286095 If sigma(n) is prime (A023194) then tau(n) is prime too. (See Crux Mathematicorum link.)
%C A286095 But the reverse is false; the numbers which verify tau(n) prime and sigma(n) not prime are in the sequence A275938.
%C A286095 All odd primes belong to the sequence A275938, but there are also in this sequence composite numbers which are all prime powers, these prime powers are here.
%H A286095 Robert Israel, <a href="/A286095/b286095.txt">Table of n, a(n) for n = 1..10000</a>
%H A286095 Peter A. Lindstrom and Andy Liu, <a href="https://cms.math.ca/crux/backfile/Crux_v6n06_Jun.pdf">Problem 465</a>, Crux Mathematicorum, page 188, Vol.6 , Jun. 80.
%e A286095 tau(49) = 3 and sigma(49) = 57 = 3 * 19.
%p A286095 for n from 2 to 550000 do p(n):=tau(n);
%p A286095 if not isprime(n) and is prime(p(n)) and not isprime(sigma(n)) then print (n,p(n),sigma(n)) else fi; od:
%p A286095 # alternative
%p A286095 N:= 10^5: # to get all terms <= N
%p A286095 P:= select(isprime, [2,seq(i,i=3..isqrt(N),2)]):
%p A286095 S:= {}:
%p A286095 for p in P do
%p A286095   k:= 1:
%p A286095   do
%p A286095     k:= nextprime(k+1)-1;
%p A286095     if p^k > N then break fi;
%p A286095     if not isprime((p^(k+1)-1)/(p-1)) then S:= S union {p^k} fi
%p A286095   od
%p A286095 od:
%p A286095 sort(convert(S,list)); # _Robert Israel_, Jun 05 2017
%t A286095 Select[Range[10^5], Function[n, And[CompositeQ@ n, Map[PrimeQ@ DivisorSigma[#, n] &, {0, 1}] == {True, False}]]] (* _Michael De Vlieger_, May 24 2017 *)
%o A286095 (PARI) lista(nn) = {forcomposite(n=1, nn, if (isprime(numdiv(n)) && !isprime(sigma(n)), print1(n, ", ")););} \\ _Michel Marcus_, May 24 2017
%Y A286095 Cf. A000005, A000203, A023194, A275938.
%K A286095 nonn
%O A286095 1,1
%A A286095 _Bernard Schott_, May 22 2017