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.

A328817 Numbers k such that at least 7 of k, k+1, ..., k+9 are divisible by their least prime factor squared.

This page as a plain text file.
%I A328817 #38 Jan 03 2021 01:33:02
%S A328817 3475,18271,25524,25623,45616,55772,72471,72472,104419,121667,133223,
%T A328817 133224,149220,164975,165568,165571,172916,180167,180168,203979,
%U A328817 203980,219123,260424,261472,261475,334516,334519,364216,381267,393320,393323,402723,412524,420467,420468
%N A328817 Numbers k such that at least 7 of k, k+1, ..., k+9 are divisible by their least prime factor squared.
%H A328817 David A. Corneth, <a href="/A328817/b328817.txt">Table of n, a(n) for n = 1..10000</a> (first 600 terms from Harvey P. Dale)
%e A328817 203980 is a member since 2^2 | 203980, 37^2 | 203981, 13^2 | 203983, 2^2 | 203984, 3^2 | 203985, 7^2 | 203987, and 2^2 | 203988.
%t A328817 Flatten[Position[Partition[Table[If[Divisible[n,FactorInteger[n][[1,1]]^2],1,0],{n,420000}],10,1],_?(Total[#]>6&)]] (* _Harvey P. Dale_, Jan 02 2021 *)
%o A328817 (PARI) upto(n) = {my(l, c, res); l = List([0,0,0,1,0,0,0,1,1,0]); c = 3; res = List(); for(i = 11, n, f = factor(i)[,2]; c -= l[1]; listpop(l, 1); if(f[1] >= 2, c++; listput(l, 1) , listput(l, 0) ); if(c >= 7, listput(res, i-9); ) ); res } \\ _David A. Corneth_, Jan 02 2021
%Y A328817 Cf. A283050.
%K A328817 nonn
%O A328817 1,1
%A A328817 _Sean A. Irvine_ (on behalf of _Donald S. McDonald_), Nov 07 2019