f := proc(n,x,y) if (n < 1) then if (n mod 2 = 1) then return (f(n+1,x,y)+1)/f(n+2,x,y) else return (f(n+1,x,y)^4+1)/f(n+2,x,y); end if; else if (n = 1) then return x; else if (n = 2) then return y; else if (n mod 2 = 1) then return (f(n-1,x,y)+1)/f(n-2,x,y) else return (f(n-1,x,y)^4+1)/f(n-2,x,y); end if; end if; end if; end if; end proc; seq(f(i,1,1),i=-10..10); q := proc(n) local i,j,k,g; g := networks[graph]({},{}); for i from 0 to n-1 do for j from 0 to n-1 do for k from 0 to n-1 do if ((i^2+j^2+k^2) mod n = (3*i*j*k) mod n and not (i+j+k = 0)) then networks[addvertex](2^(i+1)*3^(j+1)*5^(k+1),g); end if; end do; end do; end do; for i from 0 to n-1 do for j from 0 to n-1 do for k from 0 to n-1 do if ((i^2+j^2+k^2) mod n = (3*i*j*k) mod n and not (i+j+k = 0)) then if (k > 0) then networks[addedge]({2^(i+1)*3^(j+1)*5^(k+1),2^(i+1)*3^(j+1)*5^((simplify((i^2+j^2)/k) mod n)+1)},g); end if; if (j > 0) then networks[addedge]({2^(i+1)*3^(j+1)*5^(k+1),2^(i+1)*3^((simplify((i^2+k^2)/j) mod n)+1)*5^(k+1)},g); end if; if (i > 0) then networks[addedge]({2^(i+1)*3^(j+1)*5^(k+1),2^((simplify((j^2+k^2)/i) mod n)+1)*3^(j+1)*5^(k+1)},g); end if; end if; end do; end do; end do; networks[connectivity](g); end proc; size := proc(n) if (n = {}) then return(0) else return (1+size(n minus {n[1]})); end if; end proc; size({1,2,3}); q(2);