A135205 Numbers m for which Sum_digits(m!!) is a multiple of Sum_digits(m).
1, 2, 3, 4, 6, 9, 10, 11, 12, 15, 18, 20, 21, 24, 25, 27, 30, 32, 33, 36, 42, 45, 46, 54, 55, 63, 72, 75, 81, 88, 90, 91, 93, 100, 101, 102, 105, 108, 111, 112, 117, 120, 121, 122, 123, 124, 126, 127, 135, 141, 144, 153, 154, 156, 162, 171, 176, 180, 182, 189, 198
Offset: 1
Examples
11 -> 11*9*7*5*3*1=10395 -> (1+0+3+9+5)/(1+1) = 9.
Links
- G. C. Greubel, Table of n, a(n) for n = 1..1000
Programs
-
Maple
P:=proc(n) local i,j,k,w,x; for i from 1 by 1 to n do w:=0; k:=i; while k>0 do w:=w+k-(trunc(k/10)*10); k:=trunc(k/10); od; x:=i; j:=i-2; while j >0 do x:=x*j; j:=j-2; od: k:=x; x:=0; while k>0 do x:=x+k-(trunc(k/10)*10); k:=trunc(k/10); od; if trunc(x/w)=x/w then print(i); fi; od; end: P(1000);
-
Mathematica
Select[Range[100], Divisible[Total[IntegerDigits[#!!, 10]], Total[IntegerDigits[#, 10]]] &] (* G. C. Greubel, Sep 30 2016 *)
Extensions
Offset 1 and b-file adapted by Paolo P. Lava, Jun 17 2024