【算法题】最大的奇约数

小易是一个数论爱好者，并且对于一个数的奇数约数十分感兴趣。一天小易遇到这样一个问题： 定义函数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)

输入例子:
7

输出例子:
21

---

链接：https://www.nowcoder.com/questionTerminal/49cb3d0b28954deca7565b8db92c5296
来源：牛客网

因为奇数的最大奇数约数就是自己
对于偶数我们只能一直除2直到得到一个奇数即为最大奇数约数

注意：数字较大，使用long long防止溢出

---

#include <iostream>
#include <vector>
#include <string>
#include <numeric>
#include<algorithm>
using namespace std;
#define debug 1

long long func(int N)
{
    long long result(0);
    for (auto i = 1; i <= N; ++i)
    {
        if ((i & 1) != 0)//i奇数
        {
            result += i;
            continue;
        }
        else//i偶数
        {
            int tmp = i/2;
            while ((tmp&1)==0)
            {
                tmp = tmp / 2;
            }
            result += tmp;
        }
    }
    return result;
}

int main()
{
    int N;
    if (debug)
    {
        N = 7;
    }
    else
    {
        cin >> N;
    }
    cout << func(N) << endl;

    return 0;
}

---

但是，上述代码会超时，时间复杂度还是过高

需要进一步改进：
例如： 1 2 3 4 5 6 7 8 9 10

即n=10

此时奇数有1 3 5 7 9 我们把这几个奇数相加
然后: n = n/2
得到第二轮序列序列 1 2 3 4 5 分别对应上次的2 4 6 8 10 五个偶数，这是我们再加1 3 5
依次类推

• 公式如下：

  • 当n为偶数，就有n/2个奇数，根据等差数列求和公式 即(（首项+末项）*项数)/2,我们知道n/2个奇数和为

    $$((1+n-1)*n/2)/2$$ $$= (n/2) * (n/2)$$ ,此时n为偶数，因此 $$(n/2) * (n/2) = ((n+1)/2) * ((n+1)/2)$$

  • 当n为奇数，有(n+1)/2个奇数，此时奇数和为

    $$((n+1)/2) * ((n+1)/2)$$
    因此两种情况可以用一个等式来总结

---

#include <iostream>
#include <vector>
#include <string>
#include <numeric>
#include<algorithm>
using namespace std;
#define debug 1

long long func(int N)
{
    long long result(0);
    long long tmp;
    for (auto i = N; i > 0; i = i / 2)
    {
        tmp = (i + 1) / 2;
        result += tmp*tmp;
    }
    return result;
}

int main()
{

    int N;
    cin >> N;
    cout << func(N) << endl;

    return 0;
}