Hdu 5750 Dertouzos

BC的时候听说这个题卡时限以为会FST，没想到居然过了
赛后看了一下我的时间复杂度居然是$O(T \times min(\;\lfloor \frac{n}{p} \rfloor,rank(p)+1)\;)$
($rank(p)$是 p <script type="math/tex" id="MathJax-Element-219">p</script>的最小素因子在素数表的下标,即rank(9) =2 ,rank(4) = 1,rank(35)=3
(这个函数是我现造的。。

---

感觉也是十分的神奇

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

const int maxn = 1123456;
bool vis[maxn];
int pri[maxn],plen;

int sieve(){
    memset(vis,0,sizeof(vis));
    plen = 0;
    for(int i=2;i<maxn;i++){
        if(vis[i] == false){
            pri[plen++] = i;
            for(int j=i*2;j<maxn;j+=i)
                vis[j] = true;
        }
    }
    return plen;
}

int main(){
    int v = sieve();
    int T;
    scanf("%d",&T);
    int n,p;
    while(T-- && ~scanf("%d %d",&n,&p)){
        n--;
        int t = n / p;
        int pos = 1;
        if(t <= 1){
            puts("0");
            continue;
        }
        if(p%2==0) pos = p*2<=n;
        else
            while(pri[pos] <= t && p%pri[pos-1]!=0 && pos < v)
                pos++;
        printf("%d\n",pos);
    }
    return 0;
}