1、题目描述
小易是一个数论爱好者,并且对于一个数的奇数约数十分感兴趣。一天小易遇到这样一个问题: 定义函数f(x)为x最大的奇数约数,x为正整数。 例如:f(44) = 11.
现在给出一个N,需要求出 f(1) + f(2) + f(3).......f(N)
例如: N = 7
f(1) + f(2) + f(3) + f(4) + f(5) + f(6) + f(7) = 1 + 1 + 3 + 1 + 5 + 3 + 7 = 21
小易计算这个问题遇到了困难,需要你来设计一个算法帮助他。
输入描述:
输入一个整数N (1 ≤ N ≤ 1000000000)
输出描述:
输出一个整数,即为f(1) + f(2) + f(3).......f(N)
首先一看,不难,动手编起来。
#include <iostream>
#include <math.h>
using namespace std;
int main()
{
int MaxOddNum(int x);
int N,sum=0;
cin>>N;
for (int j=1;j<=N;j++)
{
sum+=MaxOddNum(j);
}
cout<<sum<<endl;
return 0;
}
int MaxOddNum(int x)
{
if(0!=x%2)//除以2的余数不为0,说明被除数是奇数
return x;
int MaxOdd=1;
int a;
for (int i=1;i<=sqrt(double(x));i++)
{
if (x%i==0)//说明此时i是x的一个约数
{
a=x/i;//计算i的另一个约数
if((0!=i%2)&&(i>MaxOdd))
MaxOdd=i;
if((0!=a%2)&&(a>MaxOdd))//如果这个约数又是奇数
return a;
}
}
return MaxOdd;
}
却运行时间太长,循环太多。
修改如下
#include <iostream>
using namespace std;
int main()
{
long long N,sum=0;
cin>>N;
while (N)
{
if (1==N%2)
{
sum+=(1+N)*(1+N)/4;
N=N/2;
}
else
{
sum+=N*N/4;
N=N/2;
}
}
cout<<sum<<endl;
return 0;
}