#include<bits/stdc++.h>
using namespace std;
using ll = long long;

#include<atcoder/modint>
using mint = atcoder::modint1000000007;
int n;

map<pair<int,int>,mint> memo;
mint calc(int k,int mask){
    if(memo.find(make_pair(k,mask))!=memo.end()) return memo[make_pair(k,mask)];
    mint ans = 0;
    if(mask==0){
        if(k==0) ans = n;
        else{
            ans = n * calc(k-1,1).pow(3);
        }
    }else if(mask==1){
        if(k==0) ans = n-1;
        else{
            ans += (n-1) * calc(k-1,2) * calc(k-1,2) * calc(k-1,1);;
            ans += calc(k-1,4) * calc(k-1,4) * calc(k-1,1);
        }
    }else if(mask==2){
        if(k==0) ans = n-2;
        else{
            ans += (n-2) * calc(k-1,3) * calc(k-1,2) * calc(k-1,2);
            ans += 2 * calc(k-1,2) * calc(k-1,4) * calc(k-1,5);
        }
    }else if(mask==3){
        if(k==0) ans = n - 3;
        else{
            ans += (n-3) * (calc(k-1,3).pow(3));
            ans += 3 * calc(k-1,3) * calc(k-1,5) * calc(k-1,5);
        }
    }else if(mask==4){
        if(k==0) ans = n-1;
        else{
            ans += (n-1) * calc(k-1,5) * calc(k-1,2) * calc(k-1,2);
            ans += calc(k-1,6) * calc(k-1,4) * calc(k-1,4);
        }
    }else if(mask==5){
        if(k==0) ans = n - 2;
        else{
            ans += (n-2) * calc(k-1,5) * calc(k-1,3) * calc(k-1,3);
            ans += calc(k-1,5) * calc(k-1,5) * calc(k-1,6);
            ans += calc(k-1,5) * calc(k-1,5) * calc(k-1,5);
        }
    }else if(mask==6){
        if(k==0) ans = n - 1;
        else{
            ans += calc(k-1,6).pow(3);
            ans += (n-1) * (calc(k-1,5).pow(3));
        }
    }
    return memo[make_pair(k,mask)] = ans;
}

int main(){
    cin.tie(nullptr);
    ios::sync_with_stdio(false); 
    int k;
    cin>>k>>n;
    cout<<calc(k,0).val()<<endl;
}