素数的线性筛 && 欧拉函数
O(n) 筛选素数
#include<bits/stdc++.h>
using namespace std;
const int M = 1e6 + 10 ; int mindiv[M] ;//每个数的最小质因数
int prim[M] , pnum ;//存素数
bool vis[M] ; void prim () {
for (int i = 2 ; i < M ; i ++) {
if (!vis[i]) {
mindiv[i] = i ;
prim[ pnum++ ] = i ;
}
for (int j = 0 ; j < pnum ; j ++) {
if ( i*prim[j] >= M ) break ;
vis[ i*prim[j] ] = 1 ;
mindiv[i] = prim[j] ;
if (i % prim[j] == 0) break ;
}
}
} int main () {
prim () ;
return 0 ;
}
欧拉函数:phi[i] 为<= i 的范围内与i互质的数的数量
欧拉埃筛,写起来简单,复杂度O(log(log(N)))(zstu 幻神):
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int M = 1e6 + 10 ; int n, m, T; int euler[M]; void Euler () {
for(int i = 0; i < M ; i ++) euler[i] = i;
for(int i = 2; i < M ; i ++){
if(euler[i] == i) {
for (int j = i; j < M ; j += i) {
euler[j] = euler[j] - euler[j]/i;
}
}
}
} int main(){
Euler ();
int n ;
while (~ scanf ("%d" , &n)) printf ("%d\n" , euler[n]) ;
return 0;
}
欧拉线筛,写起来复杂点,(墨迹了我半天)复杂度O(N):
#include<bits/stdc++.h>
using namespace std;
const int M = 1e6 + 10 ;
int prim[M] , pnum ;
bool vis[M] ;
int phi[M] ; void Euler () {
for (int i = 2 ; i < M ; i ++) {
if (!vis[i]) {
prim[ pnum++ ] = i ;
phi[i] = i - 1;
}
for (int j = 0 ; j < pnum ; j ++) {
int x = i * prim[j] ;
if (x >= M ) break ;
vis[x] = 1 ;
if (i % prim[j] == 0) {
int y = i , cnt = 0 , z = prim[j] ;
while (y % prim[j] == 0) cnt ++ , y /= prim[j] , z *= prim[j] ;
if (y == 1) phi[x] = x - x/prim[j] ;
else phi[x] = phi[y] * phi[z] ;
break ;
}
else phi[x] = phi[i] * phi[ prim[j] ] ;
}
}
} int main () {
Euler () ;
int n ;
while (~ scanf ("%d" , &n)) printf ("%d\n" , phi[n]) ;
return 0 ;
}
线性欧拉跟新:
#include<cstdio>
#include<iostream>
using namespace std;
int prime[100005],phi[1000005];
int main(){
int i,j;
for(i=2;i<1000002;++i){
if(!phi[i]){
phi[i]=i-1;
prime[++prime[0]]=i;
}
for(j=1;j<=prime[0]&&(long long)i*prime[j]<1000002;++j)
if(i%prime[j])phi[i*prime[j]]=phi[i]*(prime[j]-1);
else{
phi[i*prime[j]]=phi[i]*prime[j];
break;
}
}
int T,n;
scanf("%d",&T);
while(T--){
scanf("%d",&n);
printf("%d\n",phi[n+1]);
}
}
素数的线性筛 && 欧拉函数的更多相关文章
- The Euler function(线性筛欧拉函数)
/* 题意:(n)表示小于n与n互质的数有多少个,给你两个数a,b让你计算a+(a+1)+(a+2)+......+b; 初步思路:暴力搞一下,打表 #放弃:打了十几分钟没打完 #改进:欧拉函数:具体 ...
- 积性函数&线性筛&欧拉函数&莫比乌斯函数&因数个数&约数个数和
只会搬运YL巨巨的博客 积性函数 定义 积性函数:对于任意互质的整数a和b有性质f(ab)=f(a)f(b)的数论函数. 完全积性函数:对于任意整数a和b有性质f(ab)=f(a)f(b)的数论函数 ...
- [bzoj 2190][SDOI2008]仪仗队(线性筛欧拉函数)
题目:http://www.lydsy.com/JudgeOnline/problem.php?id=2190 分析:就是要线性筛出欧拉函数... 直接贴代码了: memset(ans,,sizeof ...
- Bzoj 2186: [Sdoi2008]沙拉公主的困惑 乘法逆元,线性筛,欧拉函数,数论
2186: [Sdoi2008]沙拉公主的困惑 Time Limit: 10 Sec Memory Limit: 259 MBSubmit: 2560 Solved: 857[Submit][St ...
- BZOJ 2190 仪仗队(线性筛欧拉函数)
简化题意可知,实际上题目求得是gcd(i,j)=1(i,j<=n)的数对数目. 线性筛出n大小的欧拉表,求和*2+1即可.需要特判1. # include <cstdio> # in ...
- poj1248 (线性筛欧拉函数)(原根)
强烈鸣谢wddwjlss 题目大意:给出一个奇素数,求出他的原根的个数,多组数据. 这里先介绍一些基本性质 阶 设\((a,m)=1\),满足\(a^r \equiv 1 \pmod m\)的最小正整 ...
- BZOJ 2818 GCD 素数筛+欧拉函数+前缀和
题目链接:https://www.lydsy.com/JudgeOnline/problem.php?id=2818 题意:给定整数N,求1<=x,y<=n且Gcd(x,y)为素数的数对( ...
- noip复习——线性筛(欧拉筛)
整数的唯一分解定理: \(\forall A\in \mathbb {N} ,\,A>1\quad \exists \prod\limits _{i=1}^{s}p_{i}^{a_{i}}=A\ ...
- Farey Sequence (素筛欧拉函数/水)题解
The Farey Sequence Fn for any integer n with n >= 2 is the set of irreducible rational numbers a/ ...
随机推荐
- [U3D 导出Xcode工程包,用Xcode给U3D脚本传递参数]
1.导出Xcode工程 File->Building and setting,导出IOS工程(有错误会失败) 2.运行Xcode工程,在Classes文件夹的UI文件夹里,早到UnityAppC ...
- python时间模块-time和datetime
时间模块 python 中时间表示方法有:时间戳,即从1975年1月1日00:00:00到现在的秒数:格式化后的时间字符串:时间struct_time 元组. struct_time元组中元素主要包括 ...
- Visual Studio原生开发的10个调试技巧
这篇文章只介绍了一些有关Visual Studio的基本调试技巧,但是还有其他一些同样有用的技巧.我整理了一些Visual Studio(至少在VS 2008下)原生开发的调试技巧.(如果你是工作在托 ...
- RNN 入门教程 Part 4 – 实现 RNN-LSTM 和 GRU 模型
转载 - Recurrent Neural Network Tutorial, Part 4 – Implementing a GRU/LSTM RNN with Python and Theano ...
- linux学习笔记-dump命令的使用
http://blog.chinaunix.net/uid-29797586-id-4458302.html
- Yocto开发笔记之《Makefile编写》(QQ交流群:519230208)
开了一个交流群,欢迎爱好者和开发者一起交流,转载请注明出处. QQ群:519230208,为避免广告骚扰,申请时请注明 “开发者” 字样 =============================== ...
- 【原】javascript数组操作
继续我的第二遍<javascript高级程序设计第三版>,今天要做的笔记是array 一.数组的操作 1.数组的创建: var colors= new Array(); //创建一个数组 ...
- 之前总结的今天给大分享一下iOS
退回输入键盘 苹果 ios 开发一年的工作笔记 - (BOOL) textFieldShouldReturn:(id)textField{ [textField resignFirstResponde ...
- select例子
好长时间没有写了,其实一直在坚持学习. #include <sys/types.h> #include <sys/socket.h> #include <stdio.h& ...
- JDK的安装和配置
JDK8 是JDK的最新版本,加入了很多新特性,如果我们要使用,需要下载安装: JDK8在windows xp下安装有点问题,所以在WIN7下安装 WIN7操作系统有32位和64位,分别要下载对应的J ...