#include<vector>
#include<iostream>

using namespace std;
#define MOD 1000000007

int main(){
    long long int n,p; cin >> n >> p;
    vector<long long int>  v(n), ac(n);
    v[0]=0;v[1]=1;
    for(int i=2;i<n;i++){
        v[i] = ((p * v[i-1]) % MOD + (v[i-2])) % MOD;
    }
    ac[0] = 0;
    for(int i=1;i<n;i++) ac[i] = (ac[i-1] + v[i]) % MOD;
    long long int ans = 0;
    for(int i=0;i<n;i++) ans = (ans + ac[i] * v[i]) % MOD;
    cout << ans << endl;
    
    return 0;
}