D. Lucky Chains.Educational Codeforces Round 139 (Rated for Div. 2)

Lucky Chains

题意：给定x和y，执行下面的代码，求cnt

while(gcd(x,y)==1) x++,y++,cnt++;

思路：考虑gcd的性质，gcd(x,y)=gcd(x,y-x).那么gcd(x+k,y+k)=gcd(x+k,y-x)，我们发现y-x是一个定值，所以我们可以枚举每一个y-x的质因子p，求(x+k)%p==0时的最小的k，也就是求p-x%p的最小值.筛质数时只需要筛到3200就可以，因为3200*3200>1e7

#include <bits/stdc++.h>
#define PII pair<int, int>
typedef long long ll;
const int N = 1e7 + 10;
const int inf = 0x3f3f3f3f;

using namespace std;
int x, y;
bool st[N];
int primes[N], cnt;
void get_primes(int n)
{
    for (int i = 2; i <= n; i++)
    {
        if (!st[i])
            primes[cnt++] = i;
        for (int j = 0; primes[j] <= n / i; j++)
        {
            st[primes[j] * i] = true;
            if (i % primes[j] == 0)
                break;
        }
    }
}

void solve()
{
    cin >> x >> y;
    if (x > y)
        swap(x, y);
    int d = y - x;
    if(gcd(x,y)!=1)
    {
        cout<<0<<'\n';
        return ;
    }
    int ans=inf;
    for(int i=0;i<cnt&&primes[i]<=sqrt(d);i++)
    {
        if(d%primes[i]==0)
        {
            while(d%primes[i]==0) d/=primes[i];
            ans=min(ans,primes[i]-x%primes[i]);
        }
    }
    if(d>1) ans=min(ans,d-x%d);

    if(ans==inf) ans=-1;
    cout<<ans<<'\n';
}
signed main()
{
    ios::sync_with_stdio(0),cin.tie(0);
    get_primes(3200);
    int _t = 1;
    cin >> _t;
    while (_t--)
        solve();
    system("pause");
    return 0;
}