A298759 Numbers k such that bphi(k) = k/2, where bphi is the bi-unitary analog of Euler's totient function (A116550).

2, 6, 30, 42, 1722, 1806, 19977474

Offset

1,1

Comments

With Euler's totient function, phi(k) = k/2 only for powers of 2 (A000079, except for 1). With the unitary totient function (A047994) the corresponding sequence is A030163.

a(8) > 2*10^9, if it exists. - Amiram Eldar, Jul 16 2022

Links

Table of n, a(n) for n=1..7.

Example

42 is in the sequence since bphi(42) = 21 = 42/2.

Mathematica

bphi[1] = 1; bphi[n_] :=  With[{pp = Power @@@ FactorInteger[n]},   Count[Range[n], m_ /; Intersection[pp, Power @@@ FactorInteger[m]] == {}]]; aQ[n_] := bphi[n] == n/2; Select[Range[10000], aQ]

Prog

(PARI) udivs(n) = {my(d = divisors(n)); select(x->(gcd(x, n/x)==1), d); }
gcud(n, m) = vecmax(setintersect(udivs(n), udivs(m)));
bphi(n) = if (n==1, 1, sum(k=1, n-1, gcud(n, k) == 1));
isok(n) = bphi(n) == n/2; \\ Michel Marcus, Jan 26 2018

Crossrefs

Cf. A047994, A030163, A116550.

Sequence in context: A100194 A229882 A325986 * A360680 A369685 A127517

Adjacent sequences: A298756 A298757 A298758 * A298760 A298761 A298762

Keyword

nonn,more

Author

Amiram Eldar, Jan 26 2018

Extensions

a(7) from Amiram Eldar, Jul 16 2022

Status

approved