題目連結：http://codeforces.com/problemset/problem/900/D

數列的最大公因數為x則所有的ai均可以由x表示於是若有這樣的陣列則y%x==0。陣列可以分成多個x，於是可以把此題轉化成將y/x個x放進m個箱子裡（不允許有空箱子），於是這道題便轉化成了經典的組合數學問題方案數為C(y/x-1,m-1)把m從2一直疊加到y/x-1於是總共的方案數變化成了2^(y/x-1)-1。但是這樣子會存在非法的序列例如當x=3,y=12相當於4個3，在這裡若直接拿上述方法做會出現6,6這樣的方案數，此序列的最大公因數將不再是3而是6所以要減去這樣的方案數。可以發現產生這樣的情況是因為m為y/x的因子於是減去當m為y/x的因子的情況以即可。在這裡注意要剪枝不然會超時！

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <vector>
#include <map>
using namespace std;
typedef long long int LL;
const int mod=1e9+7;
const int maxn=100;
LL num[maxn]={};
int x,y,cnt=0;
map<int,int>vis;   //剪枝遇到相同的y/x則直接返回即可
LL num_pow(LL a,LL k){
    LL res=1;
    for(;k;k>>=1){
        if(k&1) res=(res*a)%mod;
        a=(a*a)%mod;
    }
    return res;
}
LL dfs(LL mid){
    if(mid==1) return 1;
    LL ans=(num_pow(2,mid-1)-1+mod)%mod;
    if(vis.count(mid)) return vis[mid];
    for(int i=0;i<cnt && num[i]*num[i]<=mid;i++)if(mid%num[i]==0){
        ans=(ans-dfs(mid/num[i])+mod)%mod;
        if(num[i]*num[i]!=mid) ans=(ans-dfs(num[i])+mod)%mod;
    }
    vis[mid]=ans;
    return ans;
}
int main()
{
    while(scanf("%d%d",&x,&y)==2){
        vis.clear();cnt=0;
        memset(num,0,sizeof(num));
        if(y%x!=0){
            printf("0\n");
            continue;
        }
        for(int i=2;i*i<=y/x;i++)if(y/x%i==0)//提取全部的因子
            num[cnt++]=i;
        printf("%lld\n",dfs(y/x));
    }
    return 0;
}