LightOJ - 1336 Sigma Function

LightOJ - 1336 Sigma Function

传送门LightOJ - 1336

题意

让你求从一到n中因数和为偶数的数的个数。

题解

在vj上写题的大概都看的中文题面吧，但是换回英文题面就会发现这道题甚至给了你求一个数的因数和的公式。

$$\delta (n) = \prod_{i = 1}^s{p_i^{a_i+1} - 1 \over p_i- 1}$$
或者写成
$$\prod_{i = 1}^s{\sum_{j = 0}^{a_i} p_i^j}$$
这个公式这几天写数论题总看见他哈,我们来把它拆开看看。
$$\delta(n) = (1+ p_1 + p_1^2 +...+ p_1^{a_1})\times(1 + p_2 + p_2^2 + ...+ p_2^{a_2}) \times...\times(1+ p_j + p_j^2 +... + p_j^{a_j})$$
那我们很容易想到，当且进当每一项 $(1+ p_j + p_j^2 +... + p_j^{a_j})$都为奇数时， $\delta(n)$为奇数。
那么什么情况时 $(1+ p_j + p_j^2 +... + p_j^{a_j})$为奇数呢？

1. 如果素数为2，那么那一项的$(1+ p_j + p_j^2 +... + p_j^{a_j})$一定为奇数
2. 除了2以外，其他的素数都是奇数，设当前素数为$p_j$，那么当且仅当$a_j$为偶数的时候，$(1+ p_j + p_j^2 +... + >p_j^{a_j})$为奇数。当$n$唯一因数分解出的每一个素数的个数都为偶数时，这个数就是完全平方数

那就很好理解了，n以下的$\delta(n)$为奇数的数，只有完全平方数，或者完全平方数乘以2（ps：其他的2的几次方都已经被计算在完全平方数中了）。而小于n的完全平方数个数就是$floor(\sqrt n)$个。那么小于n的$\delta(n)$为奇数的数的个数就等于$floor(\sqrt n) + floor(\sqrt {n \over 2})$，那么本题答案就为：

$$n - (floor(\sqrt n) + floor(\sqrt {n \over 2}))$$

ac代码

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <queue>
#include <cstring>
#include <string>
#include <vector>
#include <stack>
#include <cmath>
#include <ctime>
#include <map>
#include <set>
#include <iomanip>
//#include <unordered_map>

#pragma comment(linker, "/STACK:102400000,102400000")
#define fir first
#define sec second
#define lson l, mid, rt << 1
#define rson mid + 1, r, rt << 1 | 1
#define clr(x) memset(x,0,sizeof(x))
#define cld(x) memset(x,-1,sizeof(x))
#define clx(x) memset(x,63,sizeof(x))
#define cln(x) memset(x,-64,sizeof(x))
#define rush() int T;scanf("%d",&T);for(int NUM = 1; NUM <= T ; ++NUM)
#define pi 3.1415926
#define VM 100047
#define EM 400047
#define rd(x) scanf("%d",&x);
#define seed() srand((unsigned)time(NULL))
#define random(a, b) rand() % (b - a + 1) + a
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;

const int inf = 0x3f3f3f3f;
const ll llf = 0x3f3f3f3f3f3f3f3f;
const int maxn = (int) 1e6 + 7;
const double eps = 1e-10;
const ll mod1 = (int) 1e9 + 7;
const ll mod2 = 998244353;
const ll has = 99959;
const int dx[] = {0, 1, 0, -1};
const int dy[] = {1, 0, -1, 0};



int main() {
    std::ios::sync_with_stdio(false);
#ifndef ONLINE_JUDGE
    freopen("in.txt", "r", stdin);
    freopen("out.txt", "w", stdout);
#endif
    ll n;
    rush() {
        scanf("%lld", &n);
        printf("Case %d: %lld\n",NUM,n - (ll)sqrt(n) - (ll) sqrt(n / 2));
    }
    return 0;
}