#include <bits/stdc++.h>
using namespace std;
const int mod = 1e9 + 7;
void add(int &a, int b) {
  a += b;
  if (a >= mod) a -= mod;
  if (a < 0) a += mod;
}
int mul(int a, int b) { return (a * (long long)b) % mod; }
int pw(int x, int n) {
  int res = 1;
  while (n > 0) {
    if (n % 2 == 1) {
      res = (res * (long long)x) % mod;
      n--;
    }
    x = (x * (long long)x) % mod;
    n /= 2;
  }
  return res;
}
const int N = 2304;
int cnt[N];
int res[N];
int main() {
  ios::sync_with_stdio(0);
  cin.tie(0);
  int x, y;
  cin >> x >> y;
  if (y % x) {
    cout << 0 << '\n';
    return 0;
  }
  y /= x;
  int ans = 0;
  add(ans, pw(2, y - 1));
  vector<int> d;
  for (int i = 1; i * i <= y; i++) {
    if (y % i == 0) {
      if (i != 1) {
        d.push_back(i);
      }
      if (y / i != i) {
        d.push_back(y / i);
      }
    }
  }
  sort(d.begin(), d.end());
  for (int i = 0; i < d.size(); i++) {
    cnt[i] = pw(2, y / d[i] - 1);
  }
  for (int i = (int)d.size() - 1; i >= 0; i--) {
    res[i] = cnt[i];
    for (int j = i + 1; j < d.size(); j++) {
      if (d[j] % d[i] == 0) {
        add(res[i], -res[j]);
      }
    }
  }
  for (int i = 0; i < d.size(); i++) {
    add(ans, -res[i]);
  }
  cout << ans << '\n';
  return 0;
}
