HDU2588

题目:click
题意:在这里插入图片描述设g=gcd(X,N),X= g ∗ a g*a ga , N= g ∗ b g*b gb,由于N>=X,b>=a,可以直接去枚举约数g,由于gcd(a,b)=1,a<=b,直接欧拉函数求可得。

#include<cmath>
#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<cstdlib>
#include<istream>
#include<vector>
#include<stack>
#include<set>
#include<map>
#include<algorithm>
#include<queue>
#define inf 0x3f3f3f3f
#define llinf 0x3f3f3f3f3f3f3f3f
#define MAX_len 200005*4
using namespace std;
typedef long long ll;
typedef pair<int,int> PP;
const int mod=998244353;
const int MAXlen=1e5+10;
long double eps=1e-9;
ll euler(ll n)
{
     ll i,j,res=n,temp=n;
     for(i=2;i*i<=temp;i++)
    {
         if(temp%i==0)
         {
             res=res/i*(i-1);//先进行除法是为了防止中间数据爆范围
             while(temp%i==0)
                temp/=i;
         }
     }
     if(temp>1)
        res=res/temp*(temp-1);
     return res;
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        ll ans=0;
        ll n,i,m,j;
        scanf("%lld %lld",&n,&m);
        for(i=1;i*i<=n;i++)
        {
            if(i*i==n&&i>=m)
            {
                ans+=euler(i);
                continue;
            }
            if(i>=m&&n%i==0)
            {
                ans+=euler(n/i);
            }
            if(n%i==0&&n/i>=m)
            {
                ans+=euler(i);
            }
        }
        printf("%lld\n",ans);
    }
    return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值