#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
scalar_t__ EOF;
int /*<<< orphan*/ printf(char *, long long);
scalar_t__ scanf(char *, long long *, long long *);
int main() {
  long long int i, j, k, u, n, m, a, b, ar[110][110];
  while (scanf("%lld %lld", &a, &b) != EOF) {
    for (i = 1; i <= 101; i++) {
      ar[i][0] = 1;
      ar[0][i] = 0;
    }
    ar[0][0] = 0;
    for (i = 1; i <= a; i++) {
      for (j = 1; j <= b; j++) {
        ar[i][j] = ar[i - 1][j] * (j + 1) + ar[i - 1][j - 1] * (i - j);
        ar[i][j] %= 2007;
      }
    }
    printf("%lld\n", ar[a][b] % 2007);
  }
  return 0;
}


x86 assembly (clang, -O0):
<main>:
push   %rbp
mov    %rsp,%rbp
sub    $0x17a80,%rsp
movl   $0x0,-0x4(%rbp)
lea    0xec7(%rip),%rdi
lea    -0x40(%rbp),%rsi
lea    -0x48(%rbp),%rdx
call   1040 <scanf@plt>
mov    0x2e9b(%rip),%rcx
cmp    (%rcx),%rax
je     12ea <main+0x1ca>
movq   $0x1,-0x10(%rbp)
cmpq   $0x65,-0x10(%rbp)
jg     11a3 <main+0x83>
lea    -0x17a70(%rbp),%rax
imul   $0x370,-0x10(%rbp),%rcx
add    %rcx,%rax
movq   $0x1,(%rax)
mov    -0x10(%rbp),%rax
movq   $0x0,-0x17a70(%rbp,%rax,8)
mov    -0x10(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x10(%rbp)
jmp    115e <main+0x3e>
movq   $0x0,-0x17a70(%rbp)
movq   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%rax
cmp    -0x40(%rbp),%rax
jg     12b2 <main+0x192>
movq   $0x1,-0x18(%rbp)
mov    -0x18(%rbp),%rax
cmp    -0x48(%rbp),%rax
jg     129c <main+0x17c>
mov    -0x10(%rbp),%rcx
sub    $0x1,%rcx
lea    -0x17a70(%rbp),%rax
imul   $0x370,%rcx,%rcx
add    %rcx,%rax
mov    -0x18(%rbp),%rcx
mov    (%rax,%rcx,8),%rdx
mov    -0x18(%rbp),%rax
add    $0x1,%rax
imul   %rax,%rdx
mov    -0x10(%rbp),%rcx
sub    $0x1,%rcx
lea    -0x17a70(%rbp),%rax
imul   $0x370,%rcx,%rcx
add    %rcx,%rax
mov    -0x18(%rbp),%rcx
sub    $0x1,%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x10(%rbp),%rcx
sub    -0x18(%rbp),%rcx
imul   %rcx,%rax
add    %rax,%rdx
lea    -0x17a70(%rbp),%rax
imul   $0x370,-0x10(%rbp),%rcx
add    %rcx,%rax
mov    -0x18(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
lea    -0x17a70(%rbp),%rax
imul   $0x370,-0x10(%rbp),%rcx
add    %rcx,%rax
mov    %rax,-0x17a78(%rbp)
mov    -0x18(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    $0x7d7,%esi
cqto
idiv   %rsi
mov    -0x17a78(%rbp),%rax
mov    %rdx,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x18(%rbp)
jmp    11cc <main+0xac>
jmp    12a1 <main+0x181>
mov    -0x10(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x10(%rbp)
jmp    11b6 <main+0x96>
lea    -0x17a70(%rbp),%rax
imul   $0x370,-0x40(%rbp),%rcx
add    %rcx,%rax
mov    -0x48(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    $0x7d7,%ecx
cqto
idiv   %rcx
mov    %rdx,%rsi
lea    0xd2a(%rip),%rdi
call   1030 <printf@plt>
jmp    1132 <main+0x12>
xor    %eax,%eax
add    $0x17a80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<main>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x17a38,%rsp
lea    0xeb8(%rip),%rdi
lea    0x8(%rsp),%rsi
mov    %rsp,%rdx
call   1050 <scanf@plt>
mov    0x2e8c(%rip),%r12
cmp    (%r12),%rax
jne    1176 <main+0x46>
xor    %eax,%eax
add    $0x17a38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
movabs $0x829d66ae3cf8796f,%r14
66 9d 82
mov    %rsp,%r13
lea    0x10(%rsp),%r15
jmp    11ed <main+0xbd>
nopw   0x0(%rax,%rax,1)
imul   $0x370,%r8,%rax
add    %rsp,%rax
add    $0x10,%rax
mov    (%rax,%r11,8),%rsi
mov    %rsi,%rax
imul   %r14
add    %rsi,%rdx
mov    %rdx,%rax
shr    $0x3f,%rax
sar    $0xa,%rdx
add    %rax,%rdx
imul   $0x7d7,%rdx,%rax
sub    %rax,%rsi
lea    0xe40(%rip),%rdi
call   1030 <printf@plt>
lea    0xe2a(%rip),%rdi
lea    0x8(%rsp),%rsi
mov    %r13,%rdx
call   1050 <scanf@plt>
cmp    (%r12),%rax
je     1162 <main+0x32>
mov    $0x328,%edx
lea    0x18(%rsp),%rdi
xor    %esi,%esi
call   1040 <memset@plt>
mov    $0xdc0,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movq   $0x1,-0xa40(%rsp,%rax,1)
cmp    $0x16580,%rax
je     1250 <main+0x120>
movq   $0x1,-0x6d0(%rsp,%rax,1)
movq   $0x1,-0x360(%rsp,%rax,1)
movq   $0x1,0x10(%rsp,%rax,1)
add    $0xdc0,%rax
jmp    1210 <main+0xe0>
nopl   (%rax)
movq   $0x0,0x10(%rsp)
mov    0x8(%rsp),%r8
mov    (%rsp),%r11
test   %r8,%r8
jle    1190 <main+0x60>
mov    $0x1,%r9d
lea    0x388(%rsp),%rbp
xor    %r10d,%r10d
jmp    129b <main+0x16b>
xchg   %ax,%ax
lea    0x1(%r9),%rax
add    $0x1,%r10
add    $0x370,%rbp
cmp    %r8,%r9
mov    %rax,%r9
je     1190 <main+0x60>
test   %r11,%r11
jle    1280 <main+0x150>
imul   $0x370,%r9,%rax
mov    -0x370(%rax,%r15,1),%rcx
ff
mov    %r10,%rsi
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    -0x370(%rbp,%rbx,8),%rdi
ff
lea    0x2(%rbx),%rax
imul   %rdi,%rax
imul   %rsi,%rcx
add    %rax,%rcx
mov    %rcx,%rax
imul   %r14
add    %rcx,%rdx
mov    %rdx,%rax
shr    $0x3f,%rax
sar    $0xa,%rdx
add    %rax,%rdx
imul   $0x7d7,%rdx,%rax
sub    %rax,%rcx
mov    %rcx,0x0(%rbp,%rbx,8)
add    $0x1,%rbx
add    $0xffffffffffffffff,%rsi
mov    %rdi,%rcx
cmp    %rbx,%r11
jne    12c0 <main+0x190>
jmp    1280 <main+0x150>

LLVM IR (-O0):
@.str = [10 x i8] c"%lld %lld\00"
@EOF = global i64 0
@.str.1 = [6 x i8] c"%lld\0A\00"
define i32 @main() {
%1 = alloca i32
%2 = alloca i64
%3 = alloca i64
%4 = alloca i64
%5 = alloca i64
%6 = alloca i64
%7 = alloca i64
%8 = alloca i64
%9 = alloca i64
%10 = alloca [110 x [110 x i64]]
store i32 0, i32* %1
br label %11
11: ; preds = %80, %0
%12 = call i64 @scanf(i8* getelementptr inbounds ([10 x i8], [10 x i8]* @.str, i64 0, i64 0), i64* %8, i64* %9)
%13 = load i64, i64* @EOF
%14 = icmp ne i64 %12, %13
br i1 %14, label %15, label %88
15: ; preds = %11
store i64 1, i64* %2
br label %16
16: ; preds = %26, %15
%17 = load i64, i64* %2
%18 = icmp sle i64 %17, 101
br i1 %18, label %19, label %29
19: ; preds = %16
%20 = load i64, i64* %2
%21 = getelementptr inbounds [110 x [110 x i64]], [110 x [110 x i64]]* %10, i64 0, i64 %20
%22 = getelementptr inbounds [110 x i64], [110 x i64]* %21, i64 0, i64 0
store i64 1, i64* %22
%23 = getelementptr inbounds [110 x [110 x i64]], [110 x [110 x i64]]* %10, i64 0, i64 0
%24 = load i64, i64* %2
%25 = getelementptr inbounds [110 x i64], [110 x i64]* %23, i64 0, i64 %24
store i64 0, i64* %25
br label %26
26: ; preds = %19
%27 = load i64, i64* %2
%28 = add nsw i64 %27, 1
store i64 %28, i64* %2
br label %16
29: ; preds = %16
%30 = getelementptr inbounds [110 x [110 x i64]], [110 x [110 x i64]]* %10, i64 0, i64 0
%31 = getelementptr inbounds [110 x i64], [110 x i64]* %30, i64 0, i64 0
store i64 0, i64* %31
store i64 1, i64* %2
br label %32
32: ; preds = %77, %29
%33 = load i64, i64* %2
%34 = load i64, i64* %8
%35 = icmp sle i64 %33, %34
br i1 %35, label %36, label %80
36: ; preds = %32
store i64 1, i64* %3
br label %37
37: ; preds = %73, %36
%38 = load i64, i64* %3
%39 = load i64, i64* %9
%40 = icmp sle i64 %38, %39
br i1 %40, label %41, label %76
41: ; preds = %37
%42 = load i64, i64* %2
%43 = sub nsw i64 %42, 1
%44 = getelementptr inbounds [110 x [110 x i64]], [110 x [110 x i64]]* %10, i64 0, i64 %43
%45 = load i64, i64* %3
%46 = getelementptr inbounds [110 x i64], [110 x i64]* %44, i64 0, i64 %45
%47 = load i64, i64* %46
%48 = load i64, i64* %3
%49 = add nsw i64 %48, 1
%50 = mul nsw i64 %47, %49
%51 = load i64, i64* %2
%52 = sub nsw i64 %51, 1
%53 = getelementptr inbounds [110 x [110 x i64]], [110 x [110 x i64]]* %10, i64 0, i64 %52
%54 = load i64, i64* %3
%55 = sub nsw i64 %54, 1
%56 = getelementptr inbounds [110 x i64], [110 x i64]* %53, i64 0, i64 %55
%57 = load i64, i64* %56
%58 = load i64, i64* %2
%59 = load i64, i64* %3
%60 = sub nsw i64 %58, %59
%61 = mul nsw i64 %57, %60
%62 = add nsw i64 %50, %61
%63 = load i64, i64* %2
%64 = getelementptr inbounds [110 x [110 x i64]], [110 x [110 x i64]]* %10, i64 0, i64 %63
%65 = load i64, i64* %3
%66 = getelementptr inbounds [110 x i64], [110 x i64]* %64, i64 0, i64 %65
store i64 %62, i64* %66
%67 = load i64, i64* %2
%68 = getelementptr inbounds [110 x [110 x i64]], [110 x [110 x i64]]* %10, i64 0, i64 %67
%69 = load i64, i64* %3
%70 = getelementptr inbounds [110 x i64], [110 x i64]* %68, i64 0, i64 %69
%71 = load i64, i64* %70
%72 = srem i64 %71, 2007
store i64 %72, i64* %70
br label %73
73: ; preds = %41
%74 = load i64, i64* %3
%75 = add nsw i64 %74, 1
store i64 %75, i64* %3
br label %37
76: ; preds = %37
br label %77
77: ; preds = %76
%78 = load i64, i64* %2
%79 = add nsw i64 %78, 1
store i64 %79, i64* %2
br label %32
80: ; preds = %32
%81 = load i64, i64* %8
%82 = getelementptr inbounds [110 x [110 x i64]], [110 x [110 x i64]]* %10, i64 0, i64 %81
%83 = load i64, i64* %9
%84 = getelementptr inbounds [110 x i64], [110 x i64]* %82, i64 0, i64 %83
%85 = load i64, i64* %84
%86 = srem i64 %85, 2007
%87 = call i32 @printf(i8* getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), i64 %86)
br label %11
88: ; preds = %11
ret i32 0
}

LLVM IR (-O2):
@.str = [10 x i8] c"%lld %lld\00"
@EOF = global i64 0
@.str.1 = [6 x i8] c"%lld\0A\00"
define i32 @main() {
%1 = alloca i64
%2 = alloca i64
%3 = alloca [110 x [110 x i64]]
%4 = bitcast i64* %1 to i8*
%5 = bitcast i64* %2 to i8*
%6 = bitcast [110 x [110 x i64]]* %3 to i8*
%7 = call i64 @scanf(i8* getelementptr inbounds ([10 x i8], [10 x i8]* @.str, i64 0, i64 0), i64* %1, i64* %2)
%8 = load i64, i64* @EOF
%9 = icmp eq i64 %7, %8
br i1 %9, label %63, label %10
10: ; preds = %0
%11 = getelementptr inbounds [110 x [110 x i64]], [110 x [110 x i64]]* %3, i64 0, i64 0, i64 0
%12 = getelementptr inbounds [110 x [110 x i64]], [110 x [110 x i64]]* %3, i64 0, i64 0, i64 1
%13 = bitcast i64* %12 to i8*
br label %14
14: ; preds = %10, %55
br label %15
15: ; preds = %20, %14
%16 = phi i64 [ 1, %14 ], [ %26, %20 ]
%17 = getelementptr inbounds [110 x [110 x i64]], [110 x [110 x i64]]* %3, i64 0, i64 %16, i64 0
store i64 1, i64* %17
%18 = add nuw nsw i64 %16, 1
%19 = icmp eq i64 %18, 102
br i1 %19, label %27, label %20
20: ; preds = %15
%21 = getelementptr inbounds [110 x [110 x i64]], [110 x [110 x i64]]* %3, i64 0, i64 %18, i64 0
store i64 1, i64* %21
%22 = add nuw nsw i64 %16, 2
%23 = getelementptr inbounds [110 x [110 x i64]], [110 x [110 x i64]]* %3, i64 0, i64 %22, i64 0
store i64 1, i64* %23
%24 = add nuw nsw i64 %16, 3
%25 = getelementptr inbounds [110 x [110 x i64]], [110 x [110 x i64]]* %3, i64 0, i64 %24, i64 0
store i64 1, i64* %25
%26 = add nuw nsw i64 %16, 4
br label %15
27: ; preds = %15
store i64 0, i64* %11
%28 = load i64, i64* %1
%29 = icmp slt i64 %28, 1
%30 = load i64, i64* %2
br i1 %29, label %55, label %31
31: ; preds = %27
%32 = icmp slt i64 %30, 1
br label %33
33: ; preds = %31, %52
%34 = phi i64 [ 1, %31 ], [ %53, %52 ]
br i1 %32, label %52, label %35
35: ; preds = %33
%36 = add nsw i64 %34, -1
%37 = getelementptr inbounds [110 x [110 x i64]], [110 x [110 x i64]]* %3, i64 0, i64 %36, i64 0
%38 = load i64, i64* %37
br label %39
39: ; preds = %35, %39
%40 = phi i64 [ %38, %35 ], [ %43, %39 ]
%41 = phi i64 [ 1, %35 ], [ %44, %39 ]
%42 = getelementptr inbounds [110 x [110 x i64]], [110 x [110 x i64]]* %3, i64 0, i64 %36, i64 %41
%43 = load i64, i64* %42
%44 = add nuw i64 %41, 1
%45 = mul nsw i64 %43, %44
%46 = sub nsw i64 %34, %41
%47 = mul nsw i64 %40, %46
%48 = add nsw i64 %47, %45
%49 = getelementptr inbounds [110 x [110 x i64]], [110 x [110 x i64]]* %3, i64 0, i64 %34, i64 %41
%50 = srem i64 %48, 2007
store i64 %50, i64* %49
%51 = icmp eq i64 %41, %30
br i1 %51, label %52, label %39
52: ; preds = %39, %33
%53 = add nuw i64 %34, 1
%54 = icmp eq i64 %34, %28
br i1 %54, label %55, label %33
55: ; preds = %52, %27
%56 = getelementptr inbounds [110 x [110 x i64]], [110 x [110 x i64]]* %3, i64 0, i64 %28, i64 %30
%57 = load i64, i64* %56
%58 = srem i64 %57, 2007
%59 = call i32 @printf(i8* getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), i64 %58)
%60 = call i64 @scanf(i8* getelementptr inbounds ([10 x i8], [10 x i8]* @.str, i64 0, i64 0), i64* %1, i64* %2)
%61 = load i64, i64* @EOF
%62 = icmp eq i64 %60, %61
br i1 %62, label %63, label %14
63: ; preds = %55, %0
ret i32 0
}