n = 2
p = 1
q = 0

mod = 1000000007
jie = []
ni = []

def power(a,n,mod):
    p = 1
    while n>0:
        if n%2==1:
            p = p*a
            p = p%mod
        n = n//2
        a = a*a%mod
    return p%mod

def C(n,k):
    return jie[n]*ni[k]%mod*ni[n-k]%mod

jie.append(1)
ni.append(1)
for i in range(1,10):
    sg = jie[i-1]*i%mod
    jie.append(sg)
    nio = power(sg, mod-2, mod)
    ni.append(nio)

sum =0
for i in range(p,n-q+1):
    sum = (sum + C(n,i)*(i%mod)%mod)%mod

#print(sum)
mu = 0
for i in range(p,n-q+1):
    mu = (mu + C(n,i))%mod
#print(mu)
ans = sum * power(mu, mod-2, mod)%mod
print(ans)

AC代码