HDU 多校 6434(欧拉函数)

版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/ffgcc/article/details/81951917

解析:
其实规律比较好找
即对于每个 i, 求有多少个小于它的 a 满足 gcd(i, a) = 1 且 a 是奇数
当 i 是奇数时, 答案为 φ(i)/2
当 i 是偶数时, 答案为 φ(i).
注意 i = 1 时, 答案为 0.
记个前缀和就好了, 复杂度为 O(N + T)

#include<bits/stdc++.h>
using namespace std;
#define LL long long
const int N=2e7;
const int maxn=N+10;
int T,n;
int primes,prime[maxn],phi[maxn];
bool vis[maxn];
LL Ans[maxn];
void init()
{
    int i,j;
    phi[1]=1;
    for (i=2;i<=N;i++)
    {
        if (!vis[i])
        {
            prime[++primes]=i;
            phi[i]=i-1;
        }
        for (j=1;j<=primes&&i*prime[j]<=N;j++)
        {
            vis[i*prime[j]]=1;
            if (i%prime[j])
                phi[i*prime[j]]=phi[i]*(prime[j]-1);
            else
            {
                phi[i*prime[j]]=phi[i]*prime[j];
                break;
            }
        }
    }
    for (i=2;i<=N;i++)
        Ans[i]=Ans[i-1]+((i%2!=0)?phi[i]/2:phi[i]);
}
int main()
{
    #ifdef local
        freopen("D://r.txt","r",stdin);
    #endif
    init();
    scanf("%d",&T);
    while (T--)
        scanf("%d",&n),printf("%lld\n",Ans[n]);
}

猜你喜欢

转载自blog.csdn.net/ffgcc/article/details/81951917
今日推荐