SIEVE 线性筛
今天来玩玩筛
英文:Sieve
有什么筛?
这里介绍:素数筛,欧拉筛,约数个数筛,约数和筛
为什么要用筛?
顾名思义,筛就是要漏掉没用的,留下有用的。最终筛出来1~n的数的一些信息。
为什么要用线性筛?
考虑最基础的线性筛素数,是O(n)的。
而一般的做法是:
1.对于每个m暴力枚举1~sqrt(m)看能否被整除。O(nsqrt(n))
2.对于每个找到的素数,用它去将所有它的倍数的数都干掉。O(nlogn)
但是,即使是第二种,也有一个log
这是因为一个合数会被它的所有质因子筛一次。要重复质因子个数次,除第一次之外都没用。
所以用线性筛
线性筛原理:
一个算法,使得每个合数只被它的最小质因子筛一次。
怎么保证呢?
素数线性筛:
先看代码:
#include<bits/stdc++.h>
using namespace std;
const int N=+;
int ps[N],cnt;
bool v[N];
int n,m;
void sieve(){
for(int i=;i<=n;i++){
if(v[i]==){
ps[++cnt]=i;
}
for(int j=;j<=cnt;j++){
if(i*ps[j]>n) break;
v[i*ps[j]]=;
if(i%ps[j]==) break;
}
}
}
int main()
{
scanf("%d%d",&n,&m);
sieve();
v[]=;
int t;
for(int i=;i<=m;i++){
scanf("%d",&t);
if(v[t]) printf("No\n");
else puts("Yes");
}
return ;
}
看不懂...................
解释:
对于一个质数,之前没有被标记。肯定只会有一次查到。把它记录到素数集合里去。
然后,不论这个数是否为质数,都将已有的质数从1~cnt循环一遍,把所有的i*ps[j]标记。
当i*ps>n break,可以理解。
当i%ps==0 break.???
这个时候,ps和i不互质了,而ps第一次整除i,所以ps就是i的最小质因数。叫他ps0
而之后,ps更大,ps*i的最小质因数就不是ps了。因为i里有ps0,这个合数就等着以后i更大了,通过ps0筛掉。
ps再大,后面的ps*i的最小质因数都不是ps,所以之前直接break掉就好。
由于每个合数只会被i*ps的形式找到一次(那一次的ps就是这个合数的最小质因子)。而内层循环每一次都对应一个将v[ps*i]=1的操作。
所以内层循环均摊O(1),总共O(n)
完毕。
欧拉线性筛:
代码:fai(i) 1~i中和i互质的数的个数。
可以容斥推出公式:假设:i=p1^q1*p2^q2*....pn^qn
那么,fai(i)=p1^(q1-1)*(p1-1) * p2^(q2-1)*(p2-1) * ......pn^(qn-1) * (pn-1)
证明不是本篇想讲的。
void sieve(){
fai[]=;
for(int i=;i<=n;i++){
if(v[i]==){
fai[i]=i-;
pri[++cnt]=i;
}
for(int j=;j<=cnt;j++){
if(i*pri[j]>n) break;
v[i*pri[j]]=;
if(i%pri[j]==) {
fai[i*pri[j]]=fai[i]*pri[j];break;
}
else{
fai[i*pri[j]]=fai[i]*(pri[j]-);
}
}
}
}
并不想从积性函数性质入手解释。
显然的,当处理到fai[i]的时候,i的值应该就知道了。i是质数就现成赋值。
考虑公式。
当i%ps==0 时,i的质因子中有ps,那么i*ps的质因子ps的次数就大于一,那么,就是fai[i]*ps了
否则,i*ps 的 ps的次数就是1,那么,ps^(1-1)*(ps-1)=(ps-1) 所以是fai[i]*(ps-1)
之后的各种操作基于线性筛的要求和特质。(即每个数只被它的最小质因子筛一次)
例题:SDOI2008 仪仗队
约数个数线性筛:
推荐:线性筛约数个数和、约数和
设x=p1^q1*p2^q2*....pn^qn
要知道公式:个数=(q1+1)*(q2+1)*...*(qn+1) 乘法原理就可以知道。
设t[i]表示i的约数个数
设e[i]表示i的最小素因子个数
①i是质数:t[i]=2,e[i]=1;
②i%pj!=0 这个时候,pj里面没有i,根据积性函数,或者乘法原理,t[i*pj]=t[i]*t[pj]=2t[i];
而 e[i*pj]=1
③i%pj==0 这个时候,pj里面至少有一个i,i也是pj的最小质因子。
t[i*pj]=t[i]/(e[i]+1)*(e[i]+2) 考虑公式,i*pj只在pj的位置上加了1,所以先除掉,再乘上去。
e[i*pj]=e[i]+1 最小素因子个数多了一个。
约数和的线性筛:
(很详细的解释)
设x=p1^q1*p2^q2*....pn^qn
首先还是要知道公式:和=(1+p1^1+...+p1^q1)*(1+p2^1+...+p2^q2)*...*(1+pn^1+...+pn^qn)
证明很简单,加数的个数显然就是约数个数,每次选择就是这个约数个数的质因数分解形式,数值就是这个约数的数值。
设t[i]表示i的约数和
设e[i]表示i的最小素因子对约数和的答案的贡献,即:(1+p1^1+...+p1^q1)(假设p1是最小质因子)
①当i是质数的时候,t[i]=i+1;e[i]=i+1;
②i%pj!=0 根据公式、积性函数性质 : t[i*pj]=t[i]*t[pj]
e[i*pj]=1+pj;
③i%pj==0
t[i*pj]=t[i]/e[i]*(pj*e[i]+1)
证明:考虑公式,i里面有pj的贡献,乘了一个pj,相当于多了一个pj^(qj+1)所以除掉后,乘上错位,再加一
而 e[i*pj]=e[i]*pj+1
就这样。代码参考上面的写就是了,没什么难度。
莫比乌斯函数筛:
知道定义就好说:
μ(i)={
0 i有平方因子
1 i的质因子个数为偶数
-1 i的质因子个数为奇数
}
根据定义直接筛就好了。
void sieve(){
u[]=;
for(int i=;i<=N;i++){
if(!vis[i]){
u[i]=-;
pr[++cnt]=i;
}
for(int j=;j<=cnt;j++){
if(pr[j]*i>N) break;
vis[pr[j]*i]=;
if(i%pr[j]==) {
u[pr[j]*i]=;break;
}
else u[pr[j]*i]=-u[i];
}
}
}
例题:bzoj2440 完全平方数
SIEVE 线性筛的更多相关文章
- BZOJ 2693: jzptab [莫比乌斯反演 线性筛]
2693: jzptab Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 1194 Solved: 455[Submit][Status][Discu ...
- BZOJ 2818: Gcd [欧拉函数 质数 线性筛]【学习笔记】
2818: Gcd Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 4436 Solved: 1957[Submit][Status][Discuss ...
- BZOJ 3309: DZY Loves Math [莫比乌斯反演 线性筛]
题意:\(f(n)\)为n的质因子分解中的最大幂指数,求\(\sum_{i=1}^n \sum_{j=1}^m f(gcd(i,j))\) 套路推♂倒 \[ \sum_{D=1}^n \sum_{d| ...
- Codeforces 893E Counting Arrays:dp + 线性筛 + 分解质因数 + 组合数结论
题目链接:http://codeforces.com/problemset/problem/893/E 题意: 共q组数据(q <= 10^5),每组数据给定x,y(x,y <= 10^6 ...
- 洛谷 - P1891 - 疯狂LCM - 线性筛
另一道数据范围不一样的题:https://www.cnblogs.com/Yinku/p/10987912.html $F(n)=\sum\limits_{i=1}^{n} lcm(i,n) $ $\ ...
- Codeforces 1047C (线性筛+因数分解)
题面 传送门 分析 1.暴力做法 首先先把每个数除以gcd(a1,a2-,an)gcd(a_1,a_2 \dots,a_n )gcd(a1,a2-,an) 可以O(namax)O(n\sqrt ...
- 莫比乌斯反演/线性筛/积性函数/杜教筛/min25筛 学习笔记
最近重新系统地学了下这几个知识点,以前没发现他们的联系,这次总结一下. 莫比乌斯反演入门:https://blog.csdn.net/litble/article/details/72804050 线 ...
- bzoj2693--莫比乌斯反演+积性函数线性筛
推导: 设d=gcd(i,j) 利用莫比乌斯函数的性质 令sum(x,y)=(x*(x+1)/2)*(y*(y+1)/2) 令T=d*t 设f(T)= T可以分块.又由于μ是积性函数,积性函数的约束和 ...
- 【BZOJ-4514】数字配对 最大费用最大流 + 质因数分解 + 二分图 + 贪心 + 线性筛
4514: [Sdoi2016]数字配对 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 726 Solved: 309[Submit][Status ...
随机推荐
- 编写自己的dapper lambda扩展-使用篇
前言 这是针对dapper的一个扩展,支持lambda表达式的写法,链式风格让开发者使用起来更加优雅.直观.现在暂时只有MsSql的扩展,也没有实现事务的写法,将会在后续的版本补充. 这是个人业余的开 ...
- Linux服务器性能压力测试
对于新采购的服务器,需要进行有必要的性能测试.这里选择UnixBench工具进行性能测试.记录如下: 1)安装使用下面的脚本使用了最新版UnixBench5.1.3来测试,注释了关于graphic的测 ...
- linux-IO重定向-文本流重定向
输出重定向的追加和覆盖 标准输出就这两种: 覆盖和追加 >> 是重定向操作符 1 是 命令的文件描述符 重定向操作符合文件描述符之间不能存在空白符 否则1会被当做是文件被读取 将正确和错误 ...
- javaScript——DOM1级,DOM2级,DOM3级
DOM0,DOM2,DOM3事件处理方式区别:http://www.qdfuns.com/notes/11861/e21736a0b15bceca0dc7f76d77c2fb5a.html JS中do ...
- 《Linux内核设计与分析》第四章读书笔记
<内核设计与实现>第四章读书笔记 第四章:进程调度 进程(操作系统)程序的运行态表现形式. 进程调度程序,它是确保进程能有效工作的一个内核子系统. 调度程序负责决定将哪个进程投入运行,何时 ...
- b总结
Beta 答辩总结 评审表 组名 格式 内容 ppt 演讲 答辩 总计 天机组 15 15 13 15 14 72 PMS 16 16 15 16 16 79 日不落战队 16 17 17 17 17 ...
- QT下opencv的编译和使用
需要的文件 qt-opensource-windows-x86-mingw491_opengl-5.4.0.exe cmake-3.12.0-rc1-win64-x64.msi opencv-2.4. ...
- SQL Server 递归查询上级或下级组织数据(上下级数据通用查询语法)
查询上级组织数据: WITH OCTE AS ( AS LVL FROM IOV_Users U LEFT JOIN IOV_Organization O ON U.OrgId=O.ID UNION ...
- PAT 甲级 1045 Favorite Color Stripe
https://pintia.cn/problem-sets/994805342720868352/problems/994805437411475456 Eva is trying to make ...
- [转帖] linux下面 vim 数字键无法插入的解决办法
感谢原作者: https://blog.csdn.net/guoyuqi0554/article/details/11477597 这个问题困扰自己好久了.. 刚才解决了 rlwrap的问题 这会儿 ...