%I #44 Jun 14 2018 17:51:29
%S 0,0,0,1,1,1,2,2,2,3,3,3,4,4,4,4,4,4,5,5,5,5,5,5,6,6,6,7,7,7,7,7,7,7,
%T 7,7,8,8,8,8,8,8,8,8,8,8,8,8,9,9,9,9,9,9,10,10,10,10,10,10,10,10,10,
%U 10,10,10,10,10,10,10,10,10,11,11,11,11,11,11,11
%N Number of integers less than or equal to n whose largest prime factor is 3.
%C a(n) increases when n has the form 2^a*3^b, with a >= 0 and b > 0.
%C A distinct sequence can be generated for each prime number; this sequence is for the prime number 3. For an example using another prime number see A301506.
%F From _David A. Corneth_, Mar 27 2018 (Start)
%F a(n) - a(n - 1) = 1 if and only if n is in 3 * A003586. If n isn't in that sequence then a(n) = a(n - 1).
%F a(3 * n + b) = A071521(n), n > 0, 0 <= b < 3. (End)
%e a(12) = a(2^2 * 3^1); 3 is the largest prime factor, so a(12) exceeds the previous term by 1. For a(13), 13 is a prime, so there is no increase from the previous term.
%t Accumulate@ Array[Boole[FactorInteger[#][[-1, 1]] == 3] &, 80, 0] (* _Michael De Vlieger_, Apr 21 2018 *)
%o (MATLAB)
%o clear;clc;
%o prime = 3;
%o limit = 10000;
%o largest_divisor = ones(1,limit+1);
%o for k = 0:limit
%o f = factor(k);
%o largest_divisor(k+1) = f(end);
%o end
%o for i = 1:limit+1
%o FQN(i) = sum(largest_divisor(1:i)==prime);
%o end
%o output = [0:limit;FQN]'
%o (PARI) gpf(n) = if (n<=1, n, vecmax(factor(n)[,1]));
%o a(n) = sum(k=1, n, gpf(k)==3); \\ _Michel Marcus_, Mar 27 2018
%Y Cf. A003586, A065119, A071521.
%Y Cf. A301506.
%K nonn
%O 0,7
%A _Ralph-Joseph Tatt_, Mar 21 2018
|