f:=proc(k, n) if k+1<=n then RETURN(0) else RETURN(k!/(k - n)!) fi: end; a:=n->(1/6!)*(f(64, n) - 30*f(48, n) + 120*f(40, n) + 60*f(36, n) + 60*f(34, n)- 12*f(33, n) - 345*f(32, n) - 720*f(30, n) + 810*f(28, n) + 120*f(27, n) + 480*f(26, n) + 360*f(25, n) - 480*f(24, n) - 720*f(23, n) - 240*f(22, n) - 540*f(21, n) + 1380*f(20, n) + 750*f(19, n) + 60*f(18, n) - 210*f(17, n) - 1535*f(16, n) - 1820*f(15, n) + 2250*f(14, n) + 1800*f(13, n) - 2820*f(12, n) + 300*f(11, n) + 2040*f(10, n) + 340*f(9, n) - 1815*f(8, n) + 510*f(7, n) - 1350*f(6, n) + 1350*f(5, n) + 274*f(4, n) - 548*f(3, n) + 120*f(2, n)); seq(a(n), n=0..20); # Pab Ter (pabrlos2(AT)yahoo.com), Nov 06 2005