OFFSET
1,1
LINKS
David A. Corneth, Table of n, a(n) for n = 1..10000 (first 2262 terms from M. F. Hasler)
FORMULA
EXAMPLE
a(10) = 26 as the 10th prime divides the 26th partition number, i.e., prime(10) = 29 | 2436 = A000041(26) and no smaller partition number is a multiple of 29. - David A. Corneth, Jan 05 2023
PROG
(PARI) forprime(p=2, 400, n=2:while(numbpart(n)%p, n=n+1):print1(n", "))
(PARI) first(n) = { my(res = vector(n, i, oo), todo = [1..n], pr = Set(primes(n)), prmap = Map()); for(i = 1, n, mapput(prmap, pr[i], todo[i]) ); for(i = 1, oo, c = numbpart(i); j = 1; while(j <= #pr && pr[j] <= c, if(c % pr[j] == 0, res[mapget(prmap, pr[j])] = i; pr = setminus(pr, Set(pr[j])); ); j++ ); if(#pr == 0, return(res) ) ); res } \\ David A. Corneth, Jan 05 2023
(sh) i=2; c=1; while [ $i -le 20000 ]; do if [ `factor $i | wc -w` = 2 ]; then awk '/^'$i' / {print "'$c' " $2; exit}' b046641.txt; let c+=1; if [ $i = 2 ]; then i=1; fi; fi; let i+=2; done # M. F. Hasler, Oct 20 2008
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Jan 29 2004
EXTENSIONS
Corrected by Ralf Stephan, Jul 27 2004
STATUS
approved