The number of non-negative integral solutions of the given equation is equal to the coefficient of x20 in (1−x)−1(1−x)−1(1−x)−1×(1−x4)−1, which is given by
Coefficient of x20 in (1−x)−3(1−x4)−1
= Coefficient of x20 in (1+3C1x+4C2x2+5C3x3+6C4x4+...)(1+x4+x8+...)
=1+6C4+10C8+14C12+18C16+22C20=536