login
Number of solutions m >= 2 to m - gpf(m) = n, gpf = A006530.
7

%I #25 Aug 17 2022 05:07:13

%S 0,1,1,0,1,2,1,0,1,1,1,0,1,2,2,0,1,0,1,1,2,1,1,1,1,1,0,1,1,1,1,0,2,1,

%T 2,0,1,1,1,1,1,1,1,1,2,1,1,0,1,0,2,1,1,0,2,1,1,1,1,0,1,2,1,0,1,1,1,1,

%U 2,1,1,0,1,1,1,1,2,2,1,0,0,1,1,0,2,1,1,1,1,0,2

%N Number of solutions m >= 2 to m - gpf(m) = n, gpf = A006530.

%C Number of primes p such that gpf(n+p) = p (such p must be prime factors of n).

%C Number of distinct prime factors p of n such that n+p is p-smooth.

%C Clearly we have a(n) <= omega(n) for all n, omega = A001221. The differences are given by A354527.

%C Is this sequence unbounded? Note that 4 does not appear until a(1660577).

%H Jianing Song, <a href="/A354512/b354512.txt">Table of n, a(n) for n = 1..10000</a>

%e a(78) = 2 since the prime factors of 78 are 2,3,13, and we have gpf(78+3) = 3 and gpf(78+13) = 13, so the solutions to m - gpf(m) = 78 are m = 78+3 = 81 or m = 78+13 = 91. Note that gpf(78+2) != 2.

%e a(12) = 0 since the prime factors of 12 are 2,3, and we have gpf(12+2) != 2 and gpf(12+3) != 3.

%o (PARI) gpf(n) = vecmax(factor(n)[, 1]);

%o a(n) = my(f=factor(n)[, 1]); sum(i=1, #f, gpf(n+f[i])==f[i])

%Y Cf. A006530, A076563, A001221, A354516 (indices of first occurrence of each number), A354527.

%Y Cf. A354514 (0 together with indices of positive terms), A354515 (indices of 0), A354516, A354525 (indices n for which a(n) reaches omega(n)), A354526 (indices n for which a(n) is smaller than omega(n)).

%K nonn,easy

%O 1,6

%A _Jianing Song_, Aug 16 2022