var('x,y,z')
r=x+y*(-1+sqrt(-3))/2+z*(-1-sqrt(-3))/2
s=x+y*(-1-sqrt(-3))/2+z*(-1+sqrt(-3))/2
g(x,y,z)=r^3+s^3+r*s
q(a,b,c)=0
for n in range(8):
f(x,y)=g(x,y,0)
p(a,b)=0
for m in range(11):
p(a,b)=p(a,b)+f(a,0)*b^m
f(x,y)=expand((f(x,y)-f(x+y,0))/x/y)
q(a,b,c)=q(a,b,c)+p(a,b)*c^n
g(x,y,z)=expand((g(x,y,z)-p(x+y+z,x*y+x*z+y*z))/x/y/z)
print (q(a,b,c))