J——sum【尤拉篩選+思維規律】
阿新 • • 發佈:2019-02-15
題目大意:定義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;
}