login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A321827
a(n) = Sum_{d|n, d==1 (mod 4)} d^11 - Sum_{d|n, d==3 (mod 4)} d^11.
3
1, 1, -177146, 1, 48828126, -177146, -1977326742, 1, 31380882463, 48828126, -285311670610, -177146, 1792160394038, -1977326742, -8649707208396, 1, 34271896307634, 31380882463, -116490258898218, 48828126, 350275523038332, -285311670610
OFFSET
1,3
FORMULA
a(n) = a(A000265(n)). - M. F. Hasler, Nov 26 2018
G.f.: Sum_{k>=1} (-1)^(k-1)*(2*k - 1)^11*x^(2*k-1)/(1 - x^(2*k-1)). - Ilya Gutkovskiy, Dec 06 2018
Multiplicative with a(2^e) = 1, and for an odd prime p, ((p^11)^(e+1)-1)/(p^11-1) if p == 1 (mod 4) and ((-p^11)^(e+1)-1)/(-p^11-1) if p == 3 (mod 4). - Amiram Eldar, Sep 27 2023
a(n) = Sum_{d|n} d^11*sin(d*Pi/2). - Ridouane Oudra, Sep 08 2024
MATHEMATICA
s[n_, r_] := DivisorSum[n, #^11 &, Mod[#, 4] == r &]; a[n_] := s[n, 1] - s[n, 3]; Array[a, 30] (* Amiram Eldar, Nov 26 2018 *)
f[p_, e_] := If[Mod[p, 4] == 1, ((p^11)^(e+1)-1)/(p^11-1), ((-p^11)^(e+1)-1)/(-p^11-1)]; f[2, e_] := 1; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 60] (* Amiram Eldar, Sep 27 2023 *)
PROG
(PARI) apply( A321828(n)=sumdiv(n>>valuation(n, 2), d, (2-d%4)*d^11), [1..40]) \\ M. F. Hasler, Nov 26 2018
CROSSREFS
Column k=11 of A322143.
Cf. A321543 - A321565, A321807 - A321836 for similar sequences.
Cf. A000265.
Sequence in context: A280447 A205177 A205276 * A016763 A016775 A016847
KEYWORD
sign,easy,mult
AUTHOR
N. J. A. Sloane, Nov 24 2018
STATUS
approved