Same GCDs

2020-01-31 11:03:23 阅读：278 来源： 互联网

参考：

题意很明显要求出当\(k\in [a,a+m),gcd=gcd(a,m)\)时，满足\(gcd(k,m)=gcd\)的\(k\)的个数，由欧拉函数可以转换为\(gcd(k/gcd,m/gcd)=1,k\in [a,a+m)\)，\(a\)和\(m\)肯定是\(gcd\)的倍数，那么假设\(a=x*gcd,m=y*gcd\)，令\(i=k/gcd\)，那么\(i\)的取值范围为\([x,x+y)\)，那么最终要求的东西就是当\(i\in [x,x+y)\)时，满足\(gcd(i,y)=1\)的\(i\)的个数。

把\(i \in [x,x+y)\)，分成两个区域\([x,y]\)和\((y,x+y)\)

我们先讨论\(i\in(y,x+y)\)时，由欧几里得定理可以得出\(gcd(i,y)=gcd(y,i\%y)=gcd(i\%y,y)\)，再由\(i\in (y,x+y)\)可以化简所求东西为求\(i\in (0,x)\)时，满足\(gcd(i,y)\)的\(i\)的个数，而我们把这个区间与\([x,y]\)进行合并，就可以发现，我们要求的就是\(i\in [1,y]\)时满足\(gcd(i,y)=1\)的个数，即\(i\)与\(y\)互质，可以发现求的就是\(y\)的欧拉函数值。

代码：

// Created by CAD on 2020/1/30.
#include <bits/stdc++.h>
#define ll long long
using namespace std;

ll phi(ll x){
    ll ans=x;
    for(ll i=2;i*i<=x;++i)
        if(x%i==0){
            ans-=ans/i;
            while(x%i==0) x/=i;
        }
    if(x>1) ans-=ans/x;
    return ans;
}
int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    int t;  cin>>t;
    while(t--){
        ll a,m;
        cin>>a>>m;
        cout<<phi(m/__gcd(a,m))<<'\n';
    }
    return 0;
}

标签：GCDs,ll,cin,Same,ans,gcd
来源： https://www.cnblogs.com/CADCADCAD/p/12244773.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9

Same GCDs