%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