In this paper, we consider the relations between Bernoulli polynomials, Legendre polynomials and combinatoric convolution sums of divisor functions. In addition, we give examples of approximate normal distribution derived from combinatoric convolution sums of divisor functions.