#include<bits/stdc++.h>
using namespace std;
using LL = long long;
using ULL = unsigned long long;
#define rep(i,n) for(int i=0; i<(n); i++)

const ULL M = 1000000007;

const int matrix_size = 3;

struct Matrix {
	ULL X[matrix_size][matrix_size] = {};
	static Matrix id() {
		Matrix res;
		rep(i, matrix_size) res.X[i][i] = 1;
		return res;
	}
};

Matrix operator*(Matrix l, Matrix r) {
	Matrix res;
	rep(i, matrix_size) rep(j, matrix_size) rep(k, matrix_size)
		res.X[i][j] = (res.X[i][j] + l.X[i][k] * r.X[k][j]) % M;
	return res;
}

Matrix operator^(Matrix l, ULL i) {
	if (i == 0) return Matrix::id();
	auto res = (l * l) ^ (i / 2);
	if (i % 2 == 1) res = res * l;
	return res;
}

int main() {
	ULL N, I; cin >> N >> I;

	Matrix X = { { {1,1,0},{1,0,0},{0,0,1} } };
	X = X ^ I;

	X.X[0][2] = 1;
	X = X ^ (N + 1);

	cout << X.X[1][2] << endl;

	return 0;
}