%I #6 Mar 01 2019 23:33:27
%S 37,48,111,147,148,525,832,1036,1443,2457,3367,3468,4477,5887,6591,
%T 6993,7696,11011,12025,12096,12432,12493,12636,12691,12943,12987,
%U 13357,13377,13467,13468,333667,334668
%N Numbers k which are concatenations k=x//y such that x^2 + y^2 - x*y = k.
%e 147 is a term because 147 = 14//7 = 14^2 + 7^2 - 14*7.
%p Lton := proc(L) add( op(i,L)*10^(i-1),i=1..nops(L)) ; end:
%p for n from 10 do dgs := convert(n,base,10) ; for spli from 1 to nops(dgs)-1 do ydgs := [op(1..spli,dgs)] ; xdgs := [op(spli+1..nops(dgs),dgs)] ; if op(-1,ydgs) <> 0 then x := Lton(xdgs) ; y := Lton(ydgs) ; if y^2+x^2-x*y = n then print(n) ; fi; fi; od: od: # _R. J. Mathar_, Jul 16 2009
%K nonn,base
%O 1,1
%A _Claudio Meller_, Jul 06 2009
%E Definition simplified and 37, 48, 333667, 334668 added by _R. J. Mathar_, Jul 16 2009
|