HDU 6134(2017 多校训练:Battlestation Operational(莫比乌斯反演))

本文详细解析了HDU 6134题目的解题思路,利用欧拉函数与莫比乌斯反演将原问题转化为易于解决的形式,并给出了具体的实现代码。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题目链接:

http://acm.hdu.edu.cn/showproblem.php?pid=6134

这题就是求


考虑当Gcd(i, j)==1时,除了j为1的情况,其它时候i/j一定是小数,所以i/j向上取整相当于向下取整的结果+1。这里注意的是,题目要求的是向上取整,这里转化成了向下取整(因为出现向下取整这种情况代表i和j互质,而欧拉函数就是求互质个数的,这里还有减一代表j=1时的情况,因为这里确实是满足的,虽然欧拉函数算上1了)。

那么有:(其中φ(i)为小于i与i互质的对数,即欧拉函数)


欧拉函数因为是积性函数,可以线性求出,令


为什么上面等式成立?

对于所有的n/i,当n和i不互质时,很显然它们除掉它们的Gcd之后就互质了,而等式左边正是在枚举n的约数

也就是两种情况,一种是n和i互质,另一种是n和i互质。我们发现只有d为n时才能出现互质,因为不互质那么必然有公因子,相除以后就成了n的因子和比它小的数互质。后者n和i互质就是代表左式中的d=n的情况。

考虑莫比乌斯反演


那么有:


其中莫比乌斯函数和原公式中的欧拉函数都可以O(n)预处理

这样只要枚举后半部分就可以得出答案了

但很可惜这样复杂度仍然是O(nsqrt(n))超时!!


上面公式还可以转


其中τ(n)表示n的约数个数,是积性函数也可以线性求出

那么答案有


其中先预处理
τ(n)的前缀和,再预处理u(i)与t(n)前缀和积的前缀和

最后带上欧拉函数再算一波前缀和即可


代码:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<queue>
using namespace std;
#define LL long long
#define MAXN 1001000
const int mod=1e9+7;
int mu[MAXN];
LL ans[MAXN];
LL d[MAXN];
int phi[MAXN];
void getPhi()
{
    for(int i=1;i<MAXN;++i)
        phi[i]=i;
    for(int i=2;i<MAXN;++i)
        if(phi[i]==i)
            for(int j=i;j<MAXN;j+=i)
                phi[j]=phi[j]/i*(i-1);
}
void getMu()
{
    for(int i=1;i<MAXN;++i)
    {
        int target=i == 1?1:0;
        int delta=target-mu[i];
        mu[i]=delta;
        for(int j=i+i;j<MAXN;j+=i)
            mu[j]+=delta;
    }
}
void init()
{
    getPhi();
    getMu();
    for(int i=1;i<MAXN;++i)
        for(int j=i;j<MAXN;j+=i)
            d[j]++;
    for(int i=1;i<MAXN;++i)
        d[i]=(d[i]+d[i-1])%mod;
    for(int i=1;i<MAXN;++i)
        for(int j=i;j<MAXN;j+=i)
            ans[j]=(ans[j]+mu[j/i]*d[i])%mod;
    for(int i=1;i<MAXN;++i)
        ans[i]=(ans[i]+phi[i]-1)%mod;
    for(int i=2;i<MAXN;++i)
        ans[i]=(ans[i]+ans[i-1])%mod;
}
int main()
{
    init();
    int n;
    while(~scanf("%d",&n))
    {
        printf("%d\n",ans[n]);
    }
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值