BZOJ 3930 容斥原理

思路：

移至iwtwiioi    http://www.cnblogs.com/iwtwiioi/p/4986316.html

//By SiriusRen
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define int long long
const int mod=;
int n,k,l,h,f[];
int pow(int x,int y){
    int res=;
    while(y){
        if(y&)res=res*x%mod;
        x=x*x%mod,y>>=;
    }return res;
}
signed main(){
    scanf("%lld%lld%lld%lld",&n,&k,&l,&h);
    for(int i=h-l;i;i--){
        int L=(l-)/(k*i),R=h/(k*i);
        f[i]+=(pow(R-L,n)-R+L+mod)%mod;
        for(int j=;i*j<=h-l;j++)(f[i]-=f[i*j])%=mod;
        f[i]=(f[i]+mod)%mod;
    }printf("%lld\n",(f[]+(k>=l&&k<=h))%mod);
}