import primes

def find(f):
	factors = []
	t = 1
	s = primes.sieve()
	for p in s.sieve():
		if f % p == 0:
			if not p in factors:
				factors.append(p)
				print "Factor found: %d" % (p)
				t *= p
				if t == f:
					return factors

print "Factors: %s" % (find(600851475143))