%I#10 2020年8月16日19:09:38

%S 1,4,9,10,14,21,27,33,40,46,52,65,71,79,84,90,98100104108111133，

%电话137140162166191195210212225241249254270278283301323，

%电话：32733035235638138540040241543143944460468473497513

%N N使d（1）*1！+d（2）*2！+…+d（k）*k！是一个正方形，其中d（i）是n的十进制数字。

%C A198095包含在该序列中。

%C如果n在序列中，那么10*n.-Robert Israel_，2020年8月16日

%H Robert Israel，n的表，a（n）表示n=1..10000</a>

%e 137在序列中是因为1*1！+3*2！+7*3! = 1 + 6 + 42 = 7^2.

%p从1到520 do:l:=长度（n）：n0：=n:s:=0：对于m从l到-1 do:q:=n0:u:=irem（q，10）：v:=iquo（q，10）：n0：=v:s:=s+u*m！：od:z:=sqrt（s）：如果z=floor（z），则打印f（`%d，`，n）：否则fi:od:

%p#更简单：

%p过滤器：=程序（n）

%p局部L，k；

%p L：=换算（n，基数，10）；

%p issqr（添加（L[-k]*k！，k=1..nops（L））

%p端程序：

%p选择（过滤器，[1..1000]）；#_罗伯特·伊斯雷尔，2020年8月16日

%Y参考A198095、A198044。

%K nonn，基础

%O 1,2号机组

%2011年10月21日，拉格瑙市