|
%I #11 Jun 12 2018 21:14:51
%S 0,0,0,2,0,0,0,2,3,0,0,0,0,0,0,6,0,0,0,0,0,0,0,0,5,0,3,0,0,0,0,2,0,0,
%T 0,6,0,0,0,0,0,0,0,0,0,0,0,0,7,0,0,0,0,0,0,0,0,0,0,0,0,0,0,14,0,0,0,0,
%U 0,0,0,0,0,0,0,0,0,0,0,0,12,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,10
%N a(n) is the sum of perfect divisors of n - n, where a perfect divisor of n is a divisor d such that d^k = n for some k >= 1.
%C a(1) = 0, for n >=2: a(n) = sum of perfect divisors of n less than n.
%C a(n) > 0 for perfect powers n = A001597(m) for m > 2.
%H Antti Karttunen, <a href="/A175070/b175070.txt">Table of n, a(n) for n = 1..65537</a>
%H <a href="/index/Su#sums_of_divisors">Index entries for sequences related to sums of divisors</a>
%F a(n) = A175067(n) - n.
%o (PARI) A175070(n) = if(!ispower(n),0,sumdiv(n,d,if((d>1)&&(d<n)&&((d^valuation(n,d))==n),d,0))); \\ _Antti Karttunen_, Jun 12 2018
%o (PARI) first(n) = {my(res = vector(n)); for(i = 2, sqrtint(n), for(j = 2, logint(n, i), res[i^j] += i)); res} \\ _David A. Corneth_, Jun 12 2018
%Y Cf. A175067.
%K nonn
%O 1,4
%A _Jaroslav Krizek_, Jan 23 2010
|
