cp's OEIS Frontend

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.

A351947 Numbers k for which A327500(k) <> A351946(k).

This page as a plain text file.
%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