LightOJ1336 Sigma Function
标签
- 约数相关问题
前言
- 我的csdn和博客园是同步的,欢迎来访danzh-博客园~
简明题意
- 给定n(n<=1e12),求
∑ i = 1 n [ σ ( i ) % 2 = 0 ] \sum_{i=1}^n[\sigma(i)\%2=0] i=1∑n[σ(i)%2=0]
σ \sigma σ是约数和的意思,上面的式子用文字形式表示就是:1-n中约数和是偶数的个数。
思路
-
首先看 σ \sigma σ函数的计算式:
∏ p ( 1 + p 1 + p 2 + . . . + p 1 k ) \prod_{p}(1+p^1+p^2+...+p^k_1) p∏(1+p1+p2+...+p1k) -
想让这个累乘式是偶数,那么显然这个是一点也不显然的。因为奇数 * 奇数=奇数,奇 * 偶=偶,偶数 * 偶数=偶数。所以实际上奇数的更好求一些,所以我们转化为n-约数和为奇数的个数。
-
累乘式的结果为奇数,那么每一个被乘式都应该是奇数。所以现在考虑每一个被乘式 1 + p 1 + p 2 + . . . + p k 1 1+p^1+p^2+...+p^{k_1} 1+p1+p2+...+pk1.大家可以自行验证一下,想让这个式子是奇数,必须满足下面两条中的一条:
- p是偶数
- p是奇数且k是偶数
-
我们先只考虑第二条。比如 n = 3 2 ∗ 5 4 ∗ 7 2 n=3^2*5^4*7^2 n=32∗54∗72,这样的话是满足第二条的,所有的p是奇数,且所有的k是偶数。这样的话,是不是发现此时n是一个完全平方数呢?所以对于第二条,我们统计一下1-n中范围内完全平方数的个数就可以了。
-
现在再考虑第一条。p为偶数,其实只有2这一种情况。这样的n也是符合的:
n = 2 3 ∗ 3 2 ∗ 5 3 ∗ 7 2 n=2^3*3^2*5^3*7^2 n=23∗32∗53∗72 n = 2 4 ∗ 3 2 ∗ 5 3 ∗ 7 2 n=2^4*3^2*5^3*7^2 n=24∗32∗53∗72 -
第二个仍然是一个完全平方数,我们统计完全平方数的时候就把2的指数是偶数的情况考虑进去了。现在只有2的指数是奇数时没有考虑。这种情况也很简单,对于一个完全平方数,他的两倍,如果<=n,就可以了。
-
所以最后我们可以暴力统计,也可以O(1)计算。代码贴的是暴力统计的。而O(1)的算法,是直接输出 n + n 2 \sqrt{n}+\sqrt{\frac n2} n+2n。
注意事项
- 无
总结
- 奇数*奇数=奇数
- 奇数*偶数=偶数
- 偶数*偶数=偶数
- 只有奇数*奇数是奇数,其他相乘全部是偶数
- 累乘式的奇偶性可以通过上式判断
- n以内完全平方数的个数: n \sqrt{n} n
- n以内除以k仍然是完全平方数的个数: n k \sqrt{\frac nk} kn
AC代码
#include<cstdio>
void solve()
{
int t;
scanf("%d", &t);
for (int i = 1; i <= t; i++)
{
long long n;
scanf("%lld", &n);
int ans = 0;
for (int i = 1; 1ll * i * i <= n; i++)
{
ans++;
if (2ll * i * i <= n)
ans++;
}
printf("Case %d: %lld\n", i, n - ans);
}
}
int main()
{
solve();
return 0;
}