|
%I #11 Jan 24 2022 17:28:08
%S 1,1,4,36,144,3600,3600,176400,2822400,25401600,25401600,3073593600,
%T 110649369600,18699743462400,74798973849600,1869974346240000,
%U 29919589539840000,8646761377013760000,77820852393123840000,28093327713917706240000,112373310855670824960000
%N Smallest square divisible by n!.
%H Robert Israel, <a href="/A065886/b065886.txt">Table of n, a(n) for n = 0..407</a>
%F a(n) = A053143(A000142(n)) = A065887(n)^2 = A000142(n)*A055204(n) = A001044(n)/A055071(n)
%e a(10) = 25401600 since 10! = 3628800 and the smallest square divisible by this is 25401600 = 3628800*7 = 5040^2
%p N:= 50: # to get a(0)..a(N)
%p P:= select(isprime, [$2..N]):
%p nP:= nops(P):
%p V:= Vector(nP):
%p A[0]:= 1:
%p for n from 1 to N do
%p for i from 1 to nP do V[i]:= V[i] + padic:-ordp(n,P[i]) od;
%p A[n]:= mul(P[i]^(2*ceil(V[i]/2)),i=1..nP)
%p od:
%p seq(A[n],n=0..N); # _Robert Israel_, Jan 30 2017
%t ssd[n_]:=Module[{nf=n!,k=1},While[!IntegerQ[Sqrt[k*nf]],k++];k*nf]; Array[ssd,20,0] (* _Harvey P. Dale_, Apr 29 2012 *)
%K nonn
%O 0,3
%A _Henry Bottomley_, Nov 27 2001
|
