Q72上

import tools
import math
c = tools.clock()
c.start()

max =1000000
list_primes = tools.gen_list_primes(max)

print len(list_primes)
dict_phi = {}

for prime in list_primes:
    dict_phi[prime] = prime - 1

 

list_answer = []
for i in range(max+1):
    list_answer.append(0)

for prime1 in list_primes:
    n1 = 1
    n1 *= prime1
    if n1 <= max:
        a1 = n1*(prime1-1)/prime1
        list_answer[n1] = a1
        for prime2 in list_primes:
            n2 = n1*prime2
            if n2 <= max:
                if prime2!=prime1:
                    a2 = a1*(prime2-1)
                    #/prime2*prime2
                else:
                    a2 = a1*prime2
                list_answer[n2] = a2
            else:
                break