东软杯-进制转换（求质因子的乘积）

题目

这是一道简单的进制转换题，给出一个十进制的分数，求出能将它化为有限小数形式的最小进制数。 为了更好的理解，我们给出一些例子：
$\frac{1}{3}$ 在十进制下为0.3333…，在三进制下为0.1 （题目此处写错了）
$\frac{1}{2}$ ​ 在十进制下为0.5，但是在二进制下为0.1
$\frac{7}{8}$ ​ 在十进制下为0.875，但是在二进制下为0.111

数据范围及复杂度分析

$T$组数据，分子分母不超过$10^9$，使用$O(T\sqrt{n})$(n为分母的最大值)的算法即可通过本题。

思路

一个k进制数可以表示成这个形式$\sum bas{e}_i\ast k^i(0\le bas{e}_i<k)我们把bas{e}_i称为基数，k^i称为位权，i称为位$
例如10进制下的996.12可以表示为$9\ast 10^2+9\ast 10^1+6\ast 10^0+1\ast 10^{-1}+2\ast 10^{-2}$
再比如2进制下的11.101可以表示为$1\ast 2^1+1\ast 2^0+1\ast 2^{-1}+1\ast 2^{-2}$

那么我们发现，若一个分数可以表示成一个k进制的小数，那么我们将其通分之后，其分母必定是$k^{|i|}，i$为其最小位，例如11.101最小位为-3。我再用$\frac{7}{8}$这个例子来解释一下，8为$2^3$,我们求出7的二进制$111$后，再乘$2^{-3}$即可，此时的乘相当于小数点向左移动三位。由于分子是一个整数，整数可以用任意k进制表示成整数，当分母向左移动时一定是变成有限小数的，此时就达到了题目的要求。

我们再来凝练一下做法，我们先对分子分母化简，再将分母表示成最小的$k^i$，也就是对于一个整数$n$，你要求出最小的$k$，使得存在一个$k^i$能够整除$n$，此时的$k$即为答案。那么应该很容易想到，最小的$k$就是$n$的质因子的乘积。我们可以很容易证明这个结论，若$k$中不包含$n$的某一个质因子，那么$k^i$不可能整除$n$，若$k$包含了非$n$的质因子，那么多的这个质因子完全可以去掉。

那么最后我们的难点就落在了求质因子的乘积上，我们可以用求欧拉函数的方式来求得。思路也很简单，看代码或者再去百度一下欧拉函数即可，此处不再说明。
不用怀疑，有一进制的存在。
直接看euler函数部分即可
标程丢失，这份代码可能有问题，但是思路没有问题

#include<bits/stdc++.h>
using namespace std;
using namespace std;
typedef long long LL;
const int INF = 0x3f3f3f3f;
const double Pi = acos(-1);
namespace io{
  template <typename T> inline void read(T &x) {
    x = 0; T f = 1;char s = getchar();
    for(; !isdigit(s); s = getchar()) if(s == '-') f = -1;
    for(;  isdigit(s); s = getchar()) x = (x << 3) + (x << 1) + (s ^ 48);
    x *= f;
  }
  template <typename T> inline void write(T x) {
    if(x < 0) putchar('-'), x = -x;
    if(x > 9) write(x/10);
    putchar(x%10+'0');
  }
};
#define fio ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#define _for(n,m,i) for (int i = (n); i <  (m); ++i)
#define _rep(n,m,i) for (int i = (n); i <= (m); ++i)
#define _srep(n,m,i)for (int i = (n); i >= (m); i--)
#define _sfor(n,m,i)for (int i = (n); i >  (m); i--)
#define ef(u, i) for(int i = head[u]; i; i = e[i].next)
#define lson rt << 1, l, mid
#define rson rt << 1 | 1, mid + 1, r
#define lowbit(x) x & (-x)
#define pii pair<int,int>
#define MS0(X) memset((X), 0, sizeof((X)))
#define MS1(X) memset((X), -1, sizeof((X)))
#define ALL(X) (X).begin(), (X).end()
#define SORT_UNIQUE(c) (sort(c.begin(),c.end()), c.resize(distance(c.begin(),unique(c.begin(),c.end()))))
#define GET_POS(c,x) (lower_bound(c.begin(),c.end(),x)-c.begin())
#define CASET int ___T; scanf("%d", &___T); for(int cs=1;cs<=___T;cs++)
#define fi first
#define se second
#define dbg(x...) do { cout << "|" << #x <<" -> "; err(x); } while (0)
void err() { cout << "|" << endl; }
template<class T, class... Ts> void err(const T& arg,const Ts&... args) { cout << arg << " "; err(args...); }
template<class T> void _R(T &x) { cin >> x; }
void _R(int &x) { scanf("%d", &x); }
void _R(LL &x) { scanf("%lld", &x); }
void _R(double &x) { scanf("%lf", &x); }
void _R(char &x) { scanf(" %c", &x); }
void _R(char *x) { scanf("%s", x); }
void R() {}
template<class T, class... U> void R(T &head, U &... tail) { _R(head); R(tail...); }
template<class T> void _W(const T &x) { cout << x; }
void _W(const int &x) { printf("%d", x); }
void _W(const LL &x) { printf("%lld", x); }
void _W(const double &x) { printf("%.16f", x); }
void _W(const char &x) { putchar(x); }
void _W(const char *x) { printf("%s", x); }
template<class T> void _W(const vector<T> &x) { for (auto i = x.begin(); i != x.end(); _W(*i++)) if (i != x.cbegin()) putchar(' '); }
void W() {}
template<class T, class... U> void W(const T &head, const U &... tail) { _W(head); putchar(sizeof...(tail) ? ' ' : '\n'); W(tail...); }
int euler(int n) {
  int ans = 1;
  for(int i = 2; i * i <= n; ++i) {
    if(n % i == 0) {
      ans *= i; n /= i;
      while(n % i == 0) n /= i;
    }
  }
  if(n != 1) ans *= n;
  return ans;
}
void solve() {
  int n, m; R(n, m);
  int gcd = __gcd(n, m);
  m /= gcd;
  W(euler(m));
}
int main() {
//   #define LOCAL
  #ifdef LOCAL
  freopen("test.in", "r", stdin);
  freopen("test.out", "w", stdout);
  #endif
  CASET
  solve();
  return 0;
}