1.定义
对大于1的自然数,只包含1和它本身两个约数(只能被1和它本身整除)
2.质数的判定——试除法
试除法判定质数
给定 n 个正整数 ai,判定每个数是否是质数。
输入格式
第一行包含整数 n。
接下来 n 行,每行包含一个正整数 ai。
输出格式
共 n 行,其中第 i 行输出第 i个正整数 ai 是否为质数,是则输出 Yes
,否则输出 No
。
数据范围
1≤n≤1001≤n≤100,
1≤ai≤231−11≤ai≤231−1
输入样例:
2
2
6
输出样例:
Yes
No
#include<iostream>
using namespace std;
bool is_prime(int n)
{
if(n<=1) return false;
for(int i=2;i<=n/i;i++)
{
if(n%i==0) return false;
}
return true;
}
int main()
{
int n;
cin>>n;
while(n--)
{
int x;
cin>>x;
if(is_prime(x)) puts("Yes");
else puts("No");
}
return 0;
}
3.分解质因数——试除法
分解质因数
给定 nn 个正整数 aiai,将每个数分解质因数,并按照质因数从小到大的顺序输出每个质因数的底数和指数。
输入格式
第一行包含整数 nn。
接下来 nn 行,每行包含一个正整数 aiai。
输出格式
对于每个正整数 aiai,按照从小到大的顺序输出其分解质因数后,每个质因数的底数和指数,每个底数和指数占一行。
每个正整数的质因数全部输出完毕后,输出一个空行。
数据范围
1≤n≤1001≤n≤100,
2≤ai≤2×1092≤ai≤2×109
输入样例:
2
6
8
输出样例:
2 1
3 1
2 3
#include<iostream>
using namespace std;
void divide(int n)
{
for(int i=2;i<=n/i;i++)
{
/*假如 i 是一个合数,那么它一定可以分解成多
个质因子相乘的形式,这多个质因子同时也是 a
的质因子且比 i 要小,而比 i 小的数在之前的
循环过程中一定是被条件除完了的,所以 i 不可
能是合数,只可能是质数*/
if(n%i==0)
{
int s=0;
while(n%i==0)
{
n=n/i;
s++;
}
printf("%d %d\n",i,s);
}
}
if(n>1) printf("%d %d\n",n,1);
puts("");
}
int main()
{
int n;
cin>>n;
while(n--)
{
int x;
cin>>x;
divide(x);
}
}
4.筛质数
筛质数
给定一个正整数 nn,请你求出 1∼n1∼n 中质数的个数。
输入格式
共一行,包含整数 nn。
输出格式
共一行,包含一个整数,表示 1∼n1∼n 中质数的个数。
数据范围
1≤n≤1e6
输入样例:
8
输出样例:
4
朴素筛法
#include<iostream>
using namespace std;
int n;
const int N=1e6+10;
int primes[N],cnt;
bool st[N];
void get_prime()
{
for(int i=2;i<=n;i++)
{
if(!st[i]) primes[++cnt]=i;
for(int j=i;j<=n;j+=i)
{
st[j]=true;
}
}
}
int main()
{
cin>>n;
get_prime();
cout<<cnt;
}
埃氏筛法
#include<iostream>
using namespace std;
int n;
const int N=1e6+10;
int primes[N],cnt;
bool st[N];
void get_prime()
{
for(int i=2;i<=n;i++)
{
if(!st[i])
{
primes[++cnt]=i;
for(int j=i;j<=n;j+=i) st[j]=true;
}
}
}
int main()
{
cin>>n;
get_prime();
cout<<cnt;
}
线性筛法
#include<iostream>
using namespace std;
int n;
const int N=1e6+10;
int primes[N],cnt;
bool st[N];
void get_prime(int n)
{
for(int i=2;i<=n;i++)
{
if(!st[i]) primes[cnt++]=i;
for(int j=0;primes[j]<=n/i;j++)
{
st[primes[j]*i]=true;
if(i%primes[j]==0) break;
}
}
}
int main()
{
cin>>n;
get_prime(n);
cout<<cnt;
}