%I #29 Dec 03 2022 05:46:27
%S 1,2,3,4,5,6,7,8,9,10,12,13,14,15,16,17,18,19,20,21,24,26,27,28,30,32,
%T 34,35,36,37,38,39,40,42,45,48,51,52,54,56,57,60,63,64,65,68,70,72,73,
%U 74,76,78,80,81,84,85,90,91,95,96,97,102,104,105,108,109,111,112,114
%N Numbers with 3-smooth Euler's totient (A000010).
%C An integer n>=3 belongs to this sequence if and only if a regular n-gon can be constructed using straightedge and conic sections (details in Gibbins and Smolinsky, see below). - _Austin Shapiro_, Nov 14 2021
%C Products of 3-smooth numbers (A003586) and squarefree numbers whose prime factors are all Pierpont primes (A005109). - _Amiram Eldar_, Dec 03 2022
%H Charles R Greathouse IV, <a href="/A122254/b122254.txt">Table of n, a(n) for n = 1..10000</a>
%H Aliska Gibbins and Lawrence Smolinsky, <a href="https://www.math.lsu.edu/~smolinsk/GeometricConstructions_5_2.pdf">Geometric Constructions with Ellipses</a>, The Mathematical Intelligencer 31(1) (2009), 57-62.
%F a(n) = A048135(n-2) for n>2.
%F a(n) = A122260(n) = A048737(n) for n < 22.
%F Sum_{n>=1} 1/a(n) = 3 * Product_{p > 3 in A005109} (1 + 1/p) = 5.38288865867495675807... . - _Amiram Eldar_, Dec 03 2022
%t Select[Range@115, Max[FactorInteger[EulerPhi[#]][[All, 1]]] < 5 &] (* _Ivan Neretin_, Jul 28 2015 *)
%o (PARI) is(n)=n=eulerphi(n);n>>=valuation(n,2);n==3^valuation(n,3) \\ _Charles R Greathouse IV_, Feb 21 2013
%o (PARI) list(lim)=my(v=List(),u,t);for(i=0,log(lim--+1.5)\log(3),t=3^i;while(t<=lim,if(isprime(t+1),listput(v,t+1));t<<=1));v=vecsort(Vec(v));u=List([1]);for(i=3,#v,for(j=1,#u,t=v[i]*u[j];if(t>lim,next(2));listput(u,t)));u=vecsort(Vec(u));v=List(u);for(i=1,#u,t=u[i];while((t*=3)<=lim,listput(v,t)));u=Vec(v);v=List(u);for(i=1,#u,t=u[i];while((t<<=1)<=lim,listput(v,t)));vecsort(Vec(v)) \\ _Charles R Greathouse IV_, Feb 22 2013
%Y Cf. A000010, A003586 (3-smooth), A005109.
%Y Subsequence of A122260.
%K nonn,easy
%O 1,2
%A _Reinhard Zumkeller_, Aug 29 2006
|