Educational Codeforces Round 139 1766D- Lucky Chains(div2-思维/欧拉筛/GCD)

题目链接：https://codeforces.com/problemset/problem/1766/D

题意：

如果一个数对 $(x,y)$ 是幸运的，当且仅当 $gcd(x,y)=1$ ，一条链可以由以下规律的数对组成，$(x,y),(x+1,y+1),(x+2,y+2),\dots ,(x+k,y+k)$ ，如果说一条链是幸运的，当且仅当这条链上的每一个数对都是幸运的，每次给出一个数对，问这个以这个数对为一条链的起点，所能构成的幸运链的长度最长是多少？如果说这条链的长度是无穷的，输出 $-1$，否则输出链的长度，如果一开始就不幸运输出 $0$ 。

思路：

首先 $(1\le n\le 10^6)$ 如果 $O(n^2)$ 会超时。可以发现当 $y-x=1$ 答案就是 $-1$ 。然后要知道一个 $gcd$ 的一个结论，当 $y>x,gcd(x,y)=gcd(x,y-x)$ 于是问题转化为 $gcd(x+i,y-x)$ 下幸运链的长度。然后枚举 $y-x$ 的质因子，此时 $O(n\ast \sqrt{n})$ 仍无法通过题目，于是就需要用到筛法预处理质因子，之后枚举 $y-x$ 的最小质因子，看 $x$ 加几倍 $i$ 到达即为答案，此时 $O(n\ast \ln n)$。

代码：

#include <bits/setc++.h>
#define endl '\n'
using namespace std;
typedef long long ll;
const int N = 1e7+10;
const ll INF_ = 0x3f3f3f3f3f3f3f3f;
const int INF = 0x3f3f3f3f;
ll n, m, k, x, y, res = 0, sum = 0, ans = 0, ma = 0, mi = INF;
ll minp[N];
void euler(ll n)
{
    bitset<N> vis;
    vector<ll> primes;
    vis[0] = vis[1] = true;
    minp[1] = 1;
    for (ll i = 2; i <= n; i++) {
        if (!vis[i]) primes.push_back(i), minp[i] = i;
        for (int j = 0; j < primes.size() && i * primes[j] <= n; j++) {
            vis[i * primes[j]] = true;
            minp[i * primes[j]] = primes[j];
            if (i % primes[j] == 0) break;
        }
    }
}
void solve()
{
    cin >> n;
    euler(N);
    for (int i = 1; i <= n; i++) {
        cin >> x >> y;
        y -= x;
        if (__gcd(x, y) != 1) {
            cout << 0 << endl;
            continue;
        }
        if (y == 1) {
            cout << -1 << endl;
            continue;
        }
        ans = INF_;
        while (y > 1) {
            ll p = minp[y];
            ans = min(ans, ((-x) % p + p) % p);
            while (y % p == 0)
                y /= p;
        }
        cout << ans << endl;
    }
}
int main()
{
	ios::sync_with_stdio(0);cin.tie(0);
	// cout.tie(0);
//	ll _;cin>>_;
//	while(_--)
		solve();
	return 0;
}