первый вопрос здесь. Я пытаюсь изучить Python, пройдя через проект euler, и наткнулся на препятствие. Следующий метод (возвращает список простых множителей) отлично работает для одного вызова:
def findPrimeFactors(num, primeFactors = []):
'''Find the prime factors of an arbitrary positive integer
input: num to factorize
returns: a list containing the prime factors of the number
'''
pIndex = 2
while (num >= pIndex):
if num % pIndex == 0:
num /= pIndex
primeFactors.append(pIndex)
return FindPrimes.findPrimeFactors(num, primeFactors)
else:
pIndex += 1
return primeFactors
Однако, когда я использую его в цикле, например (этот метод еще не завершен, в настоящее время он приводит к бесконечному циклу, так как больше простых чисел не может быть найдено):
def countPrimes(n = 1001):
'''find n amount of unique primes ascending
input: number of primes to find
returns: list of n primes starting from 2 '''
primes = []
i = 2
while len(primes) < n:
primeFactors = FindPrimes.findPrimeFactors(i)
print(primeFactors) #verify method behavior
if len(primeFactors) is 1:
primes.append(primeFactors[0])
i += 1
return primes
В результате первый цикл возвращает [2], следующий возвращает [2, 3] и так далее, добавляя новые результаты к списку, который я хотел, чтобы он был пустым при первом рекурсивном вызове. Кажется, что мой список сохраняется, но я не совсем понимаю, почему? Я также прочитал Объем и списки классов Python, который дает мне некоторые подсказки, но рекурсия еще больше усложняет .
Рекурсивный также означает, что я не могу просто назначить ему пустой набор. Исходя из опыта работы с C ++, я ожидал, что переменная primeFactors должна повторно инициализироваться каждый раз, когда функция вызывается из моей программы. Все еще здесь змея.
РЕДАКТИРОВАТЬ: это итеративная версия findPrimeFactors, которую я написал. Я знаю, что это не оптимально, но я хотел бы, по крайней мере, сделать его достаточно эффективным, чтобы соответствовать правилу 1 минуты проекта Эйлера. Любые предложения по улучшению или ясности приветствуются.
PRIMES = [2,3,5,7,11,13,17,19]
import math
class FindPrimes():
'''V2 iterative'''
def findPrimeFactors(n, primeFactors = None):
'''Find the prime factors of an arbitrary positive integer
input: num to factorize
returns: a list containing the prime factors of the number
'''
if primeFactors is None:
primeFactors = []
num = n
ceil = math.sqrt(n) #currently unused
global PRIMES
knownPrimes = PRIMES
#check known primes for divisors first, then continue searching for primes by brute force
while True:
factorFound = False
for prime in knownPrimes:
if num % prime == 0:
primeFactors.append(prime)
num /= prime
factorFound = True
break #ensure that the list returned has ascending primes
if not factorFound:
break
#once attempts have been made to reduce using known primes
#search for new primes if the number is not fully reduced
i = knownPrimes[-1] + 2
while num != 1:
if num % i == 0:
knownPrimes.append(i)
primeFactors.append(i)
num /= i
i += 2
return primeFactors
def countPrimes(n = 10001):
'''find n amount of unique primes ascending
input: number of primes to find
returns: list of n primes starting from 2 '''
primes = []
i = 2
while len(primes) < n:
primeFactors = FindPrimes.findPrimeFactors(i)
if len(primeFactors) == 1:
primes.append(primeFactors[0])
#print(primeFactors[-1])
i += 1
print(len(primes))
return primes
nth = 10001
print(FindPrimes.countPrimes(nth)[nth-1]) #print the largest prime found