ACdream 1084 寒假安排(阶乘素因子分解)

题意：将A(n , m)转化成k进制后末尾0的个数。A(n , m)=n!/(n-m)!

这个题数据很大，需要预处理==

而且自己第二天做想偷懒，不分解k，果然wa了T^T

关于不拆k为什么不对：

分析例子 n=9,m=6,k=4 错的结果是2 对的结果是3 拆一下就懂了

然后还是

int getNum(int n,int p){
    if(n<p) return 0;
    return getNum(n/p,p)+n/p;
}求素因子

求A(n, m)转化成k进制以后末尾0的个数。对k素因子分解，第i个因子为fac[i]，

第i个因子的指数为num[i],然后再对A(n, m)进行素因子分解，设count[i]

代表fac[i]对应的指数，ans = min{count[i]/num[i]} 

#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
using namespace std;

const int maxn = 1e6+10;
int pri[maxn],cnt;
bool vis[maxn];
int fac[1000],num;
int Count[1000];
int sum[maxn];
void getFactor(int x){
    num=0;
    memset(Count,0,sizeof(Count));
    for(int i=0;i<cnt&&pri[i]*pri[i]<=x;i++){
        if(x%pri[i]==0){
            fac[num]=pri[i];
            while(x%pri[i]==0) Count[num]++,x/=pri[i];
            num++;
        }
    }
    if(x>1) fac[num]=x,Count[num++]=1;
}

int getNum(int n,int p){
    if(n<p) return 0;
    return getNum(n/p,p)+n/p;
}

int main()
{
  //  freopen("cin.txt","r",stdin);
    memset(vis,0,sizeof(vis));
    cnt=0;
    for(int i=2;i<maxn;i++){
        if(!vis[i]){
            pri[cnt++]=i;
            for(int j=i+i;j<maxn;j+=i) vis[j]=1;
        }
    }
    int t,n,m,k;
    scanf("%d",&t);
    while(t--){
        scanf("%d%d%d",&n,&m,&k);
        getFactor(k);
        int ans = 1000000000;
        for(int i=0;i<num;i++){
            //cout<<getNum(n,fac[i])<<" "<<getNum(m,fac[i])<<endl;
            sum[i]=getNum(n,fac[i])-getNum(n-m,fac[i]);
            ans=min(sum[i]/Count[i],ans);
        }
        printf("%d\n",ans);
    }
    return 0;
}