|
%I #42 Jun 24 2022 20:26:17
%S 1,3,4,11,6,42,8,75,31,108,12,1778,14,206,234,1099,18,5901,20,8116,
%T 452,498,24,333618,131,692,760,22166,30,810372,32,33867,1104,1176,
%U 1238,10085333,38,1466,1538,2568180,42,3112382,44,85690,91386,2142,48,255138610
%N a(n) = Sum_{d|n} A007955(d) where A007955(m) = product of divisors of m.
%H Seiichi Manyama, <a href="/A175317/b175317.txt">Table of n, a(n) for n = 1..10000</a>
%F From _Bernard Schott_, Oct 26 2021: (Start)
%F a(1) = 1 (the only fixed point).
%F a(p) = p+1 for prime p only.
%F a(2^k) = A181388(k+1). (End)
%e For n = 4, with b(n) = A007955(n), a(4) = b(1) + b(2) + b(4) = 1 + 2 + 8 = 11.
%t a[n_] := DivisorSum[n, #^(DivisorSigma[0, #]/2) &]; Array[a, 50] (* _Amiram Eldar_, Oct 23 2021 *)
%o (PARI) a(n) = sumdiv(n, d, vecprod(divisors(d))); \\ _Michel Marcus_, Dec 09 2014 and Oct 23 2021
%o (Python)
%o from math import isqrt
%o from sympy import divisor_count, divisors
%o def A175317(n): return sum(isqrt(d)**c if (c:=divisor_count(d)) & 1 else d**(c//2) for d in divisors(n,generator=True)) # _Chai Wah Wu_, Jun 24 2022
%Y Cf. A007429, A007955, A206032, A266265.
%Y Subsequences: A008864, A181388 \ {0}.
%K nonn
%O 1,2
%A _Jaroslav Krizek_, Apr 01 2010
%E Corrected by _Jaroslav Krizek_, Apr 02 2010
%E Edited and more terms from _Michel Marcus_, Dec 09 2014
|
