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 A301867 #12 Aug 19 2019 03:16:54 %S A301867 1,170,679,5770,9154,9809,14322,22413,37114,51455,82615,100821,101153, %T A301867 115430,139954,171069,198462,222893,233358,270485,349894,389045, %U A301867 391090,514294,561782,595122,608685,664853,701338,815538,1009527,1231230,1290993,1397091 %N A301867 Numbers n such that phi(n), psi(n) and sigma(n) are simultaneously perfect squares. %C A301867 Subsequence of A067781. %C A301867 The squarefree terms of A067781 are in this sequence (since if n is squarefree then psi(n) = sigma(n)). The nonsquarefree terms of this sequence are: 1, 25264008, 2171889720, 2659240584, 5893922664, 16962447600,... %H A301867 Amiram Eldar, <a href="/A301867/b301867.txt">Table of n, a(n) for n = 1..1000</a> %e A301867 170 is in the sequence since phi(170) = 64 = 8^2, sigma(170) = psi(170) = 324 = 18^2. %e A301867 25264008, the first nonsquarefree term ( > 1 ) in the sequence, has phi = 2520^2, sigma = 9600^2 and psi = 8064^2. %t A301867 a = {}; psi[n_] := DirichletConvolve[j, MoebiusMu[j]^2, j, n]; %t A301867 aQ[n_] := AllTrue[Sqrt[{DivisorSigma[1, n], EulerPhi[n], psi[n]}], IntegerQ]; %t A301867 Do[If[aQ[k], AppendTo[a, k]], {k,1,10^7}]; a %Y A301867 Cf. A000010 (phi), A000203 (sigma), A001615 (psi), A067781, A291549. %K A301867 nonn %O A301867 1,2 %A A301867 _Amiram Eldar_, Mar 28 2018