import sys
input=lambda: sys.stdin.readline().rstrip()
n,p=map(int,input().split())
mod=10**9+7
if n==1:
  print(0)
elif n==2:
  print(1)
else:
  A=[0,1]
  B=[0,1]
  for _ in range(n-2):
    a=p*A[-1]+A[-2]
    if a>=mod:
      a%=mod
    b=pow(a,2,mod)
    A.append(a)
    B.append(b)
  a,b=0,0
  for i in range(n):
    a+=A[i]
    if a>=mod:
      a%=mod
    b+=B[i]
    if b>=mod:
      b%=mod
  print(((a**2+b)*pow(2,mod-2,mod))%mod)