A306950 -id:A306950 - cp's OEIS frontend

A356747 Numbers m that divide A306070(m) = Sum_{k=1..m} bphi(k), where bphi is the bi-unitary totient function (A116550).

1, 2, 141, 1035, 2388, 3973, 5157, 14160, 37023, 68861, 99889, 116106, 117939, 627400, 1561944, 1626983, 5901444, 10054091, 12260525, 32619981, 49775099

Offset: 1

Amiram Eldar, Aug 25 2022

The corresponding quotients A306070(m)/m are 1, 1, 57, 418, ... (see the link for more values).

a(22) > 6.5*10^8, if it exists.

Amiram Eldar, Table of n, a(n), A306070(a(n))/a(n) for n = 1..21.

Cf. A116550, A306070.

Similar sequences: A048290, A306950.

phi[x_, n_] := DivisorSum[n, MoebiusMu[#]*Floor[x/#] &]; bphi[n_] := DivisorSum[n, (-1)^PrimeNu[#]*phi[n/#, #] &, CoprimeQ[#, n/#] &]; seq = {}; s = 0; Do[s = s + bphi[n]; If[Divisible[s, n], AppendTo[seq, n]], {n, 1, 10^6}]; seq

A309272 Numbers m such that m divides A173290(m) = Sum_{k=1..m} psi(k), where psi is the Dedekind psi function (A001615).

Original entry on oeis.org

1, 2, 5, 15, 31, 40, 66, 81, 315, 966, 1398, 1768, 30166, 32335, 98734, 388033, 591597, 1375056, 14966304, 15160528, 50793208, 51302236, 99253376, 110994356, 230465053, 402340268, 497982399, 2027319577, 2879855394, 18450762682, 29922126368, 31711273834, 40583934786

Offset: 1

Amiram Eldar, Oct 23 2019

nonn

The corresponding quotients are 1, 2, 4, 12, 24, 31, 51, 62, 240, 735, 1063, 1344, 22924, 24572, 75029, 294870, 449560, 1044918, 11373028, 11520620, 38598210, 38985025, 75423522, 84345597, 175132440, 305741942, 378421246, 1540578144, 2188427680, 14020898356, 22738089456, 24097678498, 30840092321, ...

			2 is in the sequence since psi(1) + psi(2) = 1 + 3 = 4 is divisible by 2.
5 is in the sequence since psi(1) + psi(2) + ... + psi(5) = 1 + 3 + 4 + 6 + 6 = 20 is divisible by 5.

Cf. A001615, A173290.

Cf. A045345, A048290, A050226, A056550, A063971, A064605, A064606, A064607, A064610, A064611, A064612, A306950, A307043, A307161, A326488.

psi[1] = 1; psi[n_] := n * Times @@ (1 + 1/Transpose[FactorInteger[n]][[1]]); seq = {}; s = 0; Do[s += psi[n]; If[Divisible[s, n], AppendTo[seq, n]], {n, 1, 10^4}]; seq

a(31)-a(33) from Giovanni Resta, Oct 24 2019

cp's OEIS Frontend

A356747 Numbers m that divide A306070(m) = Sum_{k=1..m} bphi(k), where bphi is the bi-unitary totient function (A116550).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica