Codeforces Round #680 C-Division（数论+分解质因子）

题意

给你两个数字 p , q ,让你求满足（ p % x == 0 && x % q != 0 ）的最大的 x

思路

这里我们用n,m来代替p,q
如果 n % m != 0 那么最大的 x 就是 n
如果 n % m == 0 ，
我们知道
假设

如果 n % x == 0
t1 <= e1 , t2 <= e2 … tk <= ek;
又因为 n % m == 0
f1 <= e1 , f2 <= e2 … fk <= ek
如果想让 x % m != 0
我们可以让 x 对应的质因子的幂次 ti 变为 fi - 1 ( 即 ti = fi - 1 ) , 因为我们是求最大的 x ，所以我们每次只改变一个质因子的幂次即可，其他质因子的幂次与 n 相同 , 不断更新答案，最后得到的答案就是最大的 x
（最后感谢刘爷抬一手，要不然又木得了，/(ㄒoㄒ)/~~）

代码

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#define INF 0x3f3f3f3f
using namespace std;
typedef long long ll;
const int maxn=1e7+10;
const int mod=10007;
ll qpow(ll a,ll b)
{
	ll res=1;
	while(b)
	{
		if(b&1) res=res*a;
		b>>=1;
		a=a*a; 
	}
	return res;
}
int main()
{
	ios::sync_with_stdio(false);
	ll t,n,m; 
	cin>>t;
	while(t--)
	{
		cin>>n>>m;
		if(n%m) cout<<n<<endl;
		else
		{
			ll ans=0,sum1=0,sum2=0;
			for(int i=2;i*i<=m;i++)
			{
				if(m%i==0)
				{
					if(n%i==0)
					{
						sum1=sum2=0;
						while(m%i==0)
						{
							sum2++;
							m/=i;
						}
						ll nn=n;
						while(nn%i==0)
						{
							sum1++;
							nn/=i;
						}
						ans=max(ans,n/qpow(i,sum1-sum2+1));
					}
				}
			} 
			if(m>1)
			{
				ll mm=m;
				sum1=sum2=0;
				while(mm%m==0)
				{
				    sum2++;
					mm/=m;
				}
				ll nn=n;
				while(nn%m==0)
				{
					sum1++;
					nn/=m;
				}
				ans=max(ans,n/qpow(m,sum1-sum2+1));
			}
			cout<<ans<<endl;
		}
	} 
	return 0;
}