Codeforcrs 1295 D

传送门

题目大意

给定两个整数 a , m   ( 1 ≤ a < m ≤ 1 0 10 ) a,m\ (1\le a<m\le 10^{10}) a,m (1a<m1010)。计算满足条件 0 ≤ x < m 0 \le x < m 0x<m gcd ⁡ ( a , m ) = gcd ⁡ ( a + x , m ) \gcd(a,m)=\gcd(a+x,m) gcd(a,m)=gcd(a+x,m)的整数 x x x的个数。

思路

d = g c d ( a , m ) d=gcd(a,m) d=gcd(a,m),那么 g c d ( a d , m d ) = g c d ( a + x d , m d ) = 1 gcd(\frac{a}{d},\frac{m}{d})=gcd(\frac{a+x}{d},\frac{m}{d})=1 gcd(da,dm)=gcd(da+x,dm)=1
那么答案就是 φ ( m / d ) \varphi(m/d) φ(m/d)

代码

ll euler(ll n){
	ll res=n;
	for(ll i=2;i*i<=n;i++){
		if(n%i==0){//i一定是素数
			res=res/i*(i-1);//根据公式求得
			while(n%i==0) n/=i;
		}
	}
	if(n>1)
	res=res/n*(n-1);
	return res;

}
ll gcd(ll a,ll b){
	return !b?a:gcd(b,a%b);
}
int main(){
	int t;
	cin>>t;
	while(t--){
		ll a,m;
		scanf("%lld%lld",&a,&m);
		ll d=gcd(a,m);
		ll ans=euler(m/d);
		printf("%lld\n",ans);
	}
}
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值