BZOJ 3930: [CQOI2015]选数 莫比乌斯反演

https://www.lydsy.com/JudgeOnline/problem.php?id=3930

https://blog.csdn.net/ws_yzy/article/details/50966171

求从区间[L,H]中可重复的选出n个数使其gcd=k的方案数 

就是，莫比乌斯反演，我抄的代码所以没有提前求莫比乌斯函数。

自乘的倍数个方案要去掉。现在想想我最后自己想出来的代码好像是没问题的但是我忘了加上唯一的一个自乘特判情况了，wa了太多次最后没忍住读（抄）了一份ac代码，还是意志不够坚定也不够细心。

emmmm现在发现似乎好久没有写博客了呢。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
#define LL long long
const LL maxn=100010;
const LL p=(int)1e9+7;
LL n,k,l,r;
LL ans[maxn]={};
LL mpow(LL x,LL y){
    LL z=1;
    while(y){
        if(y&1)z=(z*x)%p;
        x=(x*x)%p;
        y>>=1;
    }
    return z;
}
int main(){
    scanf("%lld%lld%lld%lld",&n,&k,&l,&r);
    LL a=r/k,b=l/k;if(l%k)++b;
    for(LL i=r-l;i>=1;--i){
        LL ww=a/i-b/i; if(b%i==0)++ww;
        if(ww<=0)continue;
        ans[i]=(mpow(ww,n)-(ww%p)+p)%p;
        for(int j=2*i;j<=r-l;j+=i)ans[i]=(ans[i]-ans[j]+p)%p;
    }
    if(b==1)ans[1]=(ans[1]+1)%p;
    printf("%lld\n",ans[1]);
    return 0;
}