login
a(n) is the smallest number k with prime sum of divisors such that tau(k) = n-th prime.
2

%I #25 Sep 08 2022 08:46:18

%S 2,4,16,64,9765625,4096,65536,262144,

%T 1471383076677527699142172838322885948765175969,

%U 10264895304762966931257013446474591264089923314972889033759201,1073741824,18701397461209715023927088008788055619800417991632621566284510161

%N a(n) is the smallest number k with prime sum of divisors such that tau(k) = n-th prime.

%C tau(n) = A000005(n) = the number of divisors of n.

%C a(11) = 1073741824; a(n) > A023194(10000) = 5896704025969 for n = 9, 10 and n >= 12.

%H Davin Park, <a href="/A278913/b278913.txt">Table of n, a(n) for n = 1..50</a>

%F a(n) = A123487(n)^(prime(n)-1). - _Davin Park_, Dec 10 2016

%e a(3) = 16 because 16 is the smallest number with prime values of sum of divisors (sigma(16) = 31) such that tau(16) = 5 = 3rd prime.

%t A278913[n_] := NestWhile[NextPrime, 2, ! PrimeQ[Cyclotomic[Prime[n], #]] &]^(Prime[n] - 1) (* _Davin Park_, Dec 28 2016 *)

%o (Magma) A278913:=func<n|exists(r){k:k in[1..10000000] | IsPrime(SumOfDivisors(k)) and NumberOfDivisors(k) eq NthPrime(n)} select r else 0>; [A278913(n): n in[1..8]]

%o (PARI) a(n) = {my(k=1); while(! (isprime(sigma(k)) && isprime(p=numdiv(k)) && (primepi(p) == n)), k++); k;} \\ _Michel Marcus_, Dec 03 2016

%Y Cf. A000005, A000203, A023194, A278914.

%K nonn

%O 1,1

%A _Jaroslav Krizek_, Nov 30 2016

%E More terms from _Davin Park_, Dec 08 2016