#include <bits/stdc++.h>
#include <math.h>
#define rep(i,n) for (int i = 0; i < (n); ++i)
#define PI 3.14159265359
using namespace std;
using ll = long long;
using P = pair<int, int>;

const int mod = 1000000007;
struct mint {
  ll x;
  mint(ll x=0):x(x%mod){}
  mint operator-() const { return mint(-x);}
  mint& operator+=(const mint a) {
    if ((x += a.x) >= mod) x -= mod;
    return *this;
  }
  mint& operator-=(const mint a) {
    if ((x += mod-a.x) >= mod) x -= mod;
    return *this;
  }
  mint& operator*=(const mint a) {
    (x *= a.x) %= mod;
    return *this;
  }
  mint operator+(const mint a) const {
    mint res(*this);
    return res+=a;
  }
  mint operator-(const mint a) const {
    mint res(*this);
    return res-=a;
  }
  mint operator*(const mint a) const {
    mint res(*this);
    return res*=a;
  }
  mint pow(ll t) const {
    if (!t) return 1;
    mint a = pow(t>>1);
    a *= a;
    if (t&1) a *= *this;
    return a;
  }
  
  mint inv() const {
    return pow(mod-2);
  }
  mint& operator/=(const mint a) {
    return (*this) *= a.inv();
  }
  mint operator/(const mint a) const {
    mint res(*this);
    return res/=a;
  }
};

mint dp[5000][5000];

int main() {
  int n, m;
  cin >> n >> m;

  mint cur = 2;
  rep(i, m) {
    dp[0][i] = cur;
    cur *= 2;
  }
  cur = 2;
  rep(i, n) {
    dp[i][0] = cur;
    cur *= 2;
  }

  for (int i=1; i<n; i++) {
    for (int j=1; j<m; j++) {
      dp[i][j] = dp[i][j-1] * 2 + (dp[i-1][j] - dp[i-1][j-1]) * (j + 2);
    }
  }

  cout << dp[n-1][m-1].x << endl;

  return 0;
}

/*
*/
