print(`The list below are the first 12 terms in the sequence`):
print(`number of 3-deck sets of cards with n cards each that are "sucker's bets" i.e.`):
print(`deck A beats deck B, deck B beats deck C and deck C beats deck A. See the`):
print(`article "a Treatise on Sucker's bets"`):
print(`http://sites.math.rutgers.edu/~zeilberg/mamarim/mamarimhtml/suckerbets.html`):

gu:=[0, 0, 15, 39, 5196, 32115, 2093199, 19618353, 960165789, 11272949151, 479538890271, 6504453085104]:

print(gu):

print(`Up to cyclic synmmetry, (i.e. dividing the above by 3):`):

print(gu/3):

mu:=[seq(gu[i]/(3*i)!*i!^3,i=1..nops(gu))]:

print(`The sequence of probabilities is`):

print(``):
print(mu):
print(``):

print(`and in floating point`):

print(``):
print(evalf(mu)):
print(``):

quit: