I 有毒的玻璃球（积性函数 + 快速幂）

最新推荐文章于 2020-08-17 20:32:32 发布

少侠，慢点走

最新推荐文章于 2020-08-17 20:32:32 发布

阅读量126

点赞数

分类专栏：比赛栏

本文链接：https://blog.csdn.net/gtgym321/article/details/90183805

版权

比赛栏专栏收录该内容

7 篇文章 0 订阅

订阅专栏

积性函数 + 快速幂

#include<bits/stdc++.h>
using namespace std;
typedef long long ll ;
double e = (sqrt(5) + 1) / 2;
const int mod = 1e9 + 7;
const int N = 1e7;
bool isprime[N + 7];
int prime[N + 7],ans, cnt, n, k, sum[N + 5];
 
ll qpow(ll s, ll k)
{
    ll ans = 1 ;
    while(k)
    {
        if(k & 1)
            ans = ans * s % mod;
        k >>= 1;
        s = s * s % mod;
    }
    return ans ;
}
 
void init()
{
    sum[1] = 1;
    for(int i  = 2 ; i <= n + 1; i ++)
    {
 
 
        if(!isprime[i])
            prime[++cnt] = i,sum[i] = qpow(i,k);
 
 
        for(int j = 1; j <= cnt && prime[j] * i <= n + 1; j++)
        {
            isprime[prime[j] * i] = 1;
            sum[prime[j] * i] = sum[i] * 1ll * sum[prime[j]] % mod;
            if(i % prime[j] == 0)
                break;
        }
    }
}
int main()
{
 
    cin >> n;
    cin >> k;
    init();
    for(int i = 1 ; i <= n ; i++)
    {
        ans = (ans + 1LL * (n / i) *sum[i] ) % mod;
 
    }
    cout << ans << endl;
    return 0;
}