z(a,b):=a+b*%i;

f(a,b):=z(a,b)^60;

f(x,y);

factor(realpart(f(x,y)));

factor(imagpart(f(x,y)));

realpart(f(-1,sqrt(3)));

imagpart(f(-1,sqrt(3)));

2^60;