杜教筛模板

#define _CRT_SECURE_NO_WARNINGS

#include<iostream>

#include<string>

#include<cstring>

#include<cmath>

#include<ctime>

#include<algorithm>

#include<utility>

#include<stack>

#include<queue>

#include<vector>

#include<set>

#include<map>

#include<unordered_map>

using namespace std;

#define LL long long

const int N = 2e6 + 10;

const long long mod = 998244353;

const double PI = acos(-1);

LL su[N], a[N], mu[N], phi[N], cn;

LL n;

map<LL, LL> mpu, mpp;

void init()

{

a[0] = a[1] = mu[1] = phi[1] = 1;

for (int i = 2; i < N; i++)

{

if (!a[i]) su[++cn] = i, mu[i] = -1, phi[i] = i - 1;

for (int j = 1; j <= cn&&(LL)i*su[j]<N; j++)

{

int pm = su[j];

a[pm * i] = 1;

if (i % pm == 0)

{

phi[i * pm] = phi[i] * pm;

break;

}

mu[i * pm] = -1 * mu[i];

phi[i * pm] = phi[i] * (pm - 1);

}

}

for (int i = 1; i < N; i++)

mu[i] += mu[i - 1], phi[i] += phi[i - 1];

}

LL seekphi(LL n)

{

if (n < N) return phi[n];

if (mpp[n]) return mpp[n];

LL ans = n * (n + 1) / 2;

for (LL l = 2, r; l <= n; l = r + 1)

{

if (n / l == 0) r = n;

else

r = min(n / (n / l),n);

ans -= (r - l + 1) * seekphi(n / l);

}

mpp[n] = ans;

return ans;

}

LL seekmu(LL n)

{

if (n < N) return mu[n];

if (mpu[n]) return mpu[n];

LL ans = 1;

for (LL l = 2, r; l <= n; l = r + 1)

{

if (n / l == 0) r = n;

else

r = min(n / (n / l), n);

ans -= (r - l + 1) * seekmu(n / l);

}

mpu[n] = ans;

return ans;

}

int main() {

init();

int T;

cin >> T;

while (T--)

{

cin >> n;

LL ans1 =seekphi(n),ans2=seekmu(n);

cout << ans1 <<" "<<ans2<< endl;

}

return 0;

}

相关推荐
巫师不要去魔法部乱说1 小时前
PyCharm专项训练4 最小生成树算法
算法·pycharm
IT猿手1 小时前
最新高性能多目标优化算法:多目标麋鹿优化算法(MOEHO)求解GLSMOP1-GLSMOP9及工程应用---盘式制动器设计,提供完整MATLAB代码
开发语言·算法·机器学习·matlab·强化学习
阿七想学习1 小时前
数据结构《排序》
java·数据结构·学习·算法·排序算法
王老师青少年编程1 小时前
gesp(二级)(12)洛谷:B3955:[GESP202403 二级] 小杨的日字矩阵
c++·算法·矩阵·gesp·csp·信奥赛
Kenneth風车2 小时前
【机器学习(九)】分类和回归任务-多层感知机(Multilayer Perceptron,MLP)算法-Sentosa_DSML社区版 (1)111
算法·机器学习·分类
eternal__day2 小时前
数据结构(哈希表(中)纯概念版)
java·数据结构·算法·哈希算法·推荐算法
APP 肖提莫2 小时前
MyBatis-Plus分页拦截器,源码的重构(重构total总数的计算逻辑)
java·前端·算法
OTWOL2 小时前
两道数组有关的OJ练习题
c语言·开发语言·数据结构·c++·算法
QQ同步助手2 小时前
C++ 指针进阶:动态内存与复杂应用
开发语言·c++
qq_433554543 小时前
C++ 面向对象编程:递增重载
开发语言·c++·算法