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 A354179 #20 Nov 04 2023 12:26:45 %S A354179 1,8,27,32,64,125,216,243,256,343,512,729,864,1000,1331,1728,1944, %T A354179 2048,2197,2744,3125,3375,4000,4913,5832,6561,6859,6912,7776,8000, %U A354179 9261,10648,10976,12167,13824,15552,15625,16384,16807,17576,19683,21952,23328,24389,25000 %N A354179 Numbers whose square has a number of divisors that is coprime to 30. %C A354179 Numbers k such that gcd(d(k^2), 30) = 1, where d(k) is the number of divisors of k (A000005). %H A354179 Michael De Vlieger, <a href="/A354179/b354179.txt">Table of n, a(n) for n = 1..10000</a> (first 1709 terms from Amiram Eldar) %H A354179 <a href="/index/Eu#epf">Index entries for sequences computed from exponents in factorization of n</a>. %F A354179 a(n) = sqrt(A354178(n)). %F A354179 Sum_{n>=1} 1/a(n) = Product_{p prime} (p + p^4 + p^6 + p^7 + p^9 + p^10 + p^12 + p^15)/(p^15 - 1) = 1.2449394393... %e A354179 8 is a term since A000005(8^2) = 7 and gcd(7, 30) = 1. %t A354179 Select[Range[25000], CoprimeQ[DivisorSigma[0, #^2], 30] &] %o A354179 (PARI) isok(m) = gcd(numdiv(m^2), 30) == 1; \\ _Michel Marcus_, May 19 2022 %Y A354179 Cf. A000005, A354178. %Y A354179 Subsequence of A350014. %K A354179 nonn %O A354179 1,2 %A A354179 _Amiram Eldar_, May 18 2022