login
First prime greater than or equal to phi(n^2).
1

%I #13 Oct 26 2017 01:51:03

%S 2,2,7,11,23,13,43,37,59,41,113,53,157,89,127,131,277,109,347,163,257,

%T 223,509,193,503,313,487,337,821,241,937,521,661,547,853,433,1361,691,

%U 937,641,1657,509,1811,881,1087,1013,2179,769,2063,1009,1637,1249,2767

%N First prime greater than or equal to phi(n^2).

%C The sequence generally goes up,down,up,down...

%H Robert Israel, <a href="/A079476/b079476.txt">Table of n, a(n) for n = 1..10000</a>

%F a(n) = A007918(A002618(n)). - _Robert Israel_, Oct 25 2017

%e phi(3^2)=phi(9)=6, therefore a(3)=7.

%p map(t -> nextprime(numtheory:-phi(t^2)-1), [$1..100]); # _Robert Israel_, Oct 25 2017

%t Array[If[PrimeQ@ #, #, NextPrime@ #] &@ EulerPhi[#^2] &, 53] (* _Michael De Vlieger_, Oct 25 2017 *)

%o (PARI) for (n=1,100, print1(nextprime(eulerphi(n^2))","))

%Y Cf. A002618, A007918.

%K nonn,look

%O 1,1

%A _Jon Perry_, Jan 15 2003

%E Corrected by _Robert Israel_, Oct 25 2017