1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:
19:

def calcNextProbVector(prev):
    '''berechnet 6^(k+1)*Wkt-Vektor [k+1,k+2,...,6k+6] für die Augenzahlsumme von (k+1) Würfen'''
    '''auf der Grundlage des 6^k*Vektors [k,k+1,...,6k] für die Augenzahlsumme von k Würfen'''
    l = len(prev)   # Es ist l = 5k+1
    return [sum(prev[j] for j in range(max(0,i-5),min(l,i+1))) for i in range(l+5)]

def calcProb(n,a,b):
    '''berechnet die 6^n*Wkt, dass die Augenzahlsumme von n Würfen im Intervall [a,b] liegt'''
    at = max(a-n,0)
    bt = min(b-n,5*n)
    if at > bt:
        return 0
    pVector = [1]
    for k in range(n):
        pVector = calcNextProbVector(pVector)  # man hangelt sich hoch von 0 bis n
    return sum(pVector[k] for k in range(at,bt+1))

if __name__ == '__main__':
    print(f"{calcProb(60,241,360)}/{6**60}")

1:
2:

p:=(x+x^2+x^3+x^4+x^5+x^6)^60:
_plus(coeff(p,x,k) $ k=241..360)