Ακολουθεί ο κώδικας εισόδου (input) που παραθέσαμε στο παραπάνω έγγραφο. Όσοι έχουν το λογισμικό Mathematica μπορούν να τον αντιγράψουν και να τον επικολλήσουν. Πιθανώς να χρειαστεί να μετατραπεί το format από text σε input (αφού έχει μαρκαριστεί πρώτα το αντίστοιχο «κελί»). Στη συνέχεια κάνουμε Shift +Enter, όπου συναντάμε κενή γραμμή κατά την αντιγραφή.
ath[m_Integer,n_Integer]:=Module[{a,i,d,diair},
diair=Divisors[Factorial[n]];d=Length[diair];
If[m<=Factorial[n],Print[diair];
a=Table[diair.IntegerDigits[k,2,d],{k,1,2^d}];
For[i=1,i<=2^d,
If[Part[a,i]==m&&Apply[Plus,IntegerDigits[i,2]]<=n,
Print[IntegerDigits[i,2,d]]];i++]]]
ath [11,4]
<<DiscreteMath`Combinatorica`
ath1[m_,n_]:=Module[{a,b,diair,d,e,s,i,j},
diair=Divisors[Factorial[n]];d=Length[diair];e={};
For[s=1,s<=d,If[Part[diair,s]<=m,
AppendTo[e,Part[diair,s]]];s++];
a={};For[i=1,i<=n,a=Join[a,KSubsets[e,i]];i++];
b={};For[j=1,j<=Length[a],
If[Apply[Plus,Part[a,j]]==m,
AppendTo[b,Part[a,j]]];j++];b]
ath1[13,5]
KS[a_List,n_Integer]:=Module[{b,c,i,j},b={a};
While[Length[Part[b,1]]>n,c={};
For[i=1,i<=Length[b],
For[j=1,j<=Length[Part[b,1]],
AppendTo[c,Delete[Part[b,i],j]];j++];i++];b=c];
