码迷,mamicode.com
首页 > 其他好文 > 详细

HDU2588

时间:2018-01-10 16:52:42      阅读:137      评论:0      收藏:0      [点我收藏+]

标签:clu   约数   turn   bit   its   namespace   body   down   pos   

结论:如果p是n的约数,那么满足gcd(i,n)==p的i的个数是Φ(n/p)

#include<bits/stdc++.h>
using namespace std;
const int maxn = 1e5+11;
typedef long long ll;
ll euler(int n){
    ll ans=n;
    for(ll i = 2; i*i <= n; i++){
        if(n%i==0){
            ans=ans/i*(i-1);
            while(n%i==0) n/=i;
        }
    }
    if(n>1) ans=ans/n*(n-1);//note
    return ans;
}
int main(){
    int T; scanf("%d",&T);
    while(T--){
        int n,m; ll ans=0;
        scanf("%d%d",&n,&m);
        for(ll i = 1; i*i<=n; i++){
            if(n%i==0){
                if(i>=m) ans+=euler(n/i);
                if(n/i>=m&&i*i!=n) ans+=euler(i);
            }
        }
        printf("%lld\n",ans);
    }
}

HDU2588

标签:clu   约数   turn   bit   its   namespace   body   down   pos   

原文地址:https://www.cnblogs.com/caturra/p/8259101.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!