#include <stdio.h>

int M, K;
struct mat{
	long long a[55][55];
	mat operator *(mat t){
		mat r;
		for (int i = 0; i < M; i++) for (int j = 0; j < M; j++){
			r.a[i][j] = 0;
			for (int k = 0; k < M; k++) r.a[i][j] = (r.a[i][j] + a[i][k] * t.a[k][j]) % 1000000007;
		}
		return r;
	}
}A,I;

int main()
{
	scanf ("%d %d", &M, &K);
	for (int i = 0; i < M; i++){
		for (int j = 0; j < M; j++){
			A.a[i][(i + j) % M]++;
			A.a[i][(i * j) % M]++;
		}
		I.a[i][i] = 1;
	}

	while (K){
		if (K & 1) I = I * A;
		A = A * A;
		K /= 2;
	}

	printf ("%lld\n", I.a[0][0]);
	return 0;
}