login
a(n) is the least k such that phi(k) + d(k) = A357916(n), where phi(k) = A000010(k) is Euler's totient function, and d(k) = A000005(k) is the number of divisors of k.
2

%I #14 Oct 25 2022 20:04:48

%S 1,2,4,16,25,81,121,256,484,1296,529,1024,1600,2116,2401,7744,11664,

%T 5041,7225,11236,20164,10201,25600,12769,30976,46656,21025,17161,

%U 44944,51076,29929,84100,73984,36481,75076,107584,54289,63001,87025,69169,101761,126025,215296,256036,252004,295936

%N a(n) is the least k such that phi(k) + d(k) = A357916(n), where phi(k) = A000010(k) is Euler's totient function, and d(k) = A000005(k) is the number of divisors of k.

%C Numbers k such that A061468(k) = phi(k) + d(k) is prime, and no smaller number gives the same value of A061468, sorted in order of the prime values.

%C All terms except 2 are squares, because if k > 2, phi(k) is even, and if d(k) is odd, k must be a square.

%C All numbers in this sequence are elements of A225983. For an example, this excludes all numbers of the form (6*m)^2 but only if m is not divisible by 6. - _Thomas Scheuerle_, Oct 20 2022

%H Robert Israel, <a href="/A357917/b357917.txt">Table of n, a(n) for n = 1..3000</a>

%F A061468(a(n)) = A000010(a(n)) + A000005(a(n)) = A357916(n).

%e a(4) = 16 because phi(16) + d(16) = 8 + 5 = 13 = A357916(4), and no smaller number than 16 works.

%p N:= 10^6:

%p pmax:= evalf(N/(exp(gamma)*log(log(N))+3/log(log(N))));

%p V:= 'V': P:= {3}: V[3]:= 2:

%p for k from 1 to sqrt(N) do

%p n:= k^2;

%p v:= numtheory:-phi(n)+numtheory:-tau(n);

%p if v <= pmax and isprime(v) and not member(v,P) then

%p P:= P union {v}; V[v]:= n;

%p fi

%p od:

%p P:= sort(convert(P,list)):

%p seq(V[p], p=P);

%Y Cf. A000005, A000010, A061468, A225983, A357916.

%K nonn

%O 1,2

%A _J. M. Bergot_ and _Robert Israel_, Oct 19 2022