CF#450 D.Unusual Sequences

题目：http://codeforces.com/contest/900/problem/D
题意：
给两个数X，Y，求有多少种可能，使得一串数a1，a2…ak，
有gcd（a1，a2…ak）= X且sum（ai）=Y
其中k任意
分析：
由题意知，Y%X!=0时，无解；
Y%X==0时，令n=Y/X；
则问题转化为：
求有多少种可能，使得一串数a1，a2…ak，
有①gcd（a1，a2…ak）= 1 且 ②sum（ai）= n
设f[i]为n=i时的答案；
若只考虑条件②：这是高中的排列组合放隔板问题，n个1，n-1个空，放入隔板，有2^(n-1)种方法；
那么有条件①时，容易知道2^(n-1)=∑f[i] (i为n的约数，包括n和1)
所以此题是快速幂+DP
代码：

#include <cstdio>
#include <algorithm>
#include <cstring>
#include <map>
#include <queue>
#include <string>
using namespace std;
typedef long long ll;
const ll MOD=1e9+7;
int x,y;
map<int,ll> MM;
ll quickpow(ll base,ll k)
{
    ll ret=1;
    while(k>0)
    {
        if((k&1)==1)
            ret=(ret*base)%MOD;
        base=(base*base)%MOD;
        k>>=1;
    }
    return ret;
}
ll dp(int n)
{
    if(MM[n]>0) return MM[n];
    int i;
    MM[n]=quickpow(2,n-1);
    for(i=1;i*i<=n;i++)
        if(n%i==0)
        {
            MM[n]=((MM[n]-dp(i))%MOD+MOD)%MOD;
            if(i*i<n&&i!=1) MM[n]=((MM[n]-dp(n/i))%MOD+MOD)%MOD;
        }
    return MM[n];
}
int main()
{
    int n;
    scanf("%d%d",&x,&y);
    if(y%x>0)
    {
        printf("0");
        return 0;
    }
    n=y/x;
    MM[1]=1;
    printf("%I64d",dp(n)%MOD);
    return 0;
}