#include using namespace std; #define REP(i,a,n) for(int i=(a); i<(int)(n); i++) #define rep(i,n) REP(i,0,n) #define FOR(it,c) for(__typeof((c).begin()) it=(c).begin(); it!=(c).end(); ++it) #define ALLOF(c) (c).begin(), (c).end() typedef long long ll; typedef unsigned long long ull; static const ll MOD = 1000000007; ll A[2000005]; ll B[2000005]; int main(){ ll N, p; cin >> N >> p; A[1] = 0; B[1] = 0; A[2] = 1; B[2] = 1; REP(i,3,N+1){ A[i] = ((p * A[i-1])%MOD) + A[i-2]; A[i] %= MOD; B[i] = B[i-1] + A[i]; B[i] %= MOD; } ll ret = 0; REP(i,1,N+1){ ret += (A[i] * B[i])%MOD; ret %= MOD; } cout << ret << endl; return 0; }