

A137818


Nonbiquadratefree "year numbers": phi(n) = 2 phi(sigma(n)) and p^4  n for some p>1.


2



295569, 1811079, 1964375, 2069469, 4473387, 5854375, 10936053, 13260625, 18029709, 21576537, 22182093, 25536875, 35595625, 46404333, 49648383, 55094375, 57044817, 58650625, 67009923, 69166467, 72681875, 76106875
(list;
graph;
refs;
listen;
history;
text;
internal format)



OFFSET

1,1


COMMENTS

See A137815 for general comments and references about "year numbers". This is the subsequence of elements of A137815 divisible by a biquadrateful number, i.e. its intersection with A046101 (numbers divisible by the 4th power of some prime). As such, it is of course also a subsequence of A137817 and a fortiori of A137816.
There are only 28 such numbers below 10^8.


LINKS

Donovan Johnson, Table of n, a(n) for n = 1..500


PROG

(PARI) for( i=1, #A137816, vecmax( factor( A137816[i])[, 2])>3 && print1(A137816[i]", "))
(PARI) for( i=1, #A046099, eulerphi(A046099[i])==2*eulerphi(sigma(A046099[i])) && print1( A046099[i] ", "))
(PARI) for( n=1, 10^9, issquarefree(n) && next; vecmax(factor(n)[, 2])>3  next; eulerphi(n)==2*eulerphi(sigma(n)) && print1(n", "))


CROSSREFS

Cf. A137815A137819, A046101.
Sequence in context: A254919 A109031 A137817 * A303357 A251387 A164946
Adjacent sequences: A137815 A137816 A137817 * A137819 A137820 A137821


KEYWORD

nonn


AUTHOR

R. K. Guy, R. J. Mathar and M. F. Hasler, Feb 11 2008


STATUS

approved



