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.

A354514 Numbers k such that m - gpf(m) = k has solutions m >= 2, gpf = A006530.

This page as a plain text file.
%I A354514 #17 Aug 17 2022 05:07:33
%S A354514 0,2,3,5,6,7,9,10,11,13,14,15,17,19,20,21,22,23,24,25,26,28,29,30,31,
%T A354514 33,34,35,37,38,39,40,41,42,43,44,45,46,47,49,51,52,53,55,56,57,58,59,
%U A354514 61,62,63,65,66,67,68,69,70,71,73,74,75,76,77,78,79,82,83,85,86,87,88
%N A354514 Numbers k such that m - gpf(m) = k has solutions m >= 2, gpf = A006530.
%C A354514 Numbers k such that there is a prime p such that gpf(k+p) = p (such p must be a prime factor of n).
%C A354514 Numbers k such that there is a prime factor p of k such that k+p is p-smooth.
%C A354514 A076563 sorted and duplicates removed.
%H A354514 Jianing Song, <a href="/A354514/b354514.txt">Table of n, a(n) for n = 1..8650</a> (all terms <= 10000)
%e A354514 0 is a term because 0 = p - gpf(p) for every prime p.
%e A354514 if k/gpf(k) <= nextprime(gpf(k)) - 2, where nextprime = A151800, then k is a term since k+gpf(k) <= gpf(k)*(nextprime(gpf(k)) - 1) implies gpf(k+gpf(k)) = gpf(k).
%o A354514 (PARI) gpf(n) = vecmax(factor(n)[, 1]);
%o A354514 isA354514(n) = if(n, my(f=factor(n)[, 1]); for(i=1, #f, if(gpf(n+f[i])==f[i], return(1))); 0, 1)
%Y A354514 Cf. A006530, A076563, A151800.
%Y A354514 0 together with indices of positive terms in A354512. Complement of A354515.
%K A354514 nonn,easy
%O A354514 1,2
%A A354514 _Jianing Song_, Aug 16 2022