51nod1040 最大公约数之和

题目

给出一个n，求1-n这n个数，同n的最大公约数的和。比如：n = 6
1,2,3,4,5,6 同6的最大公约数分别为1,2,3,2,1,6，加在一起 = 15

题解

举一个更大一点的例子，
如15：
1,2,4,7,8,11,13,14这几个数的最大公约数都是1，
3,6,9,12这几个数的最大公约数都是3，
5,10,15这几个数的最大公约数都是5。

从中就可以看出一种思路，
枚举最大公约数d，然后求出有多少个数与它的最大公约数为d。
关键是怎么求出这些数的个数，
就是求在1..n/d中与n/d互质的数的个数。
然后就可以直接用欧拉。

code

#include<queue>
#include<cstdio>
#include<iostream>
#include<algorithm>
#include <cstring>
#include <string.h>
#include <cmath>
#include <math.h>
#include <time.h>
#define ll long long
#define N 100003
#define M 103
#define db double
#define P putchar
#define G getchar
#define inf 998244353
using namespace std;
char ch;
void read(int &n)
{
    n=0;
    ch=G();
    while((ch<'0' || ch>'9') && ch!='-')ch=G();
    ll w=1;
    if(ch=='-')w=-1,ch=G();
    while('0'<=ch && ch<='9')n=(n<<3)+(n<<1)+ch-'0',ch=G();
    n*=w;
}

int max(int a,int b){return a>b?a:b;}
int min(int a,int b){return a<b?a:b;}
ll abs(ll x){return x<0?-x:x;}
ll sqr(ll x){return x*x;}
void write(ll x){if(x>9) write(x/10);P(x%10+'0');}

int n;
ll ans;

int eular(int x) 
{
    int s=x;
    for(int i=2;i*i<=x;i++)
        if(x%i==0)
        {
            s-=s/i;
            for(x/=i;x%i==0;x/=i);
        }
    if(x>1)s-=s/x;
    return s;
}

int main()
{
    //freopen(".in","r",stdin);
    //freopen(".out","w",stdout);

    read(n);

    for(int i=1;i*i<=n;i++)
    {
        if(n%i==0)
        {
            ans=ans+(ll)i*eular(n/i);
            if(i!=n/i)ans=ans+(ll)n/i*eular(i);
        }
        if(n==1)break;
    }
    write(ans);

    return 0;
}