注意事项：无。

代码：

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
bool notpr[MAXn + 10]; int cntpr, pr[MAXn + 10];
int mob[MAXn + 10];
void evamob(int n) {
    notpr[1] = 1;
    mob[1] = 1;
    for (int i = 2; i <= n; ++i) {
        if (notpr[i] == 0) {
            pr[++cntpr] = i;
            mob[i] = -1;
        }
        for (int j = 1, topprj = n / i; j <= cntpr && pr[j] <= topprj; ++j) {
            notpr[i * pr[j]] = 1;
            if (i % pr[j]) {
                mob[i * pr[j]] = -mob[i];
            } else {
                mob[i * pr[j]] = 0;
                break;
            }
        }
    }
}