def mod_inverse(x, mod):
    return pow(x, mod - 2, mod)


# return nCk % mod
# it takes O(k)
def mod_comb(n, k, mod):
    numer, denom = 1, 1
    for i in range(k):
        numer = numer * ((n - i) % mod) % mod
        denom = denom * ((i + 1) % mod) % mod

    return numer * mod_inverse(denom, mod) % mod


mod = 10**9 + 7
n, a, b = map(int, input().split())

ans = pow(2, n, mod) - 1
ans = (ans - mod_comb(n, a, mod)) % mod
ans = (ans - mod_comb(n, b, mod)) % mod
ans = (ans + mod) % mod

print(ans) def mod_inverse(x, mod):
    return pow(x, mod - 2, mod)


# return nCk % mod
# it takes O(k) mod_inverse return pow(x, mod - 2, mod) pow(x, mod - 2, mod) pow x mod - 2 mod 2 mod x x mod mod def mod_comb(n, k, mod):
    numer, denom = 1, 1
    for i in range(k):
        numer = numer * ((n - i) % mod) % mod
        denom = denom * ((i + 1) % mod) % mod

    return numer * mod_inverse(denom, mod) % mod mod_comb numer, denom = 1, 1 numer 1 denom 1 for i in range(k):
        numer = numer * ((n - i) % mod) % mod
        denom = denom * ((i + 1) % mod) % mod

     i range(k) range k numer = numer * ((n - i) % mod) % mod numer numer * ((n - i) % mod) % mod numer * ((n - i) % mod) numer (n - i) % mod n - i n i mod mod denom = denom * ((i + 1) % mod) % mod denom denom * ((i + 1) % mod) % mod denom * ((i + 1) % mod) denom (i + 1) % mod i + 1 i 1 mod mod return numer * mod_inverse(denom, mod) % mod numer * mod_inverse(denom, mod) % mod numer * mod_inverse(denom, mod) numer mod_inverse(denom, mod) mod_inverse denom mod mod n n k k mod mod mod = 10**9 + 7 mod 10**9 + 7 10**9 10 9 7 n, a, b = map(int, input().split()) n map(int, input().split()) map int input().split() ().split () input split a b ans = pow(2, n, mod) - 1 ans pow(2, n, mod) - 1 pow(2, n, mod) pow 2 n mod 1 ans = (ans - mod_comb(n, a, mod)) % mod ans (ans - mod_comb(n, a, mod)) % mod ans - mod_comb(n, a, mod) ans mod_comb(n, a, mod) mod_comb n a mod mod ans = (ans - mod_comb(n, b, mod)) % mod ans (ans - mod_comb(n, b, mod)) % mod ans - mod_comb(n, b, mod) ans mod_comb(n, b, mod) mod_comb n b mod mod ans = (ans + mod) % mod ans (ans + mod) % mod ans + mod ans mod mod print(ans) print ans ans = pow(2, n, mod) - 1 pow(2, n, mod) - 1 ans def mod_comb(n, k, mod):
    numer, denom = 1, 1
    for i in range(k):
        numer = numer * ((n - i) % mod) % mod
        denom = denom * ((i + 1) % mod) % mod

    return numer * mod_inverse(denom, mod) % mod def mod_comb(n, k, mod):
    numer, denom = 1, 1
    for i in range(k):
        numer = numer * ((n - i) % mod) % mod
        denom = denom * ((i + 1) % mod) % mod

    return numer * mod_inverse(denom, mod) % mod mod_comb def mod_inverse(x, mod):
    return pow(x, mod - 2, mod)


# return nCk % mod
# it takes O(k) def mod_inverse(x, mod):
    return pow(x, mod - 2, mod)


# return nCk % mod
# it takes O(k) mod_inverse ans = (ans + mod) % mod (ans + mod) % mod ans mod = 10**9 + 7 10**9 + 7 mod a, b = map(int, input().split()) map(int, input().split()) a ans = (ans - mod_comb(n, a, mod)) % mod (ans - mod_comb(n, a, mod)) % mod ans n, a, b = map(int, input().split()) map(int, input().split()) n b = map(int, input().split()) map(int, input().split()) b ans = (ans - mod_comb(n, b, mod)) % mod (ans - mod_comb(n, b, mod)) % mod ans
def mod_inverse(x, mod):
    return pow(x, mod - 2, mod)


# return nCk % mod
# it takes O(k)
def mod_comb(n, k, mod):
    numer, denom = 1, 1
    for i in range(k):
        numer = numer * ((n - i) % mod) % mod
        denom = denom * ((i + 1) % mod) % mod

    return numer * mod_inverse(denom, mod) % mod


mod = 10**9 + 7
n, a, b = map(int, input().split())

ans = pow(2, n, mod) - 1
ans = (ans - mod_comb(n, a, mod)) % mod
ans = (ans - mod_comb(n, b, mod)) % mod
ans = (ans + mod) % mod

print(ans)