题目链接
这题的正解其实是令B=C/A,然后不断地通过G=gcd(A,B),A/=G,B*=G;这样的操作,在保证A*B=C的情况下将A和B的公因数从A转移至B,直到A和B互质。
然而我用的是一种性能更差的算法,我将A和C分别质因数分解,对于某个质因数,如果他在A中的幂小于C中的幂,则B中该质因数的幂等于C中该质因数的幂,否则B中该质因数的幂等于1。
附上AC代码:
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<cstdio>
#include<queue>
#include<set>
#include<vector>
#include<map>
#include<string>
#include<cmath>
#define pq priority_queue
#define Pi acos(-1.0)
#define MAXX 1000000007
using namespace std;
int gcd(int a,int b)
{
if(a%b)
return gcd(b,a%b);
else
return b;
}
int prime[666666];
int v[666666],l=0;
int a,b;
int ok1(int x)
{
int ans=0;
while(a%x==0)
{
ans++;
a=a/x;
}
return ans;
}
int ok2(int x)
{
int ans=0;
while(b%x==0)
{
ans++;
b=b/x;
}
return ans;
}
void go(int x,int t)
{
for(int i=0;i<t;i++)
{
b*=x;
}
}
int main()
{
for(int i=2;1LL*i*i<=10000000;i++)
{
if(!prime[i])
{
v[l++]=i;
for(int j=2;1LL*i*j*i*j<=10000000;j++)
prime[i*j]=1;
}
}
int t,l1,l2;
cin>>t;
while(t--)
{
scanf("%d%d",&a,&b);
if(b%a || b<a)
cout<<"NO SOLUTION"<<endl;
else
{
for(int i=0;i<l;i++)
{
if(a%v[i]==0)
{
l1=ok1(v[i]);
l2=ok2(v[i]);
if(l1<l2)
go(v[i],l2);
}
}
cout<<b/a<<endl;
}
}
return 0;
}
Memory: 0 KB Time: 370 MS
Language: C++ 4.8.2 Result: Accepted