ZOJ The Sum of Unitary Totient（min

The Sum of Unitary Totient

积性函数，满足质数点是多项式，直接 min_25 了，由于单次求解，所以使用递归的 min_25 会较快。

#include <bits/stdc++.h>using namespace std;const int N = 1e5 + 10;int prime[N], a[N], id1[N], id2[N], cnt, n, T, m;long long sum1[N], sum2[N], g1[N], g2[N], s[N];bool st[N];inline int ID(int x) {return x <= T ? id1[x] : id2[n / x];
}void init() {T = sqrt(n + 0.5);m = 0;for (int l = 1, r; l <= n; l = r + 1) {r = n / (n / l);a[++m] = n / l;a[m] <= T ? id1[a[m]] = m : id2[n / a[m]] = m;g1[m] = 1ll * a[m] * (a[m] + 1) /  2 - 1;g2[m] = a[m] - 1;}for (int j = 1; j <= cnt && 1ll * prime[j] * prime[j] <= n; j++) {for (int i = 1; i <= m && 1ll * prime[j] * prime[j] <= a[i]; i++) {g1[i] -= 1ll * prime[j] * (g1[ID(a[i] / prime[j])] - sum1[j - 1]);g2[i] -= (g2[ID(a[i] / prime[j])] - sum2[j - 1]);}}/*下面是递推的 min_25。*/// for (int i = 1; i <= m; i++) {//   s[i] = g1[i] - g2[i];// }// int top = 1;// while (1ll * prime[top] * prime[top] <= n) {//   top++;// }// top--;// for (int j = top; j >= 1; j--) {//   for (int i = 1; i <= m && 1ll * prime[j] * prime[j] <= a[i]; i++) {//     for (int cur = prime[j]; 1ll * cur * prime[j] <= a[i]; cur *= prime[j]) {//       s[i] += 1ll * (cur - 1) * (s[ID(a[i] / cur)] - sum1[j] + sum2[j]) + cur * prime[j] - 1;//     }//   }// }
}long long solve(int n, int m) {if (n < prime[m]) {return 0;} long long res = g1[ID(n)] - g2[ID(n)] - sum1[m - 1] + sum2[m - 1];for(int j = m; j <= cnt && 1ll * prime[j] * prime[j] <= n; j++) {for(int i = prime[j]; 1ll * i * prime[j] <= n; i *= prime[j]) {res += (i - 1) * solve(n / i, j + 1) + i * prime[j] - 1;}} return res;
}int main() {// freopen("in.txt", "r", stdin);// freopen("out.txt", "w", stdout);for (int i = 2; i < N; i++) {if (!st[i]) {prime[++cnt] = i;sum1[cnt] = sum1[cnt - 1] + i;sum2[cnt] = sum2[cnt - 1] + 1;}for (int j = 1; j <= cnt && 1ll * i * prime[j] < N; j++) {st[i * prime[j]] = 1;if (i % prime[j] == 0) {break;}}}while (scanf("%d", &n) != EOF) {init();printf("%lld\n", solve(n, 1) + 1);}return 0;
}

ZOJ The Sum of Unitary Totient（min_25 筛）相关推荐

ZOJ The Sum of Unitary Totient (Min_25筛)
题意: 给你一个函数f(n)=(p1a1-1)(p2a2-1)...(prar-1) ,n=p1a1p2a2...prar 求\[\sum\limits_{i = 1}^n {f(i)} \],\[n ...
ZOJ The Sum of Unitary Totient
ZOJ 3808 %%%唐老师杜教筛博客课后练习题不过网上好像没题解啊,好像唐老师也是分段打表的- 但是可以Min_25筛显然这个ϕ∗\phi^*ϕ∗是积性函数, ϕ∗(p)=p−1\phi^*( ...
loj #6053 简单的函数 min_25筛
\(\color{#0066ff}{ 题目描述 }\) 某一天,你发现了一个神奇的函数\(f(x)\),它满足很多神奇的性质: \(f(1)=1\) \(f(p^c)=p \oplus c\) (\( ...
Min_25筛学习Tip+链接
前言机房里差不多都会Min_25筛了,我也赶紧补一波坑v_v 参考: txc的Min_25筛学习笔记 yx的Min_25筛学习笔记由于前面两位dalao的标题都是笔记,所以我这里就是小记了,因为这 ...
Min_25 筛小结
Min_25 筛这个东西,完全理解花了我很长的时间,所以写点东西来记录一些自己的理解. 它能做什么对于某个数论函数 \(f\),如果满足以下几个条件,那么它就可以用 Min_25 筛来快速求出这个函 ...
Min_25筛有关求解次小质因子
#188. [UR #13]Sanrd 题意化简就是求次小质因子,这一步我们可以在Min_25筛的ans计算中得到, S(n, j)表示的是最小质因子大于等于primejprime_jprimej的 ...
LOJ6053简单的函数（min_25筛）
LOJ6053简单的函数 https://loj.ac/p/6053 min_25筛模板题,但是要注意质数点有2比较特殊,需要在y==0的时候特判质数包含2的情况. #include<bits/ ...
【UOJ188】 Sanrd【类min_25筛】
题意:设f(i)f(i)f(i)表示iii的不严格次大质因子(没有为000),求∑i=lrf(i)\sum_{i=l}^rf(i)∑i=lrf(i) l≤r≤1011l\leq r\leq10^{1 ...
[51NOD1847]奇怪的数学题（杜教筛+min_25筛+第二类斯特林数）
f(x)f(x)f(x)表示xxx的次大约数,有f(x)=xx的最小质因数f(x)=\frac{x}{x的最小质因数}f(x)=x的最小质因数x,那么 ∑i=1n∑j=1nsgcd(i,j)k=∑i ...

ZOJ The Sum of Unitary Totient（min_25 筛）

The Sum of Unitary Totient

ZOJ The Sum of Unitary Totient（min_25 筛）相关推荐

最新文章

热门文章