|
%I #18 May 20 2019 06:34:25
%S 3,5,7,11,13,17,19,29,31,37,41,43,53,61,67,73,79,89,97,101,103,109,
%T 113,127,137,139,149,151,157,163,173,181,193,197,199,211,223,229,233,
%U 241,257,269,271,277,281,283,293,307,313,317,331,337,349,353,367,373
%N Primes represented by cyclotomic binary forms.
%C A cyclotomic binary form over Z is a homogeneous polynomial in two variables which has the form f(x, y) = y^EulerPhi(k)*CyclotomicPolynomial(k, x/y) where k is some integer >= 3. An integer n is represented by f if f(x, y) = n has an integer solution.
%H Peter Luschny, <a href="/A325143/b325143.txt">Table of n, a(n) for n = 1..10000</a>
%H Étienne Fouvry, Claude Levesque, Michel Waldschmidt, <a href="http://dx.doi.org/10.4064/aa171012-24-12">Representation of integers by cyclotomic binary forms</a>, Acta Arithmetica 184 (2018), 67-86; <a href="https://arxiv.org/abs/1712.09019">arXiv:1712.09019</a>, arXiv:1712.09019 [math.NT], 2017.
%o (Julia) using Nemo
%o function isA325143(n)
%o (n < 3 || !isprime(ZZ(n))) && return false
%o R, x = PolynomialRing(ZZ, "x")
%o K = floor(Int, 5.383*log(n)^1.161) # Bounds from
%o M = floor(Int, 2*sqrt(n/3)) # Fouvry & Levesque & Waldschmidt
%o N = QQ(n)
%o for k in 3:K
%o e = Int(eulerphi(ZZ(k)))
%o c = cyclotomic(k, x)
%o for m in 1:M, j in 0:M if max(j, m) > 1
%o N == m^e*subst(c, QQ(j,m)) && return true
%o end end end
%o return false
%o end
%o [n for n in 1:373 if isA325143(n)] |> println
%Y Subsequence of A296095. Complement A325145. Number of A325141.
%Y Cf. A293654, A299214, A299498, A299733, A299928, A299930, A299956, A299964.
%K nonn
%O 1,1
%A _Peter Luschny_, May 16 2019
|
