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

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

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

鬼题

如果k<=10^12是可以稳过的

k<=10^16就只剩pollard Rho(超出提高组模拟范围了)或玄学优化

首先k进制最后x位为0表示n！能整除k数x次(显然而重要的条件)

于是问题变为了n!含有多少个k的因子

k为合数不好处理，分解质因数k=p1^a1*p2^a2*p3^a3.....

对于求n!含k的数量，转变为求p1....的因子数，在除以a1...，取最小值

给出2种方法

1.

n!=(1k*2k*3k*4k*5k...mk)*a

=k^m*m!*a

a为不含k的积，m=n/k

这样就可以log递归求解

2.

含有一个k：n/k

含有二个k：n/k/k

.....

含有p个k：n/(k^p)

把它们相加

 #include<iostream>
 #include<cstdio>
 #include<algorithm>
 #include<cstring>
 #include<cmath>
 #define count COUNT
 using namespace std;
 long long prime[],mi[],tot;
 long long n,k,cnt,ans=,s;
 long long count(long long p,long long k,long long Mod)
 {
   long long num=,ss=;
   while (p)
     {
       num+=p/k;
       if (num>=Mod) ss+=,num-=Mod;
       p/=k;
     }
   return ss;
 }
 int main()
 {register long long i;
   while (scanf("%lld%lld",&n,&k)==)
    {tot=;ans=;
   long long p=sqrt(k);
   for (i=;i<=p;i++)
     {
       if (k%i==)
     {
       prime[++tot]=i;
       cnt=;
       while (k%i==)
         {
           k/=i;
           cnt++;
         }
       mi[tot]=cnt;
     }
       if (k==) break;
     }
   if (k!=)
     {tot++;
       prime[tot]=k;
       mi[tot]=;
     }
   for (i=;i<=tot;i++)
     {
       s=count(n,prime[i],mi[i]);
       if (ans>s) ans=s;
     }
   printf("%lld\n",ans);
  }
 }