【牛客】1005 Forsaken喜欢数论

题目

题目描述

        Forsaken有一个有趣的数论函数。对于任意一个数xxx,f(x)f(x)f(x)会返回xxx的最小质因子。如果这个数没有最小质因子,那么就返回0。

        现在给定任意一个nnn,Forsaken想知道∑i=1nf(i)\sum_{i = 1}^{n}{f(i)}∑i=1n​f(i)的值。

        

输入描述:

一个整数nnn。

输出描述:

一个整数代表上面的求和式的值。

示例1

输入

4

输出

7

备注:

1≤n≤3e7

题解

思路:本题本质上就是个欧拉筛,在筛数的时候加上最小质因数的值就好了。

ACcode

#include<iostream>
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll N=3e7+5;
ll st[N],prim[N];
ll sai(ll n)//欧拉筛
{
    ll cnt=1,ans=0;
    st[1]=1;
    for(ll i=2;i<=n;i++)
    {
        if(!st[i])
        {
            prim[cnt++]=i;
            ans+=i;
        }
        for(ll j=1;prim[j]*i<=n;j++)
        {
            st[i*prim[j]]=1;
            ans+=prim[j];
            if(i%prim[j]==0) break;
        }
    }
    return ans;
}
int main()
{
    
    ll x;cin>>x;
    cout<<sai(x)<<endl;
    return 0;
}

根据提供的引用内容,可以看出Forsaken喜欢数论,并且在数论方面有一定的研究。以下是一些与数论相关的Forsaken的例子: 1. 线性筛法求质数 ```python def get_primes(n): is_prime = [True] * (n + 1) primes = [] for i in range(2, n + 1): if is_prime[i]: primes.append(i) for j in range(len(primes)): if i * primes[j] > n: break is_prime[i * primes[j]] = False if i % primes[j] == 0: break return primes ``` 2. 分组算法 ```python def group_students(students, k): groups = [[] for _ in range(k)] for i in range(len(students)): groups[i % k].append(students[i]) return groups ``` 3. 正方形计数 ```python def count_squares(n): count = 0 for i in range(1, int(n ** 0.5) + 1): count += (n // i - i + 1) return count ``` 4. 三维数点 ```python def count_3d_points(points): bit = [0] * (max(points, key=lambda x: max(x)) + 1) count = 0 for i in range(len(points)): x, y, z = points[i] count += i - sum(bit[j] for j in range(x)) for j in range(x, len(bit)): bit[j] += 1 return count ``` 5. 独一无二的最小生成树 ```python def unique_minimum_spanning_tree(n, edges): edges.sort(key=lambda x: x[2]) parent = list(range(n)) rank = [0] * n def find(x): if parent[x] != x: parent[x] = find(parent[x]) return parent[x] def union(x, y): px, py = find(x), find(y) if px == py: return False if rank[px] < rank[py]: parent[px] = py elif rank[px] > rank[py]: parent[py] = px else: parent[px] = py rank[py] += 1 return True mst = [] for u, v, w in edges: if union(u, v): mst.append((u, v, w)) return mst --相关问题--:
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值