#422(div2)D. My pretty girl Noora

题意：N个人比赛，将N个人平均分为N/X组，每组对f[N]的贡献是x*(x-1)/2;N/X个晋级，然后继续操作，贡献为f[N/X].t^0*f[l]+t^1*f[l+1]......+t^(r-l)*f[r]最小

思路：我们改变X使其f[i]发生改变，当i为质数只能为1组，即贡献为i*(i-1)/2。那么合数的贡献最小为什么，易发现当每组为他的最小质因子时，最优，记得开long long ，还有看见t^0,t^1，下意识就用了快速幂，T到死，还百度快速求素数，然而我SB了

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll mod=1e9+7;
const int N=5*1e6+10;

ll a[N];
ll b[N];


int main(){
    ll t,l,r;
    cin>>t>>l>>r;
    for(ll i=2;i<=r+1;i++){
        if(b[i]) {continue;}
        for(ll j=i*i;j<=r+1;j+=i){
            if(b[j]==0)
             b[j]=i;
        }
    }
    a[2]=1;a[3]=3;
    for(ll i=4;i<=r;i++){
        if(b[i]){
            ll x=b[i];
            a[i]=(a[x]*(i/x)+a[i/x]+mod)%mod;
        }
        else a[i]=(i*(i-1)/2+mod)%mod;
    }
    ll xx=1;
    ll sum=0;
    for(ll i=l;i<=r;i++){
        sum=(sum+xx*a[i]+mod)%mod;
        xx=(xx*t+mod)%mod;
    }
    cout<<sum<<endl;
}