OFFSET
1,1
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
PROG
(PARI) valp(n, p)=my(s); while(n\=p, s+=n); s
is(n)=my(t=(12*n-1)\13*13+13, s=valp(t, 13)-n); while(s<0, s+=valuation(t+=13, 13)); s>0 \\ Charles R Greathouse IV, Oct 19 2016
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Carl R. White, Jan 21 2008
STATUS
approved