bzoj3930: [CQOI2015]选数（Dp）

题目传送门
。

解法：
跟这题思路蛮像的
。
就是作为公因数不一定作为最大公因数。
那么求作为公因数的方案很简单。只需知道有多少个数是他的倍数。然后用个数^N就行了。
f[i]表示最大公因数为i*K的方案。
那么先算出公因数的方案。
然后在减去f[i的倍数]就行了呀。

代码实现：

#include
#include
#include
#include
#include
#include
#include
using namespace std;
typedef long long ll;
ll n,k,l,r;ll f[110000];const ll mod=1000000007;
ll pow_mod(ll a,int b) {
    ll ans=1;a%=mod;
    while(b!=0) {if(b%2==1)ans=(ans*a)%mod;a=(a*a)%mod;b/=2;}
    return ans;
}
int main() {
    scanf("%lld%lld%lld%lld",&n,&k,&l,&r);
    ll L=l/k,R=r/k;if(l%k)L++;
    for(int i=R-L+1;i>=1;i--) {
        ll LL=L/i,RR=R/i;if(L%i)LL++;
        if(LL<=RR) {
            f[i]=pow_mod(RR-LL+1,n);f[i]=(f[i]-(RR-LL+1)+mod)%mod;
            for(int j=i*2;j<=R-L+1;j+=i) f[i]=(f[i]-f[j]+mod)%mod;
        }
    }if(L==1)f[1]++,f[1]%=mod;
    printf("%lld\n",f[1]);
    return 0;
}