【JZOJ】3188 找数

Description

求第$N$个最小质因子为$P$的数。
$1 <= N , P <= 10^9$
若该数大于$10^9$则输出0

Brute Force

如果要暴力，该怎么暴力？

void BruteForce()
{
    int cur = P;
    for ( ; cur <= 1e9 && n ; cur += P)
        if (MinPrime[cur] >= P) n --;
    if (n) ans = 0; else ans = cur - P;
}

就是暴力地跑即可。

Solution

观察暴力，我们发现，其实$cur$根本没必要那么大。
因为我们枚举的东西一定是$P$的倍数，只要枚举$cur/P$即可。

$cur/P$仍然可能非常大。此时我们就要分段。
对于$P>70$
$cur/P$只有$1.5 * 10^7$，用上面的暴力去跑就行了。注意此时我们枚举的是$/P$之后的结果，每次加1就行。
对于$P<=70$
其实这一段的素数非常少，只有$20+$个。
这样我们就可以容斥了。
求出比$P$小的质数是哪些。
对于$cur/P$二分答案，看看里面出去比$P$小的质数的倍数之后剩下的数是不是比$N$个数多。
得解。

比赛的时候没细想，暴力都没打。真的是非常地智障。去RUSH第一题最难打的一种算法也不来打这题的暴力。这题暴力打完观察下就能出正解，因为正解也不过是个暴力。

暴力出奇迹。

代码如下：

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

const int MAXN = 15000010;

int Miprime[MAXN],prime[MAXN / 10],n,P,ans,cnt;

void Preprocessing()
{
    for (int i = 2 ; i < MAXN ; i ++)
    {
        if (! Miprime[i]) prime[++ prime[0]] = i , Miprime[i] = i;
        for (int j = 1 ; j <= prime[0] ; j ++)
        {
            if (i * prime[j] > MAXN) break;
            Miprime[i * prime[j]] = prime[j];
            if (i % prime[j] == 0) break;
        }
    }
}

int Calc(int limit)
{
    int mxt = 1 << cnt , sum = limit;
    for (int i = 1; i < mxt ; i ++)
    {
        long long p = 1 , nt = 0;
        for (int j = 1 ; j <= cnt ; j ++)
            if (i & (1 << (j - 1))) nt ++ , p *= prime[j];
        if (nt & 1) sum -= limit / p; else sum += limit / p;
    }
    return sum;
}

int main()
{
    scanf("%d%d" , &n , &P);

    Preprocessing();

    if (P > 70)
    {
        int p = (int)(1e9) / P + 1 , i;
        for (i = 1 ; i < p && n; i ++) if (Miprime[i] >= P || i == 1) n --; i --;
        if (n) puts("0"); else printf("%d\n" , i * P);
    } else {
        cnt = 1;
        for ( ; prime[cnt] < P ; cnt ++); cnt --;
        int l = 1 , r = (int)(1e9) / P , ans = 0;
        while (l <= r)
        {
            int mid = l + r >> 1;
            if (Calc(mid) >= n) ans = mid , r = mid - 1 ; else l = mid + 1;
        }
        printf("%d\n" , ans * P);
    }

    return 0;
}