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.

A010847 Number of numbers <= n with a prime factor that does not divide n.

This page as a plain text file.
%I A010847 #20 Jun 16 2016 07:17:36
%S A010847 0,0,1,1,3,1,5,4,6,4,9,4,11,8,10,11,15,8,17,12,16,15,21,13,22,19,23,
%T A010847 20,27,12,29,26,27,26,30,22,35,30,33,29,39,23,41,35,37,38,45,33,46,38,
%U A010847 45,43,51,38,50,45,51,50,57,34,59,54,55,57,60,44,65,58,63,50,69,54,71
%N A010847 Number of numbers <= n with a prime factor that does not divide n.
%H A010847 Alois P. Heinz, <a href="/A010847/b010847.txt">Table of n, a(n) for n = 1..10000</a>
%F A010847 a(p) = p-2, for p prime; that is all numbers between 2 and p-1 inclusive. - _Michel Marcus_, May 31 2014
%F A010847 a(n) = n - A010846(n). - _Anthony Browne_, Jun 07 2016
%e A010847 For n=5, the three numbers 2,3 and 4 have a prime factor that is not found in 5. Hence a(5) = 3.
%t A010847 Table[Sum[1-Floor[n^k/k]+Floor[(n^k-1)/k], {k,n}],{n,100}] (* _Anthony Browne_, Jun 07 2016 *)
%o A010847 (PARI) a(n) = {pfn = factor(n)[,1]~; nb = 0; for (i=2, n, pfi = factor(i)[,1]~; for (j=1, #pfi, if (! vecsearch(pfn, pfi[j]), nb++; break););); nb;} \\ _Michel Marcus_, May 31 2014
%Y A010847 Cf. A010846.
%K A010847 nonn,easy
%O A010847 1,5
%A A010847 _Olivier Gérard_