%I #28 Feb 11 2014 01:15:56
%S 489,4629,296206,460029,29589106,46000029,2927272726,4045046518,
%T 21223345084,29600331295,296151515206,460000000029
%N Numbers n such that phi(n) + sigma(n) = reversal(n) - 4.
%C If p = 15*10^k+10+(10^k-1)/3 is prime then 3*p is in the sequence.
%C a(7) is greater than 1.3*10^8.
%C a(11) > 10^11. - _Donovan Johnson_, Nov 08 2013
%C If p=(1/11)*(23*100^m-1) is prime then 14*p is a term of the sequence. - _Farideh Firoozbakht_, Nov 08 2013
%C a(13) > 10^13. - _Giovanni Resta_, Feb 08 2014
%C If p = (1685*10^(2k+2)+31)/33 is prime then 58*p is in the sequence. For k = 0, 3, 9, 30, 42, 51, 120, 846, ... p is prime. - _Farideh Firoozbakht_, Feb 10 2014
%e 489 is in the sequence because phi(489)+sigma(489) = 324+656 = 984-4 = reversal(489)-4.
%t Do[If[FromDigits@Reverse@IntegerDigits@n-4 == EulerPhi[n] + DivisorSigma[1, n], Print[n]], {n, 130000000}]
%o (PARI) is(n)=subst(Polrev(digits(n)),'x,10)-4==eulerphi(n)+sigma(n) \\ _Charles R Greathouse IV_, Nov 08 2013
%Y Cf. A000010, A000203, A004086, A070272, A230004, A230006, A230019, A136544, A237521, A237522.
%K nonn,base,more
%O 1,1
%A _Farideh Firoozbakht_, Nov 07 2013
%E a(7)-a(10) from _Donovan Johnson_, Nov 08 2013
%E a(11)-a(12) from _Giovanni Resta_, Feb 06 2014