|
| |
|
|
A307161
|
|
Numbers n such that A307159(n) = Sum_{k=1..n} bsigma(k) is divisible by n, where bsigma(k) is the sum of bi-unitary divisors of k (A188999).
|
|
3
|
|
|
|
1, 2, 17, 37, 50, 56, 391, 919, 1399, 2829, 6249, 13664, 28829, 62272, 67195, 585391, 5504271, 6798541, 10763933, 866660818, 3830393407, 11044287758, 23058607363, 83159875881, 206501883259, 297734985607, 1087473543732, 1184060078117, 2789730557061, 2821551579466, 3529184155643
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
|
OFFSET
|
1,2
|
|
|
COMMENTS
|
The bi-unitary version of A056550.
The corresponding quotients are 1, 2, 13, 28, 38, 43, ... (see the link for more values).
a(32) > 10^13. - Giovanni Resta, May 28 2019
|
|
|
LINKS
|
Table of n, a(n) for n=1..31.
Amiram Eldar and Giovanni Resta, Table of n, a(n), A307159(a(n))/a(n) for n=1..31
|
|
|
MATHEMATICA
|
fun[p_, e_] := If[OddQ[e], (p^(e+1)-1)/(p-1), (p^(e+1)-1)/(p-1)-p^(e/2)]; bsigma[1] = 1; bsigma[n_] := Times @@ (fun @@@ FactorInteger[n]); seq={}; s = 0; Do[s = s + bsigma[n]; If[Divisible[s, n], AppendTo[seq, n]], {n, 1, 10^6}]; seq
|
|
|
CROSSREFS
|
Cf. A056550, A064611, A188999, A307043, A307159.
Sequence in context: A041031 A041965 A095075 * A294560 A258977 A069042
Adjacent sequences: A307158 A307159 A307160 * A307162 A307163 A307164
|
|
|
KEYWORD
|
nonn
|
|
|
AUTHOR
|
Amiram Eldar, Mar 27 2019
|
|
|
EXTENSIONS
|
a(23)-a(31) from Giovanni Resta, Apr 20 2019
|
|
|
STATUS
|
approved
|
| |
|
|
