A037872 - OEIS

%I #8 Oct 19 2015 14:30:55

%S 0,0,0,1,0,0,0,2,1,0,0,3,2,1,0,1,1,1,1,1,0,0,0,2,1,0,0,3,2,1,0,2,2,2,

%T 2,2,1,1,1,2,1,0,0,3,2,1,0,3,3,3,3,3,2,2,2,3,2,1,1,3,2,1,0,1,1,1,1,1,

%U 1,1,1,2,2,1,1,3,3,2,1,2,1,1,1,1,0,0,0,2,1,0

%N (1/2)*Sum{|d(i)-e(i)|}, where Sum{d(i)*4^i} is base 4 representation of n and e(i) are digits d(i) in nonincreasing order, for i=0,1,...,m.

%p A037872 := proc(n)

%p     local a,dgs,dgsE ;

%p     a := 0 ;

%p     dgs := convert(n,base,4);

%p     dgsE := sort(dgs,`>`) ;

%p     for i from 1 to nops(dgs) do

%p         a := a+ abs(op(i,dgs)-op(i,dgsE)) ;

%p     end do:

%p     a/2 ;

%p end proc: # _R. J. Mathar_, Oct 19 2015

%K nonn,base

%O 1,8

%A _Clark Kimberling_

%E Definition swapped with A037881. - _R. J. Mathar_, Oct 19 2015