（p+q+r）^nの係数をパスカルの三角形のように出力するやり方です。
１）フロベニウスの定理から３の分解ケースをソーティングして出す（逆順になっている）
l[n_]:=FrobeniusSolve[Table[1,{i,3}],n]
ｎは次数
２）各項の係数を算出
mul[n_]:=Table[Multinomial[l[n] k/. {x___}->x] ,{k,Length[l[n]],1,-1}]
３）剰余類による色分け。
これは係数に紐づいていればＯＫ
col1=Table[Which[ Mod[mulk,2]==0, Red,Mod[mulk,2]==1, Cyan],{k,1,Length[mul]}]
４）互いに接する円を色付けするコマンド
Graphics[Riffle[col1,Apply[Circle,Flatten[Table[{{i,-Sqrt[3]/2 j},1/2},{j,-4,4},{i,-j/2,j/2}],1],1]]]
