noip模拟赛 SAC E#1 - 一道中档题 Factorial

题目背景

数据已修改

SOL君(炉石主播)和SOL菌(完美信息教室讲师)是好朋友。

题目描述

SOL君很喜欢阶乘。而SOL菌很喜欢研究进制。

这一天,SOL君跟SOL菌炫技,随口算出了n的阶乘。

SOL菌表示不服,立刻就要算这个数在k进制表示下末尾0的个数。

但是SOL菌太菜了于是请你帮忙。

输入输出格式

输入格式:

 

每组输入仅包含一行:两个整数n,k。

 

输出格式:

 

输出一个整数:n!在k进制下后缀0的个数。

 

输入输出样例

输入样例#1:
10 40
输出样例#1:
2

说明

对于20%的数据,n <= 1000000, k = 10

对于另外20%的数据,n <= 20, k <= 36

对于100%的数据,n <= 10^12,k <= 10^12

update

1.一组数据

2.K不会==1

3.现在std没有爆long long

4.对数据有问题联系icy (建议大家不要面向数据编程)

分析:k=10就是一个非常经典的问题,将n!分解质因数,看看有多少个5就可以了.现在考虑k≠10的情况,模拟一下进制转换,就是不断地%k,/k,我们就是要求一开始k能整除n!多少次.

      n!,k过大,显然不能直接算,涉及到倍数,我们可以分解一下k,看看k的质因子是不是n!都有,并且n中的次数大于等于k中的次数.假设n!有一个因子p1,次数为b1,它在k中的次数为b2,每做一次除法就是b1-=b2,每次做除法后都必须保证k中的质因子n!都有,并且次数不能小于k中的.能做多少次除法呢?把相同质因子的次数相除,取个min就可以了.

      还有一个问题:k可以直接分解质因数,但是n!就非常大了,我们根本算不出来,怎么才能统计k中的每个质因数在n!中出现的次数呢?对于一个质因数p,1~n个数中有n / p个数包含这个质因子,有n / (p ^ 2)个数包含p^2,也就是相当于包含了p.我们求出包含了p^i的个数,最后相加一下就是答案了.因为阶乘嘛,最后乘起来等于指数的相加.

#include 
#include 
#include 
#include 
#include 

using namespace std;

long long n, k;
long long cnt, p[1000010], tot[1000010], ans = 10000000000000000;

void f(long long x)
{
    for (long long i = 2; i * i <= x; i++)
    {
        if (x % i == 0)
        {
            p[++cnt] = i;
            while (x % i == 0)
            {
                tot[i]++;
                x /= i;
            }
        }
    }
    if (x != 1)
    {
        p[++cnt] = x;
        tot[x]++;
    }
}

long long cal(long long x, long long y)
{
    //printf("%lld %lld\n", x, y);
    if (y > x)
        return 0;
    return x / y + cal(x / y, y);
}

int main()
{
    scanf("%lld%lld", &n, &k);
    f(k);
    for (int i = 1; i <= cnt; i++)
        ans = min(ans, cal(n, p[i]) / tot[p[i]]);
    printf("%lld\n", ans);

    return 0;
}

 

转载于:https://www.cnblogs.com/zbtrs/p/7653465.html

你可能感兴趣的:(noip模拟赛 SAC E#1 - 一道中档题 Factorial)