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.

A073183 Sum of divisors of n that are not greater than the cubefree kernel of n.

This page as a plain text file.
%I A073183 #20 Jul 14 2018 18:33:06
%S A073183 1,3,4,7,6,12,8,7,13,18,12,28,14,24,24,7,18,39,20,42,32,36,24,36,31,
%T A073183 42,13,56,30,72,32,7,48,54,48,91,38,60,56,50,42,96,44,84,78,72,48,36,
%U A073183 57,93,72,98,54,39,72,64,80,90,60,168,62,96,104,7,84,144,68,126,96,144,72
%N A073183 Sum of divisors of n that are not greater than the cubefree kernel of n.
%C A073183 a(n) >= A073185(n).
%H A073183 Antti Karttunen, <a href="/A073183/b073183.txt">Table of n, a(n) for n = 1..10000</a>
%e A073183 The cubefree kernel of 56 = 7 * 2^3 is 28 = 7 * 2^2 and the divisors <= 28 of 56 are {1, 2, 4, 7, 8, 14, 28}, therefore a(56) = 1 + 2 + 4 + 7 + 8 + 14 + 28 = 64.
%t A073183 sdcfk[n_]:=Module[{cf=Times@@Flatten[Table[#[[1]],#[[2]]]&/@({#[[1]],If[ #[[2]]>2,2,#[[2]]]}&/@FactorInteger[n])]},Total[Select[Divisors[n],#<= cf&]]]; Array[sdcfk,80] (* _Harvey P. Dale_, Jul 14 2018 *)
%o A073183 (PARI) a007948(n) = my(f=factor(n)); for (i=1, #f~, f[i, 2] = min(f[i, 2], 2)); factorback(f);
%o A073183 a(n) = sumdiv(n, d, d*(d<=a007948(n))); \\ _Michel Marcus_, Feb 07 2015
%Y A073183 Cf. A000203, A007948, A073182, A073181.
%K A073183 nonn,look
%O A073183 1,2
%A A073183 _Reinhard Zumkeller_, Jul 19 2002