codeforces1152C. Neko does Maths

题目链接 琪亚娜世界第一可爱

给出a和b，求
$$\underset{k\ge 0}{\max }lcm(a+k,b+k)$$

---

首先，gcd与lcm有如下关系
$$lcm(a+k,b+k)=\frac{(a+k)\cdot (b+k)}{gcd(a+k,b+k)}$$
$$gcd(a+k,b+k)=gcd(a-b,a+k)$$
所以根据上述关系，可知a+k与b+k的公因数一定是a-b的因数，先假设这个公因数为a-b的因数q，那么a+k与b+k也需要满足 $$a+k\equiv b+k(\mathrm{m}\mathrm{o}\mathrm{d}q)$$ $$a+k\mathrm{m}\mathrm{o}\mathrm{d}q=0$$
也就是
$$a\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}q)$$
$$a-b\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}q)$$
这样就证明了存在一个k，能使得q为a+k和b+k的最大公因数。
当a mod q !=0 时k为
$$k=q-a\mathrm{m}\mathrm{o}\mathrm{d}q$$
否则k为0

---

#include <stdio.h>
#include <climits>
#include <cstring>
#include <time.h>
#include <math.h>
#include <iostream>
#include <algorithm>
#include <stack>
#include <queue>
#include <set>
#include <map>
#include <utility>
#include <vector>
#include <string>

#define INF 0x3f3f3f3f
#define ll long long
#define Pair pair<int,int>
#define re return

#define getLen(name,index) name[index].size()
#define mem(a,b) memset(a,b,sizeof(a))
#define Make(a,b) make_pair(a,b)
#define Push(num) push_back(num)
#define rep(index,star,finish) for(register int index=star;index<finish;index++)
#define drep(index,finish,star) for(register int index=finish;index>=star;index--)
using namespace std;

ll gcd(ll a,ll b)
{
	if(a<b)
		return gcd(b,a);
	if(b==0)
		return a;
	if(a%2)
	{
		if(b%2)
			return gcd(b,a-b);
		else
			return gcd(a,b>>1);
	}
	else
	{
		if(b%2)
			return gcd(a>>1,b);
		else
			return 2*gcd(a>>1,b>>1);
	}
}
ll lcm(ll a,ll b){
    re a*b/gcd(a,b);
}
int main(){
    ios::sync_with_stdio(false);
    cin.tie(NULL);

    ll a,b;
    cin>>a>>b;
    if(a<b)     //maxn=a
        swap(a,b);
    int k=0;
    ll dif=a-b;
    ll num=lcm(a,b);
    for(ll i=1;i*i<=dif;i++){
        if(dif%i)
            continue;
        ll add=i-a%i;
        ll n1=(a+add)/i,n2=(b+add)/i;
        ll temp=lcm(n1,n2)*i;
        if(temp<num){
            num=temp;
            k=add;
        }

        ll n=dif/i;
        add=n-(a%n);
        n1=(a+add)/n,n2=(b+add)/n;
        temp=lcm(n1,n2)*n;
        if(temp<num){
            num=temp;
            k=add;
        }
    }

    cout<<k<<endl;

    re 0;
}