|
|
|
|
1, 5, 6, 2, 5, 4, 24, 9, 3, 10, 11, 11, 37, 12, 34, 4, 42, 33, 13, 45, 11, 15, 48, 18, 5, 16, 52, 17, 17, 48, 19, 18, 56, 18, 55, 6, 61, 59, 37, 20, 21, 65, 66, 12, 65, 64, 69, 22, 7, 50, 39, 23, 73, 71, 75, 16, 24, 72, 23, 40, 19, 25, 128, 8, 25, 26, 26, 41, 13, 52, 61, 27, 86, 28
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
LINKS
|
Michael De Vlieger, Table of n, a(n) for n = 1..10000
|
|
MATHEMATICA
|
Table[If[IntegerQ@ Sqrt@ n, Sqrt@ n, k = 1;
Function[a,
Function[t,
While[Function[b,
Function[w,
Nand[Times @@ Boole@ Map[w[[#1]] >= #2 & @@ # &, #] == 1,
Length@ SequencePosition[
Select[b,
Function[d, MemberQ[#[[All, 1]], d /. 0 -> 10]]],
a] > 0]]@ RotateLeft@ BinCounts[b, {0, 10, 1}]]@
IntegerDigits[k^2] &@
Apply[Join,
Map[Lookup[t, #] /.
d_ /; IntegerQ@ d :> If[d > 0, {d, #}, {10, #}] &,
Keys@t]], k++]]@
KeyDrop[PositionIndex@ RotateLeft@ BinCounts[a, {0, 10, 1}], 0]]@
IntegerDigits@ n; k], {n, 74}] (* Michael De Vlieger, May 05 2017, Version 10.1 *)
|
|
CROSSREFS
|
Sequence in context: A060011 A021068 A286300 * A038690 A154801 A253094
Adjacent sequences: A091870 A091871 A091872 * A091874 A091875 A091876
|
|
KEYWORD
|
base,easy,nonn
|
|
AUTHOR
|
Ray Chandler, Feb 07 2004
|
|
STATUS
|
approved
|
|
|
|