%I #15 Sep 08 2022 08:44:40

%S 0,0,1,4,21,115,750,5628,47811,453999,4765190,54782414,684597510,

%T 9239952802,133952757201,2075907567336,34247231626809,599244998013963,

%U 11084682598687312,216127622080213321

%N a(n) = floor( (n/e)^n ).

%H G. C. Greubel, <a href="/A015554/b015554.txt">Table of n, a(n) for n = 1..445</a>

%p Digits := 64:f := n->floor(evalf((n/E)^n));

%t Table[Floor[(n/E)^n],{n,20}] (* _Harvey P. Dale_, Oct 18 2012 *)

%o (PARI) for(n=1, 20, print1(floor((n*exp(-1))^n), ", ")) \\ _G. C. Greubel_, Dec 30 2017

%o (Magma) [Floor(n^n*Exp(-n)): n in [1..30]]; // _G. C. Greubel_, Dec 30 2017

%K nonn

%O 1,4

%A _Simon Plouffe_