【洛谷10月月赛R1提高组】 SAC E#1 - 一道中档题 Factorial

最新推荐文章于 2022-01-29 16:55:57 发布

愤怒的愣头青

最新推荐文章于 2022-01-29 16:55:57 发布

阅读量405

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/w_yqts/article/details/78170884

版权

数论专栏收录该内容

37 篇文章 0 订阅

订阅专栏

题目背景

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

对于60%的数据，n <= 10^15，k <= 10^12

对于100%的数据，n <= 10^18，k <= 10^16

题解
即求n！能整除几个k
sqrt（k）做法，对k分解质因数，判断n！有多少因子。

代码

#include<bits/stdc++.h>
#define mod 666623333
#define inf 10000000
typedef long long ll;
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9'){if (ch=='-')f=-1;ch=getchar();}
    while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
ll n,k,ans;
ll calc(ll n,ll p)
{
    ll ans=0LL,P=p;
    while (P<=n)
    {
        ans+=n/P;
        P*=p;
    }
    return ans;
}
int main()
{
    while (scanf("%lld%lld",&n,&k)!=EOF)
    {
        ans=-1LL+4294967295LL*2147483648LL+2147483648LL;
        for (ll i=2;i*i<=k;i++)
        {
            if (k%i==0)
            {
                ll num=0;
                while (k%i==0){k/=i;num++;}
                ans=min(ans,calc(n,i)/num);
            }
            if (k==1) break;
        }
        if (k!=1) ans=min(ans,calc(n,k));
        printf("%lld\n",ans);
    }
    return 0;
}