#include <bits/stdc++.h>
/*#include <boost/multiprecision/cpp_dec_float.hpp>
#include <boost/multiprecision/cpp_int.hpp>
*/#pragma GCC target("avx")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")    
using namespace std;
/*namespace mp = boost::multiprecision;
using Bint = mp::cpp_int;
using Real = mp::number<mp::cpp_dec_float<1024>>;
*/#define lli long long int
#define uli unsigned long long int
#define INF 999999999999
#define rep(i,m,n) for(lli i = m;i < n;i++)
#define rrep(i,m,n) for(lli i=m-1;i>=n;i--)
#define pb(n) push_back(n)
#define UE(N) N.erase(unique(N.begin(),N.end()),N.end());
#define Sort(n) sort(n.begin(), n.end())
#define Rev(n) reverse(n.begin(),n.end())
#define Out(S) cout << S << endl
#define NeOut(S) cout << S
#define HpOut(S) cout << setprecision(50) << S << endl
#define Vec(K,L,N,S) vector<L> K(N,S)
#define DV(K,L,N,M,S) vector<vector<L>> K(N,vector<L>(M,S))
#define TV(K,L,N,M,R,S) vector<vector<vector<L>>> K(N,vector<vector<L>>(M,vector<L>(R,S)))
#define pint pair<lli,lli>
#define paf(L,R) pair<L,R>
#define mod 1000000007
#define MAX 10000000
#define ALL(a)  a.begin(),a.end()
#define chmax(a, b) a = (((a)<(b)) ? (b) : (a))
#define chmin(a, b) a = (((a)>(b)) ? (b) : (a))  
long long N,K,sum=0,num=0;
lli modpow(lli n,lli k,lli m){
  if(k==0)return 1;
  else if(k%2==1)return modpow(n,k-1,m)*n%m;
  else{
    lli temp=modpow(n,k/2,m);
    return temp*temp%m;
  }
}
void dfs(int n,vector<int>&P){
  if(n)rep(i,0,K){
    P[n-1]=i+1;
    dfs(n-1,P);
  }
  else{
    rep(p,0,(1<<N)-1){
      num=1;
      rep(q,0,N)if(p&(1<<(N-q-1)))num=num*P[q]%mod;
      sum=(sum+num)%mod;
    }
  }
}
int main(){
  ios::sync_with_stdio(false);
  cin.tie(0);
  cin >> N >> K;
  /*
  Vec(P,int,N,0);
  dfs(N,P);
  Out(sum);
  */
  K%=mod;
  Out(((modpow(K*(K+3)%mod,N,mod)-modpow(K*(K+1)%mod,N,mod))*modpow(modpow(2,mod-2,mod),N,mod)%mod+10LL*mod)%mod);
}