cp's OEIS Frontend

A093811 Sum of the digital products of the divisors of n.

%I A093811 #12 Jun 02 2022 20:43:23
%S A093811 1,3,4,7,6,12,8,15,13,8,2,18,4,14,14,21,8,29,10,12,13,8,7,34,16,18,27,
%T A093811 34,19,22,4,27,14,22,28,53,22,36,34,20,5,33,13,28,43,33,29,72,44,18,
%U A093811 16,32,16,63,32,72,48,61,46,28,7,18,40,51,39,62,43,74,64,34,8,83,22,52,59
%N A093811 Sum of the digital products of the divisors of n.
%C A093811 The first few n such that a(n) = n are: 1,14,27,156,196. Are there any more?
%C A093811 Inverse Moebius transform of A007954(n). a(n) = A007954(n) * A000012(n), where operation * denotes Dirichlet convolution for n >= 1. Dirichlet convolution of functions b(n), c(n) is function a(n) = b(n) * c(n) = Sum_{d|n} b(d)*c(n/d). Simultaneously holds Dirichlet multiplication: a(n) * A008683(n) = A007954(n). [From _Jaroslav Krizek_, Mar 22 2009]
%H A093811 Harvey P. Dale, <a href="/A093811/b093811.txt">Table of n, a(n) for n = 1..1000</a>
%e A093811 a(1234)=69 because the divisors of 1234 are: [1, 2, 617, 1234] and
%e A093811 1+2+(6*1*7)+(1*2*3*4) = 69.
%p A093811 A093811 := proc(n::integer)
%p A093811     local a,d ;
%p A093811     a := 0 ;
%p A093811     for d in numtheory[divisors](n) do
%p A093811         a := a+A007954(d) ;
%p A093811     end do:
%p A093811 end proc: # _R. J. Mathar_, Oct 02 2019
%t A093811 Table[Total[Times@@IntegerDigits[#]&/@Divisors[n]],{n,100}] (* _Harvey P. Dale_, Jun 02 2022 *)
%Y A093811 Cf. A007954, A008683.
%K A093811 base,easy,nonn
%O A093811 1,2
%A A093811 _Jason Earls_, May 20 2004