Hdu 5750 Dertouzos

BC的时候听说这个题卡时限以为会FST，没想到居然过了
赛后看了一下我的时间复杂度居然是O(T×min(⌊np⌋,rank(p)+1))
(rank(p)是p的最小素因子在素数表的下标,即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;}