暑期入门练习（简单数论,质数筛选）

做了入门练习，基本都是跟数学有关的，其中大多跟质数相关。

1. 约数个数定理

对于一个大于1正整数n可以分解质因数：n=p ₁^a ₁*p ₂^a ₂*p ₃^a ₃*…*p _k^a _k,

则n的正约数的个数就是（a ₁+1）（a ₂+1）（a ₃+1）…(a _k+1) .

其中a ₁、a ₂、a ₃…ak是p ₁、p ₂、p _3，…p _k的指数。

#include <stdio.h>
#include <math.h>

int main()
{
	int i,j,k,ans,n;
	while(scanf("%d",&n) && n!=0)
	{
		ans=1;
		for(i=2;n>1;i++)
		{
			k=0;
			while(n%i==0)
				n/=i,k++;
			ans*=(k+1);
		}
		printf("%d\n",ans);
	}
    return 0;
}

2.求质数

采用筛选法，每次在过程中就不断除去该素数。将一个数分解为质数乘积。

/*

将一个数的质数筛出来。利用较小因子相乘不会超过该数本身。不要忘记最后还有一个Num判断。

*/

	__int64 i,j;
	__int64 temp,k;

	yz=0;
	for(i=2;i*i<=num;i++)
	{
		if(num%i==0)
		{
			C[yz++]=i;
			while(num%i==0)
				num/=i;
			
		}
	}
	if(num>1) C[yz++]=num;//C中存放质因子

 

/*
原始的是从筛出质数角度去筛其他数。这个是从一个数字可被分解的角度。
快速筛选质数。1 的意思是将质数放入prim中。2 是避免越界的话 3.是筛 
原来的素数筛选是发现一个素数，然后将其倍数筛掉，一直筛选到最后。提升效率，一个数字被筛，只会发生在其最大因子时候 。eg: 6 ，2的时候不筛 3筛。
既然当最大因子的时候筛，那么最小因子必然在已经得到的素数表中。
那么迭代中的i这个数怎样才算做最大的因子？这个数字与一个质数相乘为止。
*/
    for(int i=2;i<10000;i++)
    {
        if(vis[i]==0)
            prim[p++]=i;//
        for(int j=0;j<p;j++)
        {
            if(i*prim[j]>10000)break;//2
            vis[i*prim[j]]=1;//3
            if(i%prim[j]==0) break;//4
        }
    }

3. 求最大最小公约数

gcd(a,b) 辗转相除法求公约数

对于乘除都有的要先除后乘，可以避免溢出大数。注：虽然使用64位做为最后输出，但在中间过程仍然是int 会有溢出产生

 

4. 求一个数字n！后有多少位0

定理：小质数的个数必然比大质数多10=2*5；则转化为求阶乘中5的个数

转为求5， 5*5， 5*5*5，....的个数   eg:5 15 20 25 30 35........刚好是到n的5的倍数。

 

5.巧留尾数        3.2.1 Factorials 阶乘

 这个题目所用为阶乘，根本无法存储。因此只保留尾数即可，但是若将后面几个0求出，每次取模也会溢出。因此每次都模10，但是剔除在乘的过程中的2 与5。5的个数用4的方法来计算。注意：在过程中要使用循环去除，因为有的数字包含多个5.

#include <stdio.h>
#include <math.h>
int main()
{
    int sum,i,j,k,m,n;
	while(scanf("%d",&n)!=EOF)
	{
		k=5;
		sum=0;
		while(n/k!=0)
		{
			sum+=n/k;
			k*=5;
		}
		k=sum;
		j=sum;
		sum=1;
		for(i=2;i<=n;i++)
		{
			sum=sum*i;
			while(sum%2==0 &&k)
				sum/=2,k--;
			while(sum%5==0 &&j)
				sum/=5,j--;
			
			sum%=10;
			
			
		}
		printf("%d\n",sum);
			
	}
    return 0;
}