POJ 1150 (数论)

题意,求n!/(n-m)!末尾第一个非零的数是什么。

思路:这题折腾了我一晚上,终于搞明白了。以10!为例,首先,要去除末尾的0,要把因子10去除,但10不是质数,所以就考虑2*5;去除所有的2,5;例如:1 2 3 4 5 6 7 8 9 10,去除了因子2和5之后,就变为了1 1 3 1 1 1 7 1 9 1;

那么这就相当于一个子问题了,我们可以用递归来求解。

这样,整个序列就只剩下1 3 7 9 这四个数了,我们现在要找的就是序列中数字末尾为3 7 9的个数。这个是本题的重点

我们可以把这个序列分为奇偶,1 3 5 7 9, 和 2 4 6 8 10;

我们尝试统计,可以发现奇数和偶数中3 7 9的个数是一样的。因此,我们依然可以递归求解

f(n) = f(n/2)+g(n)   g(n)为奇数序列中的数目。

我们继续观察奇数的序列 1 3 5 7 9 11 13 15 17 19 21 23 25;

其中5 15 25 是5的奇数倍,那么我们要统计序列中 3 7 9 的个数g(n,x) = n/10 + (n%10>=x) + g(n/5,x);

最后就是把所有这些数相乘。同过观察,2 3 7 9这四个数都是周期为4。这样就能方便的求解了。

还有就是判断一下5和2的个数:

如果num5>num2,直接输出5就行了。

如果num5<=num2,如果不相等的话就要把多出来的2先处理掉,再继续想3 7 9;

最后上代码:

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <vector>

using namespace std;
typedef long long ll;

int get2(int n)
{
    if(n==0)
        return 0;
    return (n/2)+get2(n/2);
}

int get5(int n)
{
    if(n==0)
        return 0;
    return (n/5)+get5(n/5);
}

int g(int n,int x)
{
    if(n==0)
        return 0;
    return n/10+(n%10>=x)+g(n/5,x);//1-n末尾为x的个数
}

int getx(int n,int x)
{
    if(n==0)
        return 0;
    return getx(n/2,x)+g(n,x);
}

int table[4][4] =
{
    6,2,4,8,//2
    1,3,9,7,//3
    1,7,9,3,//7
    1,9,1,9,//9
};

int main()
{
    int n,m;
    while(cin >> n >> m)
    {
        int num2 = get2(n)-get2(n-m);
        int num5 = get5(n)-get5(n-m);
        int num3 = getx(n,3)-getx(n-m,3);
        int num7 = getx(n,7)-getx(n-m,7);
        int num9 = getx(n,9)-getx(n-m,9);
        int res = 1;
        if(num5>num2)
        {
            cout<<"5"<<endl;
            continue;
        }
        else
        {
            if(num5!=num2)
            {
                res *= table[0][(num2-num5)%4];
                res %= 10;
            }
            res *= table[1][num3%4];
            res %= 10;
            res *= table[2][num7%4];
            res %= 10;
            res *= table[3][num9%4];
            res %= 10;
            cout<<res<<endl;
        }
    }
    return 0;
}


阅读更多
想对作者说点什么?

博主推荐

换一批

没有更多推荐了,返回首页