【洛谷 P4213】 【模板】杜教筛(Sum)
时间:2021-09-16
本文章向大家介绍【洛谷 P4213】 【模板】杜教筛(Sum),主要包括【洛谷 P4213】 【模板】杜教筛(Sum)使用实例、应用技巧、基本知识点总结和需要注意事项,具有一定的参考价值,需要的朋友可以参考一下。
#include <bits/stdc++.h>
using namespace std;
#define ll long long
int cnt, n;
map<int, long long> mpu, mpphi;
int prime[10000010], mu[10000010], v[10000010], usum[10000010], phi[10000010];
ll phisum[10000010];
ll getusum(int n){
if(n <= 10000000) return usum[n];
if(mpu[n]) return mpu[n];
ll ans = 1;
for(ll l = 2, r = 0; l <= n; l = r + 1){
r = n / (n / l);
ans -= (r-l+1) * getusum(n / l);
}
return mpu[n] = ans;
}
ll getphisum(int n){
if(n <= 10000000) return phisum[n];
if(mpphi[n]) return mpphi[n];
ll ans = (ll)n * ((ll)n+1) / 2;
for(ll l = 2, r = 0; l <= n; l = r + 1){
r = n / (n / l);
ans -= (r-l+1) * getphisum(n / l);
}
return mpphi[n] = ans;
}
int T;
int main(){
mu[1] = phi[1] = usum[1] = phisum[1] = 1;
for(int i = 2; i <= 10000000; ++i){
if(!v[i]){ mu[i] = -1; phi[i] = i-1; prime[++cnt] = i; }
usum[i] = usum[i-1] + mu[i];
phisum[i] = phisum[i-1] + phi[i];
for(int j = 1; j <= cnt && (ll)i * prime[j] <= 10000000; ++j){
v[i * prime[j]] = 1;
if(i % prime[j] == 0){
phi[i * prime[j]] = phi[i] * prime[j];
break;
}
mu[i * prime[j]] = -mu[i];
phi[i * prime[j]] = phi[i] * (prime[j] - 1);
}
}
scanf("%d", &T);
while(T--){
mpu.clear(); mpphi.clear();
scanf("%d", &n);
printf("%lld %lld\n", getphisum(n), getusum(n));
}
}
原文地址:https://www.cnblogs.com/Qihoo360/p/15302328.html
- java线程间通信
- android scrollview嵌套listview计算高度的问题
- 手把手教你安装深度学习软件环境(附代码)
- IntentService简介
- android 网络通信框架volly
- 仿iphone快速导航悬浮球
- WCF Data Service QuickStart
- 仿支付宝手势密码
- android仿iphone的地区选择
- WCF Data Service 的.NET Client 的不支持原生类型服务操作的解决方法
- Elcipse安装gradle插件
- android View层的绘制流程
- RESTful API生命周期管理
- android 的android httpClient详解
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法
- 一篇文章构建你的 NodeJS 知识体系(W字长文)
- 爬取天气数据并解析温度值
- FlutterDojo设计之道—状态管理之路(一)
- 即学即用Kotlin - 协程
- redis学习(三)
- update影响行数使用不当造成的bug
- 基于React+Koa实现一个h5页面可视化编辑器-Dooring
- 使用Spring Boot DevTools优化你的开发体验
- 排序算法——一篇文章搞懂常用的排序算法
- 路径寻优
- Milvus 实战 | 基于 Milvus 的 CORD-19 论文检索引擎
- Laravel 8 新特性和功能优化速览
- Linux 搭建 我的世界(MC) 基岩版服务器
- 虚函数、析构函数、静态函数、多态
- 一日一技:pylint除了检查代码风格,还能做这件事情