Problem:
幾何問題を解こう

A君は今日も幾何の問題を解いている。
幾何の問題を解く時は浮動小数点誤差に気をつけることが大事である。

浮動小数点誤差とは、2進法の有限小数で数を表す際におこる丸めによって起きる誤差である。
例えば、10進法での 0.1 は2進法では 0.00011001100110011 ... という無限小数になるが、
これを有限の桁で丸める際に誤差が発生してしまう。

正の整数 p, q が10進法で与えられる。
有理数 p / q を有限桁数の小数で表現することができるような b 進法（b は2以上の整数）を求めよ。
複数ある場合は最も小さいものを出力せよ。

Constraints

0 < p < q < 10^9

Input Format

入力は以下の形式で標準入力から与えられる。

p q

Output Format

答えを一行に出力せよ。

Sample Input 1

1 2

Sample Output 1

2

1/2 は 2 進法で 0.1 です

Sample Input 2

21 30

Sample Output 2

10

21/30 は 10 進法で 0.7 です

Solution in C++ (judge verdict: Wrong Answer)

#include <iostream>
#include<cstdlib>
#include<queue>
#include<set>
#include<vector>
#include<string>
#include<algorithm>
#include<sstream>
#include<cmath>
#include<stack>
#include<map>
#include<cstdio>
using namespace std;
#define rep(i,a) for(int i=0;i<a;i++)
#define mp make_pair
#define pb push_back
#define P pair<int,int>
#define ll long long
int p,q;
vector<int>primes;

const int MA=30000000;
bool t[MA*2];

void eratosu(){
	t[1]=1;

	for(int i=2;i<=MA;i++){
		if(t[i]==0){
			primes.push_back(i);
			for(int j=i*2;j<=MA;j+=i)t[j]=1;
		}
	}
}

bool check(int x){
	int a=x;
	int aa=x;

	while(aa<=q){
		if(aa==q)return 1;
		aa=aa*a;
	}
	return 0;
}
int gcd(int a,int b){
	if(b==0)return a;
	return gcd(b,a%b);
}
int main(){
eratosu();
	cin>>p>>q;

/*
for(int i=0;i<primes.size();i++){
	int now=primes[i];
	while(q%now==0&&p%now==0)q/=now,p/=now;
}*/


int gc=gcd(p,q);
int np=p,nq=q;
p=np/gc;
q=nq/gc;

ll ans=1,qq=q;
int cnt=0;
for(int i=0;i<primes.size();i++){
	int now=primes[i];

	if(q%now==0)ans*=now,cnt++;
	while(q%now==0)q/=now;

	//if(now>sqrt(qq))break;
}
if(cnt==0)ans*=q;
cout<<ans<<endl;
	return 0;
}