CF1444A Division

题面传送门

题意简述:给出$a$与$b$，找到一个最大的数$x$，使得$a\%x==0\&\&x\%b\ne 0$

这样的题目考虑质因数分解。

因为$b$分解是在可以接受的范围内，所以分解$b$。

设$b=a_1^{p_1}\times a_2^{p_2}\times ...\times a_n^{p_n}$

然后将$a$分解得$a_1^{k_1}\times a_2^{k_2}\times ...\times a_n^{k_n}\times last$

使得分解后的$x\%b\ne 0$的条件是任意一个$k_i<p_i$

那么从这之中取最小值即可。

代码实现:

#include<cstdio>
#include<cmath>
#include<map>
#define min(a,b) ((a)<(b)?(a):(b))
using namespace std;
int n,m,k,z,f[100039],s[100039],d[100039],head,dh,t,sh;
map<long long,int> q;
long long x,y,now,tot,pus,ans,flag;
int main(){
//	freopen("1.in","r",stdin);
	register int i,j;
	for(i=2; i*i<=1e9; i++) {
		if(!f[i]) s[++head]=i;
		for(j=1; j<=head&&i*s[j]<=100000; j++) {
			f[i*s[j]]=1;
			if(i%s[j]==0) break;
		}
	}
	scanf("%d",&t);
	while(t--){
		scanf("%lld%lld",&x,&y);now=y;
		dh=sh=0;pus=1e18;
		for(i=1;s[i]*s[i]<=y&&i<=head;i++){
			if(now%s[i]==0){
				d[++dh]=s[i];
				while(now%s[i]==0)q[s[i]]++,now/=s[i],sh++;
			}
		} 
		if(now!=1)q[now]++,d[++dh]=now,sh++;
		now=x;
		for(i=1;i<=dh;i++){
			if(x%d[i]==0){
				ans=0;
				while(x%d[i]==0) x/=d[i],ans++;
				if(ans>=q[d[i]]){
					tot=1;
					for(j=1;j<=ans-q[d[i]]+1;j++) tot*=d[i];
					pus=min(pus,tot);
				}
				sh-=min(q[d[i]],ans);
			}
		}
		if(sh)printf("%lld\n",now);
		else printf("%lld\n",now/pus);
		q.clear();
	}
}