Problem:
Score : 800 points

Problem Statement

Sig has built his own keyboard. Designed for ultimate simplicity, this keyboard only has 3 keys on it: the 0 key, the 1 key and the backspace key.

To begin with, he is using a plain text editor with this keyboard. This editor always displays one string (possibly empty). Just after the editor is launched, this string is empty. When each key on the keyboard is pressed, the following changes occur to the string:

The 0 key: a letter 0 will be inserted to the right of the string.

The 1 key: a letter 1 will be inserted to the right of the string.

The backspace key: if the string is empty, nothing happens. Otherwise, the rightmost letter of the string is deleted.

Sig has launched the editor, and pressed these keys N times in total. As a result, the editor displays a string s. Find the number of such ways to press the keys, modulo 10^9 + 7.

Constraints

1 ≦ N ≦ 5000

1 ≦ |s| ≦ N

s consists of the letters 0 and 1.

Partial Score

400 points will be awarded for passing the test set satisfying 1 ≦ N ≦ 300.

Input

The input is given from Standard Input in the following format:

N
s

Output

Print the number of the ways to press the keys N times in total such that the editor displays the string s in the end, modulo 10^9+7.

Sample Input 1

3
0

Sample Output 1

5

We will denote the backspace key by B. The following 5 ways to press the keys will cause the editor to display the string 0 in the end: 00B, 01B, 0B0, 1B0, BB0. In the last way, nothing will happen when the backspace key is pressed.

Sample Input 2

300
1100100

Sample Output 2

519054663

Sample Input 3

5000
01000001011101000100001101101111011001000110010101110010000

Sample Output 3

500886057

Solution in C++14 (GCC 5.4.1) (judge verdict: Accepted)

#include <iostream>
#include <algorithm>
#include <array>
#include <cstdint>
#include <climits>
#include <functional>
#include <map>
#include <math.h>
#include <queue>
#include <set>
#include <stack>
#include <stdlib.h>
#include <string>
#include <time.h>
#include <type_traits>
#include <utility>
#include <vector>

using int32 = std::int_fast32_t;
using int64 = std::int_fast64_t;
using uint32 = std::uint_fast32_t;
using uint64 = std::uint_fast64_t;
using intl32 = std::int_least32_t;
using intl64 = std::int_least64_t;
using uintl32 = std::uint_least32_t;
using uintl64 = std::uint_least64_t;

template<uint64 MOD>
class modint {
public:
	uint32 a;
	modint(uint32 x = 0) :a(x < MOD ? x : x%MOD) {}
	const modint operator+(const modint &o) const{
		if (a + o.a < MOD) return modint(a + o.a);
		return modint(a + o.a - MOD);
	}
	modint &operator+=(const modint &o) {
		a += o.a;
		if (a >= MOD) a -= MOD;
		return *this;
	}
	const modint operator-(const modint &o) const{
		if (a < o.a) return modint(a + MOD - o.a);
		return modint(a - o.a);
	}
	modint &operator-=(const modint &o) {
		if (a < o.a) a += MOD;
		a -= o.a;
		return *this;
	}
	const modint operator*(const modint &o) const{
		return modint((uint64)a*(uint64)o.a%MOD);
	}
	modint &operator*=(const modint &o) {
		a = (uint64)a*(uint64)o.a%MOD;
		return *this;
	}
	const modint operator/(const modint &o) const {
		return *this*~o;
	}

	const modint operator-() const{
		if (a) return modint(MOD - a);
		return modint(0);
	}
	modint operator++() {
		++a;
		if (a == MOD) a = 0;
		return *this;
	}
	modint operator++(int) {
		if (a == MOD - 1) {
			a = 0;
			return modint(MOD - 1);
		}
		++a;
		return modint(a - 1);
	}
	modint operator--() {
		if (a) {
			--a;
		}
		else {
			a = MOD - 1;
		}
		return *this;
	}
	modint operator--(int) {
		if (a) {
			--a;
			return modint(a + 1);
		}
		a = MOD - 1;
		return modint(0);
	}
	modint operator~() const {
		int64 x = a, y = MOD, u = 1, v = 0, t;
		while (y) {
			t = x / y;
			x -= t*y;
			std::swap(x, y);
			u -= t*v;
			std::swap(u, v);
		}
		return modint(u);
	}
	bool operator!() const{
		if (a) return false;
		return true;
	}
	bool operator==(const modint o) const{
		return a == o.a;
	}
	bool operator!=(const modint o) const{
		return a != o.a;
	}
	bool operator<(const modint o) const{
		return a < o.a;
	}
	bool operator<=(const modint o) const{
		return a <= o.a;
	}
	bool operator>(const modint o) const{
		return a > o.a;
	}
	bool operator>=(const modint o) const{
		return a >= o.a;
	}
};

using mint = modint<1000000007>;
int main(void) {
	std::ios::sync_with_stdio(false);
	std::cin.tie(0);
	uint32 n;
	std::string s;
	std::cin >> n >> s;
	mint dp[5003][5003] = {};
	dp[0][0] = mint(1);
	for (uint32 i = 1;i <= n;++i) {
		dp[i][0] = dp[i - 1][1] + dp[i - 1][1] + dp[i - 1][0];
		for (uint32 j = 1;j <= i;++j) {
			dp[i][j] = dp[i - 1][j - 1] + dp[i - 1][j + 1] + dp[i - 1][j + 1];
		}
	}
	std::cout << dp[n][s.length()].a << "\n";
	return 0;
}