【BZOJ2005】【codevs1937】能量采集,数论练习之二维公约数求和

传送门1
传送门2
这里写图片描述
其实phi可以前缀和处理,(n/d)*(m/d)只有(n^0.5+m^0.5)个取值,然而我这里是O(n)的偏暴力流……
代码:

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>
using namespace std;
int phi[100010],prime[100010];
bool pd[100010];
long long n,m;
long long ans;
main()
{
    scanf("%d%d",&n,&m);
    if (n>m) swap(n,m);
    for (int i=1;i<=n;i++)
    {
        if (!pd[i]) phi[i]=i-1,prime[++prime[0]]=i;
        for (int j=1;j<=prime[0];j++)
        {
            if (prime[j]*i>n) break;
            pd[i*prime[j]]=1;
            if (i%prime[j]==0)
            phi[i*prime[j]]=phi[i]*prime[j];
            else phi[i*prime[j]]=phi[i]*phi[prime[j]];
        }
    }//欧拉筛求φ
    ans=n*m;
    for (int i=1;i<=n;i++)
    ans+=2*(n/i)*(m/i)*phi[i];
    printf("%lld",ans);
}
发布了374 篇原创文章 · 获赞 74 · 访问量 25万+
展开阅读全文

没有更多推荐了,返回首页

©️2019 CSDN 皮肤主题: 技术黑板 设计师: CSDN官方博客

分享到微信朋友圈

×

扫一扫,手机浏览