def comb_mod(n,r,m):
    ans = 1
    for i in range(1,r+1):
        ans *= (n-i+1) % m
        ans *= pow(i,m-2,m) % m
        ans = ans % m
    return ans
n,a,b = map(int,input().split())
m = 10**9 + 7
ans = pow(2,n,m) - comb_mod(n,a,m) - comb_mod(n,b,m) - 1
print(int(ans) % m) def comb_mod(n,r,m):
    ans = 1
    for i in range(1,r+1):
        ans *= (n-i+1) % m
        ans *= pow(i,m-2,m) % m
        ans = ans % m
    return ans comb_mod ans = 1 ans 1 for i in range(1,r+1):
        ans *= (n-i+1) % m
        ans *= pow(i,m-2,m) % m
        ans = ans % m
     i range(1,r+1) range 1 r+1 r 1 ans *= (n-i+1) % m ans (n-i+1) % m n-i+1 n-i n i 1 m ans *= pow(i,m-2,m) % m ans pow(i,m-2,m) % m pow(i,m-2,m) pow i m-2 m 2 m m ans = ans % m ans ans % m ans m return ans ans n n r r m m n,a,b = map(int,input().split()) n map(int,input().split()) map int input().split() ().split () input split a b m = 10**9 + 7 m 10**9 + 7 10**9 10 9 7 ans = pow(2,n,m) - comb_mod(n,a,m) - comb_mod(n,b,m) - 1 ans pow(2,n,m) - comb_mod(n,a,m) - comb_mod(n,b,m) - 1 pow(2,n,m) - comb_mod(n,a,m) - comb_mod(n,b,m) pow(2,n,m) - comb_mod(n,a,m) pow(2,n,m) pow 2 n m comb_mod(n,a,m) comb_mod n a m comb_mod(n,b,m) comb_mod n b m 1 print(int(ans) % m) print int(ans) % m int(ans) int ans m def comb_mod(n,r,m):
    ans = 1
    for i in range(1,r+1):
        ans *= (n-i+1) % m
        ans *= pow(i,m-2,m) % m
        ans = ans % m
    return ans def comb_mod(n,r,m):
    ans = 1
    for i in range(1,r+1):
        ans *= (n-i+1) % m
        ans *= pow(i,m-2,m) % m
        ans = ans % m
    return ans comb_mod m = 10**9 + 7 10**9 + 7 m a,b = map(int,input().split()) map(int,input().split()) a ans = pow(2,n,m) - comb_mod(n,a,m) - comb_mod(n,b,m) - 1 pow(2,n,m) - comb_mod(n,a,m) - comb_mod(n,b,m) - 1 ans n,a,b = map(int,input().split()) map(int,input().split()) n b = map(int,input().split()) map(int,input().split()) b
def comb_mod(n,r,m):
    ans = 1
    for i in range(1,r+1):
        ans *= (n-i+1) % m
        ans *= pow(i,m-2,m) % m
        ans = ans % m
    return ans
n,a,b = map(int,input().split())
m = 10**9 + 7
ans = pow(2,n,m) - comb_mod(n,a,m) - comb_mod(n,b,m) - 1
print(int(ans) % m)