#include <bits/stdc++.h>
using namespace std;
int main(){
    long long N,p;
    cin>>N>>p;
    long long MOD=1000000007;
    long long temp=0;
    vector<long long> a(N);
    long long ans=1,sum=1;
    a[1]=1;
    for(int i=2;i<N;++i){
        a[i]=(a[i-1]*p+a[i-2])%MOD;
        sum=(sum+a[i])%MOD;
        ans=(ans+a[i]*a[i])%MOD;
    }
    ans=(ans+sum*sum)%MOD;
    ans=(ans*500000004)%MOD;
    cout<<ans<<endl;
}