codeforces 822 My pretty girl Noora(递推)

仔细想想就能想出来，f（i）是可以根据因子转移的，f(i)=f(divisor)*i/divisor+f(i/divisor);

对于i是素数的情况，也很明显答案就是c(i,2)。

从2开始递推到r就好了。

素数删可以求出最小素因子，直接用最小素因子做。

代码：

#include <bits/stdc++.h>
#define LL long long
using namespace std;
const LL mod=1e9+7;
const int maxn=5e6+5;
LL dp[maxn];
LL prim[maxn];
LL sdiv[maxn];
LL qmod[maxn];
bool pri[maxn];
int k;


void init()
{
    pri[1]=1;
    int i, j;
    for(i=2; i<maxn-2; i++)sdiv[i]=100000000;
    for(i=2; i<maxn-2; i++)
    {
        if(pri[i]==0)
        {
            prim[k++]=i;
        }
        for(j=0; j<k && prim[j]*i<maxn-2; j++)
        {
            pri[prim[j]*i]=1;
            sdiv[prim[j]*i]=min(sdiv[prim[j]*i],prim[j]);
        }
    } 
    return;
}



int main()
{
    init();
    memset(dp, -1, sizeof dp);
    LL t, l, r;
    int i, j;
    cin>>t>>l>>r;
    for(i=0; i<k; i++)
    {
        dp[prim[i]]=prim[i]*(prim[i]-1)/2;
        dp[prim[i]]%=mod;
    }
    LL x=1;
    for(i=0; i<=r-l; i++)
    {
        qmod[i]=x;   
        x=(x*t)%mod;
    }
    dp[2]=1;
    for(i=3; i<=r; i++)
    {
        if(pri[i])dp[i]=dp[sdiv[i]]*(i/sdiv[i])+dp[i/sdiv[i]];
    }
    LL res=0;
    for(i=0; i<=r-l; i++)
    {
        res=(res+(qmod[i]*dp[i+l])%mod);
        if(res>=mod)res-=mod;
//        printf("%lld %lld %lld\n", i+l, res, f(i+l));
    }
    cout<<res<<endl;
}