题目
给出一个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;
}