标签：prod const int ll ans 例题 欧拉 left

B - Function

\(f(n)=\left\{\begin{array}{ll} 1, & n \in\{1\} \bigcup \text { Prime } \\ p f\left(p^{k-2}\right), & n=p^{k}(p \in \text { Prime; } k>1) \\ f\left(p_{1}^{e 1}\right) \prod_{i=2}^{r} p_{i}^{e_{i}}, & n=\prod_{i=1}^{r} p_{i}{ }^{e_{i}}\left(p_{1}<p_{2}<\cdots<p_{r} ; p_{i} \in \text { Prime; } r \geq 2\right) \end{array}\right.\)​

求\(\sum_{i=1}^{n} f(i)\)​

假设一个合数\(n\)，代入第三个式子，可得\(f(n)=p_1f(p_1^{e_1-2})\prod_{i=1}^{r} p_{i}{ }^{e_{i}}\)，再将\(p_1^{e_1-2}\)代入第二个式子，可得\(f(n)=p_1^2f(p_1^{e_1-4})\prod_{i=1}^{r} p_{i}{ }^{e_{i}}\)，以此类推，最后可推出\(f(n)=p_1^{\frac{e_1}{2}}\prod_{i=1}^{r} p_{i}{ }^{e_{i}}\).

Code:

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
#define gcd __gcd
#define endl '\n'
const int mod = 1e9 + 7;
const int inf = 0x3f3f3f3f;
const int maxn = 1e7 + 5;
int prime[maxn], vis[maxn] = {0}, cnt = 0;
ll power(ll a, ll b)
{
    ll ans = 1;
    while (b)
    {
        if (b & 1)
        {
            ans = ans * a % mod;
        }
        a = a * a % mod;
        b >>= 1;
    }
    return ans;
}
ll Euler(int n)
{
    ll res = 1;
    for (int i = 2; i <= n; i++)
    {
        if (!vis[i])
        {
            prime[cnt++] = i;
            res++;
        }
        for (int j = 0; prime[j] <= n / i; j++)
        {
            int tmp = prime[j] * i;
            vis[tmp] = 1;
            int count = 0;
            while (tmp % prime[j] == 0)
            {
                count++;
                tmp /= prime[j];
            }
            res += power(prime[j], count / 2) * tmp;
            if (i % prime[j] == 0)
            {
                break;
            }
        }
    }
    return res;
}

signed main()
{
    ll n;
    cin >> n;
    printf("%lld\n", Euler(n));
    return 0;
}

标签：prod,const,int,ll,ans,例题,欧拉,left
来源： https://www.cnblogs.com/rivego/p/15182467.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。