洛谷10月月赛R1·浴谷八连测R1·提高组 T1

思路很简单,先对k分解质因数,存入一个数组中去。
然后对于n的阶乘分解那个数组的质因数,除数量,取最小值!
刚开始忘特判tot==0,ans开的不够大!

更快更稳定的方法RHO大质数分解,我懒的打了。
不会的可以参考我之前的BLOG(求一个很大数的欧拉函数)

#include <cstdio>
#include <iostream>
#define ll long long
using namespace std;
const int maxm=1e7;
ll prime[maxm],cnt;
bool vis[maxm];
ll z[maxm],num[maxm],tot;
void pre()
{
    vis[1]=1;

    for(int i=2;i<=1e7;i++)
    {
        if(!vis[i])
         prime[++cnt]=i;
        for(int j=1;j<=cnt&&i*prime[j]<=maxm;j++)
         {
            vis[i*prime[j]]=1;
            if(i%prime[j]==0)
             break;
         }
    }
}
int main()
{
    pre();
    ll n,k;

    scanf("%lld%lld",&n,&k);

    for(int j=1;j<=cnt&&k!=1;j++)
    {
        if(prime[j]>n) break;
        if(k%prime[j]==0)
         {
            z[++tot]=prime[j];
            while(k%prime[j]==0)
             k/=prime[j],num[tot]++;
         }
    }
    ll ans=1e17;
    for(int i=1;i<=tot;i++)
    {
        ll t=n;
        ll sum=0;
        while(t)
        {
            t/=z[i];
            sum+=t;
        }
        //printf("%lld\n",sum);
        ans=min(ans,1ll*sum/num[i]);
    }
    ll w=1e17;
    if(ans==w) ans=0;
    /*printf("%lld\n",tot);
    for(int i=1;i<=tot;i++)
     printf("%lld %lld\n",z[i],num[i]);*/
    printf("%lld",ans);
} 
  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值