Benefit
Recently
Yaghoub
is playing a new trick to sell some more. When somebody gives him
A
Tomans, he who never has appropriate changes, asks for
B
Tomans such that lowest common multiple of
A
and
B
equals to
C
and he will pay back a round bill. Or otherwise take some snack instead of the remaining of his money. He believes that finding such a number is hard enough that dissuades students from paying that.
Benefit |
You should write a program that help poor students giving the appropriate amount of money to Yaghoub. Of course if there are several answers you go for students' benefit which is the lowest of them.
Input
The first line begin with an integer T ( T![$ \le$](http://uva.onlinejudge.org/external/118/11889img1.png)
![$ \le$](http://uva.onlinejudge.org/external/118/11889img1.png)
![$ \le$](http://uva.onlinejudge.org/external/118/11889img1.png)
Output
Print the lowest integer B such that LCM(A, B) = C in a single line. If no such integer exists, print " NO SOLUTION " instead. (Quotes for clarity)Sample Input
3 2 6 32 1760 7 16
Sample Output
3 55 NO SOLUTION
题意:
输入两个整数A和C(1<=A,C<= 10^7),求最小的整数B使得LCM(A,B)=C,若无解输出NO SOLUTION
有 A*B = LCM(A,B) * GCD(A,B)
设B' = LCM(A,B) / A
则B = B' * GCD(A,B)
一种方法是枚举B'的倍数,一直到LCM(A,B)
二是不断把LCM(A,B')吸收到B'中去
#include <cstdio>
#include <iostream>
#include <vector>
#include <algorithm>
#include <cstring>
#include <string>
#include <map>
#include <cmath>
#include <queue>
#include <set>
using namespace std;
//#define WIN
#ifdef WIN
typedef __int64 LL;
#define iform "%I64d"
#define oform "%I64d\n"
#define oform1 "%I64d"
#else
typedef long long LL;
#define iform "%lld"
#define oform "%lld\n"
#define oform1 "%lld"
#endif
#define S64I(a) scanf(iform, &(a))
#define P64I(a) printf(oform, (a))
#define P64I1(a) printf(oform1, (a))
#define REP(i, n) for(int (i)=0; (i)<n; (i)++)
#define REP1(i, n) for(int (i)=1; (i)<=(n); (i)++)
#define FOR(i, s, t) for(int (i)=(s); (i)<=(t); (i)++)
const int INF = 0x3f3f3f3f;
const double eps = 10e-9;
const double PI = (4.0*atan(1.0));
LL gcd(LL a, LL b) {
return b ? gcd(b, a%b) : a;
}
int main() {
int T;
scanf("%d", &T);
while(T--) {
int a, c;
scanf("%d%d", &a, &c);
if(c % a) { puts("NO SOLUTION"); continue; }
int b = c / a;
int g = gcd(a, b);
while(g != 1) {
b *= g;
a /= g;
g = gcd(a, b);
}
printf("%d\n", b);
}
return 0;
}