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 A351947 #4 Apr 03 2022 13:53:30 %S A351947 450,882,900,1764,2178,2450,3042,3675,4356,4900,5202,6050,6084,6498, %T A351947 7350,7560,8450,9075,9522,10404,11025,11858,11880,12100,12675,12996, %U A351947 13500,14040,14450,15138,16562,16632,16900,17298,17787,18050,18150,18360,19044,19656,20520,21000,21675,23716,24642,24840,24843,25350 %N A351947 Numbers k for which A327500(k) <> A351946(k). %o A351947 (PARI) %o A351947 A051903(n) = if((1==n),0,vecmax(factor(n)[, 2])); %o A351947 A181819(n) = factorback(apply(e->prime(e),(factor(n)[,2]))); %o A351947 A351564(n) = issquarefree(factorback(apply(e->prime(e),(factor(n)[,2])))); %o A351947 A327499(n) = fordiv(n,d,if(A351564(n/d), return(d))); %o A351947 A327500(n) = { my(u=A327499(n)); if(u==n, 0, 1+A327500(u)); }; %o A351947 A351946(n) = A051903(A181819(n)); %o A351947 isA351947(k) = (A327500(k) != A351946(k)); %Y A351947 Cf. A327500, A351946. %K A351947 nonn %O A351947 1,1 %A A351947 _Antti Karttunen_, Apr 03 2022