素数的线性筛 && 欧拉函数
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/ ...
随机推荐
- 手机卫士开发记录之handler错误
- Django_collections01
python manage.py flush python manage.py changepassword username python manage.py createsuperuser pyt ...
- CentOS7安装hive-2.1.0
环境: CentOS7 Hadoop-2.6.4,配置两个节点:master.slave1 mysql-server 过程: 下载.解压hive-2.1.0到/usr/hadoop-2.6.4/thi ...
- Axure7.0汉化方法
下载汉化包 AxureRP7CN_汉化包.rar 首先退出正在运行中的 Axure (如果您正在使用). 将 汉化包.rar 文件解压, 得到 lang 文件夹, 然后将其复制到 Axure 安装目 ...
- select()
select(),确定一个或多个套接口的状态,本函数用于确定一个或多个套接口的状态,对每一个套接口,调用者可查询它的可读性.可写性及错误状态信息,用fd_set结构来表示一组等待检查的套接口,在调用返 ...
- C# 判断 当前设备的IP地址、默认网关、子网掩码在不在同一网段内
要判断两个IP地址是不是在同一个网段,就将它们的IP地址分别与子网掩码做与运算,得到的结果一网络号,如果网络号相同,就在同一子网,否则,不在同一子网. 例:假定选择了子网掩码255.255.254.0 ...
- hdu 1873 看病要排队(优先级队列)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1873 题目大意: 三个医生看病,病人排队看病,病人有优先级,优先级高的提前看病,同样的优先级按先后.I ...
- B1/B2签证拒签
http://www.mcdvisa.com/html/News/USA_visa_news/201529/152917GE.html
- jquery 停止animate动画,并且回复最初状态
// 热门推荐悬浮效果 $("#recom_con li img").mouseenter(function(){ $(this).stop(true, true); $w = p ...
- EF DbModelBuilder
protected override void OnModelCreating(DbModelBuilder modelBuilder) { var model = modelBuilder.Buil ...