C. Neko does Maths (數論+思維題)

link:http://codeforces.com/problemset/problem/1152/C

題目大意:題目要求是求一個k使得(a+k)*(b+k)/ __gcd(a+k,b+k)最小,且在最小的公倍數中找到最小的k.

方法:重點在於:__gcd(a+k,b+k),而__gcd(a+k,b+k)=__gcd(abs(a-b),a+k)   (這也可以是__gcd( abs(a-b) ,b+k)  ),而gcd( abs(a-b) ,a+k) 一定是abs(a-b) 的因子,找出因子後通過倒推求出k,具體如下:

當(a+k)%因子  !=  (b+k)%因子 時,即a%因子 != b%因子,表示這個因子不是gcd(a+k,b+k),排除。

當....................==.....................時,表示這個因子就是gcd(a+k,b+k),判斷一下題目要求,再求出k(=因子-a%因子)

紅字來由:

//http://codeforces.com/problemset/problem/1152/C
#include<bits/stdc++.h>

using namespace std;

typedef long long ll;

#define endl '\n'

ll a,b;

void solve(){
	ll i,j;
	ll d=abs(a-b);
	vector<ll> v;
	for(i=1;i*i<=d;++i){
		if(d%i==0){
			v.push_back(i);
			if(i*i!=d){
				v.push_back(d/i);
			}
		}
	}
	ll ansg=(1ll<<62),ans=0;
	sort(v.begin(),v.end());
	for(i=0;i<v.size();++i){
		int t=v[i];
		if(a%t!=b%t) continue;
		else if(a%t==0){
			ll g=a*b/__gcd(a,b);
			if(ansg>g){
				ansg=g;
				ans=0;
			}
		}
		else{
			ll g=(a+t-a%t)*(b+t-b%t)/__gcd((a+t-a%t),(b+t-b%t));
			if(g<ansg){
				ansg=g;
				ans=t-a%t;
			}
		}
	}
	cout<<ans<<endl;
}

int main(){
	
	cin>>a>>b;
	solve();
	
	return 0;
}

 

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章