cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A377734 Number of integers less than n that have the same smallest prime factor as n.

This page as a plain text file.
%I A377734 #14 Feb 16 2025 08:34:07
%S A377734 0,0,0,1,0,2,0,3,1,4,0,5,0,6,2,7,0,8,0,9,3,10,0,11,1,12,4,13,0,14,0,
%T A377734 15,5,16,2,17,0,18,6,19,0,20,0,21,7,22,0,23,1,24,8,25,0,26,3,27,9,28,
%U A377734 0,29,0,30,10,31,4,32,0,33,11,34,0,35,0,36,12,37,2,38,0,39
%N A377734 Number of integers less than n that have the same smallest prime factor as n.
%H A377734 Amiram Eldar, <a href="/A377734/b377734.txt">Table of n, a(n) for n = 1..10000</a>
%H A377734 Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/LeastPrimeFactor.html">Least Prime Factor</a>.
%F A377734 a(n) = |{j < n : lpf(j) = lpf(n)}|.
%F A377734 a(n) = A078898(n) - 1.
%F A377734 Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = Sum_{k>=1} (A038110(k)/A038111(k))^2 = 0.2847976823663... . - _Amiram Eldar_, Nov 21 2024
%t A377734 Table[Length[Select[Range[n - 1], If[# == 1, 1, FactorInteger[#][[1, 1]]] == If[n == 1, 1, FactorInteger[n][[1, 1]]] &]], {n, 80}]
%t A377734 seq[len_] := Module[{t = Table[FactorInteger[n][[1,1]], {n, 1, len}], s = Table[0, {len}]}, Do[s[[i]] = Count[t[[1;;i-1]], t[[i]]], {i, 1, len}]; s]; seq[80] (* _Amiram Eldar_, Nov 21 2024 *)
%o A377734 (PARI) a(n) = if (n>1, my(p=vecmin(factor(n)[,1])); sum(k=2, n-1, p == vecmin(factor(k)[,1])), 0); \\ _Michel Marcus_, Nov 16 2024
%Y A377734 Cf. A020639, A038110, A038111, A047983, A078898, A334655, A335097, A377730.
%K A377734 nonn
%O A377734 1,6
%A A377734 _Ilya Gutkovskiy_, Nov 05 2024