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.

A177899 Nonsquarefree numbers that are not in A177880.

This page as a plain text file.
%I A177899 #16 Aug 31 2020 02:40:34
%S A177899 8,16,24,27,40,48,54,56,80,81,88,104,112,120,125,135,136,152,162,168,
%T A177899 176,184,189,208,216,232,240,248,250,264,270,272,280,296,297,304,312,
%U A177899 328,336,343,344,351,368,375,376,378,405,408,424,432,440,456,459,464,472,488,496,512,513,520,528,536
%N A177899 Nonsquarefree numbers that are not in A177880.
%H A177899 Amiram Eldar, <a href="/A177899/b177899.txt">Table of n, a(n) for n = 1..10000</a>
%F A177899 Let B(x) be the counting function for terms not exceeding x. Then for x tends to infinity, B(x)=C*x+o(x^(0.5+eps), where C = Product_{i=p^(3^k) with prime p and k>=0}(1-1/(i^2+i+1)) - 1/zeta(2).
%p A177899 isA005836 := proc(n) convert(convert(n,base,3),set) intersect {2} ; % = {} ; end proc:
%p A177899 isA177880 := proc(n) local f; for f in ifactors(n)[2]  do if not isA005836(op(2,f)) then return true; end if;  end do: return false; end proc:
%p A177899 isA177899 := proc(n) not numtheory[issqrfree](n) and not isA177880(n) ; end proc:
%p A177899 for n from 1 to 1000 do if isA177899(n) then printf("%d,",n) ; end if; end do: # _R. J. Mathar_, Dec 20 2010
%t A177899 Select[Range[500], AnyTrue[(e = FactorInteger[#][[;; , 2]]), #1 > 1 &] && AllTrue[e, DigitCount[#1, 3, 2] == 0 &] &] (* _Amiram Eldar_, Aug 31 2020 *)
%Y A177899 Cf. A013929, A177880.
%K A177899 nonn
%O A177899 1,1
%A A177899 _Vladimir Shevelev_, Dec 15 2010