BZOJ 2705 [SDOI2012]Longge的问题题解

最新推荐文章于 2020-07-11 16:29:12 发布

csyzcyj

最新推荐文章于 2020-07-11 16:29:12 发布

阅读量1k

点赞数

分类专栏： BZOJ(大视野测评)

本文链接：https://blog.csdn.net/csyzcyj/article/details/17332237

版权

BZOJ(大视野测评) 专栏收录该内容

19 篇文章 4 订阅

订阅专栏

2705: [SDOI2012]Longge的问题

Time Limit: 3 Sec Memory Limit: 128 MB
Submit: 628 Solved: 402
[ Submit][ Status]

Description

Longge的数学成绩非常好，并且他非常乐于挑战高难度的数学问题。现在问题来了：给定一个整数N，你需要求出∑gcd(i, N)(1<=i <=N)。

Input

一个整数，为N。

Output

一个整数，为所求的答案。

Sample Input

Sample Output

HINT

【数据范围】

对于60%的数据，0<N<=2^16。

对于100%的数据，0<N<=2^32。

Source

round1 day1

【分析】：

求出∑gcd(i, N)(1<=i <=N)，用欧拉函数，可以参考

http://blog.csdn.net/csyzcyj/article/details/17332211 < BZOJ 2190 [SDOI2008]仪仗队题解 >

【代码】：

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<algorithm>
#include<iostream>
using namespace std;
#define MAX 4294967297
int N;
long long phi(int x)
{
      long long ans=x;
      for(int i=2;i*i<=x;i++)
      {
            if(!(x%i))
            {
                  ans=ans/i*(i-1);   
                  while(x%i==0){x/=i;}
            }
      }
      if(x>1) ans=ans/x*(x-1);  
      return ans;
}
int main()
{
      //freopen("input.in","r",stdin);
	  //freopen("output.out","w",stdout); 
      scanf("%lld",&N);
      long long ans=0;
      for(int i=1;i*i<=N;i++)
      {
            if(!(N%i))
            {
                ans+=i*phi(N/i);
                if(i*i!=N)
                     ans+=(N/i)*phi(N/(N/i));
            }
      }
      printf("%lld\n",ans);
	  //system("pause");
      return 0;
}

转载注明出处： http://blog.csdn.net/u011400953