codeforces1152C. Neko does Maths

最新推荐文章于 2024-05-04 16:55:46 发布

white_156

最新推荐文章于 2024-05-04 16:55:46 发布

阅读量493

点赞数 1

分类专栏： codeforces 数论文章标签： gcd 数论

本文链接：https://blog.csdn.net/white_156/article/details/89682827

版权

codeforces 同时被 2 个专栏收录

76 篇文章 0 订阅

订阅专栏

数论

13 篇文章 0 订阅

订阅专栏

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

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

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