#include <bits/stdc++.h>
using namespace std;
const int MOD = 1000000007;
long long po(long long a, long long b) {
  return !b ? 1 : po(a * a % MOD, b / 2) * (b & 1 ? a : 1) % MOD;
}
int mu(int n) {
  int cnt = 0;
  for (int i = 2; i * i <= n; i++) {
    if (n % i == 0) {
      cnt++;
      n /= i;
      if (n % i == 0) return 0;
    }
  }
  if (n > 1) cnt++;
  if (cnt & 1)
    return -1;
  else
    return 1;
}
int main() {
  ios::sync_with_stdio(false);
  cin.tie(0);
  int x, y;
  cin >> x >> y;
  if (y % x) {
    cout << 0 << endl;
    return 0;
  }
  y /= x;
  vector<int> divs;
  int i = 1;
  for (; i * i < y; i++) {
    if (y % i == 0) {
      divs.push_back(i);
      divs.push_back(y / i);
    }
  }
  if (i * i == y) divs.push_back(i);
  long long ans = 0;
  for (int div : divs) {
    ans = (ans + (mu(y / div) * po(2, div - 1) + MOD) % MOD) % MOD;
  }
  cout << ans << endl;
  return 0;
}
