polynomial

edit

.

p=5
f1=(p^2-1)
f2=divisors(f1)
f2
for  i in f2 :
    f=cyclotomic_polynomial(i,'z')
    print (i,f)

Result:

[1, 2, 3, 4, 6, 8, 12, 24]
(1, z - 1)
(2, z + 1)
(3, z^2 + z + 1)
(4, z^2 + 1)
(6, z^2 - z + 1)
(8, z^4 + 1)
(12, z^4 - z^2 + 1)
(24, z^8 - z^4 + 1)

i want the above result in the following form

phi1= z - 1
phi2= z + 1
phi3= z^2 + z + 1
phi4= z^2 + 1
phi6= z^2 - z + 1
phi8= z^4 + 1
phi12= z^4 - z^2 + 1
phi24= z^8 - z^4 + 1

how to obtained.