brz的函数

brz的函数

https://ac.nowcoder.com/acm/contest/8282/D

推公式题。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn = 5e4+7;

ll ans[maxn], mu_fac_sum[maxn];
bool vis[maxn];
int prime[maxn];
int Mob[maxn];
void init(){
     int cnt = 0;
     vis[1] = 1;
     Mob[1] = 1;
     for(int i = 2; i <= maxn; i++){
        if(!vis[i])
            prime[cnt++] = i, Mob[i] = - 1;
        for(int j = 0; j < cnt && 1LL * prime[j] * i <= maxn; j++){
            vis[prime[j] * i] = 1;
            Mob[i * prime[j]] = (i % prime[j] ? -Mob[i]: 0);
            if(i % prime[j] == 0)
                break;
        }
    }

    ans[1] = 1;
    mu_fac_sum[1] = 1;
    for(int i = 2; i < maxn; i++)
    {
        ans[i] = ans[i-1];
        for(int j = 1; j * j <= i; j++){
            if(i % j == 0)
            {
                ans[i] += Mob[j]*(Mob[i]*Mob[i] + 2*Mob[i]*mu_fac_sum[j]);
                mu_fac_sum[j] += Mob[i];

                if(j * j != i)
                {
                    ans[i] += Mob[i/j]*(Mob[i]*Mob[i] + 2*Mob[i]*mu_fac_sum[i/j]);
                    mu_fac_sum[i/j] += Mob[i];
                }
            }
        }
    }
}

int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0); cout.tie(0);
    init();
    int t;
    cin >> t;
    while(t--)
    {
        int n;
        cin >> n;
        cout << ans[n] << endl;
    }
}
全部评论

相关推荐

06-26 15:33
青岛工学院 Java
积极的秋田犬要冲国企:他现在邀请我明天面试
点赞 评论 收藏
分享
不要停下啊:大二打开牛客,你有机会开卷了,卷起来,去找课程学习,在牛客上看看大家面试笔试都需要会什么,岗位有什么需求就去学什么,努力的人就一定会有收获,这句话从来都经得起考验,像我现在大三了啥也不会,被迫强行考研,炼狱难度开局,啥也不会,找工作没希望了,考研有丝丝机会
点赞 评论 收藏
分享
评论
2
收藏
分享

创作者周榜

更多
牛客网
牛客网在线编程
牛客网题解
牛客企业服务