【反演复习计划】【COGS2433】&&【bzoj3930，CQOI2015选数】爱蜜莉雅的冰魔法

同bzoj3930。

（日常盗题图）

#include<bits/stdc++.h>
#define N 1000010
#define yql 1000000007
#define ll long long
using namespace std;
int n;
ll m,l,r,k;ll f[N];
ll Pow(ll x,int y){
    ll ans=;
    while(y){if (y&) ans=ans*x%yql;  x=x*x%yql;  y>>=;  }
    return ans;
}
ll read(){
    ll x=;int f=;char ch;
    do{ch=getchar();if(ch=='-')f=-;}while(ch<''||ch>'');
    do{x=x*+ch-'';ch=getchar();}while(ch>=''&&ch<='');
    return x*f;
}
int main(){
    freopen("aimiliyadeicemagic.in","r",stdin);
    freopen("aimiliyadeicemagic.out","w",stdout);
    scanf("%d",&n);k=read();l=read();r=read();
    for (ll i=r-l;i>=;i--){
        ll lx=(l-)/(k*i),rx=r/(k*i);
        f[i]=(Pow(rx-lx,n)-(rx-lx)+yql)%yql;
        for (int j=;i*j<=r-l;j++)  f[i]=(f[i]-f[i*j]+yql)%yql;
    }
    if (l<=k&&k<=r) f[]++;
    printf("%d\n",f[]);
}