Я правильно понял что вы пытаетесь сделать список с перебором всех подходящих бросков кубика? А потом для подсчета вероятности поделить длину списка на 20^3?

А если так? (сильно не бить, haskell давно не трогал :))
diff x a | x - а > 0 = х - а
| х - а <= 0 = 0
erfolgG (a1,a2,a3) n = [(x,y,z) | x<-[1..20], y<-[1..20], z<-[1..20], x+y<40, x+z<40, y+z<40,
(diff x a1)+(diff y a2)+(diff z a3) <= n]
Но. Про варианты с двумя выброшенными единицами вы забыли
Вам надо бы что-то вроде (с case-ом вроде понятнее)
check x y z a1 a2 a3 n = case (x,y,z) of
(1,1,_) -> True
(1,_,1) -> True
(_,1,1) -> True
(20,20,_) -> False
(20,_,20) -> False
(_,20,20) -> False
(_,_,_) -> (diff x a1)+(diff y a2)+(diff z a3) <= n
и
erfolgG (a1,a2,a3) n = [(x,y,z) | x<-[1..20], y<-[1..20], z<-[1..20], check x y z a1 a2 a3 n]