cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A086670 Sum of floor(d/2) where d is a divisor of n.

This page as a plain text file.
%I A086670 #28 Feb 13 2023 16:06:23
%S A086670 0,1,1,3,2,5,3,7,5,8,5,13,6,11,10,15,8,18,9,20,14,17,11,29,14,20,18,
%T A086670 27,14,34,15,31,22,26,22,44,18,29,26,44,20,46,21,41,36,35,23,61,27,45,
%U A086670 34,48,26,58,34,59,38,44,29,82,30,47,49,63,40,70,33,62,46,70,35,96,36,56
%N A086670 Sum of floor(d/2) where d is a divisor of n.
%C A086670 Inverse Mobius transform of A004526. - _R. J. Mathar_, Jan 19 2009
%F A086670 G.f.: Sum_{n>=1} floor(n/2)*x^n/(1-x^n). - _Joerg Arndt_, Jan 30 2011
%F A086670 a(n) = (A000203(n) - A001227(n)) / 2. - _Franklin T. Adams-Watters_, Jan 05 2012
%F A086670 G.f.: Sum_{k>=1} x^(2*k) / ((1 + x^k) * (1 - x^k)^2). - _Ilya Gutkovskiy_, Aug 02 2021
%F A086670 a(n) = Sum_{i=1..floor(n/2)} A135539(n,2*i). - _Ridouane Oudra_, Apr 15 2022
%e A086670 10 has divisors 1,2,5,10. floor(d/2) gives 0,1,2,5, therefore a(10)=8.
%t A086670 Table[Total[Floor[Divisors[n]/2]],{n,80}] (* _Harvey P. Dale_, Feb 13 2023 *)
%o A086670 (PARI) for (n=1,100,s=0; fordiv(i=n,i,s+=floor(i/2)); print1(","s))
%o A086670 (PARI) a(n) = my(f = factor(n)); (sigma(f) - (numdiv(f)/(valuation(n, 2)+1)))>>1 \\ _David A. Corneth_, Apr 15 2022 using _Franklin T. Adams-Watters_'s formula
%Y A086670 Cf. A000203, A001227, A004526, A079247.
%Y A086670 Cf. A135539.
%K A086670 nonn
%O A086670 1,4
%A A086670 _Jon Perry_, Jul 27 2003