codeforces1295D Same GCDs

https://codeforces.com/contest/1295/problem/D

比赛的时候这题做了好久，没时间做E了。

据说答案是phi(m/gcd(a,m))，感觉是那么回事但是没想得很清楚，留下了不会数学的泪水.jpg，蒟蒻只会容斥，还纠结了半天才过。

我们令d=gcd(a,m)，在a到a+m-1之间的一个数字b要是gcd(b,m)=d，一定是d的倍数才行，所以说最多(a+m-1)/d-(a-1)/d个，但是这中间我们不能让gcd(b,m)>d，那么gcd(b/d,m/d)=1，那么我们需要从之前的总量中减去gcd(b/d,m/d)>1的，那么我们可以把m/d的所有质因子求出来，然后用dfs去容斥，一个质因子的减去，两个质因子之积的加上，就这样子去容斥就行了，由于m/d的所有因子最多sqrt(m/d)个，那么我们dfs最多也是sqrt(m/d)次，而直接求phi(m/d)也是sqrt(m/d)的复杂度。

#include<bits/stdc++.h>
using namespace std;
 
const int maxl=1e5+10;
 
long long a,m,d,ans,md,sz;
int mu[maxl],sum[maxl];
long long p[maxl];
bool no[maxl],in[maxl];
vector <long long> b;
 
inline void shai()
{
	int j,t;
	no[1]=true;mu[1]=1;
	for(int i=2;i<maxl;i++)
	{
		if(!no[i]) p[++p[0]]=i,mu[i]=-1;
		j=1;t=p[1]*i;
		while(t<maxl && j<=p[0])
		{
			no[t]=true;
			if(i%p[j]==0)
			{
				mu[t]=0;
				break;
			}
			mu[t]=-mu[i];
			t=p[++j]*i;
		}
	}
	for(int i=1;i<maxl;i++)
		sum[i]=sum[i-1]+mu[i];
}
 
inline void prework()
{
	scanf("%lld%lld",&a,&m);
	d=__gcd(a,m);md=m/d;
	b.clear();
}
 
inline void dfs(int k,int last,long long num)
{
	if(k&1) ans+=(a+m-1)/(num*d)-(a-1)/(num*d);
	else	ans-=(a+m-1)/(num*d)-(a-1)/(num*d);
	for(int i=last+1;i<=sz;i++)
	if(num*b[i-1]<=a+m-1)
		dfs(k+1,i,num*b[i-1]);
}
 
inline void mainwork()
{
	//ans=(a+m-1)/d-(a-1)/d;
	ans=0;
	long long ymd=md;
	for(int i=1;i<=p[0];i++)
	if(md%p[i]==0)
	{
		b.push_back(p[i]);
		while(md%p[i]==0)
			md/=p[i];
	}
	if(md>1)
		b.push_back(md);
	sz=b.size();
	dfs(1,0,1);
}
 
inline void print()
{
	printf("%lld\n",ans);
}
 
int main()
{
	int t;
	scanf("%d",&t);
	shai();
	for(int i=1;i<=t;i++)
	{
		prework();
		mainwork();
		print();
	}
	return 0;
}