กำหนดmโดยnแท่งช็อกโกแลตm,nบวกออกจำนวนวิธีการแบ่งแถบออกเป็นmn1 โดย 1 ชิ้นที่แต่ละตัวแบ่งเกิดขึ้นบนเส้นตาราง

คำสั่งซื้อเป็นสิ่งสำคัญ ชิ้นยังมีความแตกต่างดังนั้นทั้งสองชิ้นที่ปลายด้านหนึ่งของ 1 โดย 3 ช็อกโกแลตบาร์ไม่เท่ากัน

ตัวอย่างเช่นสำหรับบล็อก 2 คูณ 2 เรามี:

 _ _            _   _            _   _            _   _
|_‖_|    ->    |‗| |_|    ->    |_| |‗|    ->    |_| |_|
|_‖_|          |_| |_|           _  |_|           _   _
                                |_|              |_| |_|


 _ _            _   _            _   _            _   _
|_‖_|    ->    |_| |‗|    ->    |‗| |_|    ->    |_| |_|
|_‖_|          |_| |_|          |_|  _            _   _
                                    |_|          |_| |_|


 _ _            _ _              _   _            _   _
|‗|‗|    ->    |_‖_|      ->    |_| |_|    ->    |_| |_|
|_|_|           _ _               _ _             _   _
               |_|_|             |_‖_|           |_| |_|


 _ _            _ _               _ _             _   _
|‗|‗|    ->    |_|_|      ->     |_‖_|    ->     |_| |_|
|_|_|           _ _              _   _            _   _
               |_‖_|            |_| |_|          |_| |_|

ดังนั้นจึงมี 4 วิธีในการเลิกช็อคโกแลตแท่ง 2 ต่อ 2

กฎระเบียบ

• อินพุตจะเป็นจำนวนเต็มสองจำนวนผ่านฟังก์ชั่นอินพุต, STDIN, บรรทัดคำสั่งหรือคล้ายกัน เอาท์พุทจำนวนเดียวจำนวนวิธีที่จะสลายช็อคโกแลตบาร์

• เนื่องจากตัวเลขเพิ่มขึ้นอย่างรวดเร็วไม่ต้องกังวลหากผลลัพธ์เกินขีด จำกัด จำนวนเต็มของภาษาของคุณ - การส่งของคุณจะถูกต้องตราบใดที่อัลกอริทึมทำงานได้ตามทฤษฎีสำหรับอินพุตที่เป็นไปได้ทั้งหมด

กรณีทดสอบ

เอาท์พุทไม่ขึ้นอยู่กับคำสั่งของดังนั้นกรณีทดสอบมีการระบุไว้เช่นนั้นm,nm <= n

1 1 -> 1
1 2 -> 1
1 3 -> 2
1 4 -> 6
1 5 -> 24
1 10 -> 362880

2 2 -> 4
2 3 -> 56
2 4 -> 1712
2 5 -> 92800
2 10 -> 11106033743298560

3 3 -> 9408
3 4 -> 4948992
3 5 -> 6085088256
3 10 -> 76209753666310470268511846400

4 4 -> 63352393728

A261964m+nเป็นตัวเลขช็อคโกแลตได้จัดในรูปสามเหลี่ยมดังกล่าวที่แต่ละแถวสอดคล้องกับผลรวม

Mathematica, 85 ไบต์

f=If[##==1,1,Tr[Sum[Binomial[1##-2,i#-1]f[i,#]f[#2-i,#],{i,#2-1}]&@@@{{##},{#2,#}}]]&

กรณีทดสอบ

f[4,4]
(* 63352393728 *)

Python 3, 168 , 156 , 147 ไบต์

ทำการปรับปรุงด้วยการแชท

f=lambda n:n<1or n*f(n-1);a=lambda m,n,c=lambda m,n:sum(f(m*n-2)/f(i*n-1)/f((m-i)*n-1)*a(i,n)*a(m-i,n)for i in range(1,m)):+(m+n<4)or c(m,n)+c(n,m)

Ungolfed:

f=lambda n:n<1or n*f(n-1) # Factorial
def a(m, n):
    if m+n < 4:
        return 1
    first = 0
    for i in range(1,m):
        first += f(m*n-2) * 1/f(i*n-1) * 1/f((m-i)*n-1) * a(i,n) * a(m-i,n)
    second = 0
    for i in range(1,n):
        second += f(m*n-2) * 1/f(i*m-1) * 1/f((n-i)*m-1) * a(i,m) * a(n-i,m)
    return first + second

อัลกอริทึมนั้นมาจากบทความนี้

ฉันอาจจะตัดมันมากขึ้นฉันไม่แน่ใจว่าอยู่ตรงไหน

R, 208 198 ไบต์

f=function(m,n){g=function(i,j){a=0;if(j>1)for(x in 2:j-1)a=a+choose(j*i-2,x*i-1)*M[x,i]*M[j-x,i];a};s=max(m,n);M=matrix(1,s,s);for(i in 1:s)for(j in 1:s)if(i*j>2)M[i,j]=M[j,i]=g(i,j)+g(j,i);M[m,n]}

เยื้องกับบรรทัดใหม่:

f = function(m,n){
    g=function(i,j){
        a = 0
        if(j>1) for(x in 2:j-1) a = a + choose(j*i-2,x*i-1) * M[x,i] * M[j-x,i]
        a
    }
    s = max(m,n)
    M = matrix(1,s,s)
    for(i in 1:s) for(j in 1:s) if(i*j>2) M[i,j] = M[j,i] = g(i,j) + g(j,i)
    M[m,n]
}

การใช้งาน:

> f(3,1)
[1] 2
> f(3,2)
[1] 56
> f(3,3)
[1] 9408
> f(4,3)
[1] 4948992
> f(5,3)
[1] 6085088256

Python 2, 135 ไบต์

C=lambda A:sum(C(A[:i]+A[i+1:]+[(c,H),(W-c,H)])for i,Q in enumerate(A)for W,H in(Q,Q[::-1])for c in range(1,W))or 1
print C([input()])

นี่คือสิ่งที่ฉันมาด้วย มันช้ามาก แต่นี่เป็นรุ่นที่เร็วกว่า (ต้องการrepoze.lru ):

from repoze.lru import lru_cache
C=lru_cache(maxsize=9999)(lambda A:sum(C(tuple(sorted(A[:i]+A[i+1:]+((c,H),(W-c,H)))))for i,Q in enumerate(A)for W,H in(Q,Q[::-1])for c in range(1,W))or 1)
print C((input(),))

ตัวอย่าง

$ time python2 chocolate.py <<< 2,5
92800

real    0m2.954s
user    0m0.000s
sys     0m0.015s

$ time python2 chocolate-fast.py <<< 3,5
6085088256

real    0m0.106s
user    0m0.000s
sys     0m0.015s

คำอธิบาย

รหัสกำหนดฟังก์ชั่นCที่ใช้อาร์เรย์ของชิ้นส่วน อัลกอริทึมเป็นเช่นนี้:

1. for i,Q in enumerate(A): วนรอบอาร์เรย์ของชิ้นส่วนต่างๆ
2. for W,H in(Q,Q[::-1]): คำนวณวิธีสองครั้งหมุน 90 องศา
3. for c in range(1,W): วนผ่านตำแหน่งที่เป็นไปได้ที่จะแยกที่
4. A[:i]+A[i+1:]+[(c,H),(W-c,H)]: รับรายการโดยไม่มีชิ้นส่วนแยกและชิ้นใหม่สองชิ้น
5. C(…): เรียกใช้ฟังก์ชันอีกครั้งในรายการนั้น
6. sum(…): รวมผลลัพธ์สำหรับแต่ละการแบ่งที่เป็นไปได้
7. or 1: หากไม่สามารถแยกได้จะมีวิธีหนึ่งในการแยกช็อคโกแลต

ในที่สุดรหัสจะถูกเรียกพร้อมกับอาร์เรย์ที่มีอินพุต

ES6, 141 ไบต์

c=(m,n)=>(f=n=>n<2||n*f(n-1),h=(m,n)=>[...Array(m-1)].reduce((t,_,i)=>t+f(m*n-2)/f(++i*n-1)/f((m-i)*n-1)*c(i,n)*c(m-i,n),0),h(m,n)+h(n,m))||1

ขึ้นอยู่กับสูตรที่พบโดย @CameronAavik Ungolfed:

function fact(n) {
    return n < 2 ? 1 : n * f(n - 1);
}
function half(m, n) {
    total = 0;
    for (i = 1; i < m; i++)
        total += fact(m * n - 2) / fact(i * n - 1) / fact((m - i) * n - 1) * choc(i, n) * choc(m - i, n)
    return total;
}
function choc(m, n) {
    total = half(m, n) + half(n, m);
    if (!total) total = 1;
    return total;
}