import math
mod = 10**9 + 7
n,a,b=map(int,input().split())
k=a
modinv_table = [-1] * (k+1)
modinv_table[1] = 1
for i in range(2, k+1):
    modinv_table[i] = (-modinv_table[mod % i] * (mod // i)) % mod

def binomial_coefficients(n, k):
    ans = 1
    for i in range(k):
        ans *= n-i
        ans *= modinv_table[i + 1]
        ans %= mod
    return ans
d=binomial_coefficients(n, a)
k=b

modinv_table = [-1] * (k+1)
modinv_table[1] = 1
for i in range(2, k+1):
    modinv_table[i] = (-modinv_table[mod % i] * (mod // i)) % mod

def binomial_coefficients(n, k):
    ans = 1
    for i in range(k):
        ans *= n-i
        ans *= modinv_table[i + 1]
        ans %= mod
    return ans
e=binomial_coefficients(n, b)
c = pow(2,n%(10**9+7), 10**9+7)
print((c-d-e-1)%(10**9+7)) import math math 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 k=a k a modinv_table = [-1] * (k+1) modinv_table [-1] * (k+1) [-1] -1 k+1 k 1 modinv_table[1] = 1 modinv_table[1] modinv_table 1 1 for i in range(2, k+1):
    modinv_table[i] = (-modinv_table[mod % i] * (mod // i)) % mod i range(2, k+1) range 2 k+1 k 1 modinv_table[i] = (-modinv_table[mod % i] * (mod // i)) % mod modinv_table[i] modinv_table i (-modinv_table[mod % i] * (mod // i)) % mod -modinv_table[mod % i] * (mod // i) -modinv_table[mod % i] modinv_table[mod % i] modinv_table mod % i mod i mod // i mod i mod def binomial_coefficients(n, k):
    ans = 1
    for i in range(k):
        ans *= n-i
        ans *= modinv_table[i + 1]
        ans %= mod
    return ans binomial_coefficients ans = 1 ans 1 for i in range(k):
        ans *= n-i
        ans *= modinv_table[i + 1]
        ans %= mod
     i range(k) range k ans *= n-i ans n-i n i ans *= modinv_table[i + 1] ans modinv_table[i + 1] modinv_table i + 1 i 1 ans %= mod ans mod return ans ans n n k k d=binomial_coefficients(n, a) d binomial_coefficients(n, a) binomial_coefficients n a k=b k b modinv_table = [-1] * (k+1) modinv_table [-1] * (k+1) [-1] -1 k+1 k 1 modinv_table[1] = 1 modinv_table[1] modinv_table 1 1 for i in range(2, k+1):
    modinv_table[i] = (-modinv_table[mod % i] * (mod // i)) % mod i range(2, k+1) range 2 k+1 k 1 modinv_table[i] = (-modinv_table[mod % i] * (mod // i)) % mod modinv_table[i] modinv_table i (-modinv_table[mod % i] * (mod // i)) % mod -modinv_table[mod % i] * (mod // i) -modinv_table[mod % i] modinv_table[mod % i] modinv_table mod % i mod i mod // i mod i mod def binomial_coefficients(n, k):
    ans = 1
    for i in range(k):
        ans *= n-i
        ans *= modinv_table[i + 1]
        ans %= mod
    return ans binomial_coefficients ans = 1 ans 1 for i in range(k):
        ans *= n-i
        ans *= modinv_table[i + 1]
        ans %= mod
     i range(k) range k ans *= n-i ans n-i n i ans *= modinv_table[i + 1] ans modinv_table[i + 1] modinv_table i + 1 i 1 ans %= mod ans mod return ans ans n n k k e=binomial_coefficients(n, b) e binomial_coefficients(n, b) binomial_coefficients n b c = pow(2,n%(10**9+7), 10**9+7) c pow(2,n%(10**9+7), 10**9+7) pow 2 n%(10**9+7) n 10**9+7 10**9 10 9 7 10**9+7 10**9 10 9 7 print((c-d-e-1)%(10**9+7)) print (c-d-e-1)%(10**9+7) c-d-e-1 c-d-e c-d c d e 1 10**9+7 10**9 10 9 7 k=b b k def binomial_coefficients(n, k):
    ans = 1
    for i in range(k):
        ans *= n-i
        ans *= modinv_table[i + 1]
        ans %= mod
    return ans def binomial_coefficients(n, k):
    ans = 1
    for i in range(k):
        ans *= n-i
        ans *= modinv_table[i + 1]
        ans %= mod
    return ans binomial_coefficients c = pow(2,n%(10**9+7), 10**9+7) pow(2,n%(10**9+7), 10**9+7) c a,b=map(int,input().split()) map(int,input().split()) a def binomial_coefficients(n, k):
    ans = 1
    for i in range(k):
        ans *= n-i
        ans *= modinv_table[i + 1]
        ans %= mod
    return ans def binomial_coefficients(n, k):
    ans = 1
    for i in range(k):
        ans *= n-i
        ans *= modinv_table[i + 1]
        ans %= mod
    return ans binomial_coefficients b=map(int,input().split()) map(int,input().split()) b n,a,b=map(int,input().split()) map(int,input().split()) n modinv_table = [-1] * (k+1) [-1] * (k+1) modinv_table e=binomial_coefficients(n, b) binomial_coefficients(n, b) e mod = 10**9 + 7 10**9 + 7 mod k=a a k modinv_table = [-1] * (k+1) [-1] * (k+1) modinv_table d=binomial_coefficients(n, a) binomial_coefficients(n, a) d
import math
mod = 10**9 + 7
n,a,b=map(int,input().split())
k=a
modinv_table = [-1] * (k+1)
modinv_table[1] = 1
for i in range(2, k+1):
    modinv_table[i] = (-modinv_table[mod % i] * (mod // i)) % mod

def binomial_coefficients(n, k):
    ans = 1
    for i in range(k):
        ans *= n-i
        ans *= modinv_table[i + 1]
        ans %= mod
    return ans
d=binomial_coefficients(n, a)
k=b

modinv_table = [-1] * (k+1)
modinv_table[1] = 1
for i in range(2, k+1):
    modinv_table[i] = (-modinv_table[mod % i] * (mod // i)) % mod

def binomial_coefficients(n, k):
    ans = 1
    for i in range(k):
        ans *= n-i
        ans *= modinv_table[i + 1]
        ans %= mod
    return ans
e=binomial_coefficients(n, b)
c = pow(2,n%(10**9+7), 10**9+7)
print((c-d-e-1)%(10**9+7))