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

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

def c(n, k):
    ue, sita = 1, 1
    for i in range(0, k):
        ue = ue * (n - i) % mod
       	sita = sita * (i + 1) % mod
    return ue * inv(sita) % mod


print(((pow(2, n, mod) - c(n, a) - c(n, b) - 1) % mod + mod) % mod) n, a, b = map(int, input().split()) n map(int, input().split()) map int input().split() ().split () input split a b mod = 10 ** 9 + 7 mod 10 ** 9 + 7 10 ** 9 10 9 7 def inv(x):
    return pow(x, mod - 2, mod) inv return pow(x, mod - 2, mod) pow(x, mod - 2, mod) pow x mod - 2 mod 2 mod x x def c(n, k):
    ue, sita = 1, 1
    for i in range(0, k):
        ue = ue * (n - i) % mod
       	sita = sita * (i + 1) % mod
    return ue * inv(sita) % mod c ue, sita = 1, 1 ue 1 sita 1 for i in range(0, k):
        ue = ue * (n - i) % mod
       	sita = sita * (i + 1) % mod
     i range(0, k) range 0 k ue = ue * (n - i) % mod ue ue * (n - i) % mod ue * (n - i) ue n - i n i mod sita = sita * (i + 1) % mod sita sita * (i + 1) % mod sita * (i + 1) sita i + 1 i 1 mod return ue * inv(sita) % mod ue * inv(sita) % mod ue * inv(sita) ue inv(sita) inv sita mod n n k k print(((pow(2, n, mod) - c(n, a) - c(n, b) - 1) % mod + mod) % mod) print ((pow(2, n, mod) - c(n, a) - c(n, b) - 1) % mod + mod) % mod (pow(2, n, mod) - c(n, a) - c(n, b) - 1) % mod + mod (pow(2, n, mod) - c(n, a) - c(n, b) - 1) % mod pow(2, n, mod) - c(n, a) - c(n, b) - 1 pow(2, n, mod) - c(n, a) - c(n, b) pow(2, n, mod) - c(n, a) pow(2, n, mod) pow 2 n mod c(n, a) c n a c(n, b) c n b 1 mod mod mod a, b = map(int, input().split()) map(int, input().split()) a n, a, b = map(int, input().split()) map(int, input().split()) n b = map(int, input().split()) map(int, input().split()) b def inv(x):
    return pow(x, mod - 2, mod) def inv(x):
    return pow(x, mod - 2, mod) inv def c(n, k):
    ue, sita = 1, 1
    for i in range(0, k):
        ue = ue * (n - i) % mod
       	sita = sita * (i + 1) % mod
    return ue * inv(sita) % mod def c(n, k):
    ue, sita = 1, 1
    for i in range(0, k):
        ue = ue * (n - i) % mod
       	sita = sita * (i + 1) % mod
    return ue * inv(sita) % mod c mod = 10 ** 9 + 7 10 ** 9 + 7 mod
n, a, b = map(int, input().split())
mod = 10 ** 9 + 7

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

def c(n, k):
    ue, sita = 1, 1
    for i in range(0, k):
        ue = ue * (n - i) % mod
       	sita = sita * (i + 1) % mod
    return ue * inv(sita) % mod


print(((pow(2, n, mod) - c(n, a) - c(n, b) - 1) % mod + mod) % mod)
