HDU 2608 0 or 1(数论)

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=2608

此题想要打表找规律其实并不容易，我也是看了题解严格的证明之后AC的。在此凭借对此题的理解用自己的话推导一遍答案。

     首先，题目规定了T[n]表示1~n的数当中能被n整除的数之和，S[n]表示(T[1]+T[2]+...+T[n])%2的值。因为最终问的是S[n]的值，所以对于某个T[i]的具体值是多少已没有意义，只需要知道T[i]的奇偶性即可。而此时，S[n]也转化成了对于1~n这n个数中，S[n] =(T[]值为奇数的个数)%2。所以，如果我们能知道所有的“T值为奇数”的数，存到一个数组中，然后对于每次询问值x，二分查找到x前有多少个满足“T值为奇数”的数，%2得出答案即可。【事实上满足条件的数只有不到8万个】

     然后，我们需要解决的就是找到什么样的数，满足其T值为奇数？因为对于每一个x值，都可以拆分成若干个素数的乘积，形如x = 2^e1*3^e2*5^e3*...*pi^ei。而x的任意一个因数y = 2^a1*3^a2*5^a3*...*pi^ai中的每个素数的幂次ai，都必然是0~ei中的某个值。所以所有的因数和，即T[x] = ( 2^0+2^1+...+2^e1 )*( 3^0+3^1+...+3^e2 )*...*( pm^0+pm^1+...+pm^em )。现在我们讨论T[x]的奇偶性，因为除2以外所有的素数均为奇数，所以只有当除2以外的所有素数幂次均为偶数的时候，T[x]%2 == 1。

AC代码：

#include <iostream>
#include <cstdio>
#include <vector>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long LL;
#define MAXN 2147483648
int t;
vector<LL> vec, vv;
void init() {
    for(LL i = 1; i * i <= MAXN; i++) {
        LL ret = i * i;
        vv.push_back(ret);
        vec.push_back(ret);
    }
    int len = vv.size();
    for(int i = 0; i < len; i++) {
        LL tmp = 2;
        while(1) {
            LL tp = tmp * vv[i];
            if(tp < MAXN) vec.push_back(tp);
            else break;
            tmp *= 2;
        }
    }
    sort(vec.begin(), vec.end());
    vec.erase(unique(vec.begin(), vec.end()), vec.end());
}
int main() {
    init();
    vec.push_back(MAXN + 10);
    while(~scanf("%d", &t)) {
        while(t--) {
            LL x;
            scanf("%I64d", &x);
            int ans = lower_bound(vec.begin(), vec.end(), x + 1) - vec.begin();
            printf("%d\n", ans % 2);
        }
    }
    return 0;
}