牛客练习赛53 B题调和级数

https://ac.nowcoder.com/acm/contest/1114/B

这题时间卡的比较死，多了一个快速幂的logn就过不了这题。

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int p=1e9+;
const int N=3e6+;
ll power[N];
int main()
{
    int n;
    cin>>n;
    ll ans=;
    for(int i=;i<=N-;i++)
    {
        power[i]=i;
    }
    for(int j=;j<=n;j++)
    {
        int i;
        for( i=j;i+j-<=n;i=i+j)
        {
            ll x=i/j;
            ans=(ans+1ll*(i+j+i-)*(j)/%p*power[x]%p)%p;
            power[x]=(power[x]*x)%p;
        }
        if(i<=n)
        {
            ll x=i/j;
            ans=(ans+(1ll)*(n+i)*(n-i+)/%p*power[n/i]%p)%p;
            power[x]=power[x]*x%p;
        }
    }
    cout<<ans<<"\n";
    return ;
}