# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/ Search: id:a003973 Showing 1-1 of 1 %I A003973 #60 Dec 08 2022 07:36:26 %S A003973 1,4,6,13,8,24,12,40,31,32,14,78,18,48,48,121,20,124,24,104,72,56,30, %T A003973 240,57,72,156,156,32,192,38,364,84,80,96,403,42,96,108,320,44,288,48, %U A003973 182,248,120,54,726,133,228,120,234,60,624,112,480,144,128,62,624,68 %N A003973 Inverse Möbius transform of A003961; a(n) = sigma(A003961(n)), where A003961 shifts the prime factorization of n one step towards the larger primes. %C A003973 Sum of the divisors of the prime shifted n, or equally, sum of the prime shifted divisors of n. - _Antti Karttunen_, Aug 17 2020 %H A003973 Charles R Greathouse IV, Table of n, a(n) for n = 1..10000 %H A003973 Index entries for sequences computed from indices in prime factorization. %H A003973 Index entries for sequences related to sigma(n). %F A003973 Multiplicative with a(p^e) = (q^(e+1)-1)/(q-1) where q = nextPrime(p). - _David W. Wilson_, Sep 01 2001 %F A003973 From _Antti Karttunen_, Aug 06-12 2020: (Start) %F A003973 a(n) = Sum_{d|n} A003961(d) = Sum_{d|A003961(n)} d. %F A003973 a(n) = A000203(A003961(n)) = A000593(A003961(n)). %F A003973 a(n) = 2*A336840(n) - A000005(n) = 2*Sum_{d|n} (A048673(d) - (1/2)). %F A003973 a(n) = A008438(A108228(n)) = A008438(A048673(n)-1). %F A003973 a(n) = A336838(n) * A336856(n). %F A003973 a(n) is odd if and only if n is a square. %F A003973 (End) %F A003973 Sum_{k=1..n} a(k) ~ c * n^3, where c = (1/2) * Product_{p prime} p^3/((p+1)*(p^2-nextprime(p))) = 3.39513795..., where nextprime is A151800. - _Amiram Eldar_, Dec 08 2022 %t A003973 b[1] = 1; b[p_?PrimeQ] := b[p] = Prime[ PrimePi[p] + 1]; b[n_] := b[n] = Times @@ (b[First[#]]^Last[#] &) /@ FactorInteger[n]; a[n_] := Sum[ b[d], {d, Divisors[n]}]; Table[a[n], {n, 1, 70}] (* _Jean-François Alcover_, Jul 18 2013 *) %o A003973 (PARI) aPrime(p,e)=my(q=nextprime(p+1));(q^(e+1)-1)/(q-1) %o A003973 a(n)=my(f=factor(n));prod(i=1,#f~,aPrime(f[i,1],f[i,2])) \\ _Charles R Greathouse IV_, Jul 18 2013 %o A003973 (PARI) A003973(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); sigma(factorback(f)); }; \\ _Antti Karttunen_, Aug 06 2020 %o A003973 (Python) %o A003973 from math import prod %o A003973 from sympy import factorint, nextprime %o A003973 def A003973(n): return prod(((q:=nextprime(p))**(e+1)-1)//(q-1) for p,e in factorint(n).items()) # _Chai Wah Wu_, Jul 05 2022 %Y A003973 Cf. A000203, A000290 (positions of odd terms), A003961, A007814, A048673, A108228, A151800, A295664, A336840. %Y A003973 Permutation of A008438. %Y A003973 Used in the definitions of the following sequences: A326042, A336838, A336841, A336844, A336846, A336847, A336848, A336849, A336850, A336851, A336852, A336856, A336931, A336932. %Y A003973 Cf. also A003972. %K A003973 nonn,easy,mult %O A003973 1,2 %A A003973 _Marc LeBrun_ %E A003973 More terms from _David W. Wilson_, Aug 29 2001 %E A003973 Secondary name added by _Antti Karttunen_, Aug 06 2020 # Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE