a = vector(10 000, n, n)

{
	tri = [0,0,0];

	print ("0 0");
	u = 1;
	for (xyz=1, oo,
		for (x=0, xyz,
			tri[1] = x;
			for (y=0, xyz-x,
				tri[2] = y;
				tri[3] = xyz-x-y;

				for (nb=4, oo,
					if (#tri<nb,
						tri = concat(tri, vector(#tri));
					);
					tri[nb] = tri[nb-3]+tri[nb-2]+tri[nb-1];

					if (tri[nb]>#a,
						for (k=1, nb-1,
							if (tri[k] && a[tri[k]] > xyz,
								a[tri[k]] = xyz;
							);
						);

						break;
					);
				);
			);
		);
		while (a[u]<=xyz,
			print (u " " a[u]);
			if (u++>#a,
				break (2);
			);
		);
	);
}

quit