Codeforces 1295D

Codeforces 1295D

题意:
gcd(a,m) == gcd(a+x,m) 0 <= x < m;
思路:
  设 a = sg, m = tg;
  gcd(a,m) ==gcd(sg+x,tg) == g
  所以令 x = k*g
  所以 gcd(s+k,t) == 1;
  因为 a <= a+x < a+m;
  (0,a] 和 (m,m+a) 与 m 互质 个数相同;
  则本题可化简为 gcd(k,t)==1 k<t, 的个数

#include<bits/stdc++.h>
using namespace std;
#define re register
typedef long long LL;
const int N = 1e5+10;
LL gcd(LL a,LL b)
{
 return b==0?a:gcd(b,a%b);
}
LL euler(LL n)
{
 LL ret = 1;
 for(re LL i=2;i*i<=n;i++)
 {
  if( n % i == 0 )
  {
   n /= i;
   ret *= i-1;
   while( n % i == 0 )
   {
    n /= i;
    ret *= i;
   }
  }
 }
 if( n > 1 )
  ret *= n-1;
 return ret;
}
int main()
{
 int T;
 scanf("%d",&T);
 while(T--)
 {
  LL a,m;
  scanf("%lld%lld",&a,&m);
  printf("%lld\n",euler(m/gcd(a,m)));
 }
 return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值