Educational Codeforces Round 81 (Rated for Div. 2) D. Same GCDs 题解（欧拉函数）

题目链接

题目大意

给两个整数$a,m(1\le a<m\le 10^{10})$有多少个$x$满足 $0\le x<m,gcd(a,m)=gcd(a+x,m),0\le x<m$

题目思路

emmm，看到这个题目完全没思路啊，但是仔细分析就能发现解法，首先观察数据$0<=x<m$ x的范围正好是m，
gcd(a+x,m)=gcd((a+x)%m,m)则0<=（a+x）%m<m,，然后设p=gcd(a,m)
那么题目变转化成${\sum }_{i=0}^{m-1}gcd(i,m)==p$ 然后同时除以p 变成了 ${\sum }_{i=0}^{m-1}gcd(\frac{i}{p},\frac{m}{p})==1$
那么结果显然就是$\frac{m}{p}$的欧拉函数

代码

#include<set>
#include<map>
#include<stack>
#include<cmath>
#include<cstdio>
#include<vector>
#include<string>
#include<cstring>
#include<iostream>
#include<algorithm>
#define fi first
#define se second
#define debug printf("I am here\n");
using namespace std;
typedef long long ll;
const int maxn=1e5+5,mod=998244353;
int t;
ll a,m;
ll gcd(ll a,ll b){
    return b==0?a:gcd(b,a%b);
}
ll euler(ll x){
    ll ans=x;
    for(ll i=2;i*i<=x;i++){
        if(x%i==0){
            ans=ans/i*(i-1);
            while(x%i==0){
                x=x/i;
            }
        }
    }
    if(x!=1){
        ans=ans/x*(x-1);
    }
    return ans;
}
signed main(){
    scanf("%d",&t);
    while(t--){
        scanf("%lld%lld",&a,&m);
        printf("%lld\n",euler(m/gcd(a,m)));
    }
    return 0;
}