Skip to content

Instantly share code, notes, and snippets.

@igorvanloo
Last active March 24, 2022 16:53
Show Gist options
  • Save igorvanloo/7b415fa80e316a2faf4b60d4127c7810 to your computer and use it in GitHub Desktop.
Save igorvanloo/7b415fa80e316a2faf4b60d4127c7810 to your computer and use it in GitHub Desktop.
admissible numbers
def admissible_numbers(p, n, primes, limit):
a_n = [n] #admissible numbers
i = primes.index(p) #index
if n > limit:
return []
a_n += admissible_numbers(p, n*p, primes, limit)
if i + 1 < len(primes):
n_p = primes[i + 1] #next prime
a_n += admissible_numbers(n_p, n*n_p, primes, limit)
return a_n
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment