#include <iostream>
#include <algorithm>
using namespace std;

    long long a[2000000],p,n,sum=0,expsum=0;
    const long long MOD=1000000007;
    bool oe=false;

long long llmod(long long x){
    oe=(oe^((x/MOD)%2));
    return x%MOD;
}
int main(void){
    // Your code here!
    cin>>n>>p;
    for(long long i=0;i<n;i++){
        if(i<2)a[i]=i;
        else a[i]=llmod(llmod(a[i-1]*p)+a[i-2]);
        sum=llmod(sum+a[i]);
        expsum=llmod(expsum+llmod(a[i]*a[i]));
    }
    cout<<llmod((llmod(llmod(sum*sum)+expsum)+(oe?MOD:0))/2)<<endl;
    return 0;
}