[HDU2588]GCD（数论）

最新推荐文章于 2019-07-09 15:09:21 发布

Clove_unique

最新推荐文章于 2019-07-09 15:09:21 发布

阅读量751

点赞数

分类专栏：题解数论文章标签： HDU 数论

本文链接：https://blog.csdn.net/Clove_unique/article/details/51037102

版权

题解同时被 2 个专栏收录

1097 篇文章 4 订阅

订阅专栏

数论

19 篇文章 0 订阅

订阅专栏

题目描述

传送门

题解

根据题目化出一坨式子：
$\sum_{i=1}^N[(i,N)>=M]$
$=\sum_{i=1}^N\sum_{d=M}^N[(i,N)=d]$
$=\sum_{d=M}^N\sum_{i=1}^N[({i\over d},{N\over d})=1]$
令i=dj
$=\sum_{d=M}^N\sum_{j}^{N\over d}[({j},{N\over d})=1]$
$=\sum_{d=M}^N\phi({N\over d})$

『转化』
枚举d=M到N转化为枚举N的约数。

代码

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
#define LL long long

LL T,n,m,ans;

inline LL phi(LL x){
    LL ans=x;
    for (LL i=2;i*i<=x;++i)
      if (x%i==0){
        ans=ans*(i-1)/i;
        while (x%i==0) x/=i;
      }
    if (x>1) ans=ans*(x-1)/x;
    return ans;
}

int main(){
    scanf("%lld",&T);
    while (T--){
        scanf("%lld%lld",&n,&m);
        ans=0;
        for (LL i=1;i*i<=n;++i)
        if (n%i==0){
            if (i>=m) ans+=phi(n/i);
            if (n/i>=m&&i!=n/i) ans+=phi(i);
        }
        printf("%lld\n",ans);
    }
}