#include <bits/stdc++.h>
using namespace std;
#define MOD 1000000007

long long o;

struct X
{
   long long p[70][70];
} pole[100], out;

bool soudelnost[70][70];

long long in, inn;
long long m;

void ini()
{
    for(long long i = 1; i < 67; i++)
    {
        for(long long j = 0; j < 67; j++)
        {
            long long  k = 1;
            while(k < 70)
            {
                k++;
                if(i%k==0 and j%k==0) soudelnost[i][j]=1;
            }
        }
    }
}


X sloz(X a , X b)
{
    X ret;
    for(long long i = 0; i<67; i++)
        for(long long j=0; j<67; j++)
         ret.p[i][j]=0;
    for(long long i = 1; i < m+1; i++)
    {
        for(long long j=1; j<m+1; j++)
        {
            for(long long k = 1; k <m+1; k++)
            {
                for(long long l = 1; l < m+1; l++)
                {
                    if(soudelnost[j][k]==0)
                    {
                        ret.p[i][l]+=a.p[i][j]*b.p[k][l];
                        ret.p[i][l] %= MOD;
                    }
                }
            }
        }
    }
    return ret;
}


int main()
{
    //cout << 70*70*70*70*70 << endl;
    ini();
    cin >> m >> in;
    for(long long i = 1; i<=m; i++ )
    {
        pole[0].p[i][i]=1;
    }

    for(long long i = 1; (1ll<<i) <= in; i++)
    {
    //long long pom = (1ll<<i);
    //cout << pom << " " << i<< endl;
        pole[i] = sloz(pole[i-1], pole[i-1]);





    }

    inn = in;

    out.p[1][1]=1;

        for(int i = 1; in; i++)
    {
        if (in%(1ll<<i) )
        {
            out = sloz(pole[i-1], out);

            in -= (1ll<<(i-1));
        }
    }

    for(int i = 1; i < m+1; i++)
    {
        for(int j = 0; j < m+1; j++)
        {
            o += out.p[i][j];
            o %= MOD;
        }
    }

    cout << o << endl;








    return 0;
}