Codeforces B. Different Divisors (Round #696 Div.2) (贪心&数学&质数)

传送门

题意：寻找一个最小正整数a，使得a至少有4个因子，且其任意两个因子之间的差至少为整数d。

思路：

• 刚开始会错了题意，其实最少4个因子且因子间差至少为d的情况下，其实往往因子都为质数。
• 对于第一个因子为1毋庸置疑，那么只需要找到另外两个质数a和b，满足它俩以及1之间的差距大于等于d即可，而答案就是a和b的最小公倍数。

代码实现：

#include<bits/stdc++.h>
#define endl '\n'
#define null NULL
#define ll long long
#define int long long
#define pii pair<int, int>
#define lowbit(x) (x &(-x))
#define ls(x) x<<1
#define rs(x) (x<<1+1)
#define me(ar) memset(ar, 0, sizeof ar)
#define mem(ar,num) memset(ar, num, sizeof ar)
#define rp(i, n) for(int i = 0, i < n; i ++)
#define rep(i, a, n) for(int i = a; i <= n; i ++)
#define pre(i, n, a) for(int i = n; i >= a; i --)
#define IOS ios::sync_with_stdio(0); cin.tie(0);cout.tie(0);
const int way[4][2] = {{1, 0}, {-1, 0}, {0, 1}, {0, -1}};
using namespace std;
const int  inf = 0x7fffffff;
const double PI = acos(-1.0);
const double eps = 1e-6;
const ll   mod = 1e9 + 7;
const int N = 1e6 + 5;

int t, d;
int pr[N], res;
bool st[N];

void get_primes(int n) //线性素数筛
{
    for(int i = 2; i <= n; i ++){
        if(!st[i]) pr[res ++ ] = i;
        for(int j = 0; pr[j] <= n / i; j ++){
            st[pr[j] * i] = 1;
            if(i % pr[j] == 0) break;
        }
    }
}

signed main()
{
    IOS;

    get_primes(N);
    cin >> t;
    while(t --){
        cin >> d;
        int cnt = 1, tmp = 1, a, b;
        for(int i = 0; i < res; i ++){
            if(pr[i]>=tmp+d){
                tmp = pr[i];
                cnt ++;
                if(cnt==2) a = pr[i];
                else{
                    b = pr[i];
                    break;
                }
            }
        }
        cout << a*b << endl;
    }

    return 0;
}