J——sum【欧拉筛选+思维规律】

题目来源：2018ACM—ICPC南京赛区网络预选赛

题目大意：定义f[i]函数代表i=a*b的对数，其中a和b都不能是平方数的倍数，a*b与b*a不相同，t组样例，给出n，求1~n的f[i]之和。
题目总结：思路没到位，逼哥倒是想到了质因子筛选，但我俩苦于处理各种因子的情况，没想到答案这种思路，就是每个质因子的倍数决定了它对f[n]函数的贡献程度。还有就是欧拉筛用的不熟练，学习一波~

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <math.h>
using namespace std;

typedef long long ll;

/*
    题解思路：筛选质因子+思维
    首先，由题意，我们可以将数分解为p1^e1*p2^e2....*pn^en,
    其中p1,p2...pn都为质因子，那么f[n]与这个有什么关系呢？
    遍历所有的ei，有三种情况：
    1.ei=1,此时说明有满足条件的质因子，f[i]=f[i]*2;
    2.ei=2,此时f[i]不变;
    3.ei>2,此时f[i]=0;
    所以我们可以一次筛选出2e7内的所有素数，由素数表筛选出数n的质因子，这样我们
    可以将数n写为n=prime[i]*j,判断j与prime[i]的关系，就可以确定这个质因子所贡献的f[n]的值

*/

const int maxn=2e7+1;

const int maxp=2*7000000;

int prime[maxp];

int vis[maxn];

ll f[maxn];//每个数的因子组合

int tot;

//筛选质因子
void sieve()
{
    tot=0;
    memset(vis,0,sizeof(vis));
    memset(f,0,sizeof(f));
    f[1]=1;
    //欧拉筛素数法
    for(int i=2;i*i<=maxn;i++)
    {
        for(int j=i*i;j<=maxn;j+=i*i)
        {
            vis[j]=1;
        }
    }
    //生成素数表，并将素数的f值置为2
    for(int i=2;i<maxn;i++)
    {
        if(!vis[i])
        {
            prime[tot++]=i;
            f[i]=2;
        }
        for(int j=0;j<tot;j++)
        {
            int k=prime[j]*i;
            if(k>=maxn)break;
            vis[k]=1;
            if(i%prime[j])f[k]=f[i]*2;
            else{
                int tp=prime[j]*prime[j];
                if(i%tp==0)f[k]=0;
                else{
                    f[k]=f[k/tp];
                }
                break;

            }
        }
    }
}


int main()
{
    int n;
    int T;
    sieve();
    scanf("%d",&T);

    for(int i=1;i<maxn;i++)
    {
        f[i]+=f[i-1];
    }
    while(T--)
    {
        scanf("%d",&n);
        printf("%lld\n",f[n]);
    }
    return 0;
}