A327875 Numbers k such that s(k) = s(k+1) where s(k) is the sum of unitary, squarefree divisors of k, including 1 (A092261).

8, 14, 288, 675, 735, 957, 1334, 1634, 2685, 2871, 5750, 8055, 9800, 12104, 12167, 20145, 33998, 42818, 71994, 74918, 79826, 79833, 84134, 111506, 122073, 138237, 144990, 147454, 166934, 201597, 235224, 274533, 289454, 324423, 332928, 347738, 383594, 400315

Offset

1,1

Links

Amiram Eldar, Table of n, a(n) for n = 1..1000

Formula

8 is in the sequence since A092261(8) = A092261(9) = 1.

Mathematica

f[p_, e_] := If[e==1, p+1, 1]; s[n_] := Times @@ f @@@ FactorInteger[n]; s1=0; seq = {}; Do[s2 = s[n]; If[s1 == s2, AppendTo[seq, n-1]]; s1 = s2, {n, 1, 10000}]; seq

Prog

(PARI) s(n)={sumdiv(n, d, d*issquarefree(d)*(gcd(d, n/d) == 1))}
{ for(k=1, 10^6, if(s(k)==s(k+1), print1(k, ", "))) } \\ Andrew Howroyd, Sep 28 2019

Crossrefs

Cf. A002961, A064115, A064125, A092261, A164522, A164557, A293183, A306985, A324295, A325175.

Sequence in context: A219778 A180756 A279708 * A158967 A253287 A357514

Adjacent sequences: A327872 A327873 A327874 * A327876 A327877 A327878

Keyword

nonn

Author

Amiram Eldar, Sep 28 2019

Status

approved