Sum 南京网络赛J题

题意：

　　统计每个数的因子的对数，如果因子能被某个平方数整除，则不统计在内，每对因子有序

解析：

　　我们对某个数n进行质因子分解，如果某个质因子的指数大于2则 f(n) = 0,

例 N = X³ * M = R * T

因为要分成两部分 所以无论怎样分 R 或 T 总有一部分X的指数大于等于2

如果指数为2 则只能是R 和 T 两部分每个部分有一个X，所以只有一种情况，贡献为1

如果指数为1 则这个X可以在R 和 T的任意一部分 ，所以有两种情况

#include <iostream>
#include <cstdio>
#include <sstream>
#include <cstring>
#include <map>
#include <cctype>
#include <set>
#include <vector>
#include <stack>
#include <queue>
#include <algorithm>
#include <cmath>
#include <bitset>
#define rap(i, a, n) for(int i=a; i<=n; i++)
#define rep(i, a, n) for(int i=a; i<n; i++)
#define lap(i, a, n) for(int i=n; i>=a; i--)
#define lep(i, a, n) for(int i=n; i>a; i--)
#define rd(a) scanf("%d", &a)
#define rlld(a) scanf("%lld", &a)
#define rc(a) scanf("%c", &a)
#define rs(a) scanf("%s", a)
#define pd(a) printf("%d
", a);
#define plld(a) printf("%lld
", a);
#define pc(a) printf("%c
", a);
#define ps(a) printf("%s
", a);
#define MOD 2018
#define LL long long
#define ULL unsigned long long
#define Pair pair<int, int>
#define mem(a, b) memset(a, b, sizeof(a))
#define _  ios_base::sync_with_stdio(0),cin.tie(0)
//freopen("1.txt", "r", stdin);
using namespace std;
const int maxn = 20000010, INF = 0x7fffffff, LL_INF = 0x7fffffffffffffff;
int prime[maxn+1];
int minprime[maxn+1];
LL a[maxn+1];

void get_prime()
{
    mem(prime, 0);
    for(int i=2; i<=maxn; i++)
    {
        if(!prime[i]) prime[++prime[0]] = i, minprime[i] = i;
        for(int j=1; j<=prime[0] && prime[j] <= maxn/i; j++)
        {
            prime[prime[j]*i] = 1;
            minprime[prime[j]*i] = prime[j];
            if(i % prime[j] == 0) break;
        }
    }
}

void init()
{
    a[0] = 0;
    a[1] = 1;
    for(int i=2; i<=maxn; i++)
    {
        int mm = minprime[i];
        if((LL)mm*mm < maxn && (LL)mm*mm*mm < maxn && i%(mm*mm*mm) == 0)
            a[i] = 0;
        else if((LL)mm*mm < maxn && i%(mm*mm) == 0) a[i] = a[i/mm/mm];
        else a[i] = 2*a[i/mm];
    }
    for(int i=1; i<=maxn; i++)
        a[i] += a[i-1];
}

int main()
{
    get_prime();
    init();
    int T;
    rd(T);
    while(T--)
    {
        int n;
        rd(n);
        pd(a[n]);
    }

    return 0;
}