This site is supported by donations to The OEIS Foundation.

User:Alexander Staunton

From OeisWiki
Jump to: navigation, search

I don't really have any credentials, nor do I see myself contributing much. I just didn't know how else to submit the following (python) code for sequences A091588 and A091839:

def valuation(k,n):

   i = 0
   while (n % (k**(i+1)) == 0):i += 1
   return i

def series(k,n): return 0 if (n == 0) else (k * records(k,n - 1) + (k + 1) * series(k,n-1)) def records(k,n): return 1 if (n == 0) else (records(k + 1,valuation(k+2,n)) + series(k,n - 1) + 2 * records(k,n - 1)) def A091588_list(n): return [records(1,i) for i in range(n)] def A091839_list(n): return [records(1,valuation(2,i)) for i in range(1,n+1)]