read MultiAlmkvistZeilberger:
print(MAZ(1,1/x/y/z/(1-x-y-z+4*x*y*z)^a,1/x/y/z,[x,y,z],n,N,{a}));

time();

quit: