Problem: In "Takahashi-ya", a ramen restaurant, basically they have one menu: "ramen", but N kinds of toppings are also offered. When a customer orders a bowl of ramen, for each kind of topping, he/she can choose whether to put it on top of his/her ramen or not. There is no limit on the number of toppings, and it is allowed to have all kinds of toppings or no topping at all. That is, considering the combination of the toppings, 2^N types of ramen can be ordered.

Akaki entered Takahashi-ya. She is thinking of ordering some bowls of ramen that satisfy both of the following two conditions:

* Do not order multiple bowls of ramen with the exactly same set of toppings.
* Each of the N kinds of toppings is on two or more bowls of ramen ordered.



You are given N and a prime number M. Find the number of the sets of bowls of ramen that satisfy these conditions, disregarding order, modulo M. Since she is in extreme hunger, ordering any number of bowls of ramen is fine.

Constraints

* 2 \leq N \leq 3000
* 10^8 \leq M \leq 10^9 + 9
* N is an integer.
* M is a prime number.

Input

Input is given from Standard Input in the following format:


N M


Output

Print the number of the sets of bowls of ramen that satisfy the conditions, disregarding order, modulo M.

Examples

Input

2 1000000007


Output

2


Input

3 1000000009


Output

118


Input

50 111111113


Output

1456748


Input

3000 123456791


Output

16369789

Solution (in CPP):

Solution:
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <queue>
#include <map>
using namespace std;
const int maxn=3e3+5;
typedef long long ll;
inline ll fpm(ll a,ll p,ll mod){
	ll res=1;
	while(p){
		if (p&1) res=res*a%mod;
		a=a*a%mod;
		p/=2;
	}
	return res;
}
int n;
int mod;
ll s[maxn][maxn];
ll fac[maxn],inv[maxn];
inline void init(){
	s[0][0]=1;
	for (int i=1;i<=n+1;i++){
		for (int j=1;j<=i;j++){
			s[i][j]=(s[i-1][j-1]+1ll*j*s[i-1][j]%mod)%mod;
		}
	}
	fac[0]=inv[0]=1;
	for (int i=1;i<=n;i++){
		fac[i]=1ll*fac[i-1]*i%mod;
		inv[i]=fpm(fac[i],mod-2,mod);
		//cout<<"print fac inv "<<fac[i]<<' '<<inv[i]<<endl;
	}
}
inline ll getc(ll n,ll m){
	if (n<m||n<0||m<0) return 0;
	return fac[n]*inv[m]%mod*inv[n-m]%mod;
}
inline ll f(int i){
	ll cur=1,a=fpm(2,n-i,mod);
	ll b=fpm(2,fpm(2,n-i,mod-1),mod);
	ll ans=0;
	for (int k=0;k<=i;k++){
		ans=(ans+s[i+1][k+1]*cur%mod*b%mod)%mod;
		cur=cur*a%mod;
	}
	return ans;
}
int main(){
	cin>>n>>mod;
	init();
	ll ans=0;
	for (int i=0;i<=n;i++){
		ll tp=getc(n,i)*f(i)%mod;
		//cout<<"print getc f "<<getc(n,i)<<' '<<f(i)<<endl;
		if (i%2==0) ans=(ans+tp)%mod;
		else ans=(ans+mod-tp)%mod;
	}
	cout<<ans<<endl;
}