This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A066111 #12 Apr 06 2020 03:15:18 %S A066111 2,3,5,13,17,31,43,61,83,109,121,125,131,229,239,257,263,269,311,313, %T A066111 343,361,443,463,503,571,593,599,619,641,647,653,659,701,797,811,853, %U A066111 953,967,1009,1031,1039,1063,1123,1373,1459,1483,1499,1663,1669,1693 %N A066111 Prime powers m such that sigma_4(m^2)/sigma_2(m^2) is prime. %C A066111 Numbers m = p^w such that A001159(m^2)/A001157(m^2) is prime, i.e., m^2 is in A066109. %C A066111 Also m is the square root of a term from A066109 (omitting the term 20). Apart from 20, up to 10000000 A066109 consists of squares of prime powers. %H A066111 Amiram Eldar, <a href="/A066111/b066111.txt">Table of n, a(n) for n = 1..10000</a> %e A066111 m=125: m^2 = 15625 = A066109(13), sigma_4(15625) = 59700165039453751, sigma_2(15625) = 254313151, sigma_4/sigma_2 = 234750601 = A066110(13) is prime. Observe also that sigma_2 is close to sigma_4/sigma_2. %o A066111 (PARI) isok(m) = isprimepower(m) && isprime(sigma(m^2, 4)/sigma(m^2, 2)); \\ _Michel Marcus_, Apr 06 2020 %Y A066111 Cf. A001157, A001159, A066109, A066110. %K A066111 nonn %O A066111 1,1 %A A066111 _Labos Elemer_, Dec 06 2001