算术基本定理

来源:互联网 发布:怎样使用多个域名 编辑:程序博客网 时间:2024/04/29 22:38

例题一 计算N!末尾0的个数

输入:第一行上有个数字,表示接下来要输入数字的个数。然后是m行,每行包含一个确定的正整数n,1<=n<=1 000 000 000

输出:对输入行中每一个数据n,输出一行,其内容是n!中末尾0的个数

分析:对于任意一个正整数,那么其末尾0必然可以分解成2*5,每一个0必然和一个因子5对应。但一个0不一定对应着一个0,因为还需要一个因子2,才能实现一一对应。

对于你n!,在因式分解中,2的因子个数要大于5的因子个数,所以本题就变成了求n!中因子5的个数。

#include<cstdio>#include<iostream>#include<cstring>using namespace std;int main(){    int m;    cin>>m;    while(m--)    {        int n;        while(scanf("%d",&n)!=EOF)        {            int five=5;            int sum=0;            while(five<=n)            {                sum=sum+n/five;                five=five*5;            }            cout<<sum<<endl;        }    }    return 0;}

例题二 组合素数

题目大意:给你两个整数N和P,求出C(2*N,N)被素数p整数的次数

思路:

由算术基本定理可得到N!被素数P整除的次数。

来看这道题,C(2*N,N) = (2*N)! / (N! * N!)。最终结果就是从(2*N)!能被素数P整除的

次数里边减去N!能被素数整除的次数*2。最终结果为:

[2*N/P] + [2*N/P^2] + … + [2*N/P^t] - 2*([N/P] + [N/P^2] + … + [N/P^t])。

其中次数t = logP(2*N),即log10(2*N) / log10(P)。

输入:第一行是一个正整数t,表示测试数据的组数,接下来每组两个数分别是n,p的值,这里1<=n<=1 000 000 000

输出:输出被素数p整除的次数,如果不能被其整除时,输出0.

#include<cstdio>#include<iostream>#include<cstring>#include<cmath>using namespace std;int main(){    int t;    cin>>t;    while(t--)    {        int n,p;        int sum=0;        int tag=1,q;///tag的用法,有点巧妙        double s;        scanf("%d%d",&n,&p);        s=log10(2.0*n)/log10(p);        q=(int)s;        for(int i=1; i<=q; i++)        {            tag=tag*p;            sum=sum+(int)(2*n/tag)-2*(int)(n/tag);        }        cout<<sum<<endl;    }    return 0;}




0 0
原创粉丝点击