Forsaken喜欢数论

题目地址：

https://ac.nowcoder.com/acm/problem/53079

基本思路：

其实就是一个素数筛，我们想一下在埃氏筛法的过程中，我们是从小到大的对于每个质数，将它的倍数都筛掉，那么这个质数其实就是这些被筛掉的非质数的质因子，并且由于我们是从小到大遍历的，所以很明显先出现的就是最小质因子，因此每次筛去非质数的时候累加这个质数就是了。

参考代码：

#pragma GCC optimize(2)
#pragma GCC optimize(3)
#include <bits/stdc++.h>
using namespace std;
#define IO std::ios::sync_with_stdio(false)
#define int long long
#define rep(i, l, r) for (int i = l; i <= r; i++)
#define per(i, l, r) for (int i = l; i >= r; i--)
#define mset(s, _) memset(s, _, sizeof(s))
#define pb push_back
#define pii pair <int, int>
#define mp(a, b) make_pair(a, b)
#define INF (int)1e18

inline int read() {
  int x = 0, neg = 1; char op = getchar();
  while (!isdigit(op)) { if (op == '-') neg = -1; op = getchar(); }
  while (isdigit(op)) { x = 10 * x + op - '0'; op = getchar(); }
  return neg * x;
}
inline void print(int x) {
  if (x < 0) { putchar('-'); x = -x; }
  if (x >= 10) print(x / 10);
  putchar(x % 10 + '0');
}

const int maxn = 3e7+ 10;
int n;
bool is_prime[maxn];
signed main() {
  IO;
  cin >> n;
  mset(is_prime,true);
  int ans = 0;
  for(int i = 2 ; i <= n ; i++){
    if(is_prime[i]){
      ans += i; //自己本身;
      for(int j = 2 * i ; j <= n ; j+= i) {
        if(is_prime[j]) is_prime[j] = false,ans += i; 
      }
    }
  }
  cout << ans << '\n';
  return 0;
}