[51NOD1244]莫比乌斯函数之和

本文介绍了一种利用杜教筛算法解决特定数学问题的方法。针对给定的整数范围,文章详细阐述了如何高效地计算莫比乌斯函数的前缀和,并提供了具体的代码实现。通过对小于特定阈值的数进行预处理,结合分块计算策略,实现了良好的时间复杂度。

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

题目大意

给定 a b,试求出

i=abμ(i)

2ab1010


题目分析

杜教筛裸题。
S(n)=ni=1μ(i) 。用恒等函数 1 卷上 μ 函数得到单位函数 ϵ ,于是有:

S(n)=1i=2nS(ni)

只需要计算形如 nx 的数的 S 值。对小于等于n23的数使用线性筛预处理,对于其它的数使用分块来计算。
积分算得时间复杂度:
On13x=1nxdx=O(n23)

更多积性函数前缀和相关请参考2016年任之洲的国家集训队论文《积性函数求和的几种方法》以及 唐老师的博客


代码实现

第一道杜教筛,写得及其丑。
关于那个 S(nx) 怎么存,对于 nx 开一个表直接存,对于 nx>n 另外开一个表存在 x 这个位置。

#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdio>
#include <cmath>

using namespace std;

typedef long long LL;

const int L=4641600;
const int N=100050;
const int M=N<<1;

int mu[L],sum[L],pri[L],f[L];
bool mark[L];
LL S[2][N];
LL num[M];
int l,s,cnt;
LL a,b;

int id(LL x,LL n){return x>s?n/x:x;}

void pre(LL n)
{
    memset(mark,0,sizeof mark),memset(pri,0,sizeof pri);
    l=pow(n,2.0/3.0),s=sqrt(n);
    mark[1]=1,f[1]=1,mu[1]=1;
    for (int i=2;i<=l;++i)
    {
        if (!mark[i]) mark[i]=1,f[pri[++pri[0]]=i]=i,mu[i]=-1;
        for (int j=1,k;j<=pri[0];++j)
        {
            if (1ll*pri[j]*i>l) break;
            mark[k=pri[j]*i]=1,f[k]=pri[j];
            mu[k]=f[k]==f[i]?0:-mu[i];
            if (!(i%pri[j])) break;
        }
    }
    for (int i=1;i<=l;++i) sum[i]=sum[i-1]+mu[i];
    cnt=0;
    for (LL st=1,en,x;st<=n;st=en) num[++cnt]=x=n/st,en=n/x+1;
}

LL sieve(LL n)
{
    pre(n);
    for (int sign,idn,i=cnt;i>=1;--i)
    {
        idn=id(num[i],n),sign=idn==num[i];
        if (num[i]<=l)
        {
            S[sign][idn]=sum[num[i]];
            continue;
        }
        LL res=1;
        for (LL st=2,en,x;st<=num[i];st=en)
        {
            x=num[i]/st,en=num[i]/x+1;
            int idn_=id(x,n),sign_=idn_==x;
            res-=(en-st)*S[sign_][idn_];
        }
        S[sign][idn]=res;
    }
    return S[0][1];
}

int main()
{
    freopen("mobius.in","r",stdin),freopen("mobius.out","w",stdout);
    scanf("%lld%lld",&a,&b);
    printf("%lld\n",sieve(b)-sieve(a-1));
    fclose(stdin),fclose(stdout);
    return 0;
}
### 关于51Nod平台上编号为1020的问题详情与解答 #### 问题描述 在51Nod平台上的第1020号问题是关于计算两个大整数相加的结果[^1]。给定两个正整数A和B,长度不超过10^6位,要求编写程序来求解这两个数的和。 #### 输入格式说明 输入数据由多组测试案例组成;每组测试案例占两行,分别表示要相加的大整数A和B。对于每一组测试案例,应当单独输出一行结果,即A+B的值。 #### 解决方案概述 解决此问题的关键在于处理超大数据类型的运算,在大多数编程语言中内置的数据类型无法直接支持如此大规模数值的操作。因此,可以采用字符串的方式来存储这些大整数,并实现逐位相加逻辑,同时考虑进位情况。 下面是一个Python版本的具体实现方法: ```python def add_large_numbers(a: str, b: str) -> str: # Reverse strings to make addition easier from least significant digit a = a[::-1] b = b[::-1] carry = 0 result = [] max_length = max(len(a), len(b)) for i in range(max_length): digit_a = int(a[i]) if i < len(a) else 0 digit_b = int(b[i]) if i < len(b) else 0 total = digit_a + digit_b + carry carry = total // 10 current_digit = total % 10 result.append(str(current_digit)) if carry != 0: result.append(str(carry)) return ''.join(reversed(result)) if __name__ == "__main__": while True: try: num1 = input().strip() num2 = input().strip() print(add_large_numbers(num1, num2)) except EOFError: break ``` 该代码片段定义了一个函数`add_large_numbers`用于接收两个作为参数传入的大整数(形式上为字符串),并返回它们之和同样作为一个字符串。通过反转输入字符串使得最低有效位位于索引位置0处从而简化了按位累加的过程。最后再将得到的结果列表反向拼接成最终答案输出。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值