Sigma Function LightOJ 1336（唯一分解定理求因数和，思维，数论）

Sigma Function LightOJ 1336

题意：

求1到n有多少数的因数和是偶数。

思路：

我们的目的是求出因子和为偶数的数的个数。

他给了2个公式

唯一分解：

因数求和：

可以看出来是等比数列求和之后和的乘积。

所以其实也可以搞成这样

F(n)=(1+p1+p1^2+……+p1^e1)*(1+p2+p2^2+……+p2^e2)*……*(1+pi+pi^2+……+pi^ei)

那么什么时候因数和会为偶数呢，只要其中一项为偶数就可以了，但是这样想无法解题。

我们需要判断1到n的数中有多少个数的因数和为偶数，想直接判断偶数的方法解题绝对超时爆空间。

所以考虑反面，我们不如求出1到n其中因数和为奇数的个数，再用总数减去即可。

那么什么时候因数和F(n)会是奇数呢？

很简单，每一项都是奇数就可以了。

那么什么时候 其中一项1+pi+pi^2+……+pi^ei 会等于奇数呢？

如果pi不等于2，那么pi绝对是奇数,既然pi是奇数，那么pi的任意非负次方都是奇数，那么这个算式的项数只要是奇数，那么结果就是奇数了。因为有奇数项，那么此时ei必定为偶数。(素数中只有2是偶数)

如果pi等于2，那么显然结果为奇数（因为对于这个算式，除了1之外的每一项都是2的倍数，即为偶数，偶数加1肯定是奇数啊）

现在我们就有了这个结论：在 1+pi+pi^2+……+pi^ei 为奇数时，如果pi不等于2，此时ei必定为偶数。

那么这个结论有什么用呢？往下看

因为n = p1^e1*p2^e2*……*pi^ei（唯一分解）

假设任意pi都不等于2，那么任意的ei都是偶数，那么存在x使n=x*x，那么n是一个完全平方数。

假设存在一个px=2，那么当px的指数ex为偶数时，依旧存在一个x使n=x*x,那么n依旧是一个完全平方数。

                                     那么当px的指数ex为奇数是，存在一个x是n=2*x*x,那么n/2是一个完全平方数。

此时可以发现，只要n或n/2是一个完全平方数，那么n的因子和为奇数。

现在问题可以转化成求1到n之中有多少个数是完全平方数或者是完全平方数的二倍。

代码：

很简单，随便说一种，设存在完全平方数m=x*x,所以，1到这个完全平方数m，就有x个完全平方数，那么对于任意一个数n你求一下sqrt(n),再对他向下取整即可。

所以 ans = n - floor (  sqrt( n ) ) - floor (  sqrt( n/2 ) )

代码：

#include <bits/stdc++.h>
using namespace std;

int main()
{
    int T;
    cin>>T;
    for(int i=1;i<=T;i++)
    {
        long long n;
        cin>>n;
        long long ans = n - (long long)sqrt(n) - (long long)sqrt(n/2);
        cout<<"Case "<<i<<": "<<ans<<endl;
    }
    return 0;
}