NOIP中数学---第2课 质数

什么是质数

思考三秒钟，注意1不是质数
质数（prime number）又称素数，质数定义为在大于1的自然数中，除了1和它本身以外不再有其它因数。在前面学过的筛质数章节中我们已经接触过质数了。2 3 5 7 是前4个质数。
合数就是非质数，每一个合数都可以以唯一形式被写成质数的乘积，即分解质因数。

试除法求质数

例题1：求质数

给定n个正整数ai，判定每个数是否是质数。

输入格式 第一行包含整数n。

接下来n行，每行包含一个正整数ai。

输出格式 共n行，其中第 i 行输出第 i 个正整数ai是否为质数，是则输出“Yes”，否则输出“No”。

数据范围 1≤n≤100, 1≤ai≤2^31−1
输入样例：
2
2
6
输出样例：
Yes
No

纯暴力就可以了，注意优化

#include<iostream>
using namespace std;
int main()
{
    int n;
    cin >> n;
    int x;
    for(int i = 1;i <= n; i ++)
    {
        cin >> x;
        if(x == 1) cout<<"No"<<endl;
        else
        {
            int flag = 1;
        for(int j = 2;j <= x / j; j ++)//不推荐 j*j<=x，因为x=2147483647会越界
        {
            if(x % j == 0) 
            {
                flag = 0;
                break;
            }
        }
        if(flag) cout << "Yes" <<endl;
        else cout <<"No"<< endl;
        }
        
    }
}

注意：
for(int j = 2;j <= x / j; j ++）此处不建议使用如下方式
1、j<=sqrt(x);,调用函数浪费时间
2、不推荐 j*j<=x，因为x=2147483647会越界未完待续

例题2：分解质因数

给定n个正整数ai，将每个数分解质因数，并按照质因数从小到大的顺序输出每个质因数的底数和指数。
输入格式
第一行包含整数n。
接下来n行，每行包含一个正整数ai。
输出格式
对于每个正整数ai,按照从小到大的顺序输出其分解质因数后，每个质因数的底数和指数，每个底数和指数占一行。
每个正整数的质因数全部输出完毕后，输出一个空行。
数据范围
1≤n≤100,
1≤ai≤2∗109
输入样例：
2
6
8
输出样例：
2 1
3 1
2 3

题解概要，暴力枚举n的所有因子，然后写成因子乘积形式。由于枚举的时候从2开始枚举，当枚举到i时，2到i-1的因子已经分解完成，所以这样分解的因子一定是质因子，不存在合数的问题。

#include<iostream>
using namespace std;
int main()
{
    int n;
    cin >> n;
    for(int i = 1;i <= n; i ++)
    {
    	int x;
    	cin >> x;
    	for(int j = 2;j <= x/j; j++)//暴力枚举x的所有因子，因为是从2开始枚举，所以当枚举j时，2到j-1的 因子都已经被枚举过了，所以得到的一定是质因子 
    	{
    		if( x % j ==0)//说明j是质因子
    		{
    			int s = 0;//质因子的指数 
    		while(x % j == 0) 
    		{
    		   
    			x /= j;
    			s++;
				
			}
		    cout << j <<" "<<s<<endl;
			}	
		}
		if(x>1) cout <<x<<" "<<1<<endl;	
		cout<<endl;	
	}
}

练习题
1、分解质因数
2、质数的和与积
3、ProjectEuler 7
4、ProjectEuler 37

素数的定义

1、只能被1和它自己整除的自然数称为素数，特别规定1不属于素数。

根据素数的定义，很明显，如果一个数是素数<==>它的因子只包含1和它本身。因此可以根据判别某个数的因子的方法来判断其是否是素数。

int isprime(int n)
{
    int i;
    for(i=2;i<=(int)sqrt((double)n);i++)
    {
        if(n%i==0)    //如果n存在其它因子,则必定不是素数
         {
            return 0;
        }
    }
    return 1;
}

但是如果要求求出1000000以内的所有素数，上面的方法效率就很低，因此通常采用筛选法去求素数。筛选法：对于一个数n，如果是素数，那么2n，3n，4*n，必定不是素数。

bool isprime[1000001];
int prime[80000];
int num=0;
void getPrime()     //用筛选法求算素数
{
    int i,j;
    for(i=0;i<1000001;i++)
    {
        isprime[i]=true;
    }
    for(i=2;i<=1000;i++)        //如果isprime[i]==true,即i是素数，那么i，2*i，3*i必定不是素数
     {
        for(j=i+i;j<=1000000;j+=i)
        {
            if(isprime[i]==true)
                isprime[j]=false;
        }
    }
    for(i=2;i<1000001;i++)
    {
        if(isprime[i]==true)
        {
            prime[num++]=i;
        }
    }
}

筛选法求素数有一个很通用的算法，就是在遍例该集合时，比方检验一个数N是否素数，用N除以2-N的开方，只要有一个能整除，就说明N不是素数。另外这道题要求用数组来计算。

线筛

所谓"筛选法"指的是"埃拉托色尼(Eratosthenes)筛法"。他是古希腊的著名数学家。他采取的方法是，在一张纸上写上1到100全部整数，然后逐个判断它们是否是素数，找出一个非素数，就把它挖掉，最后剩下的就是素数。

具体做法如下：

先将1挖掉(因为1不是素数)。
用2去除它后面的各个数，把能被2整除的数挖掉，即把2的倍数挖掉。
用3去除它后面的各数，把3的倍数挖掉。
分别用4、5…各数作为除数去除这些数以后的各数。这个过程一直进行到在除数后面的数已全被挖掉为止。例如找1～50的素数，要一直进行到除数为47为止（事实上，可以简化，如果需要找1～n范围内素数表，只需进行到除数为n²(根号n)，取其整数即可。例如对1～50，只需进行到将50²作为除数即可。）
如上算法可表示为：

挖去1;
用刚才被挖去的数的下一个数p去除p后面各数，把p的倍数挖掉;
检查p是否小于n^2的整数部分（如果n=1000, 则检查p<31?），如果是，则返回(2)继续执行，否则就结束;
纸上剩下的数就是素数。

#include <stdio.h>
#include <math.h>
int main(void)
{
    int i;
    int j;
    int a[101];                // 为直观表示，各元素与下标对应，0号元素不用
    for (i = 1; i <= 100; i++) // 数组各元素赋值
        a[i] = i;
    for (i = 2; i < sqrt(100); i++)     // 外循环使i作为除数
        for (j = i + 1; j <= 100; j++)  // 内循环检测除数i之后的数是否为i的倍数
        {
            if (a[i] != 0 && a[j] != 0) // 排除0值元素
                if (a[j] % a[i] == 0)
                    a[j] = 0;           // i后数若为i的倍数，刚将其置0(挖去)
        }
    int n = 0;    // 对输出素数计数, 以控制换行显示
    for (i = 2; i <= 100; i++)    // 输出素数
    {
        if (a[i] != 0)
        {
            printf("%-5d", a[i]); // 输出数组中非0元素(未挖去的数)
            n++;
        }
        if (n == 10)
        {
            printf("\n");         // 每行10个输出
            n = 0;
        }
    }
    printf("\n");
    return 0;
}

练习题
质数中的质数