login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A197862
Prime divisor of n which appears the fewest times previously in the sequence, with ties to the larger prime.
1
2, 3, 2, 5, 3, 7, 2, 3, 5, 11, 3, 13, 7, 5, 2, 17, 3, 19, 5, 7, 11, 23, 2, 5, 13, 3, 7, 29, 5, 31, 2, 11, 17, 7, 3, 37, 19, 13, 5, 41, 7, 43, 11, 5, 23, 47, 2, 7, 2, 17, 13, 53, 3, 11, 7, 19, 29, 59, 5, 61, 31, 7, 2, 13, 11, 67, 17, 23, 7, 71, 3, 73, 37, 5, 19
OFFSET
2,1
COMMENTS
Up to n = 100, this differs from the greatest prime factor function A006530 only at n = 24, 48, 50, 80, and 98.
EXAMPLE
The only prime divisor of 4 is 2, so a(4) = 2.
The prime divisors of 6 are 2 and 3; in the sequence to that point (2,3,2,5), there are two 2's and 1 3, we take the less common one, so a(6) = 3.
The prime divisors of 12 are 2 and 3; these occur equally often in the sequence to that point, so we take the larger one; a(12)=3.
PROG
(PARI) al(n)={local(ns=vector(primepi(n)), r=vector(n-1), ps);
for(k=1, n-1,
ps=factor(k+1)[, 1]~;
r[k]=ps[1];
for(j=2, #ps, if(ns[primepi(ps[j])]<=ns[primepi(r[k])], r[k]=ps[j]));
ns[primepi(r[k])]++);
r}
CROSSREFS
Sequence in context: A276440 A162325 A324371 * A006530 A327398 A323616
KEYWORD
nonn
AUTHOR
STATUS
approved