%I #11 Aug 25 2023 01:14:59
%S 3697,3940,4802,5845,5905,21127,21715,22182,22867,22897,23380,23640,
%T 24367,26815,28812,28910,32192,33705,33815,35000,35065,35070,35430,
%U 35977,37082,37712,40277,44535,122915,125947,126762,128350,129670,130290,133092,134397
%N Numbers k > 1 such that, in base 6, k and k^2 contain the same digits in the same proportion.
%H Ray Chandler, <a href="/A061660/b061660.txt">Table of n, a(n) for n = 1..1000</a> (First 500 terms from Alois P. Heinz)
%e 3697 = 25041_6 and 3697^2 = 1204540521_6.
%p p:= n-> add(x^i, i=convert(n, base, 6)):
%p a:= proc(n) option remember; local k;
%p for k from 1+`if`(n=1, 0, a(n-1))
%p while p(k)*2<>p(k^2) do od; k
%p end:
%p seq(a(n), n=1..50); # _Alois P. Heinz_, May 10 2015
%Y Cf. A061656, A061657, A061658, A061659, A061661, A061662, A061663, A114258, A061664.
%K base,easy,nonn
%O 1,1
%A _Erich Friedman_, Jun 16 2001
%E Offset changed to 1 by _Alois P. Heinz_, May 10 2015