#include <bits/stdc++.h>
using namespace std;
const long long MAX = 1e6 + 8;
const long long MOD = 1e9 + 7;
const long long sz = 1e5 + 5;
void f_io() {
  ios::sync_with_stdio(0);
  cin.tie(0);
  cout.tie(0);
  ;
}
long long power(long long a, long long b) {
  if (!b) return 1;
  if (b == 1) return a;
  if (b & 1) return (power(a, b - 1) * a) % MOD;
  long long q = power(a, b / 2);
  return (q * q) % MOD;
}
map<long long, long long> F;
long long f(long long N) {
  if (N == 1) return 1;
  if (F.find(N) != F.end()) return F[N];
  long long ret = power(2, N - 1) - 1;
  for (long long i = 2; i * i <= N; i++) {
    if (!(N % i)) {
      ret -= f(i);
      ret += MOD;
      ret %= MOD;
      if (i * i != N) {
        ret -= f(N / i);
        ret += MOD;
        ret %= MOD;
      }
    }
  }
  return (F[N] = (1ll * ret + MOD) % MOD);
}
int32_t main() {
  f_io();
  long long tc = 1;
  while (tc--) {
    long long n, m;
    cin >> n >> m;
    if (m % n) return 0 * printf("0\n");
    cout << f(m / n) << "\n";
  }
}
