light oj 1336 （求约数和为偶数的数的个数）

最新推荐文章于 2023-03-21 09:17:41 发布

acm er

最新推荐文章于 2023-03-21 09:17:41 发布

阅读量450

点赞数 1

分类专栏：数论文章标签：约数和为奇数

本文链接：https://blog.csdn.net/weixin_43323647/article/details/97109780

版权

数论专栏收录该内容

6 篇文章 0 订阅

订阅专栏

https://vjudge.net/problem/LightOJ-1336
这个是题目链接，题意很简单，就是求1-n中约束和为偶数的数的个数，用为n减去奇数的个数即为所求，这题我做的时候用的打表看出来的一点规律，发现为奇数的全部都有两个特点，要不就是完全平方数，要不就是完全平方数的二倍，因而在O(1)复杂度内解决就是n - sqrt(n) - sqrt(n/2);

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
typedef long long ll;
ll poww(ll a,ll b)
{
    ll ret = 1;
    while(b)
    {
        if(b & 1)
        {
            ret = ret * a;
        }
        a = a * a;
        b /= 2;
    }
    return ret;
}
ll getsum(ll n)//求约数和代码。
{
    ll res = 1;
    for(ll i = 2;i * i <= n;i++)
    {
        ll k = 0;
        while(n % i == 0)
        {
            n /= i;
            k++;
        }
        res *= ((1 - poww(i,k + 1))/(1 - i));
    }
    if(n != 1) res *= (n + 1);
    return res;
}
int main()
{
    /*for(int i = 1;i <= 1000;i++)
    {
        if(getsum(i) % 2 == 1) cout << i << " " << getsum(i) << endl;//这是当时打表过程，，
    }*/
    int T;
    scanf("%d", &T);
    for(int kase = 1; kase <= T; kase++)
    {
        ll n, sum;
        scanf("%lld", &n);
        sum = n;
        sum -= (int)sqrt(n);
        sum -= (int)sqrt(n/2);
        printf("Case %d: %lld\n", kase, sum);
    }
    return 0;
}

具体证明就得参考网上的大佬证明过程了，这里总结了一下，大致理解了：
分析：

数x的因子和 f(x)= （1+p1+p1^2+p13+…+p1^{a1）*(1+p2+p2}2+…+p2^{a2)*…*(1+pn+pn}2+…+pn^an);

因为偶数乘偶数还是偶数，奇数乘奇数还是奇数，奇数乘偶数是偶数，所有必须让每个括号内都是奇数，然后减去约数和为奇数的个数就是答案了。

1.当x有素因子2的时候，2所对应的括号内的和肯定是一个奇数，因为偶数加1一定是奇数。

2.除了2以外，所有的素数都是奇数，要使x得其他素因子对应的括号内的和为奇数，就必须保证x有偶数个该素因子，即ai必须为偶数。

3.满足上面两个条件的数，就是一个平方数，也就是说约数和为奇数的数x，它必定是一个平方数，当然这个数x乘上2也是满足2*x的约数和为奇数的。

所以只要减去用n减去sqrt(n)和sqrt(n/2)就是答案了。也可以找n以内的平方数的个数，以及2*平方数不超过n的数的所有个数和，用n减完之后就是答案。