User:Roman Czyborra

Lerch Transcendent
$\sum_{n=N}^\infty R^n =\frac{1-R}{1-R}\sum_{n=N}^\infty R^n =\frac{\sum_{n=N}^\infty R^n-R\sum_{n=N}^\infty R^n}{1-R} =\frac{\sum_{n=N}^\infty R^n-\sum_{n=N}^\infty R^{n+1}}{1-R} =\frac{\sum_{n=N}^\infty R^n-\sum_{n=N+1}^\infty R^n}{1-R} =\frac{R^N}{1-R}$

$\sum_{n=N}^\infty nR^n =\frac{\sum_{n=N}^\infty nR^n-\sum_{n=N+1}^\infty[n-1]R^n}{1-R} =\frac{NR^N+\sum_{n=N+1}^\infty 1R^n}{1-R} =\frac{NR^N+\frac{R^{N+1}}{1-R}}{1-R} =\frac{R^N}{1-R}\left[N+\frac{R}{1-R}\right]$

$\sum_{n=N}^\infty n^2R^n =\frac{\sum_{n=N}^\infty n^2R^n-\sum_{n=N+1}^\infty[n-1]^2R^n}{1-R} =\frac{\sum_{n=N}^\infty n^2R^n-\sum_{n=N+1}^\infty n^2R^n+2\sum_{n=N+1}^\infty nR^n-\sum_{n=N+1}^\infty R^n}{1-R} =\frac{N^2R^N+2\frac{R^{N+1}}{1-R}\left[N+1+\frac{R}{1-R}\right]-\frac{R^{N+1}}{1-R}}{1-R} =\frac{R^N}{1-R}\left[N^2+\frac{R}{1-R}\left[2N+1+\frac{2R}{1-R}\right]\right]$

$\sum_{n=N}^\infty n^3R^n =\frac{N^3R^N+3\frac{R^{N+1}}{1-R}\left[N^2+2N+1+\frac{R}{1-R}\left[2N+3+\frac{2R}{1-R}\right]\right]-3\frac{R^{N+1}}{1-R}\left[N+1+\frac{R}{1-R}\right]+\frac{R^{N+1}}{1-R}}{1-R} =\frac{R^N}{1-R}\left[N^3+\frac{R}{1-R}\left[3N^2+3N+1+\frac{R}{1-R}\left[6N+6+\frac{6R}{1-R}\right]\right]\right]$

$\sum_{n=N}^\infty n^4R^n =\frac{R^N}{1-R}\left[N^4+\frac{R}{1-R}\left[4N^3+6N^2+4N+1+\frac{R}{1-R}\left[12N^2+24N+14+\frac{R}{1-R}\left[24N+36+\frac{24R}{1-R}\right]\right]\right]\right]$

apparently $\sum_{n=N}^\infty n^pR^n=\frac{R^N}{1-R}\sum_{r=0}^p(\frac{R}{1-R})^r\sum_{n=0}^rN^nϤ(p,r,n)$ where $Ϥ(p,r,n)=[r=0][p=n]-[r\ne0]\sum_{s=1}^p(-1)^s{p\choose s}\sum_{m=n}^r{m\choose n}Ϥ(p-s,r-1,m)$

[//en.wikipedia.org/wiki/Exact_cover#Exact_hitting_set Exact Covers and Hittings]
$c(s) = \{ r\subseteq s : \forall x \in \cup s : \exists! y \in r: y \ni x \}$

$h(s) = \{ r\subseteq \cup s: \forall x \in s : \exists! y \in r: y \in x \}$

exacthittings sets = if null sets then else [ pivot : extension | pivot <- minimumBy (comparing length) sets, let (deletables, undeletables) = partition (elem pivot) sets, extension <- exacthittings (map (\\ unions deletables) undeletables) ] test0 = exacthittings (words "ab ef de abc cd de acef") == words "bdf" sudokus = exacthittings (concat [ [       row x >< [y],       col x >< [y],       [x] >< [y] >< [0..8],       box (divMod x 3) >< [y] ] | (x,y) <- [0..8] >< [0..8] ]) where (><) = liftM2 row i = [i] >< [0..8] col i = [0..8] >< [i] box (i,j) = map (3*i+) [0..2] >< map (3*j+) [0..2] test1 = sudokus

[//en.wikipedia.org/wiki/Sedenion Sedenion Multiplication]
tee --append PowerNum.lhs << 0 - partner (umbox - row) col where umbox = until (>row) (2*) 1 LT -> 0 - flip partner row col EQ -> row EOF

$ \begin{bmatrix} +P0&-P1&-P2&-P3&-P4&-P5&-P6&-P7&-P8&-P9&-PA&-PB&-PC&-PD&-PE&-PF \\ +P1&+P0&-P3&+P2&-P5&+P4&+P7&-P6&-P9&+P8&+PB&-PA&+PD&-PC&-PF&+PE \\ +P2&+P3&+P0&-P1&-P6&-P7&+P4&+P5&-PA&-PB&+P8&+P9&+PE&+PF&-PC&-PD \\ +P3&-P2&+P1&+P0&-P7&+P6&-P5&+P4&-PB&+PA&-P9&+P8&+PF&-PE&+PD&-PC \\ +P4&+P5&+P6&+P7&+P0&-P1&-P2&-P3&-PC&-PD&-PE&-PF&+P8&+P9&+PA&+PB \\ +P5&-P4&+P7&-P6&+P1&+P0&+P3&-P2&-PD&+PC&-PF&+PE&-P9&+P8&-PB&+PA \\ +P6&-P7&-P4&+P5&+P2&-P3&+P0&+P1&-PE&+PF&+PC&-PD&-PA&+PB&+P8&-P9 \\ +P7&+P6&-P5&-P4&+P3&+P2&-P1&+P0&-PF&-PE&+PD&+PC&-PB&-PA&+P9&+P8 \\ +P8&+P9&+PA&+PB&+PC&+PD&+PE&+PF&+P0&-P1&-P2&-P3&-P4&-P5&-P6&-P7 \\ +P9&-P8&+PB&-PA&+PD&-PC&-PF&+PE&+P1&+P0&+P3&-P2&+P5&-P4&-P7&+P6 \\ +PA&-PB&-P8&+P9&+PE&+PF&-PC&-PD&+P2&-P3&+P0&+P1&+P6&+P7&-P4&-P5 \\ +PB&+PA&-P9&-P8&+PF&-PE&+PD&-PC&+P3&+P2&-P1&+P0&+P7&-P6&+P5&-P4 \\ +PC&-PD&-PE&-PF&-P8&+P9&+PA&+PB&+P4&-P5&-P6&-P7&+P0&+P1&+P2&+P3 \\ +PD&+PC&-PF&+PE&-P9&-P8&-PB&+PA&+P5&+P4&-P7&+P6&-P1&+P0&-P3&+P2 \\ +PE&+PF&+PC&-PD&-PA&+PB&-P8&-P9&+P6&+P7&+P4&-P5&-P2&+P3&+P0&-P1 \\ +PF&-PE&+PD&+PC&-PB&-PA&+P9&-P8&+P7&-P6&+P5&+P4&-P3&-P2&+P1&+P0 \\ \end{bmatrix} \begin{bmatrix} Q0\\Q1\\Q2\\Q3\\Q4\\Q5\\Q6\\Q7\\Q8\\Q9\\QA\\QB\\QC\\QD\\QE\\QF \end{bmatrix} =: \begin{bmatrix} P0\\P1\\P2\\P3\\P4\\P5\\P6\\P7\\P8\\P9\\PA\\PB\\PC\\PD\\PE\\PF \end{bmatrix} \begin{bmatrix} Q0\\Q1\\Q2\\Q3\\Q4\\Q5\\Q6\\Q7\\Q8\\Q9\\QA\\QB\\QC\\QD\\QE\\QF \end{bmatrix}$

sudo apt-get install open-axiom ; open-axiom octon(p0,p1,p2,p3,p4,p5,p6,p7)*octon(q0,q1,q2,q3,q4,q5,q6,q7)-conjugate(octon(q8,q9,qa,qb,qc,qd,qe,qf))*octon(p8,p9,pa,pb,pc,pd,pe,pf) octon(q8,q9,qa,qb,qc,qd,qe,qf)*octon(p0,p1,p2,p3,p4,p5,p6,p7)+octon(p8,p9,pa,pb,pc,pd,pe,pf)*conjugate(octon(q0,q1,q2,q3,q4,q5,q6,q7)) )quit

Graphemica
(🏭)=map; 囪=filter; 㪭=fold