hdu-5974

最新推荐文章于 2020-07-09 23:06:09 发布

虐猫Da人薛定谔

最新推荐文章于 2020-07-09 23:06:09 发布

阅读量168

点赞数

分类专栏：数论

原文链接：https://www.cnblogs.com/kimsimple/p/6792395.html

版权

数论专栏收录该内容

4 篇文章 0 订阅

订阅专栏

转载原博客：https://www.cnblogs.com/czsharecode/p/9595577.html

题目描述：

Given two positive integers a and b,find suitable X and Y to meet the conditions:
X+Y=a
Least Common Multiple (X, Y) =b

Input

Input includes multiple sets of test data.Each test data occupies one line,including two positive integers a(1≤a≤2*10^4),b(1≤b≤10^9),and their meanings are shown in the description.Contains most of the 12W test cases.Output For each set of input data,output a line of two integers,representing X, Y.If you cannot find such X and Y,output one line of "No Solution"(without quotation).

Sample Input

6 8

798 10780

Sample Output

No Solution

308 490

题目大意：给定正整数a,b;求两个正整数 x,y，使得 x + y == a && LCM(x,y) == b，如果找不到则输出No solution.

题解：由于test case 和 a,b规模都很大，不能使用暴力，必然是通过数学方法直接求解。

不妨设x = ki, y = kj; gcd(x,y) = k

易知 i,j互质（如果不互质则gcd必然大于k）

gcd(a,b) = gcd( k*(i+j) , k*(i*j) )

由于i，j互质，则（i+j）和 (i*j)必然互质，证明如下：

对于i的任意因子p（1除外），i % p = 0, (i*j) % p = 0

(i+j) % p = (i%p + j%p) % p = j%p, 由于i,j互质则p必然不是j的因子，所以 p 不是 (i+j) 的因子

所以对于i的所有因子（1除外）i+j都没有，但i*j都有；同理对于j的所有因子（1除外），i+j也没有，但i*j都有

所以i*j的所有因子（1除外），i+j都没有即 (i+j) , (i*j) 互质

我们可以得出以下结论：

（1）如果 i,j互质，那么i 和(i+j) 互质，j和（i+j）互质

（2）如果 i,j互质，那么(i+j) 和（i*j）互质

对于此题我们推出了gcd(a,b) = gcd(x,y) = k

原方程：LCM（x,y） = x*y / gcd(x,y) = b xy = bk = b*gcd(a,b)

又有x + y = a ， a,b已知

可以把y表示成x带入解一元二次方程；

也可以用(x-y)² = (x + y)² - 4xy求出x - y进而求出x和y

#include<bits/stdc++.h>
#define ll long long

using namespace std;

int main(){
	int a,b;
	while(~scanf("%d%d",&a,&b)){
		ll c=__gcd(a,b);
		ll d=a*a-4*c*b;
		if(d<0){
			printf("No Solution\n");
		}else{
			ll i=(a-sqrt(d))/(2*c);
			ll j=a/c-i;
			ll x=i*c;
			ll y=j*c;
			if(x/__gcd(a,b)*y!=b){
				printf("No Solution\n");
			}else{
				printf("%lld %lld\n",x,y);
			}
		}
	}
	return 0;
}