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.

A387725 Number of unitary divisors d of n for which A107758(d) > 2*d, where A107758 is sigma+, multiplicative function with a(p^e) = 1+sigma(p^e).

This page as a plain text file.
%I A387725 #11 Sep 06 2025 21:40:37
%S A387725 0,0,0,0,0,1,0,0,0,1,0,1,0,1,1,0,0,1,0,1,1,1,0,1,0,1,0,1,0,4,0,0,0,1,
%T A387725 0,1,0,1,0,1,0,4,0,1,1,1,0,1,0,1,0,1,0,1,0,1,0,1,0,4,0,1,0,0,0,3,0,1,
%U A387725 0,3,0,1,0,1,1,1,0,3,0,1,0,1,0,4,0,1,0,1,0,4,0,1,0,1,0,1,0,1,0,1,0,3,0,1,3
%N A387725 Number of unitary divisors d of n for which A107758(d) > 2*d, where A107758 is sigma+, multiplicative function with a(p^e) = 1+sigma(p^e).
%H A387725 Antti Karttunen, <a href="/A387725/b387725.txt">Table of n, a(n) for n = 1..65537</a>
%F A387725 a(n) = Sum_{d|n} [gcd(d,n/d)==1 and A107758(d) > 2*d], where [ ] is the Iverson bracket.
%o A387725 (PARI)
%o A387725 A107758(n) =  { my(f = factor(n)); prod(k=1, #f~, 1+sigma(f[k, 1]^f[k, 2])); };
%o A387725 A387725(n) = sumdiv(n,d,(1==gcd(d,n/d)) && (A107758(d)>(2*d)));
%Y A387725 Cf. A107758, A387721 (positions of positive terms).
%Y A387725 Cf. also A387715, A337345.
%K A387725 nonn,new
%O A387725 1,30
%A A387725 _Antti Karttunen_, Sep 06 2025