#include <bits/stdc++.h>
 
using namespace std;

const long long int MOD = 1e9 + 7;
long long int f[2000005];
long long int sum[2000005];

int main(void)
{
	cin.tie(0);
	ios::sync_with_stdio(false);

	long long int n,p;

	cin >> n >> p;

	f[1] = 0;
	f[2] = 1;
	for(int i=3;i<=2000000;i++)
	{
		long long int val = p;
		val = ((val%MOD)*(f[i-1]%MOD)%MOD);
		val%=MOD;
		val += f[i-2];
		val%=MOD;
		f[i] = val;
	}

	sum[0] = 0;
	for(int i=1;i<=2000000;i++)
	{
		sum[i] += sum[i-1];
		sum[i] += f[i];
		sum[i]%=MOD;
	}

	long long int res = 0;

	for(int i=1;i<=n;i++)
	{
		long long int val = f[i];
		val = ((val%MOD)*(sum[i]%MOD)%MOD);
		val%=MOD;
		res += val;
		res%=MOD;
	}

	cout << res << '\n';

	return 0;
}