SPOJ LCMSUM (数论)

转载请注明出处,谢谢http://blog.csdn.net/ACM_cxlove?viewmode=contents    by---cxlove

对于 n,求出sigma (LCM (i , n))    n >=  i >=1

重要的是复习一下推导过程。。。

LCM(i , n) = i * n / gcd (i , n)

所以我们按分母分类合并,gcd (i , n)显然是n的约数,sqrt(n)级别的。

g = gcd (i , n) ,将所有分母为g的求和。

即gcd (i / g , n / g) = 1,sigma(i / g) * g。

这里有个结论,对于 m >= 2,与m的互质的数的和为m * phi (m) / 2。

如果gcd (a , m) = 1那么gcd (m - a , m) = 1。所以所有的互质的数,可以两两分组。

那么分母为g的求和便是 n * phi(n / g) * n / g / 2

做法就比较显然了,枚举所有的g,即所有的约数,预处理所有欧拉函数 + Q * sqrt(n)。

遗憾的是会超时。

如果 令n / g = d,那么对于 n ,ans[n] = sigma (n * phi(d) * d / 2) 。 n | d

所以可以预处理出所有的ans[]。枚举d,更新n。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <set>
#include <map>
#include <string>
#include <vector>
#include <queue>
#include <stack>
#define lowbit(x) (x & (-x))
#define Key_value ch[ch[root][1]][0] 
#pragma comment(linker, "/STACK:1024000000,1024000000")    
using namespace std;
typedef long long LL;
const int N = 1000005;
int t , n;
LL phi[N] , ans[N];
void Init () {
    for (int i = 2 ; i < N ; i ++) {
        if (phi[i]) continue;
        for (int j = i ; j < N ; j += i) {
            if (!phi[j]) phi[j] = j;
            phi[j] = phi[j] * 1LL / i * (i - 1);
        }
    }
    ans[1] = 1LL;
    for (int i = 2 ; i < N ; i ++)
        ans[i] = (1LL * i * i * phi[i] >> 1) + i;
    for (int i = 2 ; i * i < N ; i ++) {
        ans[i * i] += 1LL * i * i * phi[i] * i >> 1;
        for (int j = i * i + i , k = i + 1 ; j < N ; j += i , k ++) {
            ans[j] += (1LL * j * phi[i] * i >> 1) + (1LL * j * k * phi[k] >> 1);
        }
    } 
}
int main () {
    #ifndef ONLINE_JUDGE
        freopen ("input.txt" , "r" , stdin);
        // freopen ("output.txt" , "w" , stdout);
    #endif
    Init ();
    scanf ("%d" , &t);
    while (t --) {
        scanf ("%d" , &n);
        printf ("%lld\n" , ans[n]);
    }
    return 0;
}


评论 1 您还未登录,请先 登录 后发表或查看评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
©️2022 CSDN 皮肤主题:大白 设计师:CSDN官方博客 返回首页

打赏作者

ACM_cxlove

你的鼓励将是我创作的最大动力

¥2 ¥4 ¥6 ¥10 ¥20
输入1-500的整数
余额支付 (余额:-- )
扫码支付
扫码支付:¥2
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值