Codeforces 1017F The Neutral Zone（数论 埃氏筛）

题目链接：The Neutral Zone

题意

注意本题内存限制为 $16MB$。
定义函数：

$$f(x)=Ax^3+Bx^2+Cx+D\\ exlog_f(1)=0\\ exlog_f(p_1^{a_1}p_2^{a_2}\cdots p_k^{a_k})=a_1f(p_1)+a_2f(p_2)+\cdots a_kf(p_k)$$ 设 $n=p_1^{a_1}p_2^{a_2}\cdots p_k^{a_k}$，其中的每个 $p$ 是 $n$ 的质因数， $a$ 值是质因数对应的指数。
给定整数 $n$ 和参数 $A,B,C,D$ 的值，计算以下公式：
$$\sum_{i=1}^nexlog_f(i)$$

输入

输入包含 $4$ 个整数 $n,A,B,C,D~(1\leq n\leq3\times10^8,0\leq A,B,C,D\leq10^6)$。

输出

将计算结果对 $2^{32}$ 取模后输出。

样例

输入
12 0 0 1 0
输出
63

输入
4 1 2 3 4
输出
136

题解

从 $1$ 到 $n$ 内每个质数 $p$ 对答案的贡献为 $\lfloor\frac{n}{p}\rfloor+\lfloor\frac{n}{p^2}\rfloor+\cdots\lfloor\frac{n}{p^k}\rfloor$，其中 $k$ 为满足 $p^k\leq n$ 的最大值，因此最后的答案就是：

$$\sum_{i=1}^nf(p_i)(\lfloor\frac{n}{p_i}\rfloor+\lfloor\frac{n}{p_i^2}\rfloor+\cdots\lfloor\frac{n}{p_i^k}\rfloor)$$ 对 $2^{32}$ 取模，就是用 unsigned int 进行计算并让它们自然溢出。
最后是内存限制，如果我们用一个线性筛，就需要存每个数字是否是质数，如果用 bool 数组， $3\times10^8$ 的 bool 数组需要 $286MB$，用 bitset 压位需要 $\frac{286}{8}\approx35.75MB$，如果用埃氏筛，就可以在筛的过程中跳过 $2$ 和 $3$ 的倍数，我们将删去 $2$ 和 $3$ 的倍数后连续的数字用对应的 $1,2,...,n$ 表示：

$5$	$7$	$11$	$13$	$17$	$19$	$23$	$25$	$29$	$31$	$35$	$\cdots$
$1$	$2$	$3$	$4$	$5$	$6$	$7$	$8$	$9$	$10$	$11$	$\cdots$

可以发现数字与下标之间是 $Index=\lfloor\frac{n}{3}\rfloor$ 的关系，因此如果跳过 $2$ 和 $3$ 的倍数，可以减少 $\frac{2}{3}$ 的空间，最后只需要 $\frac{35.75}{3}\approx11.9MB$。

过题代码

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <climits>
#include <cstring>
#include <string>
#include <vector>
#include <list>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <bitset>
#include <algorithm>
#include <functional>
#include <iomanip>
using namespace std;

#define LL long long
unsigned n, a, b, c, d;
unsigned ans;
bitset<100000001> bit;

unsigned f(unsigned x) {
    return a * x * x * x + b * x * x + c * x + d;
}

unsigned Get(unsigned x) {
    unsigned ret = f(x);
    unsigned tmp = 0;
    for(unsigned i = 1; i <= n / x; i *= x) {
        tmp += n / (i * x);
    }
    return ret * tmp;
}

int main() {
    #ifdef LOCAL
    freopen("test.txt", "r", stdin);
//    freopen("test1.out", "w", stdout);
    #endif // LOCAL
    ios::sync_with_stdio(false);

    while(cin >> n >> a >> b >> c >> d) {
        ans = Get(2) + Get(3);
        bit.reset();
        for(unsigned i = 5; i <= n; ++i) {
            if(i % 2 == 0 || i % 3 == 0) {
                continue;
            }
            if(bit[i / 3] == 0) {
                ans += Get(i);
                for(unsigned j = i; j <= n / i; ++j) {
                    if(i * j % 2 == 0 || i * j % 3 == 0) {
                        continue;
                    }
                    bit[i * j / 3] = 1;
                }
            }
        }
        cout << ans << endl;
    }

    return 0;
}