[BZOJ2226][Spoj5971][数论]LCMSum[好题]

最新推荐文章于 2019-02-18 19:34:04 发布

BenjaminPMLee

最新推荐文章于 2019-02-18 19:34:04 发布

阅读量1.9k

点赞数

分类专栏： OI - Number Theory 文章标签： bzoj 数论欧拉phi函数 spoj

本文链接：https://blog.csdn.net/BenjaminPMLee/article/details/44947809

版权

OI - Number Theory 专栏收录该内容

1 篇文章 0 订阅

订阅专栏

题干

T组询问，每组给出一个n,求：

\sum 1 ⩽ i ⩽ n l c m (i, n), (1 ⩽ T ⩽ 300000, 1 ⩽ n ⩽ 1000000)

$\sum_{1\leqslant i\leqslant n}{lcm(i,n)},\;(1 \leqslant T \leqslant 300000,\;1\leqslant n \leqslant 1000000)$

题解

先进行变形：

\sum 1 ⩽ i ⩽ n l c m (i, n) = \sum 1 ⩽ i ⩽ n i \times n g c d ( i , n )

$\sum_{1\leqslant i\leqslant n}{lcm(i,n)} =\sum_{1\leqslant i\leqslant n}{\frac{i \times n}{gcd(i,n)}}$
枚举gcd(i, n)，有：

\sum 1 ⩽ i ⩽ n i \times n g c d ( i , n ) = n \times \sum d | n ⎛ ⎝ \sum 1 ⩽ i ⩽ n i d \times [g c d (i, n) = d] ⎞ ⎠ = n \times \sum d | n ⎛ ⎝ ⎜ \sum 1 ⩽ j ⩽ n d j \times [g c d (j, n d) = 1] ⎞ ⎠ ⎟, (j = i d)

$\sum_{1 \leqslant i \leqslant n}{\frac{i \times n}{gcd(i,n)}}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\\ = n \times \sum_{d|n}{\left( \sum_{1 \leqslant i \leqslant n}{ \frac{i }{d} \times \left[gcd(i,n)=d\right]} \right) } \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\\ = n \times \sum_{d|n}{\left( \sum_{1 \leqslant j \leqslant \frac{n}{d}}{ j \times \left[gcd(j,\frac{n}{d})=1\right]} \right)},\left( j=\frac{i}{d} \right)$
显然d与n/d是等价的，于是：

\sum 1 ⩽ i ⩽ n l c m (i, n) = n \times \sum d | n ⎛ ⎝ \sum 1 ⩽ i ⩽ n i \times [g c d (i, d) = 1] ⎞ ⎠

$\sum_{1\leqslant i\leqslant n}{lcm(i,n)}=n \times \sum_{d|n}{\left( \sum_{1 \leqslant i \leqslant n}{ i\times \left[gcd(i,d)=1\right]} \right) }$
即我们要求小于i且与i互素的所有自然数之和。有结论：

$V$

由于本人数学功底不好所以就不证明了。直接利用结论：预处理phi函数，对于n枚举因数即可。

代码

（不要随便看代码，不要随便看代码，不要随便看代码。因为很重要所以说三遍。）

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <cmath>
using namespace std;

//Global Variables & Definitions
#define LL long long
#define MAXN 1000010
int T, _n;
//End Global Variables & Definitions

//Number Theory
int p[MAXN];
int vis[MAXN];
int ph[MAXN];
int pcnt = 0;

void init_p() {
    ph[1] = 1; ph[2] = 1;

    LL temp;
    for(int i = 2;i < MAXN;++i) {
        if(!vis[i]) { p[pcnt] = i; ph[i] = i - 1; ++pcnt; }

        for(int j = 0;j < pcnt && (temp = (LL)p[j] * i) < MAXN;++j) {
            vis[temp] = 1;

            if(i % p[j] == 0) { ph[temp] = ph[i] * p[j]; break; }
            else ph[temp] = ph[i] * (p[j] - 1);
        }
    }

    //for(int i = 1;i <= 20;++i) printf("phi[%d] = %d\n", i, ph[i]);
}

inline LL phi(int n) {
    return ph[n];
}

LL solve(int n) {
    LL ans = 0ll;   
    int half = (int) (sqrt(n) + 0.01);

    if(half * half == n) { ans += phi(half) * half / 2; --half; }
    ans += 1; //for 1
    ans += phi(n) * n / 2;
    for(int i = 2;i <= half;++i)
        if(n % i == 0) {
            ans += phi(i) * i / 2;
            ans += phi(n / i) * n / i / 2;
        }

    return ans * n;
}
//End Number Theory

//Main Structure
inline void ir() {
    scanf("%d", &T);
    init_p();
}

int main() {
    ir();

    while(T--) {
        scanf("%d", &_n);
        printf("%lld\n", solve(_n));
    }
    return 0;
}