%I #8 Sep 10 2019 21:49:38
%S 1,2,3,4,5,6,7,8,9,1,3,9,2,3,6,9,3,9,4,6,9,12,5,9,6,9,7,9,8,9,12,9,1,
%T 3,9,2,3,4,6,9,3,9,6,7,9,5,9,6,8,9,12,9,9,9,10,12,15,18,2,3,6,9,11,3,
%U 9,4,6,8,9,12,5,9,6,9,13,9,9,12,14,9,10,15
%N a(n) is the digit sum of the n-th Niven number (or Harshad number).
%F a(n) = A005349(n)/A113315(n).
%F a(n) = A007953(A005349(n)).
%o (PARI) is_a005349(n) = n%sumdigits(n)==0
%o terms(n) = my(i=0); for(k=1, oo, if(i>=n, break); if(is_a005349(k), print1(sumdigits(k), ", "); i++))
%o terms(80) \\ Print initial 80 terms
%Y Cf. A005349, A007953, A113315.
%K nonn,base,easy
%O 1,2
%A _Felix Fröhlich_, Sep 06 2019