|
|
A060954
|
|
Largest prime factor of 10*n + 1.
|
|
1
|
|
|
11, 7, 31, 41, 17, 61, 71, 3, 13, 101, 37, 11, 131, 47, 151, 23, 19, 181, 191, 67, 211, 17, 11, 241, 251, 29, 271, 281, 97, 43, 311, 107, 331, 31, 13, 19, 53, 127, 23, 401, 137, 421, 431, 7, 41, 461, 157, 37, 491, 167, 73, 521, 59, 541, 29, 17, 571, 83, 197, 601, 47
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
LINKS
|
|
|
EXAMPLE
|
21 = 3*7, so a(2) = 7.
|
|
MATHEMATICA
|
Table[First[Last[FactorInteger[10n + 1]]], {n, 1, 80}]
|
|
PROG
|
(PARI) { for (n=1, 1000, f=factor(10*n + 1)~; write("b060954.txt", n, " ", f[1, length(f)]); ) } \\ Harry J. Smith, Jul 19 2009
(Python)
from sympy import factorint
def a(n): return max(factorint(10*n+1))
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|