[BZOJ 3944]Sum

Description

求 $ans1=sumlimits_{i=1}^n varphi(i)$ 和 $ans2=sumlimits_{i=1}^n mu(i)$ 。 $n<2^{31}$ ，多组询问。

Input

一共T+1行

第1行为数据组数T(T<=10)

第2~T+1行每行一个非负整数N，代表一组询问

Output

一共T行，每行两个用空格分隔的数ans1,ans2

Sample Input

6
1
2
8
13
30
2333

Sample Output

1 1
2 0
22 -2
58 -3
278 -3
1655470 2

题解

杜教筛的板子。已经有许多博客有讲，我也不凑这个热闹了。

记 $S_n$ 为积性函数的前缀。

对于 $mu$ ： $$S(n)=1-sum_{i=2}^nSleft(leftlfloorfrac{n}{i} ight floor ight)$$

对于 $varphi$ ： $$S(n)=frac{n(n+1)}{2}-sum_{i=2}^nSleft(leftlfloorfrac{n}{i} ight floor ight)$$

取 $2333333$ 作为阈值，线性筛出之前的，然后记忆化之后的。

//It is made by Awson on 2018.1.24
#include <set>
#include <map>
#include <cmath>
#include <ctime>
#include <queue>
#include <stack>
#include <cstdio>
#include <string>
#include <vector>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
#define LL long long
#define Abs(a) ((a) < 0 ? (-(a)) : (a))
#define Max(a, b) ((a) > (b) ? (a) : (b))
#define Min(a, b) ((a) < (b) ? (a) : (b))
#define Swap(a, b) ((a) ^= (b), (b) ^= (a), (a) ^= (b))
#define writeln(x) (write(x), putchar('
'))
#define lowbit(x) ((x)&(-(x)))
using namespace std;
const int N = 2333333;
void read(LL &x) {
    char ch; bool flag = 0;
    for (ch = getchar(); !isdigit(ch) && ((flag |= (ch == '-')) || 1); ch = getchar());
    for (x = 0; isdigit(ch); x = (x<<1)+(x<<3)+ch-48, ch = getchar());
    x *= 1-2*flag;
}
void write(LL x) {
    if (x > 9) write(x/10);
    putchar(x%10+48);
}

LL mu[N+5], phi[N+5], n;
struct num {
    LL ans1, ans2;
    num() {}
    num(LL _ans1, LL _ans2) {ans1 = _ans1, ans2 = _ans2; }
}ans;
map<LL, num>mp;

int prime[N+5], isprime[N+5], tot;
void get_pre() {
    memset(isprime, 1, sizeof(isprime)); isprime[1] = 0, mu[1] = phi[1] = 1;
    for (int i = 2; i <= N; i++) {
    if (isprime[i]) prime[++tot] = i, mu[i] = -1, phi[i] = i-1;
    for (int j = 1; j <= tot && i*prime[j] <= N; j++) {
        isprime[i*prime[j]] = 0;
        if (i%prime[j]) mu[i*prime[j]] = -mu[i], phi[i*prime[j]] = phi[i]*(prime[j]-1);
        else {mu[i*prime[j]] = 0, phi[i*prime[j]] = phi[i]*prime[j]; break; }
    }
    mu[i] += mu[i-1], phi[i] += phi[i-1];
    }
}
num Less (const num &a, const num &b) {num ans; ans.ans1 = a.ans1 - b.ans1, ans.ans2 = a.ans2-b.ans2; return ans; }
num Times (const num &a, const LL &x) {num ans; ans.ans1 = a.ans1*x , ans.ans2 = a.ans2*x; return ans; }
num cal(LL x) {
    if (x <= N) return num(phi[x], mu[x]);
    if (mp.count(x)) return mp[x];
    num ans = num(x*(x+1)/2, 1);
    for (LL i = 2, last; i <= x; i = last+1) {
    last = x/(x/i); ans = Less(ans, Times(cal(x/i), (last-i+1)));
    }
    return mp[x] = ans;
}
void work() {
    read(n); ans = cal(n);
    write(ans.ans1), putchar(' ');
    if (ans.ans2 < 0) putchar('-'), writeln(-ans.ans2);
    else writeln(ans.ans2);
}
int main() {
    get_pre(); LL t; read(t);
    while (t--) work();
    return 0;
}