[BZOJ 4916]神犇和蒟蒻

Description

很久很久以前，有一只神犇叫yzy;

很久很久之后，有一只蒟蒻叫lty;

Input

请你读入一个整数N;1<=N<=1E9,A、B模1E9+7;

Output

请你输出一个整数 $A=\sum_{i=1}^N{\mu (i^2)}$ ;

请你输出一个整数 $B=\sum_{i=1}^N{\varphi (i^2)}$ ;

Sample Input

1

Sample Output

1
1

题解

首先注意到 $A$ 直接输出 $1$ 得满分。因为只有 $\mu(1^2)=1$ 。至于为什么，想想莫比乌斯函数是什么。

对于第二问容易发现 $\varphi(i^2)=i\cdot\varphi(i)$ 。至于为什么，想想你是怎么线性筛欧拉函数的。

对于函数 $f(i)=i\cdot\varphi(i)$ 容易发现，这也是个积性函数。我们可以将阈值内的线性筛筛出来。对于阈值外的，考虑杜教筛。

求 $S(n)=\sum\limits_{i=1}^nf(i)$

上述式子 $$g(1)S(n)=\sum_{i=1}^n(g*f)(i)-\sum_{i=2}^ng(i)S\left(\left\lfloor\frac{n}{i}\right\rfloor\right)$$

考虑到 $\sum\limits_{d\mid n}\varphi(d)=n$ ，又由于 $(g*f)(n)=\sum\limits_{d\mid n}\varphi(d)d\cdot g\left(\frac{n}{d}\right)$ 。我们考虑让 $g(n)=id(n)$ ，那么 $(id*f)(n)=\sum\limits_{d\mid n}n\cdot\varphi(d)=n^2$ 。由于 $\sum\limits_{i=1}^ni^2=\frac{n(n+1)(2n+1)}{6}$ 。显然这个卷积的前缀为 $\sum\limits_{i=1}^n(g*f)(i)=\frac{n(n+1)(2n+1)}{6}$ 。

故对于 $f$ $$S(n)=\frac{n(n+1)(2n+1)}{6}-\sum_{i=2}^ni\cdot S\left(\left\lfloor\frac{n}{i}\right\rfloor\right)$$

 //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('\n'))
 #define lowbit(x) ((x)&(-(x)))
 using namespace std;
 const int MOD = 1e9+;
 const int N = ;
 void read(int &x) {
     char ch; bool flag = ;
     for (ch = getchar(); !isdigit(ch) && ((flag |= (ch == '-')) || ); ch = getchar());
     for (x = ; isdigit(ch); x = (x<<)+(x<<)+ch-, ch = getchar());
     x *= -*flag;
 }
 void write(int x) {
     if (x > ) write(x/);
     putchar(x%+);
 }

 int n, inv2, inv6, f[N+];
 int prime[N+], isprime[N+], tot;
 map<int, int>mp;

 int quick_pow(int a, int b) {
     int ans = ;
     while (b) {
     if (b&) ans = 1ll*ans*a%MOD;
     a = 1ll*a*a%MOD, b >>= ;
     }
     return ans;
 }
 void get_f() {
     memset(isprime, , sizeof(isprime)); isprime[] = , f[] = ;
     for (int i = ; i <= N; i++) {
     if (isprime[i]) prime[++tot] = i, f[i] = 1ll*i*(i-)%MOD;
     for (int j = ; j <= tot && i*prime[j] <= N; j++) {
         isprime[i*prime[j]] = ;
         if (i%prime[j]) f[i*prime[j]] = 1ll*f[i]*(prime[j]-)%MOD*prime[j]%MOD;
         else {f[i*prime[j]] = 1ll*f[i]*prime[j]%MOD*prime[j]%MOD; break; }
     }
     (f[i] += f[i-]) %= MOD;
     }
 }
 int cal(int x) {
     if (x <= N) return f[x];
     if (mp.count(x)) return mp[x];
     int ans = 1ll*x*(x+)%MOD*(*x+)%MOD*inv6%MOD;
     for (int i = , last; i <= x; i = last+) {
     last = x/(x/i); (ans -= 1ll*(last+i)*(last-i+)%MOD*inv2%MOD*cal(x/i)%MOD) %= MOD;
     }
     return mp[x] = ans;
 }
 void work() {
     get_f(); inv2 = quick_pow(, MOD-), inv6 = quick_pow(, MOD-); read(n);
     writeln(); writeln((cal(n)+MOD)%MOD);
 }
 int main() {
     work();
     return ;
 }