关于一道题引发的惨案

 [牛客]第三届超越杯程序设计团体赛重现赛Find And Cal

Ruri 遇到了两个函数：f(x)表示大于 x的最小的质数。g(x,y) 表示x!包含因子 y的个数。Ruri 想知道对于一个给定 (x,y) 对，g(f(x),y) 的大小。

输入描述:

第一行一个整数 T(T<=)，表示询问组数。​
后面 T 行为两个正整数 x,y(1≤x≤1e16,2≤y≤1e9)

输出描述:

 T 行，每行一个整数，g(f(x),y)g(f(x),y)g(f(x),y) 的大小。答案对 998244343 取模。

这道题的数据有问题，我已经联系过出题人了（有问题！）下面是我认为的AC代码

#include<bits/stdc++.h>
#define ll __int128
using namespace std;
inline ll read(){
    char ch=getchar();ll s=0,w=1;
    while(ch<48||ch>57){
        if(ch=='-')w=-1;ch=getchar();
    }
    while(ch>=48&&ch<=57){
        s=(s<<1)+(s<<3)+ch-48;ch=getchar();
    }
    return s*w;
}

inline void write(ll x){
    if(x<0)putchar('-'),x=-x;
    if(x>9)write(x/10);
    putchar(x%10+'0');
}

inline ll qpow(ll x,ll n,ll mod){//x^n mod mod
    x%=mod;
    ll ans=1;
    while(n){
        if (n&1)ans=ans%mod*x%mod;
        n>>=1;
        x=x%mod*x%mod;
    }
    return ans;
}

bool Rabin_Miller(ll p){//p=2^b*a;
    if (p<3 || p%2==0) return p==2;
    ll a=p-1,b=0;
    while (a&1==0)a>>=1,++b;
    // test_time 为测试次数,建议设为不小于8的整数以保证正确率,但也不宜过大,否则会影响效率
    int test_time=8;
    for (int i=1,j; i<=test_time; ++i) {
        int x=rand()%(p-2)+2,v=qpow(x,a,p);
    if(v==1)continue;
    for(j=0;j<b;++j){
      if(v==p-1)break;
      v=(ll)v*v%p;
    }
    if (j>=b)return false;
  }
  return true;
}

const int maxn=34000;
ll prime[maxn+100],cnt;
bool st[maxn+100];//0表示不是合数,即为素数
ll factor[100];
const ll mod=998244343;

int main(){
    ll T;
    T=read();
    for(int i=2;i<=maxn;i++){
        if(st[i]==0)prime[cnt++]=i;
        for(int j=0;prime[j]*i<=maxn;++j){
                
            st[i*prime[j]]=true;
            if(i%prime[j]==0)break;
        }
    }

    while(T--){
        ll x,y;
        x=read(),y=read();
        ++x;
        while(!Rabin_Miller(x))++x;

        ll ans=0;
        for(ll j=0;prime[j]*prime[j]<=y;++j){
            if(Rabin_Miller(y))break;
            if(y%prime[j]==0){
                ll cnt=0,x1=x,cmp=0;
                while(y%prime[j]==0)++cnt,y/=prime[j];
                while(x1){cmp+=x1/prime[j]; x1/=prime[j]; }
                cmp/=cnt;
                if(cmp>ans)ans=cmp;
            }
        }
        if(y>1){
            ll cmp=0;
            while(x){cmp=(cmp+x/y)%mod; x/=y;}
            if(cmp>ans)ans=cmp;
        }
        write(ans%mod);putchar('\n');
    }
    return 0;
}

虽然是错的，但是让我学到了很多东西（郁闷了好久，一直wawa）

首先是对于大素数的检测，一般情况下使用Rabin_Miller进行检测，famer也可以

然后对于一个数的质因子分解，首先可以用欧拉筛打表，得到小数据规模的质数，

然后分解过程中，y不断除它的质因子减小规模，当然保证尝试质数的最大上限为根号y

考虑到y可能含有较大的质因子，可以在过程中对y使用Rabin_Miller进行检测