login
Search: a001269 -id:a001269
     Sort: relevance | references | number | modified | created      Format: long | short | data
Number of prime factors of 2^n + 1 (counted with multiplicity).
+10
29
1, 1, 2, 1, 2, 2, 2, 1, 4, 3, 2, 2, 2, 3, 4, 1, 2, 4, 2, 2, 4, 3, 2, 3, 4, 4, 6, 2, 3, 6, 2, 2, 5, 4, 5, 4, 3, 4, 4, 2, 3, 6, 2, 3, 7, 5, 3, 3, 3, 7, 6, 3, 3, 6, 6, 3, 5, 3, 4, 4, 2, 5, 7, 2, 6, 6, 3, 4, 5, 7, 3, 5, 3, 5, 7, 4, 6, 10, 2, 3, 10, 5, 6, 5, 4, 5, 5, 4, 4, 11, 6, 2, 5, 4, 5, 3, 5, 6, 9, 6, 2, 9, 3
OFFSET
1,3
COMMENTS
The length of row n in A001269.
FORMULA
a(n) = A046051(2n) - A046051(n). - T. D. Noe, Jun 18 2003
a(n) = A001222(A000051(n)). - Amiram Eldar, Oct 04 2019
EXAMPLE
a(3) = 2 because 2^3 + 1 = 9 = 3*3.
MATHEMATICA
a[q_] := Module[{x, n}, x=FactorInteger[2^n+1]; n=Length[x]; Sum[Table[x[i]][2]], {i, n}][j]], {j, n}]]
A054992[n_Integer] := PrimeOmega[2^n + 1]; Table[A054992[n], {n, 200}] (* Vladimir Joseph Stephan Orlovsky, Jul 22 2011 *)
PROG
(PARI) a(n)=bigomega(2^n+1) \\ Charles R Greathouse IV, Apr 29 2015
CROSSREFS
bigomega(b^n+1): A057934 (b=10), A057935 (b=9), A057936 (b=8), A057937 (b=7), A057938 (b=6), A057939 (b=5), A057940 (b=4), A057941 (b=3), this sequence (b=2).
Cf. A046051 (number of prime factors of 2^n-1).
Cf. A086257 (number of primitive prime factors).
KEYWORD
nonn
AUTHOR
Arne Ring (arne.ring(AT)epost.de), May 30 2000
EXTENSIONS
Extended by Patrick De Geest, Oct 01 2000
Terms to a(500) in b-file from T. D. Noe, Nov 10 2007
Deleted duplicate (and broken) Wagstaff link. - N. J. A. Sloane, Jan 18 2019
a(500)-a(1062) in b-file from Amiram Eldar, Oct 04 2019
a(1063)-a(1122) in b-file from Max Alekseyev, Jul 15 2023
STATUS
approved
Smallest prime factor of 2^n + 1.
(Formerly M2385 N0947)
+10
15
3, 5, 3, 17, 3, 5, 3, 257, 3, 5, 3, 17, 3, 5, 3, 65537, 3, 5, 3, 17, 3, 5, 3, 97, 3, 5, 3, 17, 3, 5, 3, 641, 3, 5, 3, 17, 3, 5, 3, 257, 3, 5, 3, 17, 3, 5, 3, 193, 3, 5, 3, 17, 3, 5, 3, 257, 3, 5, 3, 17, 3, 5, 3, 274177, 3, 5, 3, 17, 3, 5, 3, 97, 3, 5, 3, 17, 3, 5, 3, 65537, 3, 5, 3, 17, 3, 5
OFFSET
1,1
COMMENTS
Conjecture: a(8+48*k) = 257 and a(40+48*k) = 257, where k is a nonnegative integer. - Thomas König, Feb 15 2017
Conjecture is true: 257 divides 2^(8+48*k)+1 and 2^(40+48*k)+1 but no prime < 257 ever does. Similarly, a(24+48*k) = 97. - Robert Israel, Feb 17 2017
From Robert Israel, Feb 17 2017: (Start)
If a(n) = p, there is some m such that a(n+m*j*n) = p for all j.
In particular, every member of the sequence occurs infinitely often.
a(k*n) <= a(n) for any odd k. (End)
REFERENCES
J. Brillhart et al., Factorizations of b^n +- 1. Contemporary Mathematics, Vol. 22, Amer. Math. Soc., Providence, RI, 2nd edition, 1985; and later supplements.
M. Kraitchik, Recherches sur la Théorie des Nombres, Gauthiers-Villars, Paris, Vol. 1, 1924, Vol. 2, 1929, see Vol. 2, p. 85.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
LINKS
Max Alekseyev, Table of n, a(n) for n = 1..3967 (first 500 terms from T. D. Noe)
J. Brillhart et al., Factorizations of b^n +- 1, Contemporary Mathematics, Vol. 22, Amer. Math. Soc., Providence, RI, 3rd edition, 2002.
E. Lucas, Théorie des Fonctions Numériques Simplement Périodiques, I", Amer. J. Math., 1 (1878), 184-240, 289-321. See pages 239 and 240.
Edouard Lucas, The Theory of Simply Periodic Numerical Functions, Fibonacci Association, 1969. English translation of article "Théorie des Fonctions Numériques Simplement Périodiques, I", Amer. J. Math., 1 (1878), 184-240.
S. S. Wagstaff, Jr., The Cunningham Project
FORMULA
a(n) = 3, 5, 3, 17, 3, 5, 3 for n == 1, 2, 3, 4, 5, 6, 7 (mod 8). (Proof. Let n = k*odd with k = 1, 2, or 4. As 2^k = 2, 4, 16 == -1 (mod 3, 5, 17), we get 2^n + 1 = 2^(k*odd) + 1 = (2^k)^odd + 1 == (-1)^odd + 1 == 0 (mod 3, 5, 17). Finally, 2^n + 1 !== 0 (mod p) for prime p < 3, 5, 17, respectively.) - Jonathan Sondow, Nov 28 2012
EXAMPLE
a(2^k) = 3, 5, 17, 257, 65537 is the k-th Fermat prime 2^(2^k) + 1 = A019434(k) for k = 0, 1, 2, 3, 4. - Jonathan Sondow, Nov 28 2012
MATHEMATICA
f[n_] := FactorInteger[2^n + 1][[1, 1]]; Array[f, 100] (* Robert G. Wilson v, Nov 28 2012 *)
FactorInteger[#][[1, 1]]&/@(2^Range[90]+1) (* Harvey P. Dale, Jul 25 2024 *)
PROG
(PARI) a(n) = my(m=n%8); if(m, [3, 5, 3, 17, 3, 5, 3][m], factor(2^n+1)[1, 1]); \\ Ruud H.G. van Tol, Feb 16 2024
KEYWORD
nonn,nice
EXTENSIONS
More terms from James A. Sellers, Jul 06 2000
Definition corrected by Jonathan Sondow, Nov 27 2012
STATUS
approved
Table T(n,k) in which n-th row lists prime factors of 2^n + 1 (n >= 0), without repetition.
+10
4
2, 3, 5, 3, 17, 3, 11, 5, 13, 3, 43, 257, 3, 19, 5, 41, 3, 683, 17, 241, 3, 2731, 5, 29, 113, 3, 11, 331, 65537, 3, 43691, 5, 13, 37, 109, 3, 174763, 17, 61681, 3, 43, 5419, 5, 397, 2113, 3, 2796203, 97, 257, 673, 3, 11, 251, 4051
OFFSET
0,1
COMMENTS
Rows have irregular lengths.
The length of row n is A046799(n).
REFERENCES
J. Brillhart et al., Factorizations of b^n +- 1. Contemporary Mathematics, Vol. 22, Amer. Math. Soc., Providence, RI, 2nd edition, 1985; and later supplements.
LINKS
T. D. Noe, Rows n = 0..500 of triangle, flattened (derived from Brillhart et al.)
J. Brillhart et al., Factorizations of b^n +- 1, Contemporary Mathematics, Vol. 22, Amer. Math. Soc., Providence, RI, 3rd edition, 2002.
S. S. Wagstaff, Jr., The Cunningham Project.
EXAMPLE
Triangle begins:
2;
3;
5;
3,17;
3,11;
5,13;
3,43;
257;
...
MATHEMATICA
Flatten[Table[Transpose[FactorInteger[2^n+1]][[1]], {n, 0, 25}]] (* Harvey P. Dale, Aug 10 2011 *)
PROG
(PARI) apply( A060444_row(n)=factor(2^n+1)[, 1]~, [0..10]) \\ M. F. Hasler, Nov 19 2018
CROSSREFS
Cf. A001269 (factors with repetition), A046799 (number of prime divisors).
KEYWORD
nonn,tabf
STATUS
approved
A list of primes written in order of their first appearance in a table of prime factorizations of 2^k+1, k=1,2,... .
+10
1
3, 5, 17, 11, 13, 43, 257, 19, 41, 683, 241, 2731, 29, 113, 331, 65537, 43691, 37, 109, 174763, 61681, 5419, 397, 2113, 2796203, 97, 673, 251, 4051, 53, 157, 1613, 87211, 15790321, 59, 3033169, 61, 1321, 715827883
OFFSET
1,1
COMMENTS
This sequence has the property that if a(n) appears first in the table as a prime factor of 2^m+1 for some m then a(n)=2*k*m+1 for some k.
When, for some m, 2^m+1 has more than one prime factor appearing in the table for the first time, we adopt the convention of entering them in ascending order. For example, the entries ..., 29, 113, ... both arise from 2^14+1.
LINKS
Harvey P. Dale and Charles R Greathouse IV, Table of n, a(n) for n = 1..4017 (first 650 terms from Dale)
EXAMPLE
2^1+1=3, 2^2+1=5, 2^3+1=3^2 and 2^4+1=17. Thus a(1)=3, a(2)=5 and a(3)=17, on noting that 2^3+1 contributes no new prime factors.
MATHEMATICA
DeleteDuplicates[Flatten[Table[Transpose[FactorInteger[2^k+1]][[1]], {k, 50}]]] (* Harvey P. Dale, Mar 30 2014 *)
PROG
(PARI) lista(n)=prs = Set(); for (k=1, n, f = factor(2^k+1); for (i=1, length(f~), onef = f[i, 1]; if (! setsearch(prs, onef), print1(onef, ", "); prs = setunion(prs, Set(onef)); ); ); ); \\ Michel Marcus, Apr 18 2013
(PARI) G=1; for(n=1, 500, g=gcd(f=2^n+1, G); while(g>1, g=gcd(g, f/=g)); f=factor(f)[, 1]; if(#f, for(i=1, #f, print1(f[i]", ")); G*=factorback(f))) \\ Charles R Greathouse IV, Jan 03 2018
CROSSREFS
Subsequence of A001269.
KEYWORD
nice,nonn
AUTHOR
Martin Griffiths, Mar 29 2009
STATUS
approved

Search completed in 0.033 seconds