BZOJ 2005 [Noi2010]能量采集 ——Dirichlet积

最新推荐文章于 2022-09-13 16:35:06 发布

weixin_34127717

最新推荐文章于 2022-09-13 16:35:06 发布

阅读量85

点赞数

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/SfailSth/p/6435474.html

本文介绍了一种利用卷积算法进行分块处理的方法，并通过O(n)的时间复杂度来解决特定数学问题。文章详细展示了如何初始化并计算phi值，再通过循环迭代逐步求得最终答案。

【题目分析】

卷积一卷。

然后分块去一段一段的求。

O(n)即可。

【代码】

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define maxn 100005 
#define ll long long
#define F(i,j,k) for (ll i=j;i<=k;++i)
ll n,m,pr[maxn],top,d;
ll phi[maxn];
 
void init()
{
    F(i,2,maxn-1)
    {
        if (!phi[i]) pr[++top]=i,phi[i]=i-1;
        for (ll j=1;j<=top&&(ll)pr[j]*i<(ll)maxn;++j)
        {
            if (i%pr[j]==0) {phi[i*pr[j]]=phi[i]*pr[j];break;}
            else phi[i*pr[j]]=phi[i]*phi[pr[j]];
        }
    }
}
 
int main()
{
    init();
    scanf("%d%d",&n,&m);
    F(i,1,maxn-1) phi[i]+=phi[i-1];
    ll ans=0;
    ll la,lb,nowa,nowb,l,r=n;
    while (r)
    {
        nowa=n/r; nowb=m/r;
        la=n/(nowa+1)+1;lb=m/(nowb+1)+1;
        l=max(la,lb);
        ans+=(ll)nowa*nowb*(phi[r]-phi[l-1]);
        r=l-1;
    }
    printf("%lld\n",n*m+2*ans);
}

转载于:https://www.cnblogs.com/SfailSth/p/6435474.html