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.

A261614 Numbers that are neither prime (A000040) nor practical (A005153).

Original entry on oeis.org

9, 10, 14, 15, 21, 22, 25, 26, 27, 33, 34, 35, 38, 39, 44, 45, 46, 49, 50, 51, 52, 55, 57, 58, 62, 63, 65, 68, 69, 70, 74, 75, 76, 77, 81, 82, 85, 86, 87, 91, 92, 93, 94, 95, 98, 99, 102, 105, 106, 110, 111
Offset: 1

Views

Author

Frank M Jackson, Nov 18 2015

Keywords

Examples

			a(5)=21 and it is neither prime nor practical. It is the 5th such occurrence.
		

Crossrefs

Programs

  • Mathematica
    PracticalQ[n_] := Module[{f, p, e, prod = 1, ok = True}, If[n<1||(n>1&&OddQ[n]), False, If[n==1, True, f=FactorInteger[n]; {p, e}=Transpose[f]; Do[If[p[[i]]>1+DivisorSigma[1, prod], ok=False; Break[]]; prod=prod*p[[i]]^e[[i]], {i, Length[p]}]; ok]]]; Select[Range[1, 1000], ! PracticalQ[#] && ! PrimeQ[#] &] (* using T. D. Noe's program A005153 *)
  • PARI
    isok(n) = bittest(n, 0) && return(n==1); my(P=1); n && !for(i=2, #n=factor(n)~, n[1, i]>1+(P*=sigma(n[1, i-1]^n[2, i-1])) && return) ;
    for(n=1, 200, if(!isok(n) && !isprime(n), print1(n, ", "))) \\ Altug Alkan, Nov 19 2015