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.

A066733 Numbers such that the product of the digits of its square is a square > 0.

This page as a plain text file.
%I A066733 #15 Dec 06 2024 20:19:13
%S A066733 1,2,3,7,12,17,21,38,88,106,107,108,109,117,122,128,129,141,146,164,
%T A066733 168,171,173,178,191,196,204,206,207,208,209,212,221,222,236,263,276,
%U A066733 278,288,306,342,359,364,367,372,377,394,432,463,478,479,518,537,538
%N A066733 Numbers such that the product of the digits of its square is a square > 0.
%H A066733 Harry J. Smith, <a href="/A066733/b066733.txt">Table of n, a(n) for n = 1..1000</a>
%e A066733 17 is in the sequence because the square of 17 is 289 and 2*8*9 = 144 = 12^2.
%t A066733 Do[a = Apply[Times, IntegerDigits[n^2]]; If[ a != 0 && IntegerQ[a^(1/2)], Print[n]], {n, 1, 10^4} ]
%t A066733 nzpQ[n_]:=Module[{prod=Times@@IntegerDigits[n^2]},prod!=0 && IntegerQ[ Sqrt[ prod]]]; Select[Range[600],nzpQ] (* _Harvey P. Dale_, May 27 2012 *)
%o A066733 (PARI) isok(k)= { my(p=vecprod(digits(k^2))); p > 0 && issquare(p) } \\ _Harry J. Smith_, Mar 20 2010
%Y A066733 Cf. A067071.
%K A066733 nonn,base
%O A066733 1,2
%A A066733 _Robert G. Wilson v_, Jan 15 2002